摘要：棲霞笏山西陡崖金礦位于蓬萊棲霞金成礦帶南側(cè)，是產(chǎn)于新太古代變質(zhì)基底與中生代燕山期花崗巖接觸帶附近的破碎蝕變巖型金礦。本次研究在礦區(qū)外圍開展了詳細的野外地質(zhì)調(diào)查，并對該礦床進行了詳細研究，在Ⅰ號主礦體深部施工DKZK03鉆孔，控制金礦體4層，累計真厚度6.85m。其中，Ⅰ1、Ⅰ2號金礦體埋深分別為1273.80～1275.45m、1072.95～1076.55m，Ⅰ1號金礦體為該區(qū)目前最深金礦體。采用趨勢外推法對區(qū)內(nèi)深部金礦體進行定位預測，為今后在棲霞地區(qū)金礦深部找礦勘查提供了依據(jù)，該地區(qū)金礦深部具有較大的找礦潛力。

關(guān)鍵詞：金礦；深部鉆探；成礦預測；笏山西陡崖；山東棲霞

中圖分類號：P618.51""" 文獻標識碼：A""" doi：10.12128/j.issn.16726979.2024.03.004

引文格式：劉強，張巖，蔡娜，等.山東棲霞笏山西陡崖金礦地質(zhì)特征及資源潛力分析［J］.山東國土資源，2024，40（3）：3341. LIU Qiang， ZHANG Yan， CAI Na， et al. Analysis on Geological Characteristics and Resource Potential of Hushan-Xidouya Gold Deposit in Qixia City in Shandong Province［J］.Shandong Land and Resources，2024，40（3）：3341.

0 引言

棲霞金礦位于棲霞蓬萊成礦帶上［13］，是膠東地區(qū)重要的金及多金屬礦成礦區(qū)［47］，近年來，大量的中外地質(zhì)、礦產(chǎn)專家來區(qū)內(nèi)進行調(diào)查研究，取得了豐富的地質(zhì)資料［844］。隨著山東省開展找礦突破戰(zhàn)略行動，在棲霞金礦田探明了笏山西陡崖大型蝕變巖型金礦床，累計查明金資源量30余噸，伴生銀125 t［31］，為棲霞金礦田的首個大型金礦床。該礦床的發(fā)現(xiàn)，改變了傳統(tǒng)的棲霞蓬萊地區(qū)沒有大型金礦的看法，開拓了在該地區(qū)主要斷裂帶深部找礦的新方向。棲霞地區(qū)以往勘查開采深度不夠，對深部控礦條件、成礦要素、成礦規(guī)律等的研究相對缺乏［45］，本文以棲霞笏山西陡崖金礦為研究對象，在總結(jié)礦床地質(zhì)特征的基礎上，進行深部成礦預測并進行鉆探驗證，為在該地區(qū)開展深部金礦找礦工作提供依據(jù)。

1 區(qū)域地質(zhì)背景

棲霞地區(qū)金礦位于膠東半島的中北部，西為臺前陡崖斷裂、東臨棲霞斷裂（圖1），大地構(gòu)造位置隸屬華北板塊的膠遼隆起區(qū)（ⅡⅠ）膠北隆起（ⅡⅠa）及膠萊盆地北東部膠北斷隆（Ⅲa1）、回里養(yǎng)馬島斷隆（Ⅲa2）及萊陽斷陷（Ⅲb4）。區(qū)內(nèi)有臺前金礦、西陡崖金礦、笏山金礦、馬家窯金礦、山城金礦等多處金礦，是棲霞地區(qū)重要的金成礦區(qū)，成礦條件優(yōu)越［2333，4647］。

研究區(qū)出露地層主要有古元古代荊山群和粉子山群、新元古代蓬萊群及新生代臨朐群和第四系。區(qū)內(nèi)構(gòu)造以NNE向和NEE向壓扭性斷裂為主，另發(fā)育少量NNW、近SN向斷裂，韌性剪切帶多發(fā)育在臺前陡崖斷裂中。研究區(qū)巖漿侵入活動頻繁，以新太古代TTG、中生代燕山期侵入巖最為發(fā)育。其中，前寒武紀侵入巖均遭受不同程度的變質(zhì)變形，經(jīng)區(qū)域變質(zhì)作用形成一套片麻巖類；中生代侵入巖受近EW、NE和NNE向斷裂制約，總體呈近EW或NE向展布的巖基、巖株、巖瘤狀產(chǎn)出。

