付海濤, 萬方來, 蔣麗麗, 徐 華, 李祎昕, 楊 磊
1)遼寧省地質勘探礦業集團, 遼寧沈陽 110032; 2)遼寧省第六地質大隊有限責任公司, 遼寧大連 116200;3)遼寧省物測勘查院有限責任公司, 遼寧沈陽 110121
遼寧省是我國重要的金伯利巖分布區, 遼南的瓦房店地區、遼東的桓仁地區、遼北的鐵嶺地區和遼西的葫蘆島地區都有金伯利巖出露(圖1a)。其中,遼南的瓦房店地區最為典型, 從 1971年在該區發現第一個金伯利巖體以來, 已發現超過 100個金伯利巖體, 既有巖管也有巖脈, 該區已提交 4個大型原生金伯利巖型金剛石礦床和3個近源小型金剛石砂礦床。金伯利巖既是金剛石礦的載體, 也是研究地球深部信息的重要對象, 其作為地球上來源最深的巖漿巖含有大量的深部物質信息, 因而具有重要的科學價值(路鳳香和桑隆, 2002)。長期以來, 專家學者對瓦房店地區的金伯利巖做了大量研究工作,認為該區金伯利巖體的巖漿來源深度約為 215 km,以18 m/s的上升速度沿著一定寬度的裂隙快速上侵而形成目前看到的金伯利巖(鄭建平等, 1989), 侵位時間約為 465 Ma(張宏福和楊岳衡, 2007), 金伯利巖受斷裂構造控制, 巖管產出于構造交匯部位(齊玉興等, 1998), 金伯利巖的分布與構造盆地有關(付海濤, 2020)。最新研究結果表明, 瓦房店地區存在逆沖推覆構造, 可能對金伯利巖體的形態造成了破壞(張國仁等, 2020)。
在金剛石礦勘查過程中對金伯利巖體地質特征進行了研究, 認為有必要對區內金伯利巖分布規律重新認識, 將原來的4個礦帶劃分為3個不同產出特征的區域可能更為合理(圖 1)。同時認為, 逆沖推覆構造是影響區內金伯利巖體形態特征的重要因素。
瓦房店金剛石礦田位于遼寧省南部, 行政區隸屬于大連市。大地構造位置位于華北陸塊遼東新元古代—古生代坳陷帶, 跨兩個四級構造單元, 東側為城子坦—莊河太古宙基底隆起, 西側為大連新元古代—古生代凹陷, 總面積約4000 km2。東側隆起區主要出露新太古界片麻巖; 西側凹陷區, 出露有古元古界、新元古界、古生界、中生界, 以及溝谷中的新生界。其中, 新元古界分布面積最大, 約占全區總面積的60%。
瓦房店礦田的金伯利巖可簡化分為三類五種巖性, 即塊狀金伯利巖、含巖球金伯利巖和角礫狀金伯利巖。其中塊狀金伯利巖分為兩種巖性, 一是斑狀金伯利巖, 呈暗綠色, 圓斑狀結構、塊狀構造,圓斑晶為蛇紋石化、滑石化、碳酸鹽化橄欖石(假像)和橢圓餅狀金云母。斑晶含量一般 15%~20%, 基質為橄欖石(假像)、鉻鐵礦、銳鈦礦及隱晶或微晶質方解石等組成; 二是斑狀金云母金伯利巖, 呈灰綠色, 地表風化后呈灰-黃褐色, 斑狀結構、塊狀構造。金云母含量大于15%。斑晶為蛇紋石化、碳酸鹽化橄欖石(假像), 一般直徑 3~5 mm, 最大可達10 mm以上, 含量一般小于30%。金云母呈淡黃色,風化后呈銀白色薄片狀, 片徑3~5 mm。含巖球斑狀金伯利巖, 巖石呈灰綠-黑綠色, 風化后表面呈黃綠色、黃褐色, 具有特殊的球狀構造。風化強烈的巖石巖球脫落, 在其殘坡積物中可見大小不等的圓形球蛋。此種金伯利巖在1號、30號和42號等巖管均有見到。其中 42號巖管, 巖球含量平均達25.6%, 球徑一般為 10~30 mm, 最大為 150 mm,巖球由核心及外殼組成, 核心一般為渾圓狀橄欖石(假像), 少數核心為蝕變片麻巖, 外殼由細粒斑狀金伯利巖組成。這類金伯利巖除巖球外, 其它則是斑狀金云母金伯利巖或斑狀金伯利巖。角礫狀金伯利巖也分為兩種巖性, 一是含圍巖角礫金伯利巖,呈灰綠、暗綠色, 發育強碳酸鹽化時呈灰白色, 地表風化后呈黃褐色。碎屑結構、角礫狀構造。角礫為圍巖角礫和早期的金伯利凝灰角礫巖, 角礫大小由幾 mm, 到幾百 mm 不等, 圍巖角礫邊部常見有滑石化和綠泥石化, 膠結物為斑狀(金云母)金伯利巖; 二是含金伯利巖物質角礫巖, 該巖性中角礫為不規則棱角狀圍巖角礫, 角礫含量 95%以上, 角礫大小不一, 大的可達2 m, 膠結物主要由小于2 mm的各種圍巖細碎屑及金伯利巖組成。
巖管以角礫狀金伯利巖為主, 其中最著名的50號巖管角礫狀金伯利巖約占 80%, 塊狀金伯利巖約占20%。而脈狀金伯利巖則基本上都是塊狀金伯利巖。區內金伯利巖的蝕變作用主要為蛇紋石化和碳酸鹽化, 地表褐鐵礦化、綠泥石化比較發育,深部逐漸變弱。強烈碳酸鹽化是礦區金伯利巖顯著的、普遍的蝕變作用。
已發現的金伯利巖體主要集中在礦田中部長約33 km、寬約15 km的范圍內(稱為礦田核心區),這個區內分布有113個金伯利巖體, 其中24個為巖管。礦田北部, 在礦田核心區西北側27 km處發現7個疑似金伯利巖體。這些金伯利巖體成群成帶分布, 局部地段巖體成帶性非常明顯, 在小比例尺圖面上可以清晰看到金伯利巖體呈線狀或孤狀展布(圖1b)。以往在礦田核心區將金伯利巖劃分為3個礦帶, 將西北側出露有可疑巖體的地段劃分為第Ⅳ礦帶。本次研究, 根據金伯利巖體相鄰關系將金伯利巖劃分為1個巖群區、1個巖帶區和1個單巖體區(圖 1)。

