黃河流域長孝巖溶水系統水文地質特征及蓄水模式分析

摘要：基于在黃河流域長孝巖溶水系統開展的1∶5萬水文地質調查工作，系統總結了該水文地質單元的邊界條件、含水巖組分布及富水性空間分布特征，分析了區內地下水補給、徑流、排泄條件與地下水動態變化特征，探討了黃河長孝單元與周邊地下水的相互作用。黃河向北側補給孔隙水，向南側補給孔隙水及巖溶水，其中對孔隙水側滲補給帶范圍至王營村—龐道口—太平莊—姚河門村一帶，對巖溶水補給主要發生在望口山村至柳河圈村一帶、孝里鎮勝利村附近。圈定了以朱砂洞組為含水巖組的雙泉鎮段店村裂隙巖溶水富水地段、以三山子組馬家溝群為含水巖組的歸德孝里馬集裂隙巖溶水富水地段，分析其富水機理，總結了層間巖溶斷層阻水型及單斜斷裂型兩種基巖蓄水構造模式，可為該區域找水定井提供參考。

關鍵詞：長孝巖溶水系統；地下水動態特征；蓄水模式；黃河流域

中圖分類號：P641.134""" 文獻標識碼：A""" doi：10.12128/j.issn.16726979.2024.03.015

0 引言

我國北方巖溶水多以相對獨立單元進行循環，構成一系列規模不等的巖溶水系統［1］，具有“巖溶含水層連續沉積巨厚，資源要素構成多、地下水資源分布不均、轉化關系復雜、環境質量脆弱”等特點［2］。巖溶地下水在魯中地區城鄉居民供水、工農業發展及能源基地建設中發揮著不可替代的支撐性作用，隨著經濟社會快速發展，地下水開采量急劇增加，人類活動影響加劇，資源型缺水及水質型缺水已成為制約區域社會經濟發展的重大問題［34］。

濟南市面臨供水和保泉的雙重壓力，水資源供需矛盾尤為嚴峻［5］。長孝水文地質單元位于黃河流域濟南段，為典型的單斜巖溶水系統，區內匯水面積廣，地下調節庫容大、水量水質穩定，建有長孝水源地（曹樓、三官廟、薛莊井組）向濟南市區供水約5萬m3/d。前人對該巖溶水系統的研究主要針對長孝水源地開展，例如1986年開展了長清孝里鋪地區供水水文地質勘探，2011年開展了長孝水源地曹樓井組階段性開采抽水試驗等，而對巖溶水上游補給區及孝里鎮以西、黃河北等下游排泄區勘查精度仍然較低，對于黃河水與周邊地下水的相互作用關系研究較少。

本文在總結前人研究基礎上，結合山東省自然資源廳地質勘查項目“山東省1∶5萬孝里、石橫幅水文地質調查”（魯勘字〔2021〕46號），系統掌握長孝巖溶水系統水文地質條件、動態特征、地下水循環及賦存規律，總結蓄水模式，尋找優質富水地段作為濟南市后備供水水源，對指導地下水資源開發利用規劃、城鎮建設、服務黃河流域生態保護和高質量發展等具有重要意義。

1 研究區概況

1.1 氣象水文

長孝巖溶水系統屬暖溫帶大陸性季風氣候，四季分明。區內多年平均氣溫14.2℃，無霜期202d，年平均降水量為657.7mm（1951—2022年）。降水多集中在7—8月，占全年降水總量的60%左右。

研究區位于黃河流域，地表水系較為發育，黃河是區內最主要河流，自平陰縣望口山村進入研究區，沿濟南市境逶迤東北，其支脈河流均從東側南岸匯入，主要有南大沙河、清水溝等［6］。

1.2 地形地貌

研究區位于魯中南構造侵蝕為主中低山丘陵區與魯西北沖積平原區的過渡地帶［7］，地形東南高，西北低。東南部屬低山丘陵區，基巖裸露良好，海拔一般80～441m，切割深度一般為50～300m，屬構造剝蝕成因。西北部黃河沖積平原呈自東南向西北緩傾的平原地形，海拔一般30～35m。

1.3 地層巖性