2 礦區(qū)特征

2.1 礦床地質(zhì)特征

礦區(qū)地層發(fā)育古元古代粉子山群祝家夼組與新生代第四系；區(qū)內(nèi)主要斷裂構(gòu)造為NNE及NE向兩組，構(gòu)造主體為臺前陡崖斷裂帶，另有少量的次級構(gòu)造；巖漿巖主要出露新太古代譚格莊序列、棲霞序列和馬連莊序列，以及中生代玲瓏序列大莊子單元，另外還有一些較小的脈巖（圖2），脈巖走向NNE，大致與構(gòu)造線分布一致。

2.2 礦體地質(zhì)特征

棲霞笏山西陡崖金礦區(qū)分為笏山礦區(qū)和西陡崖礦區(qū)。笏山礦區(qū)共有礦體21個，均為隱伏礦體，其中Ⅰ號礦體為主礦體（圖3），其他20個礦體均為小規(guī)模礦體；Ⅰ號礦體分布于532～523線之間，主要由ZK53234等26個鉆孔控制，礦體呈脈狀、透鏡狀，走向20°～60°，傾向SE，傾角20°～45°。區(qū)內(nèi)已控制走向長1400m（北部出礦區(qū)），傾向延伸平均524m，控制標高+61m～690m，礦體厚0.71～18.14m之間，平均厚度3.59m，厚度變化系數(shù)96%，屬較穩(wěn)定型。礦體品位在1.15～20.48g/t，平均品位2.68g/t。

西陡崖礦區(qū)發(fā)育礦體7個，其中Ⅰ1礦體為主礦體。Ⅰ號礦體分布于71～47 線之間，由ZK06等12個鉆孔控制，礦體呈脈狀、透鏡狀，走向23°～29°，傾向SE，傾角27°～30°。區(qū)內(nèi)已控制走向長約670m，傾向延伸約460m，控制標高+40m～223m，礦體厚0.87～25.14m 之間，平均厚度5.12m，厚度變化系數(shù)115%，屬較穩(wěn)定型，礦體品位在1.12～13.31g/t，平均品位3.28g/t，單樣最高品位13.31g/t。

兩礦區(qū)的礦體均賦存于與臺前陡崖斷裂下盤平行的含礦蝕變帶中。

2.3 礦石特征

礦體埋藏深度較深，礦石主要為破碎蝕變巖型含金黃鐵絹英巖，金屬硫化物主要為黃鐵礦，其次為黃銅礦。礦石以黃鐵絹英巖化碎裂巖為主，金屬礦物主要為黃鐵礦，其次為磁黃鐵礦、黃銅礦（圖4a、圖4b），金礦物主要為自然金及銀金礦（圖4c、圖4d）。非金屬礦物為絹云母、石英、長石等。與成礦密切相關(guān)的黃鐵礦呈細脈、浸染狀分布，受構(gòu)造控制明顯。細脈浸染狀構(gòu)造表現(xiàn)為礦石中金屬礦物呈星點狀均勻散布或局部不規(guī)則相對集中分布于石英顆粒內(nèi)，多期活動形成浸染狀構(gòu)造；同時也發(fā)育有呈細脈狀沿裂隙填充分布，構(gòu)成細脈浸染狀構(gòu)造。

含礦熱液沿前期的構(gòu)造裂隙對構(gòu)造巖和圍巖產(chǎn)生蝕變交代作用，形成了該區(qū)較典型的中低溫蝕變，即鉀化、青磐巖化、硅化、絹云母化、黃鐵礦化、碳酸鹽化。整個蝕變受構(gòu)造控制，具多期性，硅化、黃鐵礦化越強則金品位越高。

3 深部驗證

根據(jù)前期詳細地質(zhì)調(diào)查，本研究選擇在ZK53234孔以東524勘探線布設DKZK03鉆孔，終孔深度1500.85m，達到施工地質(zhì)目的，為深部找礦和成礦預測提供依據(jù)（圖5）。