圖1 遼寧金伯利巖區域分布(a)和瓦房店金剛石礦田地質簡圖(b)(根據遼寧省地質勘查院, 2017改編)Fig.1 Regional distribution of kimberlites in Liaoning (a) and a simplified geological map of the Wafangdian diamond ore field (b) (adapted from Liaoning Provincial Institute of Geological Exploration, 2017)
①巖群區金伯利巖成群分布, 僅有少量單獨的金伯利巖體, 位于原Ⅰ、Ⅱ礦帶的西南端, 分布有37個金伯利巖體, 其中 10個是巖管, 巖群區有5個巖群和3個單獨的金伯利巖體, 1、2號巖管以及50號和51-68-74號巖管兩個大型原生金剛石礦床分布在此區。
②巖帶區金伯利巖成帶分布, 范圍包括原Ⅰ礦帶中部、Ⅱ礦帶區東部和Ⅲ礦帶, 區內有 58個金伯利巖體, 其中8個是巖管, 包括30號巖管大型原生金剛石礦床和34、38、110、111號等巖管, 巖帶區內各巖體之間相距很近, 平面上構成直線狀或孤線狀4條巖帶。
③單巖體區, 雖然局部也有相距較近的金伯利巖體, 但不具備巖群或巖帶的特征, 這一區域包括原Ⅰ礦帶東端和Ⅳ礦帶, 區內有 18個金伯利巖體和7個疑似金伯利巖體, 其中6個是巖管, 42號大型原生金剛石礦床出露在此區域內。
本區金伯利巖體受斷裂構造控制特征明顯, 按產出形態分為巖脈和巖管, 巖脈形態簡單, 產狀比較穩定, 巖管形態較為復雜。
巖脈一般呈北東東至近東西向展布, 走向65°~90°, 方向性明顯, 嚴格受密集節理帶控制,相鄰脈體相互平行, 走向一般比較穩定, 傾向南東,傾角 70°~85°之間, 長度幾十米到千余米不等, 脈寬0.2~1.0 m, 個別地段可達3~5 m。一般來說, 脈體形態簡單規整, 與圍巖界線清楚, 脈幅變化較小,僅在局部有膨大或變窄現象, 單個脈沿走向具有尖滅再現分枝復合現象, 局部可見幾條相互平行的脈組成一個小脈帶。脈狀金伯利巖巖性比較簡單主要為強碳酸鹽化斑狀金云母金伯利巖, 一般均含有金剛石。
巖管的形態比較復雜, 地表出露形態有橢圓狀、棱形狀、葫蘆狀和不規則狀等(圖2)有的巖管具有明顯的長軸方向, 有的巖管則呈近等軸狀, 能看出長軸方向的基本上其長軸方向均為近東西向或北東向, 規模較大的以近東西向為主, 如30號、50號巖管。地表出露面積最大的是 42號巖管, 面積約41 200 m2, 小巖管的面積甚至小于100 m2。巖管多數向南傾斜, 傾角 75°~85°。