研究區位于華北地層區、魯西地層分區［8］。區內地層出露較為齊全，由東南向西北依次為新太古代泰山巖群變質巖；古生代寒武紀朱砂洞組白云巖，饅頭組頁巖、砂巖及薄層灰巖，張夏組灰巖，崮山組薄層灰巖夾頁巖；炒米店組竹葉狀灰巖、泥質條帶灰巖；奧陶紀三山子組白云巖，馬家溝群厚層灰巖、白云巖。石炭二疊系隱伏于黃河北趙官鎮一帶，以砂質頁巖、黏土巖為主。新生代新近系僅發育明化鎮組，分布于黃河西北部沖積平原，以鈣質黏土為主；第四系在區內主要分布于山前沖洪積平原及黃河沖積平原，以粉質黏土、粉細砂為主。

1.4 構造體系

研究區地處魯中隆起，位于泰山凸起及齊河潛凸起兩個構造單元的交界部位［9］，區內脆性斷裂發育，以張性、張扭性為主。主要發育NNE向斷裂，沿走向延伸多在10km以上，例如孝里鋪斷裂、北傅莊斷裂。NNW—NW向斷裂為區內較發育斷裂，規模較大有馬山斷裂、牛角店斷裂（圖1）。

2 水文地質特征

2.1 邊界條件

長孝單斜巖溶水系統是一個相對獨立的水文地質單元，其東邊界為馬山斷裂；南邊界為地表地下分水嶺，南邊界的東部分水嶺為晚太古代侵入巖，西部分水嶺為寒武系碳酸鹽巖夾碎屑巖；西邊界南部為黃山巖脈，北部為牛角店斷裂；北邊界為灰巖頂板埋深400m界線，西部位于黃河以北，東部在黃河南長清區的苗莊—東倉莊一帶，總面積約935.4km2。黃山巖脈具有阻水性，牛角店斷裂具有導水性；馬山斷裂南部崗辛孫莊段，斷層兩側寒武系與奧陶系接觸，導水性較差，為阻水段；老屯前隆段斷層兩側均為奧陶紀灰巖，具透水性（圖2）［1012］。因此，長孝巖溶水系統南、北邊界均為隔水邊界，東、西邊界的北部均為透水邊界，南部為隔水邊界。西邊界北部牛角店斷裂為一流入系統內的流量邊界，東邊界馬山斷裂長清老屯段具透水性，流入流出隨兩側地下水開采情況而定。

2.2 含水巖組與富水性

依據地下水賦存條件、水理性質、含水介質及其空隙特點，區內含水巖組主要可以分為松散巖類孔隙含水巖組、碳酸鹽巖類裂隙巖溶含水巖組及巖漿巖變質巖類裂隙含水巖組。受地貌單元、地層巖性、構造等影響，研究區內地下水的流動和賦存特征差異性明顯（圖3）。

（1）松散巖類孔隙含水巖組。分布在山間河谷、山前沖洪積平原及黃河沖積平原?？紫端饕患谀洗笊澈記_洪積扇附近，第四系沉積厚度從30～90m不等，含水層巖性以中粗砂，粗砂夾礫為主，在沖洪積扇頂部較粗，含卵礫較多，含水層2～3層，累計厚度10～15m。在沖洪積扇中部，含水層增至4～6層，厚度20～25m。歸德鎮莊樓村至路莊村一帶，富水性較好，單井涌水量一般500～1000m3/d，局部可達1000～3000m3/d。

（2）碳酸鹽巖裂隙巖溶含水巖組。在區內廣泛分布，寒武奧陶紀碳酸鹽巖地層呈NE—SW向展布，傾向NW。在丘陵山區，含水層巖性主要為三山子組白云巖、炒米店組灰巖、張夏組灰巖和朱砂洞組白云巖，分布地勢較高，富水性整體較差，單井涌水量一般小于500m3/d，局部受斷裂構造等影響形成小型富水地段。奧陶紀馬家溝群在山前平原區多隱伏第四系之下，地下水徑流變緩，含水層裂隙、巖溶發育，富水性增強，單井涌水量一般1000～5000m3/d，在歸德、孝里鎮以西及黃河北一帶單井涌水量可達5000～10000m3/d。

（3）巖漿巖變質巖類裂隙含水巖組。分布于研究區東南部雙泉鎮、馬山鎮東一帶基巖山區，含水巖組主要為新太古代侵入巖，巖性以花崗閃長巖、二長花崗巖類為主，風化帶厚度一般小于20m，地下水賦存在表層風化裂隙中，富水性整體較差，一般小于100m3/d。