該鉆孔上部主要為新太古代棲霞序列回龍夼單元細粒英云閃長質(zhì)片麻巖，頂部有元古代沉積變質(zhì)地層覆蓋，下部以中細粒黃鐵絹英巖化二長花崗巖為主。根據(jù)鉆孔巖性組合特征，將該孔劃分為17個大的巖性層（圖6），控制礦體4個。金礦化層主要發(fā)育在黃鐵絹英巖化碎裂巖和黃鐵絹英巖化二長花崗巖中，其中在1273.80～1275.45m新發(fā)現(xiàn)真厚度為1.53m的金礦體（Ⅰ1號礦體），品位1.09g/t；在1072.95～1076.55m處控制Ⅰ2號金礦體真厚度達3.33m，品位1.42～4.56g/t，平均品位2.90g/t。另外，在842.95～843.95m、922.30～923.35m分別見真厚度1.02m和0.97m的金礦體，品位分別為1.29g/t和1.24g/t。

DKZK03是該區(qū)目前見金礦體最深的鉆孔，較好地控制了笏山西陡崖金礦區(qū)Ⅰ號礦帶的延伸，達到地質(zhì)目的。

4 資源潛力分析

在充分收集前期礦床勘查資料的基礎上，利用礦區(qū)最新的礦床勘查資料，對勘探線剖面深部礦體的變化特征進行詳細分析研究，采用趨勢外推法對區(qū)內(nèi)深部金礦體進行定位預測，為今后棲霞地區(qū)金礦深部找礦勘查提供依據(jù)。

研究區(qū)內(nèi)圈定了1個成礦預測區(qū)：棲霞臺前艾前夼金礦成礦預測區(qū)（YCA1）。該成礦預測區(qū)位于研究區(qū)西南部的棲霞臺前—笏山—西陡崖—艾前夼一帶，呈NNE向展布，面積約175.39km2。區(qū)內(nèi)主要出露古元古代粉子山群祝家夼組變質(zhì)巖系；新太古代棲霞序列TTG、晚侏羅世玲瓏序列二長花崗巖，中酸性脈巖較發(fā)育；區(qū)內(nèi)NNE向的臺前陡崖斷裂、栆林斷裂及紫現(xiàn)頭斷裂極其發(fā)育，構(gòu)成臺前艾山金礦成礦帶。近EW向基底構(gòu)造與NNE向斷裂構(gòu)造交會部位、不同巖體接觸帶是成礦有利部位，斷裂帶下盤發(fā)育的大量花崗偉晶巖脈和巖石發(fā)生強烈的面型鉀化是礦床找礦的間接標志（表1）。

臺前陡崖斷裂帶總體上呈NE—NNE向展布，斷裂帶中多發(fā)育含金礦化蝕變，或被石英脈及煌斑巖、閃長玢巖脈充填。與該斷裂帶大致平行的斷裂帶在區(qū)域上自西向東有三山島斷裂、焦家斷裂、招平斷裂等，斷裂之間有近等間距分布特征。斷裂面往往沿各種巖體接觸帶或附近分布，顯示多期活動的特征。早期以韌性變形為主，形成灰綠色為主的糜棱巖和碎裂巖；中期為張性（張扭性）活動，構(gòu)造中常見大小不等棱角清楚的角礫，其成分主要是早期的糜棱巖及部分圍巖。晚期則是明顯的壓扭性特征，往往形成大量灰白色的糜棱巖和碎裂巖，局部有早期產(chǎn)物構(gòu)造透鏡體，斷裂帶下盤一般有大量花崗偉晶巖產(chǎn)出。斷裂在深部地段見浸染狀和細脈狀的黃鐵礦，甚至有多期次的硫化物出現(xiàn)，并出現(xiàn)低品位金礦（化）體。

該區(qū)1∶20萬區(qū)域重砂測量圈定自然金重砂異常5處，黃鐵礦重砂異常2處，銅重砂異常2處，沿臺前陡崖斷裂呈NNE向帶狀分布。而自然金、黃鐵礦、銅重砂異常主要集中在臺前—笏山—西陡崖—回兵崖一帶。

4.1 棲霞小廟前成礦潛力區(qū)

