• <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
    午夜91福利影院| 日韩av免费高清视频| 视频区欧美日本亚洲| 中文字幕亚洲精品专区| 欧美日韩视频高清一区二区三区二| 在线观看免费日韩欧美大片| 精品人妻1区二区| 久久ye,这里只有精品| 国产极品粉嫩免费观看在线| 久久国产精品大桥未久av| 亚洲欧美一区二区三区国产| 男女高潮啪啪啪动态图| 国产一区二区三区av在线| 国产免费现黄频在线看| 国产成人系列免费观看| 97人妻天天添夜夜摸| 国产国语露脸激情在线看| 色精品久久人妻99蜜桃| 中文字幕人妻丝袜制服| 1024香蕉在线观看| 久久精品久久久久久久性| 欧美人与性动交α欧美软件| 最新在线观看一区二区三区 | 日韩大码丰满熟妇| 电影成人av| 欧美黑人精品巨大| 久久久久国产精品人妻一区二区| 丁香六月欧美| 久久久久国产一级毛片高清牌| 最新的欧美精品一区二区| 91精品三级在线观看| 飞空精品影院首页| 又黄又粗又硬又大视频| 人人妻人人添人人爽欧美一区卜| 日韩av免费高清视频| 午夜激情久久久久久久| 婷婷色综合大香蕉| 久久99一区二区三区| 天堂中文最新版在线下载| 性高湖久久久久久久久免费观看| 欧美黑人欧美精品刺激| 纯流量卡能插随身wifi吗| 亚洲国产最新在线播放| 熟女少妇亚洲综合色aaa.| 久久精品久久久久久噜噜老黄| 亚洲国产欧美一区二区综合| 免费高清在线观看日韩| 欧美激情 高清一区二区三区| 亚洲国产精品国产精品| 一二三四在线观看免费中文在| 水蜜桃什么品种好| 亚洲精品一卡2卡三卡4卡5卡 | 欧美日韩综合久久久久久| 亚洲精品一二三| 国产在线观看jvid| 精品亚洲成a人片在线观看| 国产深夜福利视频在线观看| 女人高潮潮喷娇喘18禁视频| 成年人午夜在线观看视频| 黄色a级毛片大全视频| e午夜精品久久久久久久| 国产黄色免费在线视频| 看十八女毛片水多多多| 赤兔流量卡办理| 新久久久久国产一级毛片| 午夜两性在线视频| 日韩免费高清中文字幕av| 色视频在线一区二区三区| 大片电影免费在线观看免费| av福利片在线| 久久久亚洲精品成人影院| 人妻一区二区av| 国产日韩欧美在线精品| 免费黄频网站在线观看国产| 亚洲精品一区蜜桃| 秋霞在线观看毛片| 高清欧美精品videossex| 国产一区二区 视频在线| √禁漫天堂资源中文www| 十分钟在线观看高清视频www| 18禁黄网站禁片午夜丰满| 欧美精品一区二区免费开放| 男女下面插进去视频免费观看| 最近中文字幕2019免费版| 亚洲成色77777| 色94色欧美一区二区| 国产不卡av网站在线观看| 一区福利在线观看| 亚洲第一青青草原| 亚洲第一av免费看| 一区二区日韩欧美中文字幕| 男人添女人高潮全过程视频| 日韩 欧美 亚洲 中文字幕| 久久久久久人人人人人| 久久久久国产精品人妻一区二区| 欧美日韩国产mv在线观看视频| 亚洲欧洲精品一区二区精品久久久| 亚洲精品乱久久久久久| 菩萨蛮人人尽说江南好唐韦庄| 国产女主播在线喷水免费视频网站| 精品国产乱码久久久久久男人| 国产精品久久久久久人妻精品电影 | 国产免费又黄又爽又色| 女人被躁到高潮嗷嗷叫费观| 国产视频一区二区在线看| 精品少妇一区二区三区视频日本电影| 下体分泌物呈黄色| 久久久久网色| a级毛片黄视频| 夜夜骑夜夜射夜夜干| 亚洲精品av麻豆狂野| netflix在线观看网站| 欧美精品啪啪一区二区三区 | 久久性视频一级片| 夫妻性生交免费视频一级片| 久久久久网色| 在线观看免费高清a一片| 午夜福利影视在线免费观看| 9热在线视频观看99| 国产视频一区二区在线看| 国产免费视频播放在线视频| 成人18禁高潮啪啪吃奶动态图| 精品欧美一区二区三区在线| 两人在一起打扑克的视频| 国产极品粉嫩免费观看在线| 国产97色在线日韩免费| 午夜老司机福利片| 欧美人与性动交α欧美精品济南到| 精品国产乱码久久久久久小说| 日韩大片免费观看网站| 国产精品国产三级专区第一集| 亚洲av欧美aⅴ国产| 欧美在线一区亚洲| 交换朋友夫妻互换小说| 亚洲黑人精品在线| 精品福利永久在线观看| 99国产精品一区二区蜜桃av | 大香蕉久久网| 叶爱在线成人免费视频播放| 色精品久久人妻99蜜桃| 久久人妻熟女aⅴ| 大话2 男鬼变身卡| 久久免费观看电影| 国产精品 国内视频| 久久精品久久久久久噜噜老黄| 亚洲综合色网址| 黄片播放在线免费| 不卡av一区二区三区| 中文字幕精品免费在线观看视频| 99久久综合免费| 亚洲成人手机| 99国产精品99久久久久| 国产av国产精品国产| 脱女人内裤的视频| 国产99久久九九免费精品| 日本av手机在线免费观看| 两个人看的免费小视频| 国产日韩一区二区三区精品不卡| av天堂久久9| 久久久久久久久免费视频了| 欧美日韩视频高清一区二区三区二| 久久人人爽av亚洲精品天堂| 黑人欧美特级aaaaaa片| 国产一区二区 视频在线| 黑人欧美特级aaaaaa片| 精品少妇久久久久久888优播| 国产99久久九九免费精品| 亚洲精品av麻豆狂野| 亚洲成人手机| a 毛片基地| 在线观看免费高清a一片| 咕卡用的链子| 亚洲欧美一区二区三区国产| 在线观看www视频免费| 国产真人三级小视频在线观看| 国产真人三级小视频在线观看| 午夜影院在线不卡| 美女午夜性视频免费| 国产精品国产av在线观看| 91精品三级在线观看| 精品久久久久久久毛片微露脸 | 久久久久久人人人人人| 国产在线观看jvid| 国产欧美日韩一区二区三 | 婷婷色av中文字幕| 久久久国产欧美日韩av| 天天躁狠狠躁夜夜躁狠狠躁| 无遮挡黄片免费观看| 99国产精品免费福利视频| 久久鲁丝午夜福利片| 国产av精品麻豆| 久久精品国产亚洲av高清一级| 亚洲美女黄色视频免费看| 亚洲精品国产一区二区精华液| 多毛熟女@视频| 亚洲国产欧美在线一区| 国产精品亚洲av一区麻豆| 成人三级做爰电影| 亚洲精品国产色婷婷电影| 丰满迷人的少妇在线观看| 人人妻人人添人人爽欧美一区卜| 国产成人免费无遮挡视频| 好男人视频免费观看在线| 亚洲精品国产av成人精品| a级毛片黄视频| 国产一级毛片在线| 女人高潮潮喷娇喘18禁视频| 日韩视频在线欧美| 午夜福利一区二区在线看| 国产不卡av网站在线观看| 日韩,欧美,国产一区二区三区| 国产免费现黄频在线看| 亚洲成av片中文字幕在线观看| 国产淫语在线视频| 久久国产精品男人的天堂亚洲| h视频一区二区三区| 黄片小视频在线播放| 天堂中文最新版在线下载| 麻豆国产av国片精品| 视频区图区小说| 熟女少妇亚洲综合色aaa.| 青青草视频在线视频观看| 精品久久久精品久久久| 乱人伦中国视频| 18禁国产床啪视频网站| 成人免费观看视频高清| 大型av网站在线播放| 男人操女人黄网站| 国产一区亚洲一区在线观看| 国产日韩一区二区三区精品不卡| 日韩免费高清中文字幕av| 精品少妇一区二区三区视频日本电影| 精品亚洲成a人片在线观看| 国产欧美日韩精品亚洲av| 亚洲精品日韩在线中文字幕| 又粗又硬又长又爽又黄的视频| 精品一区二区三区av网在线观看 | 日韩一区二区三区影片| 一边摸一边做爽爽视频免费| 十八禁高潮呻吟视频| 精品人妻熟女毛片av久久网站| 一级,二级,三级黄色视频| 国产亚洲欧美精品永久| 一级毛片我不卡| 老司机影院成人| 在线观看www视频免费| 韩国高清视频一区二区三区| 成人国产一区最新在线观看 | 亚洲av欧美aⅴ国产| 亚洲一卡2卡3卡4卡5卡精品中文| 国产免费视频播放在线视频| 午夜视频精品福利| 国产精品二区激情视频| 欧美亚洲 丝袜 人妻 在线| 在线av久久热| 黑人猛操日本美女一级片| 国产成人a∨麻豆精品| 精品福利观看| 亚洲av成人不卡在线观看播放网 | 国产欧美日韩一区二区三区在线| 久久精品熟女亚洲av麻豆精品| 久久精品人人爽人人爽视色| 亚洲人成电影观看| 国产男人的电影天堂91| 国精品久久久久久国模美| 亚洲五月色婷婷综合| videos熟女内射| 尾随美女入室| 国产伦人伦偷精品视频| 亚洲成国产人片在线观看| 亚洲av片天天在线观看| 日韩制服丝袜自拍偷拍| www.av在线官网国产| 色网站视频免费| 欧美成狂野欧美在线观看| 男女国产视频网站| 一区二区三区四区激情视频| 夜夜骑夜夜射夜夜干| 欧美av亚洲av综合av国产av| 人体艺术视频欧美日本| 亚洲av男天堂| 亚洲精品国产色婷婷电影| 国产成人av激情在线播放| 一本综合久久免费| 亚洲精品一卡2卡三卡4卡5卡 | 多毛熟女@视频| 王馨瑶露胸无遮挡在线观看| 一级黄色大片毛片| 女人高潮潮喷娇喘18禁视频| 尾随美女入室| cao死你这个sao货| 国产成人一区二区三区免费视频网站 | 亚洲一卡2卡3卡4卡5卡精品中文| 高清av免费在线| 好男人视频免费观看在线| 另类精品久久| 精品一品国产午夜福利视频| 夫妻性生交免费视频一级片| 大香蕉久久成人网| 大码成人一级视频| 热re99久久国产66热| 国产精品国产三级国产专区5o| 少妇精品久久久久久久| 国产一区二区激情短视频 | 国产成人精品久久二区二区91| 婷婷色综合大香蕉| 18禁裸乳无遮挡动漫免费视频| av国产久精品久网站免费入址| 亚洲av成人不卡在线观看播放网 | 看免费成人av毛片| 日本a在线网址| 久久九九热精品免费| 免费久久久久久久精品成人欧美视频| 美国免费a级毛片| 国产爽快片一区二区三区| 人人妻人人爽人人添夜夜欢视频| 1024视频免费在线观看| 美女扒开内裤让男人捅视频| 又黄又粗又硬又大视频| 成人国产一区最新在线观看 | 黄色一级大片看看| 久久性视频一级片| 亚洲,一卡二卡三卡| 乱人伦中国视频| 国产亚洲精品久久久久5区| 亚洲人成电影免费在线| 亚洲欧洲精品一区二区精品久久久| a级毛片黄视频| 成人午夜精彩视频在线观看| 成人亚洲欧美一区二区av| 成人黄色视频免费在线看| av国产精品久久久久影院| 免费看av在线观看网站| 精品亚洲成国产av| 亚洲视频免费观看视频| 国产欧美日韩精品亚洲av| 成人亚洲精品一区在线观看| 亚洲七黄色美女视频| 精品免费久久久久久久清纯 | 日本欧美视频一区| 黑丝袜美女国产一区| 99精品久久久久人妻精品| 国产精品一区二区在线不卡| 国产精品久久久久成人av| 成人国产一区最新在线观看 | 国产日韩一区二区三区精品不卡| 狂野欧美激情性bbbbbb| 国产真人三级小视频在线观看| 午夜福利在线免费观看网站| 久9热在线精品视频| 久久精品久久精品一区二区三区| 久热爱精品视频在线9| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲成色77777| 美女主播在线视频| 青青草视频在线视频观看| 久久久久久免费高清国产稀缺| 一本一本久久a久久精品综合妖精| 91国产中文字幕| 亚洲欧美日韩高清在线视频 | 大话2 男鬼变身卡| 中文字幕色久视频| 午夜91福利影院| 国产日韩欧美在线精品| 日韩av不卡免费在线播放| 精品第一国产精品| 丝袜在线中文字幕| 国产成人一区二区在线| 99精国产麻豆久久婷婷| 91麻豆av在线| 黄色毛片三级朝国网站| 欧美 日韩 精品 国产| 国产亚洲欧美在线一区二区| 久久精品久久久久久噜噜老黄| 欧美性长视频在线观看| 80岁老熟妇乱子伦牲交| 婷婷色综合大香蕉| 自拍欧美九色日韩亚洲蝌蚪91| 精品高清国产在线一区| 18禁观看日本| 看十八女毛片水多多多| 欧美日韩福利视频一区二区| 国精品久久久久久国模美| 嫩草影视91久久| 国产在线免费精品| 大香蕉久久成人网| 精品少妇久久久久久888优播| 乱人伦中国视频| 大话2 男鬼变身卡| 日韩av免费高清视频| 纯流量卡能插随身wifi吗| 成年人免费黄色播放视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美日韩视频高清一区二区三区二| 丰满饥渴人妻一区二区三| 国产一卡二卡三卡精品| 欧美精品一区二区免费开放| 黄色毛片三级朝国网站| 涩涩av久久男人的天堂| 大型av网站在线播放| 欧美精品av麻豆av| 国产高清视频在线播放一区 | 国产精品香港三级国产av潘金莲 | 好男人视频免费观看在线| av在线app专区| 日韩人妻精品一区2区三区| 高潮久久久久久久久久久不卡| 国产野战对白在线观看| 只有这里有精品99| 亚洲自偷自拍图片 自拍| 日韩视频在线欧美| 交换朋友夫妻互换小说| 人人澡人人妻人| 极品少妇高潮喷水抽搐| 欧美精品一区二区大全| 欧美日韩亚洲国产一区二区在线观看 | 90打野战视频偷拍视频| 亚洲欧美成人综合另类久久久| av国产精品久久久久影院| 午夜av观看不卡| 咕卡用的链子| h视频一区二区三区| 在线观看免费日韩欧美大片| 午夜福利影视在线免费观看| 午夜福利,免费看| 日韩一卡2卡3卡4卡2021年| 一级a爱视频在线免费观看| 国产三级黄色录像| 七月丁香在线播放| 国产精品99久久99久久久不卡| 人人妻人人爽人人添夜夜欢视频| 人体艺术视频欧美日本| 国产男女内射视频| 亚洲精品中文字幕在线视频| 叶爱在线成人免费视频播放| 国产精品欧美亚洲77777| 欧美成人精品欧美一级黄| 你懂的网址亚洲精品在线观看| 国产精品三级大全| 亚洲国产成人一精品久久久| 大话2 男鬼变身卡| 丁香六月欧美| 精品人妻在线不人妻| 久久精品久久久久久噜噜老黄| 午夜激情av网站| 90打野战视频偷拍视频| 一级毛片女人18水好多 | 99香蕉大伊视频| 丝袜脚勾引网站| 午夜福利在线免费观看网站| svipshipincom国产片| 久久人人97超碰香蕉20202| 欧美精品人与动牲交sv欧美| 国产1区2区3区精品| 国产在线观看jvid| 中文欧美无线码| 亚洲欧美中文字幕日韩二区| 精品久久久久久久毛片微露脸 | 又大又爽又粗| 欧美国产精品va在线观看不卡| 大香蕉久久网| 国产精品欧美亚洲77777| 中文字幕亚洲精品专区| 亚洲七黄色美女视频| 亚洲情色 制服丝袜| 男女午夜视频在线观看| 又大又爽又粗| 丝袜美足系列| 成人午夜精彩视频在线观看| 99精品久久久久人妻精品| 免费日韩欧美在线观看| 中文字幕人妻丝袜一区二区| 最新在线观看一区二区三区 | 亚洲精品久久成人aⅴ小说| 两性夫妻黄色片| av不卡在线播放| 免费高清在线观看日韩| 黄色a级毛片大全视频| 搡老岳熟女国产| 在线观看免费午夜福利视频| 亚洲专区国产一区二区| 亚洲九九香蕉| 久久人人爽av亚洲精品天堂| 午夜影院在线不卡| 婷婷色综合大香蕉| 中文乱码字字幕精品一区二区三区| 黄色视频不卡| 夫妻午夜视频| 十八禁高潮呻吟视频| 久久九九热精品免费| 五月天丁香电影| 国产免费现黄频在线看| 国产精品一区二区在线不卡| 无遮挡黄片免费观看| 欧美日韩福利视频一区二区| 亚洲精品第二区| 你懂的网址亚洲精品在线观看| 日日摸夜夜添夜夜爱| 亚洲av美国av| 国产亚洲午夜精品一区二区久久| av在线老鸭窝| 精品一区二区三区四区五区乱码 | 欧美精品高潮呻吟av久久| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 欧美xxⅹ黑人| 久久99热这里只频精品6学生| 国产精品香港三级国产av潘金莲 | 国产一级毛片在线| 只有这里有精品99| 日本欧美国产在线视频| 最近中文字幕2019免费版| 少妇猛男粗大的猛烈进出视频| 熟女av电影| 国产男人的电影天堂91| 狠狠精品人妻久久久久久综合| av网站免费在线观看视频| 欧美日韩国产mv在线观看视频| 亚洲,欧美,日韩| 国产男女内射视频| av天堂在线播放| 后天国语完整版免费观看| 成年人午夜在线观看视频| 波野结衣二区三区在线| 久久热在线av| 久久久欧美国产精品| 色婷婷av一区二区三区视频| 国产老妇伦熟女老妇高清| 成年人免费黄色播放视频| 亚洲国产精品一区三区| 国产一卡二卡三卡精品| 国产精品 欧美亚洲| avwww免费| 久久久久视频综合| 每晚都被弄得嗷嗷叫到高潮| 国产黄色免费在线视频| 亚洲国产精品成人久久小说| 国产熟女欧美一区二区| 午夜福利视频精品| 成年人免费黄色播放视频| 人人澡人人妻人| 熟女av电影| 欧美精品av麻豆av| 欧美黑人精品巨大| 又粗又硬又长又爽又黄的视频| 91成人精品电影| 一区二区av电影网| 老司机靠b影院| 嫁个100分男人电影在线观看 | 成人亚洲精品一区在线观看| 久久国产亚洲av麻豆专区| 国产精品九九99| 91精品伊人久久大香线蕉| 成人免费观看视频高清| 真人做人爱边吃奶动态| 女人久久www免费人成看片| 久久精品久久精品一区二区三区| 亚洲,欧美精品.| 亚洲国产中文字幕在线视频| 日韩电影二区| 夫妻午夜视频| 国产又色又爽无遮挡免| 又粗又硬又长又爽又黄的视频| 一个人免费看片子| 七月丁香在线播放| 无遮挡黄片免费观看| 久久久久视频综合| 国产日韩欧美在线精品| 欧美日韩亚洲国产一区二区在线观看 | 如日韩欧美国产精品一区二区三区| 亚洲精品成人av观看孕妇| 男的添女的下面高潮视频| 在线 av 中文字幕| 亚洲欧美精品自产自拍| 欧美人与性动交α欧美软件| 人体艺术视频欧美日本| 99九九在线精品视频| 国产三级黄色录像| 免费一级毛片在线播放高清视频 | 美女高潮到喷水免费观看| 亚洲精品久久午夜乱码| 91精品三级在线观看| 精品一区二区三区av网在线观看 | 青草久久国产| 久久久精品区二区三区| 天堂中文最新版在线下载| 丰满少妇做爰视频| 成在线人永久免费视频| 咕卡用的链子| 久久人妻福利社区极品人妻图片 | 国产欧美亚洲国产| 亚洲精品av麻豆狂野| 女性生殖器流出的白浆| 91精品伊人久久大香线蕉| 日韩熟女老妇一区二区性免费视频| 国产成人免费无遮挡视频| 国产精品.久久久| 热re99久久国产66热| 天天添夜夜摸| 91九色精品人成在线观看| 国产有黄有色有爽视频| 夜夜骑夜夜射夜夜干| 成年美女黄网站色视频大全免费| 亚洲欧美中文字幕日韩二区| 国产黄频视频在线观看| www.自偷自拍.com| 狂野欧美激情性xxxx| 久久精品熟女亚洲av麻豆精品| 日本猛色少妇xxxxx猛交久久| 天天躁夜夜躁狠狠躁躁| 亚洲激情五月婷婷啪啪| 丰满饥渴人妻一区二区三| 国产黄色免费在线视频| 色婷婷久久久亚洲欧美| 中文字幕人妻丝袜制服|