2.3 地下水補給、徑流、排泄特征

（1）松散巖類孔隙水。主要接受大氣降水補給，豐水期接受南大沙河、清水溝等河道雨洪水滲漏補給及下伏巖溶水的頂托補給，在沿黃一帶接受黃河測滲補給。地下水徑流方向與地表水徑流方向一致，總體上由東南向西北方向徑流?？紫端饕判狗绞綖閺搅髋判?、農業分散開采及蒸發排泄。

（2）碳酸鹽巖類裂隙巖溶水。長孝單斜地層整體傾向NW，地表水、地下水流向與碳酸鹽巖地層傾向基本一致，形成順置式巖溶地下水系統［13］。大氣降水入滲和河流、水庫滲漏是巖溶水的主要補給來源。在山前平原的第四系較薄處，巖溶水接受孔隙水下滲補給。黃河滲漏補給巖溶水發生于平陰縣望口山至柳河圈一帶、長清區孝里鎮勝利村附近，第四系較薄，馬家溝群灰巖與黃河直接接觸，河水沿層面、裂隙和破碎帶滲入補給巖溶水。黃河北一帶巖溶水還接受牛角店斷裂以西平陰東阿巖溶水系統的側向補給。該區域巖溶地下水整體由東南丘陵山區向西北平原區徑流，向北徑流過程中在黃河北趙官鎮一帶受到石炭二疊系阻擋，又被上覆第四系、新近紀明化鎮組較厚的黏土層阻隔，巖溶水流向轉向NE或NNE。排泄方式主要為側向徑流排泄及人工開采。

（3）巖漿巖變質巖類裂隙水。該區巖漿巖出露地勢較高，直接接受大氣降水補給，在地勢低處可接受地表水補給、松散層孔隙水補給及高處裂隙水徑流補給。由于地形坡度較大，大部分降水以表流產出，少量沿裂隙發育方向滲入地下形成徑流。地下水運動受地形條件控制，隨地形坡向成散流狀態，在溝谷下游多以潛流排泄于松散層中。在溝底及構造破碎帶發育處，常呈下降泉方式排泄。裂隙水整體表現為就地補給、淺部運動、短途排泄。

2.4 地下水流場特征

地下水位是水平衡的主要反映，其主要受降水、蒸發、徑流、地形、地貌和巖性特征等自然條件及開采等人為活動影響［14］。根據2022年枯水期（6月）及豐水期（10月）統測數據分別繪制了研究區地下水等水位線。由圖4可知，巖溶水整體流向為自東南向西北徑流，在丘陵山區水力梯度較大，至山前平原一帶巖溶水徑流變緩，水力梯度減小。從年內變化特征來看，上游補給徑流區巖溶水枯豐期變幅較大，一般10～20m，下游排泄區水位年內變幅減小，一般3～5m。影響巖溶地下水流場特征的主要因素為地形及大氣降水。

孔隙水流向主要受控于地形，總體由東南向西北方向流動。黃河在研究區為地上河，其水位在枯、豐水期均高于兩側孔隙水水位。在黃河以南，孔隙水水位表現為在山前及黃河附近孔隙水水位較高，枯水期水位29～30m，豐水期水位32.7～33m左右，而孝里鎮以西至洼里村北部一帶水位較低，枯水期水位28m左右，豐水期水位32.6m左右，水位年內變幅約3～4.6m。由于孝里山前孔隙水自東南向西北徑流，黃河側滲自西北向東南徑流，在太平莊等交匯部位水流滯緩（圖5），黃河側滲補給帶范圍至王營—龐道口—太平莊—姚河門一帶。在黃河北平原，表現為黃河附近孔隙水水位最高，向西北方向水位逐漸降低，枯水期水位24～29m、豐水期水位29.5～32.5m，水位年內變幅約3.5～5m。由此可見，黃河具備側滲孔隙水的水動力場條件，枯水期水力梯度較大，豐水期水力梯度較小。影響孔隙地下水流場特征的主要因素為地形、大氣降水及黃河側滲。

2.5 地下水動態特征