對該成礦潛力區(qū)深部的Ⅰ1礦體、Ⅰ2礦體分別進行預測，其中，Ⅰ1礦體由于向南西方向沒有被鉆孔封閉，已有多個鉆孔見礦。根據(jù)該類型金礦具有礦化連續(xù)、穩(wěn)定分布的特點，在預測礦體邊界，向南西再推150m。預測范圍在65～556線之間。預測含礦地質(zhì)體鉛垂厚度3.12m，傾角18°，礦床勘查深度約2000m，預測標高550m～1750m。預測礦體投影寬度約2000m，長度約2800m，投影面積為4195031m2，體積為13088497m3。體積含礦率為0.00793，預測含礦地質(zhì)體的總體變化趨勢由上向下厚度減小，品位變化不大，礦體厚度變化較大（圖7），預測金金屬量51.896 t。

Ⅰ2礦體預測范圍在47～540 線之間。預測含礦地質(zhì)體鉛垂厚度5.83m，傾角18°，礦床勘查深度約2000m，預測標高510m～1500m。預測礦體投影寬度約1600m，長度約2660m，投影面積為2568476m2，預測含礦地質(zhì)體鉛垂厚度5.83m，傾角18°，體積為14974215m3，預測金金屬量59.373 t（圖8）。

4.2 棲霞上峴成礦潛力區(qū)

由于預測礦體向南東方向礦體沒有封閉，根據(jù)該類型金礦具有礦化連續(xù)、穩(wěn)定分布的特點，預測范圍在17～0線之間（圖9），預測含礦地質(zhì)體寬度655m，長度約833m，投影面積497846m2。預測含礦地質(zhì)體鉛垂厚度2.36m，傾角28°，體積1174917m3，預測金金屬量4.258t。

5 結(jié)論

（1）在笏山金礦床深部布設鉆孔（DKZK03）進行深部驗證，首次在1200m以下控制了金礦體。DKZK03鉆孔控制4層金礦體，其中Ⅰ1號金礦體為該區(qū)目前最深金礦體。

（2）對棲霞小廟前成礦潛力區(qū)（Ⅰ1、Ⅰ2礦體）和棲霞上峴成礦潛力區(qū)（Ⅱ號礦體）采用趨勢外推法對深部金礦體進行定位預測，為今后在棲霞地區(qū)金礦深部找礦勘查提供了依據(jù)，該地區(qū)金礦深部具有較大的找礦潛力。

參考文獻：

［1］ 田瑞聰.蓬萊棲霞金礦帶成礦作用地球化學［D］.北京：中國地質(zhì)大學（北京），2022：1120.

［2］ 田杰鵬.膠東棲蓬福礦集區(qū)中生代金多金屬礦區(qū)域成礦作用［D］.北京：中國地質(zhì)大學（北京），2020：1140.

［3］ 單偉，遲乃杰，熊玉新，等.山東省蓬萊棲霞金成礦帶地殼結(jié)構(gòu)對金成礦的約束［J］.山東國土資源，2023，39（3）：513.

［4］ 鄧軍，王慶飛，張良，等.膠東型金礦成因模型［J］.中國科學：地球科學，2023，53（10）：23232347.

［5］ 楊偉光.膠東蓬萊金礦區(qū)控礦構(gòu)造特征及構(gòu)造演化研究［D］.北京：中國地質(zhì)大學（北京），2018：176.

［6］ 侯爵，武澄瀧，俞貴平，等.膠東金礦集中區(qū)深部結(jié)構(gòu)研究進展［J］.地球物理學進展，2022（37）：952963.

［7］ 楊立強，鄧軍，王中亮，等.膠東中生代金成礦系統(tǒng)［J］.巖石學報，2014（30）：24472467.

［8］ 翟明國，范宏瑞，楊進輝，等.非造山帶型金礦—膠東型金礦的陸內(nèi)成礦作用［J］.地學前緣，2004（1）：8598.

［9］ 于學峰，田京祥，李大鵬，等.膠西北金礦集區(qū)基于超殼走滑斷裂構(gòu)造體系的金成礦機制探討：來自深反射地震探測的證據(jù)［J］.山東國土資源，2023，39（11）：19.

［10］ 朱日祥，孫衛(wèi)東.大地幔楔與克拉通破壞型金礦［J］.中國科學：地球科學，2021（51）：14441456.

［11］ 楊立強，魏瑜吉，王偲瑞，等.膠東金礦床中關(guān)鍵金屬資源儲量估算與潛力初探［J］.巖石學報，2022，38（1）：922.

