• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    巖溶區(qū)生物醫(yī)藥建設(shè)項目對地下水環(huán)境影響分析

    2017-09-03 10:31:40陸海玉李旺林
    資源環(huán)境與工程 2017年4期
    關(guān)鍵詞:巖組水文地質(zhì)含水層

    陸海玉, 李旺林

    (1.山東省水利勘測設(shè)計院,山東 濟南 250013; 2.濟南大學 資源與環(huán)境學院,山東 濟南 250022)

    巖溶區(qū)生物醫(yī)藥建設(shè)項目對地下水環(huán)境影響分析

    陸海玉1,2, 李旺林2

    (1.山東省水利勘測設(shè)計院,山東 濟南 250013; 2.濟南大學 資源與環(huán)境學院,山東 濟南 250022)

    通過等效多孔介質(zhì)模型聯(lián)合反求水文地質(zhì)參數(shù)的方法建立地質(zhì)概化后的非穩(wěn)定流巖溶區(qū)地下水數(shù)值模型,經(jīng)與模擬區(qū)地下水流場、地下水位動態(tài)曲線擬合識別、驗證,其擬合度較高。利用該模型建設(shè)項目特征因子的源、強在模型預測期對地下水的影響范圍和程度進行預測,為預防和防治地下水污染提供科學依據(jù),對巖溶區(qū)相似水文地質(zhì)條件的地下水環(huán)境影響評價具有重要的借鑒意義。

    巖溶;水文地質(zhì);等效多孔介質(zhì);數(shù)值模型;環(huán)境影響評價

    隨著地下水資源的減少、地下水污染越來越嚴重等問題的出現(xiàn),國家對地下水環(huán)境保護越來越重視。地下水環(huán)境影響評價已成為巖溶區(qū)建設(shè)項目環(huán)評的一個重要環(huán)節(jié)。巖溶含水系統(tǒng)因其高度的非均質(zhì)性和各向異性,使地下水環(huán)境影響評價時常用的類比分析法和解析法受到很大局限性;采用數(shù)值法時,數(shù)值模型的準確性又受含水層概化、水文地質(zhì)參數(shù)取值等影響較大。本文通過等效多孔介質(zhì)模型聯(lián)合反求水文地質(zhì)參數(shù)的方法建立地質(zhì)概化后的非穩(wěn)定流巖溶區(qū)地下水數(shù)值模型,經(jīng)與模擬區(qū)地下水流場、地下水位動態(tài)曲線擬合識別、驗證,其擬合度較高。并利用該模型對建設(shè)項目特征因子的源、強在模型預測期對地下水的影響范圍和程度進行預測,為預防和防治地下水污染提供科學依據(jù),對巖溶區(qū)相似水文地質(zhì)條件的地下水環(huán)境影響評價具有重要的借鑒意義。

    1 項目背景

    擬建生物醫(yī)藥項目包括生物藥真核原液、生物藥制劑等6個生產(chǎn)車間以及污水處理站、危險化學品原輔材料庫、危險固廢暫存庫、一般固廢暫存庫、事故水池等。

    研究區(qū)南部為魯中低山丘陵區(qū),北部為山前傾斜平原,地勢東南高西北低,地形高差達700 m。地貌上由南向北依次為構(gòu)造剝蝕區(qū)、剝蝕堆積區(qū)、堆積平原區(qū)[1],建設(shè)項目位于構(gòu)造剝蝕區(qū)與剝蝕堆積區(qū)的交接部位。

    2 區(qū)域水文地質(zhì)概述

    2.1 含水巖組及補徑排條件

    根據(jù)含水介質(zhì)不同及地下水在介質(zhì)中的賦存、運移特征,研究區(qū)可分為5種地下水含水巖組[2]。

    2.1.1 第四系松散巖類孔隙含水巖組

    主要分布于研究區(qū)北半部山前沖洪積平原及河流沖洪積扇等地帶,含水層巖性主要為粉細砂、中粗砂、砂礫石層,厚度2~15 m,水位埋深2~18 m,單井涌水量400~5 000 m3/d不等。地下水主要接受大氣降水、地表河水入滲補給和巖溶水頂托補給;徑流總體流向是北北西向,在局部開采量較大的地段形成孔隙水降落漏斗;其排泄方式主要有人工開采、向北徑流入小清河及局部反補給下部巖溶水。模型建立時該含水層主要考慮對碳酸巖類巖溶裂隙水的補排影響。

    2.1.2 碎屑巖類裂隙或?qū)娱g巖溶裂隙含水巖組

    分布在研究區(qū)北部,被第四系松散堆積物覆蓋,含水層巖性為石炭—二疊系砂巖、薄層灰?guī)r。地下水主要賦存在構(gòu)造裂隙中,富水性差,在構(gòu)造裂隙發(fā)育處接受上游巖溶水補給,沿構(gòu)造裂隙向下游徑流排泄。該巖組表現(xiàn)為阻擋裂隙巖溶水的徑流通道,構(gòu)成了數(shù)值模型的北邊界。

    2.1.3 碳酸巖類裂隙巖溶含水巖組

    在研究區(qū)內(nèi)分布較廣,南部出露地表;北部平原區(qū)埋藏于第四系和石炭、二疊系之下,含水層巖性為奧陶系亮甲組、馬家溝組、八陡組灰?guī)r及寒武系鳳山組、張夏組灰?guī)r。裸露區(qū)水位埋深30~200 m,單井涌水量500~1 000 m3/d或<500 m3/d;隱伏區(qū)處于地下水徑流排泄帶,水位埋深5~40 m,單井涌水量1 000~10 000 m3/d。該層巖溶水主要接受大氣降水和地表水的入滲補給;其徑流方向基本上沿巖層傾向由南向北徑流;在北部受石炭系—二疊系砂巖、薄層灰?guī)r阻擋,主要以泉群型式排泄地表,其次向第四系及石炭系排泄和人工開采。

    2.1.4 碳酸巖夾碎屑巖巖溶裂隙含水巖組

    呈北西、南東向分布在研究區(qū)南部,出露位置較高,含水層巖性主要為寒武系崮山組、長山組、毛莊組、徐莊組薄層灰?guī)r及泥灰?guī)r,由于有砂頁巖間隔各層灰?guī)r互不聯(lián)系,裂隙巖溶發(fā)育較差,單井涌水量一般<500 m3/d。該巖組主要接受大氣降水,地下水儲存于層間灰?guī)r的巖溶裂隙中,基本沿巖層傾向向北運動,受構(gòu)造及河流切割作用,常以下降泉的形式排泄轉(zhuǎn)化成地表水,沿溝谷方向向北運動,當?shù)乇硭鹘?jīng)寒武系張夏組、鳳山組及奧陶系分布區(qū)時,地表水便大量漏失補給巖溶地下水。

    2.1.5 塊狀巖類風化裂隙含水巖組

    分布于研究區(qū)外,構(gòu)成了區(qū)域南部的分水嶺,巖性以花崗片麻巖類為主,因深部巖石堅硬致密,僅在淺部發(fā)育細小的風化裂隙,富水性差,單井涌水量50~100 m3/d。以大氣降水為補給源,其地下水的運動特點是就地補給、短距離排泄,呈散流狀態(tài)向溝谷匯集。

    2.2 巖溶發(fā)育規(guī)律及特征

    2.2.1 巖溶發(fā)育具有層控性

    研究區(qū)內(nèi)碳酸巖分布廣、厚度大,由于各巖層生成時代、沉積環(huán)境不同,礦物成分、化學成分、結(jié)構(gòu)、構(gòu)造均有差別,造成巖溶發(fā)育的差異。

    2.2.2 巖溶發(fā)育的分帶性

    (1) 水平分帶特征:從補給區(qū)到徑流排泄區(qū),巖溶發(fā)育有逐漸增強的特征。補給區(qū)地表溶溝、溶槽、溶孔較發(fā)育,地下巖溶發(fā)育相對較弱。在補給徑流區(qū)為中等巖溶發(fā)育,在徑流排泄區(qū)為強烈?guī)r溶發(fā)育。另外,沿干谷附近巖溶發(fā)育強烈,而地表分水嶺附近巖溶發(fā)育相對較弱。

    (2) 垂直分帶特征:上部巖溶發(fā)育,下部巖溶相對減弱。裸露灰?guī)r區(qū)位置高,在水位變幅帶附近巖溶發(fā)育。隱伏灰?guī)r區(qū)較強巖溶發(fā)育帶底界由南向北逐漸變深,較強巖溶發(fā)育帶以下為中等巖溶發(fā)育帶。

    3 數(shù)值模型的建立

    地下水數(shù)值模型是根據(jù)模擬區(qū)水文地質(zhì)條件,先建立地下水系統(tǒng)概念模型,再通過數(shù)值建模建立地下水流數(shù)學模型,通過不斷調(diào)試模型參數(shù)進行重復計算,使數(shù)值模擬工作成為“數(shù)值模擬實驗”,以觀測結(jié)果和計算結(jié)果擬合差為標準,不斷校正模型,最后達到利用數(shù)值模擬實驗方法研究地下水運動機理的目的[3]。

    3.1 水文地質(zhì)概念模型

    3.1.1 含水層概化

    建設(shè)項目場區(qū)位于構(gòu)造剝蝕區(qū)與剝蝕堆積區(qū)的交接部位,屬碳酸巖類巖溶裂隙水的補給徑流區(qū),上部第四系松散層常年基本無水,確定第一層穩(wěn)定含水層碳酸巖類巖溶裂隙水為本次評價的目標含水層。模型建立時重點考慮目標含水層,其它含水層主要考慮對碳酸巖類巖溶裂隙水的補排影響[4]。

    根據(jù)研究區(qū)含水層的水文地質(zhì)特征及補給條件,將下部碳酸巖類裂隙巖溶含水巖組、碳酸巖夾碎屑巖巖溶裂隙含水巖組概化為一層,即碳酸巖類巖溶裂隙含水層,該含水層為區(qū)域內(nèi)大厚度含水巖組,地下水流受構(gòu)造、地質(zhì)條件控制,總體在徑流過程中很少受阻,全區(qū)地下水存在著普遍的水力聯(lián)系,有較統(tǒng)一的自由水面(承壓水頭),在南部山區(qū)表現(xiàn)為潛水,至山前地帶由于上部分布的第四系更新統(tǒng)含礫粘土層及石炭系砂頁巖,過渡為承壓水[5]。

    將上部第四系松散巖類孔隙含水巖組概化為第四系松散巖類孔隙含水層,主要分布在北部山前平原地段;中部將第四系更新統(tǒng)含礫粘土層及石炭系砂頁巖概化為相對隔水層,構(gòu)成上部第四系松散巖類孔隙含水層與下部碳酸巖類巖溶裂隙含水層之間的隔水層。

    3.1.2 邊界條件的確定與概化

    根據(jù)研究區(qū)水文地質(zhì)條件及地下水系統(tǒng)特點,西側(cè)以東塢斷裂為界,概化為隔水邊界;北部灰?guī)r埋藏在石炭系—二疊系以下,最大揭露深度在400 m左右,以灰?guī)r頂板埋深400 m為界,概化為隔水邊界;北部小清河為較大的地表河流,概化為第四系孔隙水定水頭邊界;東側(cè)巨野河河流與該區(qū)域內(nèi)巖溶地下水有水力聯(lián)系,概化為河流邊界,水頭高度根據(jù)各河流多年月平均水位按月份分別賦值,模擬河流與含水層之間的側(cè)向補排關(guān)系;南側(cè)山體為地下水分水嶺,概化為隔水邊界。模型范圍總面積約338 km2。

    垂向上邊界根據(jù)鉆孔揭露含水層深度進行刻畫;下部為太古代泰山巖群變質(zhì)巖,概化為隔水邊界。地表高程等值線根據(jù)地形圖等高線進行刻畫。

    綜上,可將研究區(qū)概化為一個統(tǒng)一的非均質(zhì)、各向同性、三維非穩(wěn)定地下水流系統(tǒng),模型三維概化示意圖見圖1。

    圖1 水文地質(zhì)概化示意圖Fig.1 Generalization sketch of hydrogeology

    3.2 數(shù)學模型的建立及求解

    3.2.1 水流模型

    研究區(qū)地下水系統(tǒng)概化成非均質(zhì)、各向同性、三維非穩(wěn)定地下水流系統(tǒng),地下水模型數(shù)學模型為:

    式中:D為滲流區(qū)域;H為地下水水頭(m);Kx、Ky、Kz為x、y、z方向上的滲透系數(shù)(m/d);S為自由水面之下的含水層單位儲水系數(shù)(L/m);W為承壓含水層源匯項(L/d);p為潛水面的蒸發(fā)量和降水的補給量等(L/d);μ為潛水含水層重力給水度;H0(x,y,z)為初始流場水頭分布值(m);H1(x,y,z,t)為第一類邊界水頭分布值(m);q(x,y,z,t)為第二類邊界單位面積流量(m3/m2·d);Kn為邊界面外線方向的滲透系數(shù)(m/d);Г1、Г2為滲流區(qū)第一、二類邊界。

    3.2.2 模型建立

    運用基于有限差分法的Visual MODFLOW軟件包建立評價區(qū)的地下水流數(shù)值模型,經(jīng)參數(shù)識別與模型檢驗后,對研究區(qū)地下水流系統(tǒng)進行模擬分析,該模型作為地下水溶質(zhì)運移模擬的基礎(chǔ)。

    根據(jù)含水層結(jié)構(gòu)和地下水徑流特征,將研究區(qū)剖分為82行×63列的單元格,垂向為三層。

    3.2.3 源匯項處理

    (1) 補給項。模型中的補給項包括降雨入滲補給、河流及水庫滲漏補給等。

    ① 降雨入滲補給。在模型中降雨入滲量按面狀補給量處理。降雨入滲補給量計算公式為:

    Q=α×F×P

    式中:α、F、P分別為降水入滲系數(shù)、降水入滲面積和降水量。

    降水入滲面積取研究區(qū)整個面積,降水入滲系數(shù)按照參數(shù)分區(qū)分別賦值,降雨量按照水文站觀測資料以月為步長賦值。

    ② 河流、水庫滲漏補給。河流與含水層的交換量可以采用下式計算求得:

    式中:QR為河流滲漏量(m3/a);B為河床寬度(m);L為河流長度(m);KR為河底滲透系數(shù);hR為河流水位(m),按月平均水位賦值;h為地下水位(m);S為河底積層厚度(m);t為河流行水時間(d)。

    研究區(qū)內(nèi)水庫年滲漏量按多年實際運行記錄滲漏量平均值賦值,即狼貓山水庫補給量62.3萬m3/a,杜張水庫補給量44萬m3/a。

    (2) 排泄項。模型中的排泄項包括蒸發(fā)、人工開采等。

    ① 蒸發(fā)量。蒸發(fā)量按照區(qū)內(nèi)水文站實測蒸發(fā)量取每月平均值,蒸發(fā)極限埋深取經(jīng)驗值4 m,蒸發(fā)量通過軟件中Evapotranspira-tion模塊自動計算并模擬。

    ② 開采量。農(nóng)業(yè)灌溉開采為季節(jié)性開采,每年3—6月份集中采水灌溉。經(jīng)調(diào)查統(tǒng)計農(nóng)業(yè)灌溉開采地下水量為1 766 m3/d。

    飲用及工業(yè)用水開采主要為巖溶裂隙水,據(jù)統(tǒng)計平均開采量約3 086.46萬m3/a。

    3.2.4 等效水文地質(zhì)參數(shù)的初始值

    根據(jù)研究區(qū)水文地質(zhì)條件,將模型劃分為四個水文地質(zhì)參數(shù)分區(qū),所用的滲透系數(shù)、給水度、蒸發(fā)強度系數(shù)以及降雨入滲系數(shù)等參數(shù)初始值采用研究區(qū)內(nèi)抽(注)水試驗及相關(guān)水文地質(zhì)資料綜合確定。

    3.3 模型驗證與識別

    地下水數(shù)學模型的識別采用試估—校正法,識別的過程同時也是等效水文地質(zhì)參數(shù)的調(diào)整過程。通過不斷的參數(shù)調(diào)整和反演,直到模擬地下流場與實測地下流場及地下水位計算值與實測水位值達到較好的擬合效果,最后得到的參數(shù)即是該模型校正后的最終參數(shù)。研究區(qū)內(nèi)上、中、下游三個碳酸巖類巖溶裂隙含水層監(jiān)測井的地下水位動態(tài)擬合曲線見圖2。

    圖2 地下水位動態(tài)擬合曲線圖Fig.2 Dynamic fitted curve of groundwater level

    4 溶質(zhì)運移模型及預測

    地下水溶質(zhì)運移模型是在地下水流場數(shù)值模型的基礎(chǔ)上建立的。本次利用Visual MODFLOW中的MT3D模塊計算地下水中污染物質(zhì)的運移情況,分析污染物在地下水系統(tǒng)中的運移擴散規(guī)律。

    4.1 溶質(zhì)運移數(shù)學模型

    4.1.1 控制方程

    本次建立的地下水溶質(zhì)運移模型是在三維水流影響下的三維彌散模型,水流主方向和坐標軸重合,溶液密度不變,存在局部平衡吸附和一級不可逆動力反應,溶解相和吸附相速率相等,即λ1=λ2。在此前提下,溶質(zhì)運移的三維水動力彌散方程的數(shù)學模型為:

    4.1.2 初始條件

    C(x,y,z,t)=C0(x,y,z)(x,y,z)∈Ω,t=0

    式中:C0(x,y,z)為已知濃度分布;Ω為模型模擬區(qū)域。

    本次模擬污染源的概化有兩種方式:一是補給濃度邊界,二是注水井邊界??蓪⒀a給濃度邊界和注水井處的初始濃度定為C0,其余地方均為0 mg/L,具體為:

    4.1.3 定解條件

    本模型將含水層邊界分為兩類:第一類邊界條件為定水頭邊界,第二類邊界為定流量邊界。穿越邊界的彌散通量為0,可表述為:

    式中:Γ2為第二類邊界。

    4.2 溶質(zhì)運移預測

    正常工況下,建設(shè)項目運行期內(nèi)污水處理站、污水管線等各部分運行正常,并采取了正確的防滲保護措施,污水無法進入土壤和地下水,正常工況下不會對地下水環(huán)境造成影響。為分析非正常工況下建設(shè)項目可能對地下水環(huán)境的影響程度,選取具有代表性的可信事故進行分析預測。

    (1) 情景一:生產(chǎn)裝置的設(shè)備或管線由于連接處(如法蘭、焊縫)開裂或腐蝕磨損等原因,發(fā)生管道的跑、冒、滴、漏等污水泄漏,污水泄漏量較小但持續(xù)滲漏,污染隱患點選取在污水處理站處。經(jīng)模型預測分析,不同模擬時段COD及氨氮在巖溶裂隙含水層中污染超標范圍和最大運移距離見表1。

    (2) 情景二:污水處理站內(nèi)調(diào)節(jié)池底部防滲局部失效,造成污水持續(xù)泄漏,污染隱患點選取在污水處理站調(diào)節(jié)池處。經(jīng)模型預測分析,不同模擬時段COD及氨氮在巖溶裂隙含水層中污染超標范圍和最大運移距離見表2。

    表1 情景一COD及氨氮對地下水環(huán)境影響預測結(jié)果Table 1 Scenario 1:Predictions on effect of COD and ammonia nitrogen on groundwater environment

    表2 情景二COD及氨氮對地下水環(huán)境影響預測結(jié)果Table 1 Scenario 2:Predictions on effect of COD and ammonia nitrogen on groundwater environment

    (3) 情景三:發(fā)生突發(fā)性事故,造成污水處理站內(nèi)調(diào)節(jié)池大面積破壞,大量污水滲漏,污染源類型為瞬時污染,污染隱患點選取在污水處理站調(diào)節(jié)池處。經(jīng)模型預測分析,不同模擬時段COD及氨氮在巖溶裂隙含水層中污染超標范圍和最大運移距離見表3。

    (4) 情景四:當發(fā)生突發(fā)事故將廢水存放于事故池內(nèi)時,恰好事故池局部防滲失效,污水滲漏并進入地下水,儲存時間假定為1天,污染隱患點選取在生物原核原液車間事故池處。經(jīng)模型預測分析,不同模擬時段COD及氨氮在巖溶裂隙含水層中污染超標范圍和最大運移距離見表4。

    表3 情景三COD及氨氮對地下水環(huán)境影響預測結(jié)果Table 3 Scenario 3:Predictions on effect of COD and ammonia nitrogen on groundwater environment

    表4 情景四COD及氨氮對地下水環(huán)境影響預測結(jié)果Table 4 Scenario 4:Predictions on effect of COD and ammonia nitrogen on groundwater environment

    5 地下水環(huán)境影響分析

    在設(shè)置的非正常工況可信事故下,根據(jù)溶質(zhì)運移計算結(jié)果可以得到如下結(jié)論:

    (1) 各污染物運移具有明顯的各向異性的特征,運移主方向基本為西北方向,因此在該方向上加強地下水動態(tài)監(jiān)測,可及時掌握建設(shè)項目對地下水環(huán)境影響程度。

    (2) 從時間變化規(guī)律上看,各污染因子的初始最大污染濃度普遍較高,但隨著時間推移最大污染濃度會逐漸降低。

    (3) 從空間變化規(guī)律上看,各污染物濃度在其運移方向上逐步降低。

    (4) 在非正常工況和事故工況下,預測時間內(nèi)污染物影響最遠距離為7.5 km,可能會對該范圍內(nèi)自備井地下水水質(zhì)產(chǎn)生影響。

    因此,應對項目場區(qū)易產(chǎn)生污染的部位采取有效防滲措施,嚴格控制非正常工況的產(chǎn)生,避免對地下水水質(zhì)產(chǎn)生污染。

    [1] 山東省地礦工程勘察院,濟南市地質(zhì)礦產(chǎn)技術(shù)監(jiān)測中心.濟南市環(huán)境地質(zhì)調(diào)查評價[R].濟南:山東省地礦工程勘察院,1999.

    [2] 山東省地質(zhì)礦產(chǎn)局八○一水文地質(zhì)工程地質(zhì)大隊.山東省濟南市白泉—武家水源地供水水文地質(zhì)勘探報告[R].濟南:山東省地質(zhì)礦產(chǎn)局,1990.

    [3] 易立新,徐鶴.地下水數(shù)值模擬:GMS應用基礎(chǔ)與實例[M].北京:化學工業(yè)出版社,2009.

    [4] 孟慶斌,邢立亭,滕朝霞.濟南泉域“三水”轉(zhuǎn)化與泉水恢復關(guān)系研究[J].山東大學學報(工學版),2008,38(5):82-87.

    [5] 汪家權(quán),吳義鋒,錢家忠,等.濟南泉域巖溶地下水三維等參有限元數(shù)值模擬[J].煤田地質(zhì)與勘探,2005,33(3):39-41.

    (責任編輯:陳姣霞)

    Analysis on Effect of Bio-pharmaceutical Construction Project onGroundwater Environment in Karst Area

    LU Haiyu1,2, LI Wanglin2

    (1.ShandongSurveyandDesignInstituteofWaterConservancy,Shandong,Jinan250013; 2.SchoolofResourcesandEnvironment,CollegeofJinanUniversity,Shandong,Jinan250002)

    Based on equivalent porous medium model united with the method of conversely calculating the hydrogeology parameters,this paper establishes karst groundwater numerical model of unsteady flow after geological generalization. Through identifying and verifying the fitted dynamic curves,the resulting model shows a high fit degree with the groundwater flow field and groundwater level in simulation area. Then,this model is used to evaluate the effect of source intensities characteristics of the project on the extent of groundwater environment at model forecasting period. This not only provides scientific basis for groundwater pollution prevention,but also shows a referential significance on the environment evaluation of karst groundwater with similar hydrogeological conditions.

    karst; hydrogeology; equivalent porous medium; numerical model; environmental impact assessment

    2017-06-15;改回日期:2017-07-13

    陸海玉(1978-),男,工程師,地質(zhì)工程專業(yè),從事水文地質(zhì)、工程地質(zhì)勘察及地下水環(huán)境科學研究工作。E-mail:Luhaiyu201@163.com

    X824

    A

    1671-1211(2017)04-0480-06

    10.16536/j.cnki.issn.1671-1211.2017.04.027

    數(shù)字出版網(wǎng)址:http://www.cnki.net/kcms/detail/42.1736.X.20170620.1320.002.html 數(shù)字出版日期:2017-06-20 13:20

    猜你喜歡
    巖組水文地質(zhì)含水層
    全球多個含水層里的水正快速流失
    銀川市地下水賦存條件及動態(tài)特征
    基于抽水試驗計算水文地質(zhì)參數(shù)的分析與確定
    基于GPRS實現(xiàn)自動化水文地質(zhì)抽水試驗
    基于轉(zhuǎn)移概率的三維水文地質(zhì)結(jié)構(gòu)劃分
    水文地質(zhì)在工程地質(zhì)勘察中的有效應用
    淮南潘集深部勘查區(qū)15-2孔工程地質(zhì)巖組劃分
    美國西部奧加拉拉含水層水位下降原因初探
    白音華煤田三號露天礦區(qū)水文地質(zhì)條件分析
    全球地下含水層下降驚人:要被抽干了
    地理教學(2015年14期)2015-03-31 20:04:53
    成人影院久久| 日本猛色少妇xxxxx猛交久久| 在线观看人妻少妇| 久久影院123| 欧美日韩视频精品一区| 国产又色又爽无遮挡免| 欧美精品亚洲一区二区| 搡女人真爽免费视频火全软件| 人妻 亚洲 视频| 女的被弄到高潮叫床怎么办| 美国免费a级毛片| av免费观看日本| 欧美av亚洲av综合av国产av | 丝袜美腿诱惑在线| 少妇熟女欧美另类| 中文字幕人妻熟女乱码| 蜜桃在线观看..| 久久久久网色| 久久午夜综合久久蜜桃| 国产 一区精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产综合精华液| 日产精品乱码卡一卡2卡三| 国产一区有黄有色的免费视频| 丝袜脚勾引网站| 中文字幕色久视频| 午夜福利在线免费观看网站| 久久久久久人人人人人| 中文字幕精品免费在线观看视频| 多毛熟女@视频| 青草久久国产| videosex国产| 涩涩av久久男人的天堂| 欧美 亚洲 国产 日韩一| 亚洲综合色网址| 在线观看免费日韩欧美大片| 久久久久久久久久人人人人人人| 亚洲国产精品国产精品| 在线 av 中文字幕| 综合色丁香网| 国产一区二区在线观看av| 午夜福利在线观看免费完整高清在| h视频一区二区三区| 免费黄色在线免费观看| 久久青草综合色| 精品午夜福利在线看| 女人久久www免费人成看片| 久久久久久免费高清国产稀缺| 一区二区三区四区激情视频| 哪个播放器可以免费观看大片| 久久精品国产亚洲av高清一级| 成人午夜精彩视频在线观看| 久久 成人 亚洲| 国产极品天堂在线| 在线天堂中文资源库| 亚洲少妇的诱惑av| 人妻少妇偷人精品九色| 波多野结衣一区麻豆| 久久av网站| 国产成人aa在线观看| 夫妻午夜视频| 搡女人真爽免费视频火全软件| 色94色欧美一区二区| 久久久久久人妻| 黄色毛片三级朝国网站| 亚洲精品国产一区二区精华液| 一区在线观看完整版| 国产熟女欧美一区二区| 在线精品无人区一区二区三| 国产爽快片一区二区三区| 99久久精品国产国产毛片| 免费久久久久久久精品成人欧美视频| 亚洲av男天堂| av网站免费在线观看视频| 久久久久久久久久久免费av| 亚洲第一av免费看| 欧美 日韩 精品 国产| 97人妻天天添夜夜摸| 久久影院123| 七月丁香在线播放| 一区二区av电影网| 国产精品人妻久久久影院| 69精品国产乱码久久久| 成人免费观看视频高清| 午夜福利一区二区在线看| 97在线人人人人妻| 久久久久久久大尺度免费视频| 久久精品国产综合久久久| 久久人人爽av亚洲精品天堂| 国产一区亚洲一区在线观看| 哪个播放器可以免费观看大片| 丝袜人妻中文字幕| 亚洲av日韩在线播放| 日韩,欧美,国产一区二区三区| 一区二区日韩欧美中文字幕| 一本大道久久a久久精品| 日日撸夜夜添| 国产麻豆69| 日产精品乱码卡一卡2卡三| 午夜福利网站1000一区二区三区| av有码第一页| 深夜精品福利| 国产av一区二区精品久久| 亚洲欧美色中文字幕在线| 涩涩av久久男人的天堂| 高清在线视频一区二区三区| 少妇熟女欧美另类| 亚洲五月色婷婷综合| 最近的中文字幕免费完整| 啦啦啦视频在线资源免费观看| 在线观看国产h片| 777米奇影视久久| 桃花免费在线播放| av在线观看视频网站免费| 亚洲av日韩在线播放| 国产白丝娇喘喷水9色精品| 十八禁网站网址无遮挡| 午夜福利,免费看| 日韩在线高清观看一区二区三区| 国产色婷婷99| 青春草视频在线免费观看| 国产亚洲精品第一综合不卡| 精品亚洲乱码少妇综合久久| 国产日韩欧美视频二区| 性色avwww在线观看| 欧美在线黄色| 少妇被粗大的猛进出69影院| 色婷婷av一区二区三区视频| 久久精品人人爽人人爽视色| 亚洲男人天堂网一区| 久久精品aⅴ一区二区三区四区 | 国产一区二区激情短视频 | 成人午夜精彩视频在线观看| 久久久久久久亚洲中文字幕| 18禁动态无遮挡网站| 在线亚洲精品国产二区图片欧美| 国产一级毛片在线| 在线亚洲精品国产二区图片欧美| 成年av动漫网址| 久久人人97超碰香蕉20202| 天天躁夜夜躁狠狠躁躁| 天天躁夜夜躁狠狠躁躁| 自线自在国产av| 一区二区av电影网| 国产黄色免费在线视频| 日韩av在线免费看完整版不卡| 亚洲精品视频女| 男的添女的下面高潮视频| 亚洲av日韩在线播放| 大香蕉久久成人网| 国产一区二区三区av在线| 街头女战士在线观看网站| 久久99蜜桃精品久久| 街头女战士在线观看网站| 黄片无遮挡物在线观看| 国产精品国产三级专区第一集| 亚洲激情五月婷婷啪啪| 亚洲一码二码三码区别大吗| 色网站视频免费| 欧美最新免费一区二区三区| 午夜福利在线免费观看网站| 免费大片黄手机在线观看| 欧美老熟妇乱子伦牲交| 青春草视频在线免费观看| 精品国产一区二区三区四区第35| 宅男免费午夜| 精品人妻偷拍中文字幕| 黄片无遮挡物在线观看| 欧美精品一区二区大全| 国产精品av久久久久免费| 久久精品国产亚洲av涩爱| 蜜桃在线观看..| 免费在线观看完整版高清| 人妻一区二区av| 在线精品无人区一区二区三| 人人妻人人爽人人添夜夜欢视频| 伦精品一区二区三区| 日韩熟女老妇一区二区性免费视频| 精品国产露脸久久av麻豆| 黄频高清免费视频| av网站免费在线观看视频| av在线观看视频网站免费| 亚洲美女视频黄频| 2022亚洲国产成人精品| 欧美日韩视频精品一区| 国产一区二区在线观看av| 欧美精品国产亚洲| 成人影院久久| 啦啦啦在线免费观看视频4| 少妇猛男粗大的猛烈进出视频| 久久精品国产综合久久久| 久久国内精品自在自线图片| 一级a爱视频在线免费观看| 另类精品久久| 亚洲图色成人| 在线观看一区二区三区激情| 日韩在线高清观看一区二区三区| 国产在线视频一区二区| 精品国产露脸久久av麻豆| 久久久久久久久久久久大奶| 青草久久国产| 人人妻人人澡人人看| 亚洲国产欧美在线一区| 国产高清国产精品国产三级| 久久国产亚洲av麻豆专区| 99久久精品国产国产毛片| 国产乱来视频区| www.自偷自拍.com| 午夜影院在线不卡| 尾随美女入室| 精品福利永久在线观看| 久久精品国产自在天天线| 在线观看三级黄色| 看免费av毛片| 国产片特级美女逼逼视频| 婷婷色av中文字幕| 国产亚洲精品第一综合不卡| 国产成人午夜福利电影在线观看| 99re6热这里在线精品视频| 一级毛片黄色毛片免费观看视频| 日韩人妻精品一区2区三区| 国产熟女午夜一区二区三区| 精品久久蜜臀av无| 天堂8中文在线网| 亚洲精品国产av成人精品| 亚洲三区欧美一区| 极品少妇高潮喷水抽搐| 宅男免费午夜| 亚洲经典国产精华液单| 国产又爽黄色视频| 久久精品国产亚洲av涩爱| 亚洲人成77777在线视频| 成人毛片a级毛片在线播放| 少妇被粗大猛烈的视频| 国产男人的电影天堂91| 人人妻人人添人人爽欧美一区卜| 热re99久久国产66热| 九色亚洲精品在线播放| av线在线观看网站| 成年女人在线观看亚洲视频| 亚洲精品美女久久av网站| 免费观看在线日韩| 国产精品无大码| 丰满乱子伦码专区| 在线观看免费高清a一片| 亚洲欧美精品综合一区二区三区 | 国产日韩欧美在线精品| 黄片小视频在线播放| 精品国产露脸久久av麻豆| 免费高清在线观看视频在线观看| 亚洲美女视频黄频| 青春草亚洲视频在线观看| 不卡av一区二区三区| 女人被躁到高潮嗷嗷叫费观| 久久久久久久亚洲中文字幕| 国产一区二区三区综合在线观看| 永久免费av网站大全| 亚洲国产av新网站| 制服人妻中文乱码| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久精品古装| 天天影视国产精品| 久久鲁丝午夜福利片| 美女大奶头黄色视频| 成人国产av品久久久| 精品国产超薄肉色丝袜足j| 女人高潮潮喷娇喘18禁视频| 熟女av电影| 国产有黄有色有爽视频| 国产精品国产三级专区第一集| av国产精品久久久久影院| 亚洲精品第二区| 欧美变态另类bdsm刘玥| 国产成人精品在线电影| 亚洲精品,欧美精品| 欧美日韩亚洲国产一区二区在线观看 | 麻豆乱淫一区二区| 一个人免费看片子| 国产在线视频一区二区| 母亲3免费完整高清在线观看 | 中文字幕最新亚洲高清| 不卡av一区二区三区| 熟妇人妻不卡中文字幕| 人妻一区二区av| 国产av国产精品国产| 一边摸一边做爽爽视频免费| 亚洲欧美色中文字幕在线| 成人亚洲精品一区在线观看| 中文欧美无线码| www.熟女人妻精品国产| 热99国产精品久久久久久7| 九九爱精品视频在线观看| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美精品自产自拍| 国产成人91sexporn| 9色porny在线观看| 下体分泌物呈黄色| 一个人免费看片子| 欧美人与善性xxx| 国产精品 国内视频| 日韩大片免费观看网站| 美女福利国产在线| 美女视频免费永久观看网站| 国产精品香港三级国产av潘金莲 | 丝袜美足系列| 蜜桃国产av成人99| 免费看不卡的av| 亚洲国产精品一区二区三区在线| 日本wwww免费看| 日本爱情动作片www.在线观看| 黄色视频在线播放观看不卡| 免费黄频网站在线观看国产| 欧美最新免费一区二区三区| 欧美日韩精品成人综合77777| 天堂8中文在线网| 少妇熟女欧美另类| 一区二区日韩欧美中文字幕| 在线天堂最新版资源| 菩萨蛮人人尽说江南好唐韦庄| 一二三四中文在线观看免费高清| 亚洲三区欧美一区| 一区二区日韩欧美中文字幕| 大话2 男鬼变身卡| 亚洲国产欧美在线一区| 一本久久精品| 亚洲伊人色综图| 人人妻人人澡人人爽人人夜夜| 波多野结衣一区麻豆| 久久综合国产亚洲精品| 一边摸一边做爽爽视频免费| 久久久精品免费免费高清| 午夜福利在线观看免费完整高清在| 亚洲av男天堂| 精品卡一卡二卡四卡免费| 久久久久久久大尺度免费视频| 日韩欧美一区视频在线观看| 国产成人精品无人区| 在线观看三级黄色| 一本大道久久a久久精品| 久久久a久久爽久久v久久| 亚洲精华国产精华液的使用体验| 欧美精品一区二区大全| 中文字幕色久视频| 我要看黄色一级片免费的| 天天躁狠狠躁夜夜躁狠狠躁| 99国产综合亚洲精品| 国产成人免费观看mmmm| 嫩草影院入口| 男女午夜视频在线观看| 久久久久久久久久人人人人人人| 久久久精品国产亚洲av高清涩受| 啦啦啦视频在线资源免费观看| 亚洲精品在线美女| 亚洲精品一二三| 亚洲精品日韩在线中文字幕| 在线观看三级黄色| 老司机亚洲免费影院| 一级爰片在线观看| 一级毛片黄色毛片免费观看视频| 女人精品久久久久毛片| 欧美 亚洲 国产 日韩一| 欧美av亚洲av综合av国产av | 亚洲第一青青草原| 日韩av在线免费看完整版不卡| 97在线视频观看| 侵犯人妻中文字幕一二三四区| 日韩精品免费视频一区二区三区| 老鸭窝网址在线观看| 久久ye,这里只有精品| 一级a爱视频在线免费观看| 宅男免费午夜| 午夜av观看不卡| 免费观看无遮挡的男女| 999久久久国产精品视频| 亚洲精品美女久久久久99蜜臀 | 久久人妻熟女aⅴ| 久热久热在线精品观看| 校园人妻丝袜中文字幕| 汤姆久久久久久久影院中文字幕| 免费人妻精品一区二区三区视频| 国产熟女午夜一区二区三区| 黑人欧美特级aaaaaa片| 999精品在线视频| 免费大片黄手机在线观看| 色哟哟·www| 晚上一个人看的免费电影| 国产免费现黄频在线看| 一区二区av电影网| 国产精品久久久久久av不卡| 久久久久久久国产电影| 久久女婷五月综合色啪小说| 国产97色在线日韩免费| 欧美人与性动交α欧美软件| 亚洲国产毛片av蜜桃av| 在线天堂最新版资源| 大码成人一级视频| av电影中文网址| 美女午夜性视频免费| 日韩av不卡免费在线播放| 男女边吃奶边做爰视频| www日本在线高清视频| 国产av一区二区精品久久| 天天躁夜夜躁狠狠久久av| 高清av免费在线| 欧美精品亚洲一区二区| 另类亚洲欧美激情| 一级黄片播放器| 亚洲国产毛片av蜜桃av| 欧美xxⅹ黑人| 一本色道久久久久久精品综合| 热99国产精品久久久久久7| 精品视频人人做人人爽| 1024视频免费在线观看| 欧美日韩视频高清一区二区三区二| 麻豆乱淫一区二区| 亚洲精品美女久久av网站| 久久久久久伊人网av| 免费日韩欧美在线观看| 赤兔流量卡办理| 久久久久精品久久久久真实原创| av国产久精品久网站免费入址| 国产一区有黄有色的免费视频| 欧美xxⅹ黑人| 激情五月婷婷亚洲| 国产一级毛片在线| 女人久久www免费人成看片| 精品国产乱码久久久久久小说| 亚洲 欧美一区二区三区| 久久人人97超碰香蕉20202| 人体艺术视频欧美日本| 亚洲国产欧美在线一区| 五月天丁香电影| 亚洲国产精品999| 成人国产麻豆网| 亚洲精品日本国产第一区| h视频一区二区三区| 亚洲精品视频女| 国产成人av激情在线播放| 亚洲国产精品国产精品| 午夜福利网站1000一区二区三区| 一级毛片电影观看| 香蕉精品网在线| 丝袜美腿诱惑在线| 黑人猛操日本美女一级片| 婷婷色av中文字幕| 人体艺术视频欧美日本| 国产视频首页在线观看| 18禁观看日本| 欧美成人精品欧美一级黄| 免费少妇av软件| 丝袜人妻中文字幕| 美女福利国产在线| 国产精品麻豆人妻色哟哟久久| 婷婷色麻豆天堂久久| 美女福利国产在线| 亚洲三级黄色毛片| 如日韩欧美国产精品一区二区三区| 99国产精品免费福利视频| 大陆偷拍与自拍| 一级片'在线观看视频| 亚洲五月色婷婷综合| 九草在线视频观看| 午夜久久久在线观看| 中文字幕亚洲精品专区| 精品福利永久在线观看| 中文字幕av电影在线播放| 日韩电影二区| 免费播放大片免费观看视频在线观看| 99久久精品国产国产毛片| 伦精品一区二区三区| 大码成人一级视频| 老汉色∧v一级毛片| 寂寞人妻少妇视频99o| 91在线精品国自产拍蜜月| av在线老鸭窝| 国产 一区精品| 亚洲欧美中文字幕日韩二区| 国产一区二区激情短视频 | 国产成人精品久久二区二区91 | 成人手机av| 成年人午夜在线观看视频| 国产探花极品一区二区| 国产男女内射视频| 18禁观看日本| 最新的欧美精品一区二区| 又黄又粗又硬又大视频| 亚洲成av片中文字幕在线观看 | 亚洲国产欧美网| 亚洲欧美成人精品一区二区| 少妇人妻久久综合中文| 一级黄片播放器| 免费观看在线日韩| 黄色怎么调成土黄色| 欧美97在线视频| 少妇人妻 视频| 日韩一区二区三区影片| 亚洲av日韩在线播放| 久久韩国三级中文字幕| 久久亚洲国产成人精品v| 国产男人的电影天堂91| 在线亚洲精品国产二区图片欧美| 成人毛片a级毛片在线播放| 国产免费一区二区三区四区乱码| 国产不卡av网站在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美亚洲日本最大视频资源| 天天躁夜夜躁狠狠久久av| 最近手机中文字幕大全| 捣出白浆h1v1| 青青草视频在线视频观看| 亚洲,欧美精品.| 亚洲av福利一区| 天天躁夜夜躁狠狠躁躁| 日产精品乱码卡一卡2卡三| 一本大道久久a久久精品| 日韩,欧美,国产一区二区三区| 伦理电影免费视频| 午夜影院在线不卡| 人人妻人人澡人人看| 欧美成人午夜精品| 成人毛片60女人毛片免费| 国产免费一区二区三区四区乱码| 中文字幕人妻丝袜一区二区 | 日韩一卡2卡3卡4卡2021年| 欧美日韩亚洲高清精品| 2022亚洲国产成人精品| 中文乱码字字幕精品一区二区三区| 人体艺术视频欧美日本| 卡戴珊不雅视频在线播放| 国产在线视频一区二区| 丰满迷人的少妇在线观看| 美女午夜性视频免费| 成人国产av品久久久| 亚洲人成电影观看| 久久久久精品性色| 99国产综合亚洲精品| 国产精品久久久久久av不卡| 不卡视频在线观看欧美| 中文字幕最新亚洲高清| 午夜免费观看性视频| 亚洲av福利一区| 精品亚洲成a人片在线观看| 日韩精品有码人妻一区| 性色av一级| 久久99精品国语久久久| 日韩 亚洲 欧美在线| 亚洲第一区二区三区不卡| 一级毛片我不卡| 波野结衣二区三区在线| 啦啦啦在线免费观看视频4| 亚洲精品美女久久久久99蜜臀 | 午夜影院在线不卡| 亚洲av.av天堂| 十八禁高潮呻吟视频| 亚洲美女视频黄频| 精品99又大又爽又粗少妇毛片| 国产精品成人在线| 国产精品免费视频内射| 中文天堂在线官网| 亚洲一区中文字幕在线| 国产成人精品久久久久久| 中文字幕精品免费在线观看视频| 黄色 视频免费看| 如日韩欧美国产精品一区二区三区| 国产在视频线精品| 男女无遮挡免费网站观看| 26uuu在线亚洲综合色| 纯流量卡能插随身wifi吗| 午夜福利在线观看免费完整高清在| 国产高清国产精品国产三级| 女人高潮潮喷娇喘18禁视频| 1024视频免费在线观看| 亚洲婷婷狠狠爱综合网| 午夜福利,免费看| 欧美日韩亚洲国产一区二区在线观看 | 91成人精品电影| 日韩伦理黄色片| 国产一级毛片在线| 香蕉国产在线看| 美女视频免费永久观看网站| 只有这里有精品99| 人妻系列 视频| freevideosex欧美| 国产精品免费视频内射| 一级片免费观看大全| 国产精品嫩草影院av在线观看| 香蕉丝袜av| 好男人视频免费观看在线| 亚洲在久久综合| 黄片播放在线免费| 国产成人一区二区在线| 国产人伦9x9x在线观看 | 免费在线观看黄色视频的| 丝袜在线中文字幕| 中文欧美无线码| 亚洲成人一二三区av| 欧美日韩视频精品一区| 亚洲精品在线美女| 18在线观看网站| 啦啦啦啦在线视频资源| 两性夫妻黄色片| 麻豆乱淫一区二区| 亚洲精品自拍成人| 国产黄色视频一区二区在线观看| 亚洲精品自拍成人| 建设人人有责人人尽责人人享有的| 亚洲av欧美aⅴ国产| 老司机亚洲免费影院| 亚洲人成77777在线视频| 卡戴珊不雅视频在线播放| 九草在线视频观看| 国产av国产精品国产| 看免费成人av毛片| av国产久精品久网站免费入址| 国产男女超爽视频在线观看|