影響地下水水位動態變化因素主要有氣象、水文、地層結構、含水層的埋藏條件、農田灌溉及人工開采等，是自然與人為作用綜合影響的結果［15］。以曹樓水源地巖溶水監測點為例（圖6），其水位動態曲線在年內表現為典型的“陡升緩降型”，反映受大氣降水影響較快，地下水具有補給迅速、補給量充沛、排泄相對緩慢等特點。2022年地下水位介于28.7～32.69m，年變幅3.99m，年均水位31.39m，較上年偏高3.0m。年水位最大值出現在10月16日，最小值出現在6月21日。年末水位較年初水位下降0.95m，為負均衡。

根據長序列監測數據（圖7），1987—2022年，曹樓附近地下水位趨于動態平衡，水位標高25～35m，埋深0～10m，水位變幅10m左右；2002年為枯水年，2003年7月6日水位標高降至多年來最低值22.88m，2003—2004年為連續豐水年，地下水位顯著上升，在2004年9月1日水位達到多年來最大值35.17m。2004—2020年，水位總體呈下降趨勢，2020年枯水期下降至近年來最低值24.86m，較2004年枯水期水位下降4.49m，平均每年降幅約28cm。2021年、2022年連續兩年降水量豐富，區域地下水得到充分補給，水位迅速回升，達到2013年來的最大值。從多年動態可以看出，區域巖溶水受大氣降水影響明顯，在降水相對平均和偏少的情況下，巖溶水位呈逐漸下降趨勢，說明區域巖溶水處于負均衡狀態。同時雨季來臨時大氣降水快速補給導致巖溶水水位迅速回升，當降水量豐富時，能夠補償多年開采形成的水位降幅，說明區內巖溶水現有的開采量情況在多年尺度下具備可調蓄性。

研究黃河與周邊地下水的水位及動態變化可以判斷其相互作用關系。以黃河水位（平陰黃金浮橋段）、安城鎮柳河圈村巖溶水監測點（XZ31）及孝里鎮宋莊村孔隙水監測點（XZ33）為例分析。該兩眼監測井均位于黃河南岸，其中XZ31距黃河約600m，井深102m，含水巖組為馬家溝群北庵莊組灰巖；XZ33距黃河500m，井深55m，含水巖組為第四系粉細砂。

由圖8可知，黃河沿岸巖溶水及孔隙水水位升降及變化幅度具有明顯的一致性。2022年9月中旬黃河水位顯著抬升，而后在9月底巖溶水及孔隙水出現響應，地下水位明顯抬升，在10月5日左右達到年內水位最高值，后隨著黃河水位緩慢下降，地下水位也隨之緩慢持續下降，10月中旬至11月底期間，柳河圈附近巖溶水水位與黃河水位高度基本一致，至2023年2月底，地下水位受農灌影響出現顯著下降。由此可見，黃河南岸柳河圈宋莊一帶，地下水主要接受大氣降水入滲補給、黃河側滲補給，巖溶水與孔隙水水力聯系密切。

2.6 巖溶發育特征

2.6.1 地表巖溶發育特征

濟南是我國北方巖溶有代表性的地區之一［16］，丘陵地帶碳酸鹽巖裸露區地表溶蝕較為發育，巖溶形態有溶痕、溶溝、溶槽、溶孔、溶洞等，有利于大氣降水及地表徑流入滲補給地下水。

2.6.2 地下巖溶發育特征

受氣候條件、地層產狀和規模、巖石的礦物成份、地質構造、水的溶蝕性和流動性等因素的影響，巖溶含水介質具有強烈的不均一性和各向異性［1718］。根據野外調查及本次施工鉆孔揭露，在水平方向上，從補給區—徑流區—排泄區，裂隙巖溶發育深度、程度逐漸增強，上游補給區多為溶蝕裂隙，排泄區多為溶孔、溶洞。在垂向上，不同巖層的生成時代、沉積環境不同，導致礦物成分及結構構造均有一定差別，因此巖溶發育表現出明顯差異。

寒武紀碳酸鹽巖由于多含頁巖、泥質條帶灰巖、砂巖等碎屑巖夾層，各組灰巖受頁巖隔水相互水力聯系差［19］，水位呈階梯狀，無統一水面，地下水主要沿構造裂隙運移，以垂向入滲作用為主，溶蝕裂隙面常見黃褐色鐵質浸染，為地下水活動痕跡，地層溶蝕作用整體較差。