［12］ 于學峰，李大鵬，田京祥，等.山東金礦深部勘查進展與成礦理論創(chuàng)新［J］.山東國土資源，2018，34（5）：113.

［13］ 李敏，劉傳娥，遲乃杰，等.膠東蓬萊棲霞金礦帶侯格莊金礦流體包裹體特征及其地質(zhì)意義［J］.山東國土資源，2022，38（10）：1018.

［14］ 宋明春，林少一，楊立強，等.膠東金礦成礦模式［J］.礦床地質(zhì)，2020，39（2）：215236.

［15］ 李洪奎，陳國棟，耿科，等.山東膠東金礦成礦地質(zhì)背景初論［J］.山東國土資源，2018，34（5）：3441.

［16］ 李愛民，呂軍陽，張朋，等.山東蓬萊地區(qū)金礦成礦背景與成礦規(guī)律探討［J］.山東國土資源，2014，30（7）：1417.

［17］ 王立功，王英鵬，朱繼托，等.山東棲霞地區(qū)水系沉積物地球化學特征及找礦預測［J］.山東國土資源，2017，33（5）：16.

［18］ 趙磊，張儒誠，孫偉清，等.華北克拉通東部膠北地體早前寒武紀表殼巖系：研究進展與綜述［J］.巖石學報，2023，39（8）：22112237.

［19］ 朱學強，秦杰，郭晶，等.膠東福山—棲霞多金屬成礦帶成礦巖體成因機制及構(gòu)造背景［J］.地質(zhì)論評，2022，68（6）：20332052.

［20］ 時偉.膠東金礦幔源成礦流體深部演化與深部找礦［J］.山東國土資源，2015，31（10）：1216.

［21］ 侯爵，徐濤，艾印雙，等.膠東成礦省巨量金成礦模型：來自地殼速度結(jié)構(gòu)的約束［J］.中國科學（地球科學），2023，53（12）：29372952.

［22］ 宋明春，周明嶺，鮑中義，等.膠東深部金礦三維成礦預測及巨大的資源潛力［J］.金屬礦山，2023（1）：246258.

［23］ 戴長國.成礦深度研究及其對深部找礦的指示意義：以膠東典型焦家式金礦為例［J］.山東國土資源，2021，37（5）：916.

［24］ 霍光，宋國政，閆春明，等.膠東紗嶺金礦床構(gòu)造疊加暈特征及深部找礦預測［J］.山東國土資源，2019，35（9）：1622.

［25］ 徐方.膠東地區(qū)中生代金礦床成礦規(guī)律與成礦模式［D］.北京：中國礦業(yè)大學（北京），2019：120.

［26］ 王瑞良，張招崇，曾慶棟，等.膠東棲霞金礦集區(qū)早白堊世花崗巖形成時代及地質(zhì)意義［J］.大地構(gòu)造與成礦學，2019，43（1）：186198.

［27］ 李理想，李洪奎，韓學林，等.膠東棲霞山城金礦流體包裹體特征及其成礦意義［J］.山東國土資源，2023，39（2）：916.

［28］ 嚴育通.膠東石英脈型和蝕變巖型金礦關(guān)系的成因礦物學研究［D］.北京：中國地質(zhì)大學（北京），2012：1129.

［29］ 史宏江，王翠彭，田少斌，等.山東棲霞呂家金礦地質(zhì)特征及礦床成因分析［J］.山東國土資源，2021，37（6）：1927.

［30］ 李洪奎，卜文峰，禚傳源，等.山東膠東礦集區(qū)中生代構(gòu)造事件與金礦成礦作用［J］.山東國土資源，2017，33（1）：110.

［31］ 廖明偉，陳玉菡，謝軍民.山東省棲霞市笏山西陡崖金礦地質(zhì)特征及深部找礦遠景［J］.山東國土資源，2014，30（4）：2126.

［32］ 鄒鍵.膠東臧家莊盆地邊緣金及多金屬礦礦床特征與成礦預測［D］.北京：中國地質(zhì)大學（北京），2018：177.

［33］ 劉冬生，李彬，許凱磊，等.棲霞市馬耳夼石英脈型金礦地質(zhì)特征及找礦方向［J］.山東國土資源，2016，32（6）：1418.

［34］ 史宏江.山東棲霞馬家疃金礦區(qū)地質(zhì)特征及找礦前景分析［J］.山東國土資源，2016，32（10）：511.