圖2 典型金伯利巖管形態示意圖Fig.2 Schematic diagram of typical kimberlite pipe morphology
自從發現金伯利巖以來, 本區的金剛石勘查工作基本上沒有中斷, 對 4個大型原生礦床施工了大量的鉆探工作量, 對巖體形態的控制也比較好。
50號巖管是區內最著名的巖管, 位于巖群區,已開發利用。據勘探報告記載, 該巖管出露標高130~200 m, 主礦體部分地表呈不規則菱形, 東西較長、南北較窄, 長寬比為 4.8:1, 東部為脈狀。巖管縱向變化較大, 在+60 ~ –20 m標高部位膨大, 地表與+60 m標高巖管面積的比值為 1:1.9, 地表與–20 m標高巖管面積的比值為1:2.3, 巖管傾向南南西, 傾角70°~78°。在–70 m標高左右急劇收縮尖滅。圖3示意了50號巖管不同標高的形態。根據礦區鉆孔資料建立的三維礦體模型顯示, 50號巖管主礦體部分呈靴子狀, 底部西端翹起似鞋跟、東端有拖尾似鞋尖, 底界平直, 總體上為一個厚大板狀體, 向南傾斜, 傾角約85°。
30號巖管是目前區內工程控制最深的巖管, 位于巖帶區。地表形態為長橢圓形, 長軸方向 70°, 長212 m, 平均寬80 m。詳查報告認為, 巖管分為30-1和 30-2兩個礦體。30-1出露于地表, 傾向南東, 傾角50°左右, 延深120 m, 在其南東側存在隱伏礦體,編號為30-2。圖4兩條勘探線剖面示意了30號巖管的深部形態。根據礦區鉆孔資料, 認為30號巖管是近直立的金伯利巖管, 在被近水平斷層錯斷后, 鉆孔揭露到的礦體向南東方向逐漸加深從而顯示側伏的特征(付海濤等, 2017)。總體來看, 30號巖管還是一個相對完整的巖管, 只是局部有些錯動。30號巖管是區內面積第二大的巖管, 邊部出露有含金伯利巖物質的角礫巖, 鉆孔中在巖管南側見有該類型角礫巖, 但與地表見到的角礫巖不具有上下對應關系。

圖4 30號巖管深部形態示意圖Fig.4 Schematic diagram of the deep morphology of No.30 kimberlite pipe
42號巖管位于單巖體區, 地表南北兩側出露有含金伯利巖物質角礫巖, 中間為金伯利巖管(圖 5),呈筒狀, 不同標高形態變化不大, 是區內唯一一個上下近等粗的筒狀巖管, 也是目前我國面積最大的金伯利巖體。地表出露標高約 150 m, 鉆探控制巖體延深超過560 m, 由于是區內4個大型原生金剛石礦床中品位最低的, 所以工作程度相對較低, 有部分鉆孔終孔在巖體之中, 也就是說這個巖管向下還有延深。

