本發(fā)明涉及礦產(chǎn)勘查，尤其涉及基于多源地質(zhì)數(shù)據(jù)融合的卡林型金礦找礦預(yù)測方法。
背景技術(shù)：
1、卡林型金礦的找礦預(yù)測長期依賴于對地質(zhì)、地球物理和地球化學信息的綜合解譯。現(xiàn)有方法通常對這些多源數(shù)據(jù)進行疊加或定性分析，以圈定異常區(qū)域。在構(gòu)造分析方面，常規(guī)技術(shù)側(cè)重于對斷裂、褶皺等構(gòu)造形跡的形態(tài)描述與分級，或利用地球物理資料反演深部構(gòu)造輪廓。對于三維建模，現(xiàn)有技術(shù)多基于解譯成果構(gòu)建靜態(tài)的實體模型，用以可視化地質(zhì)體與礦化的空間關(guān)系。
2、這些現(xiàn)有方法存在局限。傳統(tǒng)的構(gòu)造分析未能有效區(qū)分不同成因和功能的構(gòu)造要素，即深部導(dǎo)礦構(gòu)造與淺部容礦構(gòu)造往往被混為一談或僅作簡單關(guān)聯(lián)，導(dǎo)致對流體運移通道與沉淀場所空間配置關(guān)系的理解停留在概念層面，無法定量刻畫其耦合機制。同時，常規(guī)的三維地質(zhì)建模成果是一個確定性模型，缺乏對模型自身可靠性的量化評價。模型構(gòu)建所依賴的參數(shù)的選擇多基于經(jīng)驗，其變化對預(yù)測結(jié)果的影響程度不明，使得模型預(yù)測能力存在不確定性，難以指導(dǎo)勘探?jīng)Q策的風險評估。
技術(shù)實現(xiàn)思路
1、本發(fā)明的目的是解決現(xiàn)有技術(shù)中存在的缺點，而提出的基于多源地質(zhì)數(shù)據(jù)融合的卡林型金礦找礦預(yù)測方法。
2、為了實現(xiàn)上述目的，本發(fā)明采用了如下技術(shù)方案：基于多源地質(zhì)數(shù)據(jù)融合的卡林型金礦找礦預(yù)測方法，包括：
3、構(gòu)建針對導(dǎo)礦構(gòu)造與容礦構(gòu)造的多維度識別標志體系，所述多維度識別標志體系涵蓋地質(zhì)維度、地球物理維度和地球化學維度的識別要素；
4、基于所述多維度識別標志體系，對目標區(qū)域進行多源地質(zhì)數(shù)據(jù)的采集與整合，形成初始融合數(shù)據(jù)體；
5、對所述初始融合數(shù)據(jù)體實施構(gòu)造格架解析，提取與導(dǎo)礦構(gòu)造相關(guān)的深部構(gòu)造信息要素和與容礦構(gòu)造相關(guān)的淺部構(gòu)造信息要素，對所述深部構(gòu)造信息要素和所述淺部構(gòu)造信息要素進行耦合性分析，建立反映導(dǎo)礦構(gòu)造與容礦構(gòu)造空間配置關(guān)系的構(gòu)造耦合模型；
6、將地球化學異常信息與所述構(gòu)造耦合模型進行疊置分析，識別出與構(gòu)造耦合模型高吻合度的金礦化異常靶區(qū)，利用三維地質(zhì)建模技術(shù)對所述金礦化異常靶區(qū)進行立體表征，生成包含構(gòu)造控礦要素與礦化顯示要素的三維預(yù)測模型；
7、對所述三維預(yù)測模型進行不確定性量化評估，確定模型關(guān)鍵參數(shù)的敏感度與置信區(qū)間，基于所述不確定性量化評估結(jié)果，對所述三維預(yù)測模型進行迭代優(yōu)化，形成最終的找礦預(yù)測成果。
8、作為本發(fā)明的進一步方案，所述構(gòu)建針對導(dǎo)礦構(gòu)造與容礦構(gòu)造的多維度識別標志體系，包括：
9、系統(tǒng)收集區(qū)域地質(zhì)圖、地球物理勘探數(shù)據(jù)、地球化學采樣數(shù)據(jù)以及已知礦床地質(zhì)資料；
10、從所述區(qū)域地質(zhì)圖中解譯區(qū)域性滑脫構(gòu)造、擠壓背斜伴生斷裂作為導(dǎo)礦構(gòu)造的地質(zhì)維度識別標志；
11、從所述地球物理勘探數(shù)據(jù)中提取反映深部構(gòu)造特征的重力異常梯度帶、磁異常線性帶作為導(dǎo)礦構(gòu)造的地球物理維度識別標志；
12、從所述地球化學采樣數(shù)據(jù)中篩選出與深部流體活動相關(guān)的特定元素組合異常作為導(dǎo)礦構(gòu)造的地球化學維度識別標志；
13、從所述已知礦床地質(zhì)資料中歸納構(gòu)造蝕變體的空間展布特征、特定幾何參數(shù)的背斜軸部形態(tài)、穹隆構(gòu)造的分布范圍作為容礦構(gòu)造的地質(zhì)維度識別標志；
14、將從不同維度獲取的識別標志進行系統(tǒng)化歸類與編碼，形成結(jié)構(gòu)化的多維度識別標志體系數(shù)據(jù)庫。
15、作為本發(fā)明的進一步方案，所述對所述初始融合數(shù)據(jù)體實施構(gòu)造格架解析，包括：
16、應(yīng)用邊緣檢測與紋理分析算法處理地球物理數(shù)據(jù)，識別出潛在的線性構(gòu)造與邊界信息；
17、對地質(zhì)圖進行構(gòu)造要素數(shù)字化與空間拓撲關(guān)系構(gòu)建，明確不同級別斷裂的交接關(guān)系與序次；
18、將地球化學數(shù)據(jù)網(wǎng)格化并計算元素空間分布的各向異性，圈定元素富集中心與遷移軌跡；
19、綜合所述線性構(gòu)造與邊界信息、所述斷裂的交接關(guān)系與序次以及所述元素富集中心與遷移軌跡，構(gòu)建反映目標區(qū)域構(gòu)造格架的空間幾何模型與成因聯(lián)系模型。
20、作為本發(fā)明的進一步方案，所述提取與導(dǎo)礦構(gòu)造相關(guān)的深部構(gòu)造信息要素和與容礦構(gòu)造相關(guān)的淺部構(gòu)造信息要素，包括：
21、在所述構(gòu)造格架的空間幾何模型與成因聯(lián)系模型中，設(shè)定深度閾值以區(qū)分深部構(gòu)造域與淺部構(gòu)造域；
22、在深部構(gòu)造域內(nèi)，篩選具有區(qū)域延伸規(guī)模、切割深地殼層次的斷裂帶與韌性剪切帶，將其幾何形態(tài)、產(chǎn)狀參數(shù)、活動期次信息提取為深部構(gòu)造信息要素；
23、在淺部構(gòu)造域內(nèi)，識別與所述深部構(gòu)造信息要素具有成因聯(lián)系或空間交切的次級斷裂、裂隙密集帶、褶皺虛脫部位，將其空間形態(tài)、密度、組合關(guān)系信息提取為淺部構(gòu)造信息要素。
24、作為本發(fā)明的進一步方案，所述對所述深部構(gòu)造信息要素和所述淺部構(gòu)造信息要素進行耦合性分析，包括：
25、計算所述深部構(gòu)造信息要素的空間展布方向與所述淺部構(gòu)造信息要素的優(yōu)勢發(fā)育方向之間的夾角關(guān)系；
26、分析所述深部構(gòu)造信息要素中斷裂帶的多期活動歷史與所述淺部構(gòu)造信息要素中構(gòu)造的形成時序匹配關(guān)系；
27、評估所述深部構(gòu)造信息要素提供的成礦流體輸運潛力與所述淺部構(gòu)造信息要素提供的成礦空間存儲容量之間的配置效率；
28、基于所述夾角關(guān)系、所述形成時序匹配關(guān)系以及所述配置效率，建立定量描述導(dǎo)礦構(gòu)造與容礦構(gòu)造空間耦合強度的構(gòu)造耦合模型。
29、作為本發(fā)明的進一步方案，所述將地球化學異常信息與所述構(gòu)造耦合模型進行疊置分析，包括：
30、對所述地球化學采樣數(shù)據(jù)進行趨勢面分析與剩余異常計算，分離出與局部礦化作用相關(guān)的金元素異常及伴生元素異常；
31、將金元素異常及伴生元素異常的濃集中心空間坐標、異常強度、異常規(guī)模參數(shù)與所述構(gòu)造耦合模型的空間位置、耦合強度區(qū)域進行空間疊加；
32、計算每個疊加單元內(nèi)地球化學異常與構(gòu)造耦合強度的空間統(tǒng)計相關(guān)性指數(shù)；
33、依據(jù)預(yù)設(shè)的相關(guān)性指數(shù)閾值，圈定出與構(gòu)造耦合模型高吻合度的金礦化異常靶區(qū)，并記錄每個靶區(qū)的吻合度評分。
34、作為本發(fā)明的進一步方案，所述利用三維地質(zhì)建模技術(shù)對所述金礦化異常靶區(qū)進行立體表征，包括：
35、基于鉆孔數(shù)據(jù)、地球物理反演結(jié)果、地質(zhì)剖面圖，構(gòu)建目標靶區(qū)的三維地質(zhì)結(jié)構(gòu)框架模型；
36、將所述構(gòu)造耦合模型中的導(dǎo)礦構(gòu)造與容礦構(gòu)造空間關(guān)系集成到所述三維地質(zhì)結(jié)構(gòu)框架模型中，形成含構(gòu)造控礦要素的三維地質(zhì)模型；
37、將地球化學異常的三維空間分布數(shù)據(jù)、蝕變礦物組合分布數(shù)據(jù)作為屬性體賦予所述含構(gòu)造控礦要素的三維地質(zhì)模型，生成同時包含構(gòu)造控礦要素與礦化顯示要素的三維預(yù)測模型。
38、作為本發(fā)明的進一步方案，所述對所述三維預(yù)測模型進行不確定性量化評估，包括：
39、識別所述三維預(yù)測模型構(gòu)建過程中所依賴的數(shù)據(jù)來源，包括地質(zhì)解釋數(shù)據(jù)、地球物理反演數(shù)據(jù)、地球化學插值數(shù)據(jù)；
40、評估每種數(shù)據(jù)來源的數(shù)據(jù)精度、分辨率及覆蓋范圍，并分析數(shù)據(jù)之間的匹配一致性；
41、采用蒙特卡洛模擬方法，對所述三維預(yù)測模型的關(guān)鍵界面形態(tài)、構(gòu)造空間位置、礦化體邊界進行多次隨機擾動模擬；
42、統(tǒng)計多次模擬結(jié)果中模型關(guān)鍵參數(shù)的變化范圍，確定所述模型關(guān)鍵參數(shù)的敏感度與置信區(qū)間。
43、作為本發(fā)明的進一步方案，所述基于所述不確定性量化評估結(jié)果，對所述三維預(yù)測模型進行迭代優(yōu)化，包括：
44、根據(jù)所述模型關(guān)鍵參數(shù)的置信區(qū)間，調(diào)整三維地質(zhì)建模過程中地質(zhì)界面的控制點權(quán)重與插值算法參數(shù)；
45、針對敏感度高的模型關(guān)鍵參數(shù)，補充采集高精度數(shù)據(jù)或進行針對性地質(zhì)驗證，以約束模型優(yōu)化方向；
46、利用優(yōu)化后的參數(shù)重新運行三維地質(zhì)建模流程，生成更新后的三維預(yù)測模型；
47、比較更新前后的三維預(yù)測模型在已知礦點處的擬合程度差異，直至模型變化趨于穩(wěn)定且符合地質(zhì)規(guī)律，形成最終的找礦預(yù)測成果。
48、作為本發(fā)明的進一步方案，所述形成最終的找礦預(yù)測成果之后，還包括：
49、將所述最終的找礦預(yù)測成果與新一輪地質(zhì)勘查工程部署進行動態(tài)關(guān)聯(lián)；
50、根據(jù)勘查工程反饋的實見礦化信息，對所述多維度識別標志體系、所述構(gòu)造耦合模型以及三維建模參數(shù)進行校準；
51、利用校準后的參數(shù)和模型，對未經(jīng)驗證區(qū)域的預(yù)測結(jié)果進行滾動更新，實現(xiàn)找礦預(yù)測能力的持續(xù)提升。
52、與現(xiàn)有技術(shù)相比，本發(fā)明的優(yōu)點和積極效果在于：
53、通過實施構(gòu)造格架解析，專門提取深部導(dǎo)礦構(gòu)造信息要素與淺部容礦構(gòu)造信息要素，并對這兩類要素進行耦合性分析，建立了反映其空間配置關(guān)系的構(gòu)造耦合模型。這一技術(shù)路徑將構(gòu)造研究從形態(tài)學推進到功能系統(tǒng)學，明確了成礦流體從源區(qū)到沉淀場所的完整輸運網(wǎng)絡(luò)與控礦格架。基于此模型篩選的靶區(qū)，其地質(zhì)依據(jù)從“存在構(gòu)造”深化為“存在有效耦合的構(gòu)造系統(tǒng)”，提升了從復(fù)雜構(gòu)造背景中識別出真正具備成礦潛力空間的精準度。
54、對生成的三維預(yù)測模型進行不確定性量化評估，確定其關(guān)鍵參數(shù)的敏感度與置信區(qū)間，并依據(jù)該評估結(jié)果對模型進行迭代優(yōu)化。這一過程將靜態(tài)的地質(zhì)模型轉(zhuǎn)化為一個動態(tài)的、可評估的預(yù)測系統(tǒng)。敏感度分析揭示了哪些參數(shù)對預(yù)測結(jié)果影響最大，從而指導(dǎo)數(shù)據(jù)采集與解釋工作的優(yōu)化方向；置信區(qū)間則給出了預(yù)測結(jié)論的可信范圍。通過迭代反饋，模型能夠不斷融合新的認識與數(shù)據(jù)，其預(yù)測結(jié)果不再是單一結(jié)論，而是附帶可靠性度量與優(yōu)化路徑的決策支持信息，增強了預(yù)測成果在實際勘探應(yīng)用中的實用價值與風險可控性。