［35］ 魯楠，周明磊，沈立軍，等.膠東馬家窯金礦床流體包裹體及S同位素特征研究［J］.山東國土資源，2023，39（5）：916.

［36］ 張國剛.山東棲霞區(qū)域地質(zhì)特征及找礦標志和方向［J］.黃金科學技術(shù)，2008（4）：5457.

［37］ 王佳良，孫豐月，王力，等.山東棲霞馬家窯金礦床地質(zhì)及成因探討［J］.黃金，2013，34（6）：1420.

［38］ 智云寶，孫海瑞，李風華.山東棲霞笏山金礦床成因：元素地球化學與流體包裹體證據(jù)［J］.吉林大學學報（地球科學版），2020，50（5）：15521569.

［39］ 史宏江.山東棲霞馬家窯金礦區(qū)地質(zhì)特征及找礦前景分析［J］.山東國土資源，2016，32（10）：511.

［40］ 戴廣凱，李秀章，朱學強，等.棲霞福山多金屬成礦帶地質(zhì)特征及其成礦作用［J］.山東國土資源，2021，37（12）：2434.

［41］ 付晨晨.膠西北棲霞孫疃地區(qū)金成礦特征及成礦預測［D］.南昌：東華理工大學，2014：149.

［42］ 黃鑫，宋倩，王勇軍，等.山東棲霞—大柳行地區(qū)金礦資源潛力評價［J］.黃金科學技術(shù)，2020，28（1）：2131.

［43］ 彭雪峰.膠東棲霞大霞址斷裂南段金礦地質(zhì)特征及成因研究［J］.山東國土資源，2022，38（11）：2130.

［44］ 王佳良.山東棲霞馬家窯金礦床地質(zhì)特征及礦化富集規(guī)律研究［D］.長春：吉林大學，2013：176.

［45］ 徐萌.山東省棲霞市笏山金礦構(gòu)造特征及其成礦作用分析［D］.石家莊：河北地質(zhì)大學，2016：157.

［46］ 智云寶，孫海瑞，李風華.山東棲霞笏山金礦床成因：元素地球化學與流體包裹體證據(jù)［J］.吉林大學學報（地球科學版），2020，50（5）：15521569.

［47］ 彭雪峰.膠東棲霞大霞址斷裂南段金礦地質(zhì)特征及成因研究［J］.山東國土資源，2022，38（11）：2130.

［48］ 廖小明，史宏江，孫曉東，等.山東省棲霞臺前地區(qū)控礦構(gòu)造與找礦預測［J］.黃金，2020，41（6）：713.

Analysis on Geological Characteristics and Resource Potential of Hushan-Xidouya Gold Deposit in Qixia City in Shandong Province

LIU Qiang1，ZHANG Yan1，CAI Na1， ZHANG Chao1，LI Dapeng1，SHI Hongjiang2，WEI Pengfei1，GENG Ke1，XIE Wei1

（1.Shandong Institute of Geological Sciences， Key Laboratory of Gold Mineralization Processes and Resource Utilization，MNR，Shandong Provincial Key Laboratory of Metallogenic Geological Process and Resource Utilization， Shandong Ji'nan 250013， China；2.No." 273 Geological Brigade of Shandong Nuclear Industry， Shandong Yantai 264006， China）

Abstract： Hushan-Xidouya gold deposit in Qixia city in Jiaodong area is located in the south of Penglai-Qixia gold metallogenic belt. It is a crushing altered rock type gold deposit occured in the contact zone between Neoarchean metamorphic basement and Mesozoic granite. A detailed field geological survey has been carried out in the surrounding areas， and the deposit has been studied as well. DKZK03 borehole has beem laid in the deep part of No.1 main ore body， and 4 layers of gold ore bodies have been controlled，with a cumulative true thickness of 6.85m. Among them， the buried depths of No. I-1 and I-2 gold ore bodies are 1273.80～1275.45m and 1072.95～1076.55m respectively. The I-1 gold ore body is currently the deepest in the area. The trend extrapolation method is used to locate and predict the deep gold deposits in the area. It will provide a basis for future exploration and exploration of gold deposits in Qixia area. The deep gold deposits in the area have great potential for exploration

Key words： Gold deposits;deep drilling;metallogenic predication;Hushan-Xidouya;Qixia city in Shandong province