圖5 42號巖管不同標高水平投影斷面示意圖Fig.5 Schematic diagram of the morphology of No.42 kimberlite pipe at different elevation levels
相對來說, 本區金伯利巖管的工程控制程度要好于巖脈, 尤其是 4個大型原生金剛石礦床的控制程度更高。位于50號巖管北側1 km處由51號、68號和74號三個巖管組成的大型原生金剛石礦床,雖然地表出露面積相差不大, 但鉆孔控制的范圍內,51號巖管明顯比另外兩個巖管大很多, 至少向下延深了近 400 m, 不同標高處巖管的水平斷面面積略大于地表出露面積; 74號巖管是向下延深較小的一個, 僅在少數幾個鉆孔中見到; 68號巖管不同標高水平斷面的面積變化也不大。根據鉆孔資料, 西端的68號巖管在–200 m左右、東端的51號在–240 m標高處巖管突然中斷, 這兩點的連線角度與 50號巖管底部界面的連線角度基本相同。
本區金伯利巖脈一般控制程度較低, 鉆孔中見到的脈體與地表見到的脈體不論是厚度還是巖性變化都不大, 但未發現向下延深較大的脈體。
前人研究成果認為, 遼寧金伯利巖的化學組成與世界金伯利巖平均成分相比, 堿度偏低; 與南非、前蘇聯和我國山東金伯利巖成分相比, Fe質偏高; 與上地幔巖石相比, 鐵族元素降低, 稀有及稀土元素、放射性元素、鎢鉬族元素、親銅元素等相對富集(李星云等, 1989)。
有學者認為, 用混染指數C.I——(SiO2+Al2O3+Na2O)/(MgO+2K2O), 可以較好地反映金伯利巖中殼源巖石的混染程度, 受混染的金伯利巖C.I>1.5, 我國部分金伯利巖受到了混染作用的影響(董振信, 1991)。本次研究對不同分區典型巖管的樣品分析結果進行統計, 結果顯示, 本區金伯利巖的混染指數具有一定的規律性, 巖群區的 1號、50號和51號巖管的C.I<1.5, 巖帶區和單巖體區的金伯利巖的C.I≥1.5(見表1)。說明瓦房店礦田巖群區的金伯利巖受混染的影響可能較小, 更適合于開展深部物質信息的研究。同時, 對典型金伯利巖體與世界金伯利巖成分進行了對比(見表 2), 結果顯示, 瓦房店礦田的金伯利巖總體上富 Ca、Al, 特別是位于巖帶區的巖體 Ca含量更高, 甚至高出世界平均值的3~4倍; 貧Mg、Na, 巖帶區的巖體甚至只有世界平均值的 30%左右; 本次研究的巖體中除51號巖管外, 其它巖體呈現貧K特征; 位于單巖體區的 42號巖管, 貧 Mg、K、Na、P, 富 Si、Al、Ca、Ti、Fe, 是區內唯一個富Si、Fe的巖體。

表1 瓦房店礦田典型金伯利巖體混染指數統計表Table 1 Statistical table of contamination index of typical kimberlites in Wafangdian ore field

表2 瓦房店礦田典型金伯利巖體元素含量與世界平均值的比值Table 2 Ratio of element content of typical kimberlites in Wafangdian ore field compared to the world average
瓦房店礦田金伯利巖最重要的地球物理特征是其具有明顯而不均勻的磁性, 表 3記錄了礦田典型金伯利巖的磁化率和剩余磁化強度。從表3上可以看到, 各巖體的磁性變化較大, 有的巖體磁化率變化可達百倍以上, 有的巖體磁化率則變化較小。50號巖管的剩磁明顯高于其它巖體。碳酸鹽化金伯利巖(如11號脈)的磁性明顯降低, 甚至達到無磁性狀態。

表3 瓦房店礦田典型金伯利巖磁性參數統計表Table 3 Statistical table of magnetic parameters of typical kimberlites in Wafangdian ore field
根據航磁和地磁測量結果, 金伯利巖體具有磁異常。面積較大的30號和42號巖管, 1:5萬航磁測量有明顯的異常顯示, 異常強度超過100 nT, 而面積相對較小的50號巖管只有在降低測線高度至幾十米時才能觀測到很弱的異常(安戰鋒等,2017)。地磁測量是本區發現金伯利巖的重要手段,1:2000地磁ΔZ測量結果, 50號巖管, 50 nT等值線大致是巖管的出露邊界; 30號巖管, 100 nT等值線揭示了隱伏礦體的范圍。42號巖管, 1:5000地磁ΔZ測量結果 200 nT等值線大致是巖管的出露邊界。51號巖管是通過 1:2000地磁測量而發現的,而68號和74號巖管在1:2000地磁測量結果中卻沒有異常。不同的巖體磁測比例尺不同, 有些小巖管需要用 1:1000地磁測量才能用磁異常圈出巖體的范圍。
總體上看, 本內地層基本上屬無磁性地質體,金伯利巖屬中等磁性地質體, 與區內輝綠巖、橄欖玄武巖等磁性相當。
表4是30號和50號兩個巖管的電阻率統計結果。總體上, 金伯利巖屬中低電阻率地質體, 50號巖管電阻率略高于30號巖管。當圍巖是高阻的泥灰巖、片麻巖時, 金伯利巖顯示為低阻體; 當圍巖是粉砂巖、頁巖時, 金伯利巖與圍巖之間幾乎無電阻率差異, 很難在這些巖性中判斷出金伯利巖的位置。