奧陶紀碳酸鹽巖一般為質純的厚層灰巖、白云巖，多隱伏于山前平原，地下水徑流以水平運動為主，流速相對滯緩，巖溶發育多呈層狀分布，巖心常見大量溶蝕裂隙、蜂窩狀孔洞等巖溶現象，溶孔直徑一般0.5～3cm，溶洞直徑最大可達10～15cm。

3 富水地段及地下水蓄水模式

3.1 富水地段

3.1.1 雙泉鎮段店村裂隙巖溶水富水地段

段店富水地段位于雙泉鎮長孝巖溶水系統上游補給區（圖9），其范圍西至北傅莊斷裂，東北至段店白樓斷裂，南至西坦村，面積約0.79km2，該區域主要含水巖組為朱砂洞組中厚層含燧石結核白云巖，受構造斷裂帶影響，地層巖溶裂隙較為發育，為地下水賦存提供了良好的蓄水空間與導水通道，主要接受大氣降水和南大沙河滲漏補給，西南、東南2個方向為地下水主要補給區，地形上呈開口向南的簸箕狀展布，匯水條件好，地下水順地層由南向北徑流，在北傅莊斷裂及段店白樓斷裂交會楔形地帶前緣處，受到斷裂帶及侵入巖體的阻擋，巖溶地下水在斷裂帶前緣滯留、儲存，形成“楔形”富水地段，單井涌水量達5000～10000m3/d。經XLZK01鉆孔驗證，該井成井深度84.4m，水位降深1.74m時，單井涌水量達6171m3/d。

3.1.2 歸德孝里馬集裂隙巖溶水富水地段

該富水地段位于巖溶水系統下游排泄區，南部大致以東黃山組與三山子組地層界線為界、東北邊界為馬山斷裂，西北邊界至黃河北趙官一帶灰巖頂板埋深400m界線，面積約249.5km2。

收集并統計了長孝水源地、孝里洼地及東阿縣下碼頭水源地附近的鉆孔資料，該富水地段主要含水巖組為奧陶紀北庵莊組灰巖、東黃山組白云巖及三山子組細晶白云巖，巖溶發育在50～500m深度均有分布，大段蜂窩狀孔洞多集中在100～300m深度（圖10），組成脈狀網絡蓄水體系［20］，單井涌水量可達5000～12000m3/d，降深一般0.4～5m不等，富水性好，具有巨大的儲水空間，是一個天然形成的巖溶地下水庫。

例如本次在黃河以南孝里洼地施工的XLZK02鉆孔，孔深399.7m，蜂窩狀孔洞等巖溶發育段在土峪組、北庵莊組、東黃山組、三山子組均有分布，降深1.53m時，單井涌水量達6970m3/d。在黃河以北馬集鎮齊莊村施工的XLZK04鉆孔，含水層為奧陶系五陽山組，孔深332.1m，巖溶發育段位于193.4～233.3m，該孔鉆進過程中自流，自流高度0.48m，自流量達2688m3/d，降深2.71m時，單井涌水量達到7020m3/d。

3.2 蓄水構造模式

蓄水構造是由透水層（帶）、隔水層（體）及地下水補給排泄條件組成的能在水交替循環過程中不斷富集和儲存地下水的開放性地質構造［21］。賈德旺等研究了魯南地區巖溶蓄水構造模式，認為主要包括阻水型、斷塊型、接觸型和層間巖溶型等四類［22］。結合工作區的地層分布、地質構造及本次開展的水文地質調查、物探、鉆探、抽水試驗等研究成果，總結了區內兩種蓄水構造模式：層間巖溶斷層阻水型蓄水構造、單斜斷裂型基巖蓄水構造模式。

3.2.1 層間巖溶斷層阻水型蓄水構造模式

本次針對基巖山區廣泛分布的寒武系朱砂洞組提出了層間巖溶斷層阻水型蓄水構造模式。朱砂洞組埋藏于地下，巖性以白云巖、泥質白云巖為主，在該區域厚度15～20m，與寒武系基底侵入巖不整合接觸，受多期構造活動及熱液蝕變等因素影響，朱砂洞組底部順層形成不同規模的溶蝕空隙、層間虛脫、碎裂巖帶或揉皺帶等［23］，形成地下水蓄存與運移的空間環境。地下水在巖層傾伏端及埋藏區段匯流，如遇到斷層阻水發生滯留，則形成一定范圍的富水段［24］。地下水富集帶常位于阻水體附近的地形低洼處，特別是來水方向的一側。本次勘探發現的雙泉鎮段店村富水地段屬于該種蓄水構造模式（圖11）。

