DOI: https://doi.org/10.1038/s41598-026-36077-1
PMID: https://pubmed.ncbi.nlm.nih.gov/41547957
تاريخ النشر: 2026-01-17
المؤلف: Swarnendu Bhuyan وآخرون
الموضوع الرئيسي: التحليل الجيولوجي والجيوكيميائي
نظرة عامة
تبحث الدراسة في تكوين رواسب الذهب من نوع اللود في أواخر العصر الأركي، والتي تُنسب تقليديًا إلى ضخ الزلازل للسوائل الحاملة للذهب عبر آلية “صمام العطل”. باستخدام نموذج عددي ثنائي الأبعاد تم تنفيذه في COMSOL Multiphysics، تحلل الدراسة بشكل كمي تراكم ضغط السوائل وديناميات التدفق من خلال حل معادلات مرتبطة تأخذ في الاعتبار الخصائص المعتمدة على درجة الحرارة للصخور والسوائل. تشير النتائج إلى أن آلية صمام العطل قد لا تكون مستدامة لدورات تدفق السوائل المستمرة، وأن خصائص الغطاء الزلزالي، بما في ذلك وضعه ونفاذيته، تؤثر بشكل كبير على ديناميات ضغط السوائل.
علاوة على ذلك، تشير الدراسة إلى أن العيوب ذات الزاوية العالية تعمل كقنوات أكثر فعالية لحركة السوائل مقارنة بالعيوب المائلة برفق. تفترض أن الضغط الفوقي يمكن تحقيقه من خلال الضغط البعيد بدلاً من الاعتماد فقط على تدفق الحرارة الذي يتجاوز التدرج الحراري. من المهم أن تستنتج الدراسة أن الغطاء الزلزالي ليس شرطًا لعمليات التمعدن التي لوحظت في بعض معسكرات الذهب، مثل سيغما وهوتي، مما يستدعي الحاجة لاستكشاف آليات بديلة تتجاوز نموذج صمام العطل لفهم هذه الأنظمة المعدنية بشكل أفضل.
مقدمة
تناقش المقدمة أهمية رواسب الذهب من نوع اللود المتوسطة الحرارة الموجودة في تيرانات الجرانيت-الحجر الأخضر من أواخر العصر الأركي، والتي تنتشر عبر الكتل القارية وتلعب دورًا حيويًا في إمدادات الذهب العالمية. تُعرف هذه الرواسب باسم رواسب الذهب “الأوروجينية”، وترتبط بإعدادات جيولوجية معقدة تتميز بتشوهات متعددة وتحولات، ناتجة عن تقارب البروتوقارات الأركية المتأخرة وعمليات البراكين البحرية الرسوبية المبكرة. يتم تصور نظام رواسب الذهب الأوروجينية كنظام ديناميكي حراري يتحكم فيه التدفق، حيث يسهل التشوه، وتدفق السوائل، ونقل الحرارة، والتفاعلات الكيميائية نقل وت deposition الذهب.
تم اقتراح نموذج متعدد المراحل لشرح أصل هذه الرواسب، مرتبطًا بأنظمة الهيدروحرارية في قاع البحر (VMS) وإعادة التحريك اللاحقة أثناء التشوه والتحول. يفترض هذا النموذج المتحولي، الذي اكتسب زخمًا في أواخر الثمانينيات، أن السوائل التي يتم إطلاقها من التحلل المتحولي المتقدم لتسلسلات الحجر الأخضر تنقل الذهب إلى مواقع هيكلية في القشرة. المركز في هذا النموذج هو وجود طبقة غير نفاذة على أعماق تتراوح بين 10-15 كم، مما يشير إلى الانتقال من الأنظمة الهشة إلى الأنظمة اللدنة. يقترح النموذج أن ضغط السوائل يتراكم تحت هذه الطبقة، مما يؤدي إلى تمزقات يمكن أن تسبب نشاطًا زلزاليًا، وبالتالي ربط عمليات تمعدن الذهب بالديناميات التكتونية.
طرق البحث
في هذه الدراسة، استخدمنا COMSOL Multiphysics v6 لتحليل نظام كاديلاك-لاردير ليك، المعروف برواسبه العديدة من الذهب الأوروجيني. تضمنت المنهجية استخدام هذا البرنامج المتقدم لمحاكاة الخصائص الجيولوجية والجيوفيزيائية للنظام، مما سمح بإجراء تحقيق مفصل في العوامل التي تؤثر على تمعدن الذهب. مكنتنا هذه الطريقة من دمج مصادر بيانات متنوعة وإجراء محاكاة معقدة، مما ساهم في فهم أفضل للتوزيع المكاني وعمليات التكوين لرواسب الذهب داخل المنطقة.
النتائج
توضح النتائج المقدمة في هذا القسم (الأشكال 4-11) تأثير المعلمات الفيزيائية والليثولوجية المختلفة على تطور ضغط السوائل ضمن تكوينات مختلفة لنماذج صمام العطل. تبدأ التحليل بغطاء زلزالي أفقي، مستكشفة ثلاثة سيناريوهات متميزة: (أ) حالة ما قبل التمزق، (ب) عيب مائل برفق بزاوية 30° يتقاطع مع الختم، و(ج) عيب مائل بشدة بزاوية 60° (الشكل 4). بالإضافة إلى ذلك، تم تعديل النماذج لأخذ في الاعتبار غطاء زلزالي منحني بنصف قطر انحناء يبلغ 1 كم (الأشكال 4د-4و). تم تقديم تعديل للحد من نشاط العيوب إلى الحدود العليا للغطاء الزلزالي (الشكل 5)، وتمت دراسة تأثيرات ختم شبه نفاذ بنفاذية قدرها $10^{-23} \, \text{m}^2$ أيضًا (الشكل 6).
تشير النتائج إلى أن حتى انخفاضًا طفيفًا في ضغط السوائل يمكن أن يؤدي إلى ترسيب الذهب بسبب تغييرات طفيفة في الرقم الهيدروجيني، حيث يؤدي زيادة بمقدار 0.5-1.0 وحدة في الرقم الهيدروجيني إلى تقليل تركيزات الذهب المذاب (Au) بحوالي 3-10 مرات. لوحظ انخفاض كبير في ضغط السوائل بعد حدث تمزق، مع تصنيف قيم الضغط الزائد المحسوبة على أنها “فوق ضغط الليثوستاتي”، متجاوزة تراكم ضغط السوائل المبلغ عنه سابقًا والذي بلغ 320 ميجا باسكال في معسكر الذهب هوتي. كما تسلط النتائج الضوء على ظاهرة تدفق السوائل المركزة على العيوب، والتي تحدث عندما يتقاطع عيب مع نظام سوائل مضغوط، مما يسهل هجرة السوائل إلى مستويات قشرية أكثر ضحالة، كما يتضح من أسهم مجال سرعة دارسي في الشكل 4.
المناقشة
في هذه المناقشة، يقوم المؤلفون بتقييم نقدي لنموذج الضخ الزلزالي المتحول (MSP)، الذي تم قبوله على نطاق واسع في سياق رواسب الذهب الأوروجينية. يبرزون عدة قضايا غير محلولة تتعلق بالنموذج، خاصة فيما يتعلق بدور “الغطاء الزلزالي” – طبقة ذات نفاذية قريبة من الصفر يُفترض وجودها على أعماق تتراوح بين 12-15 كم. يجادل المؤلفون بأن العلاقة المكانية بين هذا الغطاء والتمعدن لا تزال غير واضحة، وأن الافتراض بأن الذهب يتم نقله بكفاءة بواسطة سوائل عالية الضغط بالقرب من الغطاء الزلزالي هو موضع تساؤل. يؤكدون على الحاجة إلى تقييم كمي للآليات وراء تراكم ضغط السوائل الفوقي وتأثيرات اتجاه العيب على هجرة السوائل.
كما ينتقد المؤلفون النماذج الحالية، بما في ذلك تلك التي اقترحها شيلدون وأورد، والتي تقترح أن العيوب المائلة بشدة أكثر فعالية لتدفق السوائل. يقترحون تعديلات على نموذج صمام العطل، مؤكدين أن نفاذية الختم ضرورية لتراكم ضغط السوائل، وأن العيوب المائلة بشدة قد تسهل بالفعل هجرة أفضل للسوائل مقارنة بالعيوب ذات الزاوية المنخفضة. علاوة على ذلك، يقدمون مفهوم آلية “تبديل الوضع”، والتي يمكن أن تفسر ديناميات السوائل دون الاعتماد على وجود غطاء زلزالي. بشكل عام، يدعو المؤلفون إلى إعادة تقييم نموذج MSP وفرضياته، داعين إلى فهم أكثر دقة للتفاعل بين ميكانيكا العيوب، وديناميات السوائل، وعمليات التمعدن في أنظمة الذهب الأوروجينية.
DOI: https://doi.org/10.1038/s41598-026-36077-1
PMID: https://pubmed.ncbi.nlm.nih.gov/41547957
Publication Date: 2026-01-17
Author(s): Swarnendu Bhuyan et al.
Primary Topic: Geological and Geochemical Analysis
Overview
The research investigates the formation of late-Archean lode-type gold deposits, traditionally attributed to seismic pumping of auriferous fluids via a ‘fault-valve’ mechanism. Utilizing a 2-D numerical model implemented in COMSOL Multiphysics, the study quantitatively analyzes fluid pressure buildup and flow dynamics by solving coupled equations that account for temperature-dependent properties of rocks and fluids. The findings suggest that the fault-valve mechanism may not be sustainable for continuous fluid flow cycles, and the characteristics of the seismic lid, including its disposition and permeability, significantly influence fluid pressure dynamics.
Moreover, the research indicates that high-angle faults serve as more effective conduits for fluid movement compared to gently dipping reverse faults. It posits that supralithostatic pressure can be achieved through far-field stress rather than solely relying on heat flux exceeding the geothermal gradient. Importantly, the study concludes that a seismic lid is not a requisite for the mineralizing processes observed in certain gold camps, such as Sigma and Hutti, prompting the need to explore alternative mechanisms beyond the fault-valve model to better understand these mineralizing systems.
Introduction
The introduction discusses the significance of mesothermal lode-type gold deposits found in late Archean granite-greenstone terranes, which are prevalent across cratonic blocks and play a crucial role in the global gold supply. These deposits, termed ‘orogenic’ gold deposits, are associated with complex geological settings characterized by multiple deformations and metamorphosis, arising from the convergence of late-Archean protocontinents and early seafloor volcanosedimentary processes. The orogenic gold deposit system is conceptualized as a flow-controlled thermodynamic system, where deformation, fluid flow, heat transfer, and chemical reactions facilitate the transport and deposition of gold.
A multi-stage model has been proposed to explain the origin of these deposits, linking them to sea-floor hydrothermal systems (VMS) and subsequent remobilization during deformation and metamorphism. This metamorphic model, which gained traction in the late 1980s, posits that fluids released from the prograde metamorphic devolatilization of greenstone sequences transport gold to structural sites in the crust. Central to this model is the existence of an impermeable layer at depths of 10-15 km, marking the transition from brittle to ductile regimes. The model suggests that fluid pressure builds up beneath this layer, leading to ruptures that can cause seismic activity, thereby linking gold mineralization processes to tectonic dynamics.
Methods
In this study, we employed COMSOL Multiphysics v6 to analyze the Cadillac-Larder Lake system, which is known for its numerous orogenic gold deposits. The methodology involved utilizing this advanced simulation software to model the geological and geophysical characteristics of the system, allowing for a detailed investigation of the factors influencing gold mineralization. This approach enabled us to integrate various data sources and perform complex simulations, ultimately contributing to a better understanding of the spatial distribution and formation processes of the gold deposits within the region.
Results
The results presented in this section (Figs. 4-11) elucidate the impact of various physical and lithological parameters on the evolution of fluid pressure within different configurations of fault-valve models. The analysis begins with a horizontal seismic lid, exploring three distinct scenarios: (a) a pre-rupture condition, (b) a gently dipping fault at a 30° angle intersecting the seal, and (c) a steeply dipping fault at 60° (Fig. 4). Additionally, models were adapted to consider a curved seismic lid with a 1 km radius of curvature (Figs. 4d-4f). A modification was introduced to limit fault activity to the upper boundary of the seismic lid (Fig. 5), and the effects of a semipermeable seal with a permeability of $10^{-23} \, \text{m}^2$ were also examined (Fig. 6).
The findings indicate that even a slight drop in fluid pressure can trigger gold precipitation due to minor pH changes, with a 0.5-1.0 pH unit increase reducing dissolved gold (Au) concentrations by approximately 3-10 times. A significant reduction in fluid pressure is observed following a rupture event, with computed fluid overpressure values classified as ‘supralithostatic’, surpassing the previously reported fluid pressure build-up of 320 MPa in the Hutti gold camp. The results also highlight the phenomenon of fault-focused fluid flow, which occurs when a fault intersects a pressurized fluid regime, facilitating fluid migration to shallower crustal levels, as illustrated by the Darcy velocity field arrows in Figure 4.
Discussion
In this discussion, the authors critically evaluate the metamorphogenic cum seismic pumping (MSP) model, which has been widely accepted in the context of orogenic gold deposits. They highlight several unresolved issues regarding the model, particularly concerning the role of the ‘seismic lid’—a near-zero permeability layer hypothesized to exist at depths of 12-15 km. The authors argue that the spatial relationship between this lid and mineralization remains unclear, and the assumption that gold is efficiently transported by high-pressure fluids near the seismic lid is questionable. They emphasize the need for a quantitative assessment of the mechanisms behind supralithostatic fluid pressure build-up and the implications of fault orientation on fluid migration.
The authors also critique the existing models, including those by Sheldon and Ord, which suggest that steeply dipping faults are more effective for fluid flow. They propose modifications to the fault-valve model, asserting that the permeability of the seal is crucial for fluid pressure build-up, and that steeply dipping faults may indeed facilitate better fluid migration than low-angle faults. Furthermore, they introduce the concept of a ‘mode-switching’ mechanism, which could explain fluid dynamics without relying on the presence of a seismic lid. Overall, the authors call for a reevaluation of the MSP model and its assumptions, advocating for a more nuanced understanding of the interplay between fault mechanics, fluid dynamics, and mineralization processes in orogenic gold systems.