表4 瓦房店礦田30號、50號金伯利巖管電阻率統計表Table 4 Resistivity of No.30 and No.50 kimberlite pipes in Wafangdian ore field
巨(大)型構造是控制金伯利巖產區的重要因素(李開善, 1981); 金伯利巖漿沿裂隙快速上侵(鄭建平等, 1989); 就位于寒武系毛莊組以下層位的地層之中, 礦田見有金伯利巖被同位素年齡為 95~144 Ma的橄欖玄武巖所穿切(趙春強等, 2018); 測年結果顯示, 30號巖管432 Ma, 42號巖管341 Ma、356 Ma、455 Ma, 50號巖管422 Ma、462 Ma, 51號巖管 463 Ma(董振信, 1994), 最老的年齡與張宏福和楊岳衡(2007)的結果相近。根據南京大學(2016)研究結果, 本區中生代存在兩期(165~166 Ma和138 Ma)對金剛石礦床起破壞作用的巖漿事件(樣品取自30號巖管鉆孔中和42號巖管旁的閃長玢巖)。在礦田東部的早侏羅統砂礫巖中發現有金剛石砂礦。雖然山東蒙陰坡里地區金伯利巖帶形成于中生代燕山晚期(褚志遠等, 2019), 但目前遼寧還沒有相似的研究成果見報。因此認為, 本區金伯利巖侵位于中寒武世之后侏羅紀之前。
本區金伯利巖體, 除42號巖管以外, 其它被稱為巖管的金伯利巖體絕大部分仍具有脈的特征, 具有明顯的拉長形態(以往將長寬比小于10的稱為巖管)。以規模較大的50號和30號巖管為例, 地表形態以及鉆孔控制的深部形態, 這兩個巖體均具有厚大板狀特征。平面上近等軸狀的有 1號、34號、38號和75號等小巖體(圖2), 但從鉆孔控制結果看,這些巖管深部有的與脈體相連, 比如38號與111號巖管通過巖墻連成一體。區內脈狀金伯利巖體的邊界與圍巖清晰, 侵位時受斷裂構造控制的特征明顯。厚板狀的巖管可能侵位于斷裂構造的復合部位,或者是不同方向斷裂構造的交叉部位。由于脈狀體未見角礫狀金伯利巖或含金伯利物質角礫巖, 說明脈狀體是沿壓(扭)性斷裂侵入的。
根據已有資料, 本區金伯利巖體在侵位后基本上未遭受太多的改造, 對巖體形態影響最大的是逆沖推覆構造。在進行三維建模過程中發現30號巖有被近水平斷層錯斷的現象, 經過野外工作發現礦田存在推覆構造體系, 初步判斷該推覆構造形成于中晚侏羅世, 晚于金伯利巖侵位時間,對區內金伯利造成了不同程度的破壞, 從而很好地解釋了區內有些金伯利巖體向下延深不大或突然尖滅的現象。這套推覆構造體系, 從南東東向北西西方向運動, 在 50號巖管附近可見該推覆構造線理產狀顯示構造面為傾角 15°左右向南東傾、由110°向 290°方向運動。位于礦田核心區東北部的42號巖管在鉆孔控制的深度內未見明顯錯斷的現象; 30號巖管位于中部偏北的位置, 雖然被錯斷了但巖管基本上還是連續的, 鉆孔資料和音頻大地電磁測量結果未見明顯的地層重疊現象; 南部的 50號巖管則被徹底錯斷, 礦區地表可見不同巖性層多次重復, 鉆孔巖心可見地層重疊現象,音頻大地電磁測量結果有重復出現的高低阻異常。重力剖面測量反演結果, 30號巖管所在位置的地層受推覆構造影響而產生的折曲明顯弱于50號巖管所在的位置, 30號巖管受破壞的程度也弱于50號巖管。礦田南部金伯利巖體受逆沖推覆構造的影響明顯大于北部。在逆沖推覆構造鋒帶局部可見規模不大的平臥褶皺導致地層發生倒轉, 目前尚未發現金伯利巖被錯斷后發生倒轉現象。
總體來看, 本區深大斷裂提供了金伯利巖巖漿的上侵通道, 淺部壓(扭)性斷裂提供了就位空間。金伯利巖侵位后受到推覆構造的影響, 巖體形態發生了變化, 顯現出西淺東深向東側伏或向下突然尖滅的形態特征。也就是說, 本區影響金伯利巖體形態的主要因素是斷裂構造, 巖漿侵位時斷裂構造控制了巖體形態, 成巖后斷裂構造對巖體形態進行了改造, 特別是逆沖推覆構造使部分金伯利巖體變成多個斷塊, 在開展深部找礦時需要對巖體的不連續性加以注意。
有研究認為, 瓦房店地區與郯廬斷裂活動相關的北北東向斷裂構造控制了金伯利巖的空間分布范圍(趙建軍等, 2011); 郯廬斷裂的雛形出現在中生代的中央造山帶東段(梁光河, 2018), 斷裂開始形成于中、晚三疊世(萬天豐等, 1996), 郯廬斷裂北段左行走滑斷裂形成于100 Ma左右(竇立榮等, 1996); 有研究認為, 瓦房店地區的金伯利巖與晚侏羅世—早白堊世的大火山巖省的有關(曾普勝等, 2021); 瓦房店地區的金伯利巖可能在早古生代之前形成于華北地臺之巖石圈底部, 而在中生代白堊紀時期才在大規模巖石圈拆沉的地質背景下快速上升(王十安等,2021)。但目前尚無支持這些認識的測年數據。也有研究認為, 郯廬斷裂帶是一條長壽斷裂, 具有長期復雜的活動歷史。在《郯廬斷裂》(國家地震局地質研究所, 1987)一書中記載: 李四光認為該斷裂控制了遼東與遼西古生代沉積盆地; 徐嘉煒認為該斷裂東支可能在太古代末產生, 長期繼承活動; 黃汲清認為形成于前震旦紀, 中、新生代重復活動; 馬杏垣認為元古代以前控制東北元古代褶皺帶的構造發展, 主要活動期在中生代; 于丕休等認為太古代以來長期活動。另外, 有研究認為, 古郯廬斷裂是震旦紀開始形成的(喬秀夫和張安棣, 2002)。本文認為,區內金伯利巖受古郯廬斷裂帶的次級斷裂控制, 郯廬斷裂在本區至少在古生代就已經存在。
本區金伯利巖成群、成帶分布, 南西部以巖群為主, 中部以巖帶為主, 東北部則以單巖體分布為主。根據巖體的化學成分, 巖群區的金伯利巖體受圍巖成分混染的程度小于巖帶區或單巖體區的金伯利巖。4個大型原生礦床中, 位于巖群區的50號和51號巖管混染指數分別為1.36和1.26, 小于1.5的判別標準; 巖帶區的 30號巖管混染指數為 1.5, 是巖帶區混染指數最小的巖體; 單巖體區的 42號巖管混染指數為1.61, 高于1.5的判別標準, 因此認為巖群區的巖體更適合用于研究深部信息。
本區的金伯利巖與其它地區的金伯利巖相比,地球化學特征和地球物理特征既有相似性也有獨特性。影響巖體形態的主要因素是斷裂構造, 斷裂構造提供了巖漿上侵通道和就位空間, 成巖后的逆沖推覆構造對巖體形態造成了不同程度的影響。
致謝: 本文為遼寧省地質勘查項目“瓦房店金剛石整裝勘查區增儲專項”(編號: LNZC2019-0078-20)資助的成果。感謝項目組全體成員和遼寧金剛石勘查技術攻關團隊各位專家、同事付出的辛勤勞動和大力支持, 感謝遼寧省第六地質大隊有限責任公司總經理韓雙、總工程師許洪斌給予的大力支持, 感謝審稿專家的悉心指導。
Acknowledgements:
This study was supported by Liaoning Province Geological Exploration Project (No.LNZC2019-0078-20).