3.2.2 單斜斷裂型蓄水構造模式

單斜構造巖層或巖塊受地質應力作用常會產生明顯的破裂位移，斷層兩側巖石破碎，構造裂隙發育，易于接受水流侵蝕，巖溶不斷發育，不僅為地下水提供儲存空間，更為地下水的徑流、排泄提供了運移通道。孝里洼馬集富水地段屬于該種蓄水構造模式，該區域主要發育一系列NE向斷裂，孝里鋪斷裂、黃河下部隱伏斷裂等將碳酸鹽巖地層切割，形成階梯狀單斜斷陷構造（圖12），斷層兩盤均為可溶性巖石，在地下水作用下形成裂隙、溶孔、溶洞發育的層狀巖溶網絡。在斷層兩側接觸帶附近，往往富水性強，單井涌水量一般在5000～10000m3/d。

4 結論

在山東省1∶5萬孝里、石橫幅水文地質調查工作的基礎之上，分析了長孝巖溶水系統水文地質條件、地下水動態特征、三水轉化關系、地下水循環及賦存規律。

（1）長孝巖溶水系統基巖地層呈單斜構造產出，傾向NW。主要發育松散巖類孔隙含水巖組、碳酸鹽巖類裂隙巖溶含水巖組及巖漿巖變質巖類裂隙含水巖組。巖溶水主要富水巖層為奧陶紀馬家溝群、三山子組、寒武紀朱砂洞組。區內地下水主要補給來源為大氣降水、地表水滲漏補給，地下水水流向整體由東南向西北徑流，受構造、地層等阻水條件影響后可形成富水區域，地下水排泄方式主要為人工開采、側向徑流。

（2）地下水動態特征表明，區域地下水位受大氣降水影響明顯，巖溶水多年水位呈逐漸下降趨勢，豐水年時水位迅速回升，能夠補償多年開采形成的水位降幅，區內巖溶水現有的開采量在多年尺度下具備可調蓄性。

（3）黃河向北側補給孔隙水；向南側補給孔隙水及巖溶水，其中對孔隙水側滲補給帶范圍至王營村—龐道口—太平莊—姚河門村一帶，對巖溶水補給主要發生在望口山至柳河圈一帶、孝里鎮勝利村附近。黃河南岸柳河圈—宋莊一帶巖溶水與孔隙水水力聯系密切。

（4）圈定了以朱砂洞組為含水巖組的雙泉鎮段店村裂隙巖溶水富水地段、以三山子組馬家溝群為含水巖組的歸德孝里洼馬集裂隙巖溶水富水地段，總結了層間巖溶斷層阻水型蓄水構造模式及單斜斷裂型蓄水構造模式，可為該區域找水定井提供參考。

參考文獻：

［1］ 王振興，李向全，侯新偉，等.三姑泉域巖溶地下水分布特征及子系統識別［J］.地球與環境，2020，48（2）：228239.

［2］ 梁永平，申豪勇，高旭波.中國北方巖溶地下水的研究進展［J］.地質科技通報，2022，41（5）：199219.

［3］ 梁永平，申豪勇，趙春紅，等.對中國北方巖溶水研究方向的思考與實踐［J］.中國巖溶，2021，40（3）：363380.

［4］ 馬祖陸，楊黃成．國外巖溶地區水資源開發利用技術及其對我國的啟示［J］．廣西科學院學報，2022，38（4）：352363.

［5］ 丁冠濤，李常鎖，魏善明，等.玉符河外源回灌對巖溶地下水影響研究［J］.中國巖溶，2023，42（5）：907916.

［6］ 張赟，王鵬，王斌.濟南市長清區水資源開發利用與保護探討［J］.山東水利，2022（9）：4849.

［7］ 齊歡.濟南長清巖溶地下水位與降水量、黃河水位的相關性［J］.長江科學院院報，2023，40（7）：2731.

［8］ 宋明春.山東省大地構造單元組成、背景和演化［J］.地質調查與研究，2008（3）：165175.

［9］ 崔華偉，鄭建常，柴光斌，等.2020年2月18日濟南長清M4.1地震震源區發震構造分析［J］.地球物理學進展，2022，37（1）：110.

［10］ 薄克庭，蔡有兄.長清孝里鋪水文地質單元巖溶地下水資源量計算及開采潛力分析［J］.山東國土資源，2015，31（12）：3742.

［11］ 韓琳，顏翠翠，趙振偉，等.基于VisualMODFLOW的地下水數值模擬研究及水位預測［J］.山東國土資源，2021，37（3）：6774.

［12］ 束龍倉，張曼琦，李虎，等.兩種典型邊界條件下的趵突泉泉域水均衡分析［J］.南水北調與水利科技（中英文），2021，19（3）：427436.

［13］ 潘曉東，唐建生，蘇春田.魯中南地區巖溶水系統模式：以泰安市岱岳區為例［J］.桂林理工大學學報，2014，34（1）：3944.

［14］ 劉小婧.地下水流場的時空演變特征及驅動力分析［J］.地下水，2023，45（5）：8588.

［15］ 李波，王金曉，吳璇，等.山東萊蕪盆地東部水文地質條件及富水塊段特征［J］.中國巖溶，2020，39（5）：637649.

［16］ 奚德蔭.魯中南地區巖溶水文地質條件及其特征［J］.中國巖溶，1988，7（3）：4348.

［17］ 申豪勇，梁永平，徐永新，等.中國北方巖溶地下水補給研究進展［J］.水文，2019，39（3）：1521.

［18］ 李嚴，王家樂，靳孟貴，等.運用水文時間序列分析識別濟南泉域巖溶發育特征［J］地球科學，2021，46（7）：25832593.

［19］ 劉元晴，文冬光，呂琳，等.沂蒙山區典型斷陷盆地巖溶地下水系統特征：以萊蕪盆地為例［J］.地質科技通報，2022，41（1）：157167.

［20］ 武東強，邢立亭，蘭曉荀，等.濟南巖溶含水介質孔隙結構特征［J］.中國巖溶，2021，40（4）：680688.

［21］ 張福存，王新峰，李偉，等.缺水基巖山區蓄水構造類型劃分及其屬性分析［J］.水文地質工程地質，2022（2）：716.

［22］ 賈德旺.魯南山區基巖蓄水構造類型及找水定井方法［J］.地質學刊，2020，44（3）：318325.

［23］ 劉元晴，周樂，李偉，等.魯中山區下寒武統朱砂洞組似層狀含水層成因分析［J］.地質論評，2019，65（3）：653663.

［24］ 關琴，徐聰聰，葉倩，等.魯中南新汶盆地地下水蓄水模式分析［J］.山東國土資源，2021，37（5）：4148.

Analysis on Hydrogeological Characteristics and Water Storage Mode of Changxiao Karst Water System in the Yellow River Basin

ZHANG Wenqiang， TENG Yue， XU Qingyu， LIU Xiaotian， LIU Haoran， ZHANG Hailin

（No.801Hydrogeological Engineering Geological Brigade of Shandong Provincial Exploration Bureau of Geology and Mineral Resources，Research Center of Groundwater Environment Protection and Restoration Engineering Technology，Shandong Ji'nan 250014，China）

Abstract：Based on hydrogeological survey with the scale of 1∶50000 conducted in Changxiao karst water system in the Yellow River basin， the boundary conditions， distribution of water bearing rock formations， and spatial distribution characteristics of water richness of the hydrogeological unit have been summarized， groundwater recharge， runoff， discharge conditions， and dynamic changes in the area have been analyzed， and the interaction between the Yellow River and surrounding groundwater has been studied. The Yellow River supplies pore water to the north and pore water and karst water to the south. The range of lateral seepage supply zone for pore water is Wangying village PangdaokouTaipingzhuangYaohemen village. The supply of karst water mainly occurs in the area from Wangkoushan village to Liuhequan village and near Shengli village in Xiaoli town. The fractured karst water rich area are Duandian village in Shuangquan town， with Zhushadong formation as the water bearing rock group， and GuideXiaoliMaji fractured karst water rich area with Sanshanzi formation" Majiagou group as the water bearing rock group. Water rich mechanism has been analyzed. Two types of bedrock water storage structures have been summarized， they are interlayer karst" fault water blocking type and monoclinic" fault type. It will provide some references for water exploration and well location in this area.

Key words：Changxiao karst water system；dynamic characteristics of groundwater；rich groundwater area；groundwater storage model