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

    北京基巖水開采對地面沉降的影響

    2017-06-09 13:14:33董佩李陽劉殷孫穎李志萍
    南水北調(diào)與水利科技 2017年3期
    關(guān)鍵詞:小湯山開采量基巖

    董佩+李陽+劉殷+孫穎+李志萍

    摘要:過量開采地下水造成的地面沉降問題已經(jīng)是北京市平原區(qū)主要的地質(zhì)災(zāi)害之一。論文選擇小湯山和后沙峪地區(qū)為典型研究區(qū),通過分析研究區(qū)基巖水開發(fā)利用狀況,開展基巖水開采引起地面沉降的研究。運用簡化的彈性模型評估基巖水開采造成的彈性變形量,結(jié)果表明,從1974年截至2010年,小湯山地區(qū)地熱田總凈開采量約9 400萬m3,引起的地面彈性變形量約為1.02 cm;小湯山地區(qū)巖溶水年均開采3 000萬m3,年均形變量約0.32 cm。后沙峪地區(qū)通過估算獲得近似值,1961年-2011年期間由于第四系地下水的開采造成的基巖變形量約為1.09 cm。

    關(guān)鍵詞:地面沉降;基巖水;小湯山;地熱田;彈性模型

    中圖分類號:P641.8 文獻標識碼:A 文章編號:1672-1683(2017)03-0120-06

    Abstract:Land subsidence induced by excessive groundwater withdrawal has constituted a major environmental geological problem in Beijing′s plain areas.This paper selected Xiaotangshan and Houshayu areas as typical research areas.By analyzing the development and utilization of bedrock water,we carried out research on the land subsidence induced by bedrock water exploitation.A simplified elastic model was used to assess the amount of elastic deformation caused by bedrock water exploitation.The results showed that from 1974 to 2010,the total exploitation of Xiaotangshan geothermal field was about 94 million cubic meters and the exploitation led to about 1.02 cm elastic deformation;the exploitation of Xiaotangshan karst water reached 30 million cubic meters annually,causing an average annual deformation of about 0.32 cm.Approximate values of the Houshayu area was obtained through estimation,and the accumulated elastic deformation caused by Quaternary groundwater exploitation was about 1.09 cm in the period from 1961 to 2011.

    Key words:land subsidence;bedrock water;Xiaotangshan;geothermal field;elastic model

    地面沉降是在自然和人為因素的作用下,由于地表松散土體的壓縮而導(dǎo)致區(qū)域性地面高程下降的一種環(huán)境地質(zhì)現(xiàn)象,是孔隙水承擔的孔隙水壓力和土骨架承擔的有效應(yīng)力發(fā)生變化的結(jié)果。隨著全球氣候變暖以及城市化進程的加速,人為因素在地面沉降中的作用日益凸顯,特別是大規(guī)模持續(xù)地開發(fā)利用地下水導(dǎo)致區(qū)域性的地面沉降迅速發(fā)展,成為地面沉降的主要影響因素[1-3]。

    目前地下水在北京城市供水中占2/3強。地下水的高強度開采,支撐著北京市的繁榮和穩(wěn)定,而其帶來的地面沉降問題日趨嚴重。2009年-2012年,北京年度平均降水分別為484.96 mm、497.10 mm、640.19 mm和653.20 mm,呈逐年增多的態(tài)勢,但是大范圍的地面沉降卻沒有出現(xiàn)預(yù)期的減緩。降水量的增加對淺部含水層(尤其是潛水)的補給是明顯的,在淺層地面沉降量較大的地區(qū),當淺部地層水位上升時,淺部地層沉降的減緩必然帶動總沉降量得以減緩。但應(yīng)該看到,這種減緩的幅度是有限的,并沒有改變深部地層地下水持續(xù)下降的趨勢。這既和深部地下水的補給困難有關(guān)系,更主要的是由于目前地下水開采逐步加深,引起深部地層沉降加速。

    本文將京北小湯山和后沙峪地區(qū)作為典型研究區(qū),受地下水開采影響,研究區(qū)地下水水位逐年下降,地下水水水位變化引起了地層的彈性、黏彈性變形,并受開采情況和地層巖性、厚度等因素影響,造成水平上的差異性地面沉降。而在垂向上,該區(qū)域內(nèi)深層地下水的開采量較大,以往對于深層水,特別是基巖水開采引起含水層變形造成地面沉降的監(jiān)測和研究工作都開展較少。本次工作嘗試通過分析研究區(qū)基巖水開發(fā)利用狀況,開展基巖水開采引起地層變形的研究。

    1 研究區(qū)自然條件

    研究區(qū)位于北京市北部平原地區(qū),包括小湯山和后沙峪地區(qū)兩個區(qū)域,行政區(qū)劃上橫跨順義區(qū)、昌平區(qū),西南邊界平原區(qū)以南口-孫河斷裂為界,山區(qū)以四橋子西梁-居庸關(guān)為界,北部邊界為石泉門-松樹溝-沙嶺-大羊山地表分水嶺,東南以順義斷裂為界(圖1),屬于潮白河沖洪積扇和溫榆河沖洪積扇交匯地帶,地勢特征大致為西北高,東南低。地貌類型主要為沖洪積平原,平原以昌平北為頂點,向南呈扇形展開。

    研究區(qū)屬于溫帶大陸性季風氣候,多年平均氣溫約12 ℃,昌平氣象站監(jiān)測1961年-2014年多年平均降水量530.9 mm。降水量多集中于5月-10月,占全年總降水量的95.87%。1979年-2013年多年平均蒸發(fā)量為1 691.64 mm,蒸發(fā)量主要集中在3月-6月,約占全年的50%以上。

    1.1 地質(zhì)條件

    研究區(qū)北部為北山隆起,東南部為順義凹陷,西部為京西北隆起,地質(zhì)構(gòu)造上表現(xiàn)為扭動構(gòu)造行跡。斷裂主要有北東向和北西向兩種構(gòu)造體系,北東向構(gòu)造包括黃莊-高麗營斷裂、良鄉(xiāng)-前門斷裂等;北西向構(gòu)造則主要是南口-孫河斷裂(圖1)。其中,北東向黃莊-高麗營斷裂和北西向南口-孫河斷裂具有較強的活動性,對新生界沉積厚度控制作用明顯,是研究區(qū)的最主要斷裂。集中在研究區(qū)中部的諸多斷裂一方面是深部地下水的通道,同時形成的斷裂破碎帶也提供了良好的地下水賦存環(huán)境。

    區(qū)內(nèi)已揭露的地層自上而下有:第四系、侏羅系,石炭系、奧陶系、寒武系,青白口系和薊縣系霧迷山組。此外,在小湯山鎮(zhèn)北側(cè)阿蘇衛(wèi)-小湯山一線分布有燕山期形成的酸性花崗巖[4-9]。

    1.2 水文地質(zhì)條件

    研究區(qū)地下水類型主要有兩大類:第四系松散巖類孔隙潛水和基巖水。第四系孔隙水主要分布于各大河流河道中,主要含水層為中粗砂、含砂黏性土、碎石土,屬沖洪積物,地下水豐富,滲透性較好;基巖水(微承壓水)大部分位于穩(wěn)定的基巖中,主要類型為基巖裂隙水和巖溶水。主要巖溶含水巖組有長城系、薊縣系、寒武系和奧陶系巖溶含水巖組。山前基巖地下水補給以大氣降水入滲補給為主,另外當十三陵水庫水位達到一定高度時,庫水滲漏補給巖溶地下水,平原區(qū)主要補給來源為側(cè)向徑流補給。

    2 地下水開發(fā)利用及地下水位監(jiān)測

    2.1 地下水開發(fā)利用概況

    目前,研究區(qū)地下水開采包括引用水源地的開采和地熱井的開采,平原區(qū)地下水開采深度主要集中在300 m以上。其中,埋深80~100 m的地下水以農(nóng)業(yè)開采為主,埋深150~180 m的地下水開采主要用于生活及工業(yè)用水,埋深200~300 m的地下水主要是生活用水。主要巖溶水水源地有鄧莊、化莊、真順、興壽等集中水源地,巖溶水年開采量達2 385萬 m3;另外研究區(qū)有豐富的地熱資源,地熱田開采利用自1974年,目前已有40年的開采歷史,小湯山更有“溫泉之鄉(xiāng)”的美譽。目前研究區(qū)有熱水開采井70多眼,年開采量110萬 m3。主要開采三個地層,分別為薊縣系霧迷山組、薊縣系鐵嶺組及寒武系地層,已有回灌井11眼,實施回灌地層主要為薊縣系霧迷山組和薊縣系鐵嶺組[10-13]。

    目前,為滿足日益增長的生活和生長需要,地下水正處于集中和超量開采狀態(tài)。由于工作區(qū)位于北京市地下水資源相對豐富的地段,區(qū)內(nèi)建設(shè)了多座水廠和用水量大的企業(yè),且隨著地熱溫泉優(yōu)勢越來越突出,一些單位的地熱用水規(guī)模也迅速擴大,在沒有其它新增替代水源的情況下,地下水的超量開采將隨著經(jīng)濟的快速發(fā)展而加劇。2000年以后,由于連續(xù)多年干旱,造成地下水位持續(xù)下降,出水量降低,同時地表水和淺層地下水的污染導(dǎo)致地下水開采越來越深,為了滿足生活用水的需求,300 m以上深層承壓水水位下降明顯,工作區(qū)內(nèi)不少地區(qū)承壓水水位下降速度都超過2 m/a,最大下降速率接近3 m/a,見圖2。

    地下水水位的持續(xù)下降,造成了工作區(qū)地下水位下降漏斗擴展、單井涌水量減少、局部地區(qū)含水層發(fā)生季節(jié)性疏干,致使包氣帶變厚,土壤含水量減少。

    2.2 基巖地下水水位監(jiān)測

    研究區(qū)現(xiàn)有地熱井長期水位動態(tài)監(jiān)測點2個,分別為湯熱灌-1井和苗圃觀測井,分別觀測霧迷山組和鐵嶺組水位變化。

    霧迷山組壓力的多年變化情況通過湯熱灌-1井的月最低水位動態(tài)曲線可以進行分析(圖3)。

    薊縣系霧迷山組地層的開采量、回灌量和水位變化整體呈現(xiàn)以下幾個階段(圖3)的特征:20世紀80年代初至20世紀90年代中期,開采量從不到100萬 m3/a逐年增加至200萬 m3/a,熱儲水位迅速下降,平均降幅約2.08 m/a,這一階段稱為快速開采期;20世紀90年代中期至2000年回灌之前,霧迷山組開采量基本保持在200~250萬 m3/a,水位下降速度變緩,平均降幅約1.29 m/a,這個階段稱為穩(wěn)定開采期;2001年-2007年,霧迷山組開采量逐漸減少到200萬 m3/a以內(nèi),同時開始接受回灌補給,年均回灌量從最初的不到10萬 m3/a,迅速增加到120萬 m3/a左右,水位下降速度大大減緩,并且2005年-2007年開始水位出現(xiàn)小幅抬升,這一個階段稱為開采恢復(fù)期;2008年以后,由于霧迷山組開采量又迅速增大至250~350萬 m3/a,回灌量僅維持在130~140萬 m3/a,水位出現(xiàn)小幅下降[14]。

    鐵嶺組地層同樣是從1974年開始進行開采,開采規(guī)模逐年增加,到1984年僅為77.25萬 m3/a,到1990年達到200.85萬 m3/a,之后開采量逐漸得到控制,20世紀90年代年開采量穩(wěn)定在100~200萬m3/a。21世紀后,隨著一些國企改制,很多老井停止開采,開采量迅速減小到100萬 m3/a以內(nèi),之后開采量一直降低,到2005年更是減小到23.17萬m3/a(圖4)。

    鐵嶺組長期水位監(jiān)測點是從2009年才開始實施監(jiān)測的,由于缺少多年長期監(jiān)測資料,目前無法對該地層的壓力變化情況進行詳細分析。2010年鐵嶺組年平均水位比2009年下降2.47 m。

    寒武系地層從1982年進行開采,長期以來開采規(guī)模比較穩(wěn)定,除了1998年開采量達到73.58萬 m3/a,其余年份基本都在60萬 m3/a以下(圖5)。由于該地層沒有水位長期監(jiān)測點,無法對該地層的壓力變化情況進行分析。

    綜上所述,由于長期大規(guī)模開發(fā),地下水資源開始枯竭,表現(xiàn)為井水位大幅度下降。在1985年-1999年間,霧迷山組地層水位由3.44 m降到30.65 m,降幅超過27 m;鐵嶺組地層水位由4.04 m降到31.90 m,降幅近28 m,寒武系地層水位由4.34 m降到24.24 m,降幅近20 m[15]。

    3 基巖形變的模型與計算

    3.1 研究區(qū)地面沉降狀況

    從圖6可以看出,研究區(qū)內(nèi)地面沉降量從山區(qū)至平原逐漸增加,且地面沉降量較大的區(qū)域大部處于昌平區(qū)和順義區(qū)的城區(qū)地段,尤以昌平八仙莊附近最為嚴重,根據(jù)《北京市地面沉降調(diào)查報告》昌平沙河-八仙莊沉降區(qū)最大累計沉降量已達到1 233 mm。2006年-2009年,昌平沙河-八仙莊地區(qū)沉降區(qū)的中心沉降速率約為75.14 mm/a。據(jù)北京市水文隊2012年度地面沉降報告分析,研究區(qū)是北京市平原區(qū)北部的主要沉降區(qū),其中八仙莊沉降量同期增加量大于15 mm,往年此地的沉降量已經(jīng)較大,在沉降量持續(xù)增大后,年度沉降量均超過100 mm,成為平原區(qū)北部的主要沉降中心。此外,在順義的西馬各莊地區(qū),近年沉降量有較大幅度增加,且在持續(xù)增加中,增幅超過20 mm,成為平原區(qū)北部2012年沉降量增加最快的地區(qū)之一,最大年沉降量已經(jīng)超過60 mm。此外研究區(qū)內(nèi)的北七家地區(qū),多年來一直延續(xù)著快速沉降的勢頭,其累計沉降量仍然是昌平沙河-八仙莊沉降區(qū)乃至全北京平原區(qū)的最大值。

    3.2 模型與實例

    Thomas L.Holzer(1979)在分析美國加利福尼亞州等地的地面沉降現(xiàn)象時發(fā)現(xiàn),長時間超采地下水造成含水層的彈性變形可能是造成加州地區(qū)地面沉降的原因之一[16]。文章指出,在基巖水開采條件下,由于承壓水的減少,引起含水層中孔隙水壓力的降低,含水層的應(yīng)力增加,使得含水層的孔隙水體積縮小,壓縮量應(yīng)等于相應(yīng)的開采量,從而引起了地面

    沉降的發(fā)生。因此對于地面沉降量已知或者可測的地區(qū),根據(jù)地面沉降量可獲得相應(yīng)的含水層開采量。承壓含水層開采量可由水位、給水度與開采區(qū)面積的乘積獲得。Thomas L.Holzer(1979)使用Timoshenko和Goodier的彈性力學理論對研究區(qū)由于地下水開采造成的地面沉降量進行計算,即對于均質(zhì)、各向同性的半無限圓形區(qū)域,當其表面存在均布載荷時引起的最大的表面變形可由下面的簡化模型獲得[17]:

    就基底變形而言,泊松比可取常量0.25[18],楊氏模量為0.7 Mbar[19],則從式(2)可知,研究區(qū)基巖彈性變形量的大小受開采量與區(qū)域面積的影響。

    研究區(qū)地面沉降的發(fā)生以地下水的長期持續(xù)過量開采為主要誘因。抽水引起的地面沉降的機理可以用有效應(yīng)力原理來解釋。在孔隙承壓含水層中抽取地下水引起承壓水水頭降低,必然要使含水層本身及其上、下相對隔水層中的孔隙水壓力隨之減小。在含水層上覆地層壓力不變的情況下,孔隙水壓力的減小導(dǎo)致土中有效應(yīng)力等量增大,結(jié)果引起孔隙體積減小,從而使土層壓縮。

    近年,隨著地下水的開發(fā)利用,導(dǎo)致了基巖水位逐年下降,相應(yīng)基巖水水位變化會引起基巖的彈性變形。嘗試利用Thomas L.Holzer(1979)的簡化模型,根據(jù)工作區(qū)地下水開發(fā)利用情況評估可能引起的彈性變形。

    從小湯山霧迷山組凈開采量與地層最大變形量關(guān)系圖(圖7)可以看出:從20世紀80年代初至90年代中期,隨著年凈開采量從不到100萬m3逐年增加至200萬m3,霧迷山組地層的年最大變形量也相應(yīng)的從0.008 9 mm增大到0.021 7 mm,平均增幅約0.001 mm/a;90年代中期至2000年回灌之前,霧迷山組凈開采量基本保持在200~260萬m3/a,地層形變速度變緩,平均增幅約0.000 78 mm/a;隨著回灌的進行,2001年-2007年,霧迷山組開采量逐漸減少到200萬m3/a以內(nèi),同時開始接受回灌補給,年均回灌量從最初的不到10萬 m3,迅速增加到120萬m3左右,地層形變增加速度大大減緩,并且2003年-2007年開始年最大變形量出現(xiàn)直線下降,從0.021 9 mm/a降到0.013 5 mm/a;2006年以后,由于霧迷山組年開采量又迅速增大至250~350萬m3,年回灌量僅維持在130~140萬 m3,地層形變又出現(xiàn)了回升。

    鐵嶺組地層從1974年開始進行開采,開采規(guī)模逐年增加,從1985年101.54萬 m3,到1990年達到205.11萬 m3,相應(yīng)的鐵嶺組地層年最大變形量從0.011 mm逐漸增加到0.022 3 mm,之后開采量逐漸得到控制,20世紀90年代年開采量穩(wěn)定在100~200萬 m3。21世紀后,隨著一些國企改制,很多老井停止開采,年開采量迅速減小到100萬 m3以內(nèi),從1990年-2004年,鐵嶺組地層年開采量從205.11萬 m3直降到48.04萬 m3,隨之而來的是鐵嶺組地層的最大彈性形變量從0.022 3 mm/a減小到0.005 2 mm/a;2004年之后年開采量一直降低,2005年-2008年年開采量更是小于回灌量,凈開采量出現(xiàn)了負值,地層得到一定的回彈(圖8)。

    從1974年截至2010年,小湯山地區(qū)地熱田總凈開采量約9 400萬m3,用簡化模型計算得在此期間由于地熱水的開采所造成的地面沉降量的貢獻約為1.02 cm;小湯山地區(qū)巖溶水年均開采3 000萬m3,年均形變量約0.32 cm/a。

    后沙峪地區(qū)目前沒有確切的開采量數(shù)據(jù),只能通過估算獲得近似值。1961年-2011年全北京平原區(qū)第四系地下水貯存損失量為104.3億m3,平原區(qū)面積6 400 km2,后沙峪面積240 km2,為平原區(qū)面積的1/29;后沙峪凹陷寬11 000 m,半徑5 500 m,由此可估算出后沙峪地區(qū)1961年-2011年由于第四系地下水的開采造成的基巖變形量總量約為1.09 cm。

    從計算結(jié)果看來,開采量較小的地熱水造成的地面沉降量遠遠小于開采量較大的巖溶水造成的地面沉降。本文使用簡化模型實際上是計算基巖的形變,即由于地熱水的開采使得荷載(重力)變化引起的形變;InSAR和水準獲取的測量結(jié)果都為地表總的形變量,是不同形變驅(qū)動因素累積形變之和,如第四系沉積地層形變和本文計算的基巖形變。在北京,抽水引起第四系沉積地層形變繼而誘發(fā)地面沉降是主因,但是沉降的產(chǎn)生還有其它驅(qū)動力。其它驅(qū)動力類型,所占沉降貢獻多少,即是本文的研究目的。此外,局部地區(qū)(地塊)的微震,除了受區(qū)域構(gòu)應(yīng)力影響外,地下水開采也是引起微震的一個重要因素。

    特別需要指出的是,深層承壓水的補給周期相對淺層地下水而言更長,有些甚至很難得到恢復(fù),研究區(qū)深層承壓水水頭持續(xù)下降的趨勢難以逆轉(zhuǎn),尤其是在北京水資源現(xiàn)狀條件下,深層承壓水超采現(xiàn)象短期內(nèi)難以得到根本改變。根據(jù)研究區(qū)附近的天竺地面沉降監(jiān)測站的分層監(jiān)測資料顯示,深部地層壓縮量占地面沉降的比例越來越大,這種現(xiàn)象與深層承壓水的過度開采且不易恢復(fù)密切相關(guān)。深層承壓水水頭的大幅下降使得深部地層以較快的速度持續(xù)壓縮,這種壓縮基本是不可逆轉(zhuǎn)的,且通常情況下由于土層變形的塑性和蠕變性,即使深層承壓水水頭停止下降,深部土層的壓縮還會持續(xù)相當長的時間。

    4 結(jié)論和建議

    4.1 結(jié)論

    論文在收集整理以往研究資料和成果報告的基礎(chǔ)上,結(jié)合地下水監(jiān)測與地面沉降監(jiān)測的實際情況,通過調(diào)查研究區(qū)(小湯山地區(qū)和后沙峪地區(qū))基巖水的開發(fā)利用現(xiàn)狀、收集整理開采條件下地下水水位相應(yīng)特征資料,重點分析由于基巖水的開發(fā)可能導(dǎo)致的形變貢獻,得出以下結(jié)論。

    (1)由于工作區(qū)地下水的開發(fā)利用,導(dǎo)致了地下水位逐年下降,相應(yīng)地下水水位變化引起地層的彈性變形。

    (2)運用簡化的彈性模型評估基巖水開采造成的彈性變形量,選取小湯山地熱田為典型研究區(qū),評估結(jié)果表明從1974年截至2010年,小湯山地區(qū)地熱田總凈開采量約9 400萬 m3,計算得引起的地面彈性變形量約為1.02 cm;小湯山地區(qū)巖溶水年均開采3 000萬 m3,年均形變量約0.32 cm。

    (3)后沙峪地區(qū)通過估算獲得近似值。1961年-2011年期間由于第四系地下水的開采造成的基巖變形量約為1.09 cm。

    4.2 建議

    (1)針對巖溶水開發(fā)利用基礎(chǔ)工作薄弱、盲目開采等問題,北京市積極布置了全面的、系統(tǒng)的巖溶水勘查。該項勘查的順利完成,不僅為保護巖溶水制定一系列的監(jiān)測和保護措施,也能夠進一步發(fā)揮巖溶水的供水潛力,將優(yōu)質(zhì)的巖溶水開采出來作為飲用水源。由于巖溶水主要埋藏在巖層中,注意采取適當?shù)拇胧┖?,其開采基本不會引發(fā)新的地面沉降問題。

    (2)北京市目前回灌井遠遠低于生產(chǎn)井的數(shù)量。為合理利用資源,應(yīng)持續(xù)推行開采-回灌制度,確保資源的可持續(xù)利用。

    參考文獻(References):

    [1] Biot M A.Theory of elasticity and consolidation for a porous anisotropic soid[J].Journal of Applied Physics,1955,26(2):182-185.

    [2] Wu J C,Shi X Q,Ye S J,et al.Numerical simulation of land subsidence induced by groundwater overexploitation in Su-Xi-Chang area,China[J].Environmental Geology,2009,57 (6):1409-1421.

    [3] Terzaghi K,Peck R B.Soil mechanics in engineering practice[M].New York:John Wiley and Sons,1967.

    [4] 吳建國,張玉貞.北京市順義區(qū)后沙峪鎮(zhèn)金泰花園別墅地熱地質(zhì)條件分析[J].中國煤田地質(zhì),2007,19 (2):69-71.(WU Jian-guo,ZHANG Yu-zhen.Analysis of geotherm geological condition in Jintai garden villa,Houshayu township,Shunyi district,Beijing[J].Coal Geology of China,2007,19(2):69-71.(in Chinese))

    [5] Sun Y.Conceptual model and potential assessment for the Xiaotangshan geothermal field,Beijing,China[J].Geothermal Training Program,2005:373-396.

    [6] Zhang Y Q,Gong H L.InSAR analysis of land subsidence caused by groundwater exploitation in Changping,Beijing,China[J].International Geoscience & Remote Sensing Symposium,2008.

    [7] 龔士良.上海地面沉降層次分析法研究[J].系統(tǒng)工程,1996,14(3):30-34.(GONG Shi-liang.A study on Shanghai land subsidence by the analytic hierarchy process[J].Syste ms Engineering,1996,14(3):30-34.(in Chinese))

    [8] 北京市地質(zhì)礦產(chǎn)勘查開發(fā)局,北京市地熱研究院.北京地熱[M].北京:中國大地出版社,2010.(Beijing Geology Prospecting & Developing Bureau,Beijing Geothermal Research Institute.Beijing geothermal energy[M].Beijing:China Land Press,2010.(in Chinese))

    [9] 北京市水文地質(zhì)工程地質(zhì)大隊.北京地下水[M].北京:中國大地出版社,2008.(Beijing Institute of Hydrogeology and Engineering Geology.Beijing groundwater[M].Beijing:China Land Press,2008.(in Chinese))

    [10] 劉成龍,車用太,呂金波,等.京北地熱田開發(fā)對地下流體動態(tài)的影響[J].地震地質(zhì),2006,28 (1):142-149.(LIU Cheng-long,CHE Yong-tai,LYU Jin-bo,et al.Influence of Jingbei geothermal field exploitation on subsurface fluid behaviors[J].Seismology and Geology,2006,28(1):142-149.(in Chinese))

    [11] 柯柏林.北京市平原區(qū)北部孫河斷裂的地熱地質(zhì)特征[J].現(xiàn)代地質(zhì),2009,23 (16):43-48.(KE Bai-lin.Geothermal and geological features of Sunhe fault in the northern part of Beijing plain[J].Geoscience,2009,23(1):43-48.(in Chinese))

    [12] 柯柏林.北京城區(qū)地熱田西北部地熱地質(zhì)特征[J].現(xiàn)代地質(zhì),2009,23 (1):49-56.(KE Bai-lin.Characteristics of geothermal geology in the northwestern part of the urban Beijing geothermal field[J].Geoscience,2009,23(1):49-56.(in Chinese))

    [13] 孫穎,劉久榮,韓征,等.北京市地熱資源開發(fā)利用狀況[J].安徽農(nóng)業(yè)科學,2009,37 (16):7564-7566.(SUN Ying,LIU Jiu-rong,HAN Zheng,et al.Exploration and utilization status of geothermal resources in Beijing city[J].Journal of Anhui Agricultural Sciences,2009,37(16):7564-7566.(in Chinese))

    [14] 朱加藍.北京市小湯山地區(qū)基巖熱儲數(shù)值模擬研究[D].北京:中國地質(zhì)大學(北京),2014.(ZHU Jia-lan.Numerical modeling of the bedrock geothermal reservoir in the Xiaotangshan area of Beijing,China[D].Beijing:China University of Geosciences,Beijing,2014.(in Chinese))

    [15] 王澤龍.北京市小湯山地區(qū)地溫場特征及地下熱水成因模式分析[D].北京:中國地質(zhì)大學(北京),2014.(WANG Ze-long.Numerical characteristics of the geothermal field and formation of the thermal groundwater in the Xiaotangshan area of Beijing,China[D].Beijing:China University of Geosciences,Beijing,2014.(in Chinese))

    [16] Holzer T L.Elastic expansion of the lithosphere caused by groundwater depletion[J].Journal of Geophysical Research,1979,84 (89):4689-4698.

    [17] Timoshenko S,Goodier J N.Theory of elasticity[M].New York:McGraw-Hill,1951.

    [18] Birch F S.Compressibility elastic constants in hand book of physical constants[M].Coloradan:Geological Society of America,1966.

    [19] Gutenberg B.Internal constitution of the Earth[M].New York:John Wiley and Sons,1951.

    猜你喜歡
    小湯山開采量基巖
    再談河北省灤平縣馬營子鄉(xiāng)高鍶天然礦泉水特征與開采量估算
    在小湯山致火神山
    草堂(2020年3期)2020-11-18 01:02:19
    北京小湯山啟用AI黑科技再“戰(zhàn)疫”
    伙伴(2020年3期)2020-04-14 04:48:35
    當上古名醫(yī)遇上新冠病毒阿里云智能化系統(tǒng)馳援鄭州“小湯山”醫(yī)院
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    基于改進物元的大壩基巖安全評價
    河北省基巖熱儲開發(fā)利用前景
    利用統(tǒng)計分析法預(yù)測地熱水可開采量應(yīng)注意的問題
    中國新疆石油開采量總額增長
    中亞信息(2016年1期)2016-03-07 01:00:42
    基載比對厚松散層薄基巖煤層安全開采的數(shù)值模擬研究
    河南科技(2014年23期)2014-02-27 14:18:47
    不卡一级毛片| 在线观看免费视频日本深夜| 亚洲五月色婷婷综合| 亚洲伊人色综图| 亚洲久久久国产精品| 天天操日日干夜夜撸| 在线av久久热| 99精品在免费线老司机午夜| 国产乱人伦免费视频| 波多野结衣av一区二区av| 国产精品九九99| 一边摸一边抽搐一进一小说 | 一边摸一边做爽爽视频免费| av中文乱码字幕在线| 亚洲久久久国产精品| 这个男人来自地球电影免费观看| 欧美日韩亚洲国产一区二区在线观看 | 老司机亚洲免费影院| 十八禁高潮呻吟视频| 国产一区二区三区综合在线观看| 桃红色精品国产亚洲av| 欧洲精品卡2卡3卡4卡5卡区| 女性被躁到高潮视频| 日韩欧美国产一区二区入口| 成人亚洲精品一区在线观看| 日韩有码中文字幕| 看免费av毛片| 欧美日韩黄片免| 18禁国产床啪视频网站| av超薄肉色丝袜交足视频| 自拍欧美九色日韩亚洲蝌蚪91| 香蕉久久夜色| av中文乱码字幕在线| 亚洲人成77777在线视频| 国产极品粉嫩免费观看在线| 99精品在免费线老司机午夜| 少妇粗大呻吟视频| 亚洲av成人av| 后天国语完整版免费观看| 黑人欧美特级aaaaaa片| 999久久久精品免费观看国产| 国产欧美日韩一区二区三| 夜夜爽天天搞| 日韩欧美在线二视频 | tocl精华| 50天的宝宝边吃奶边哭怎么回事| 一区二区三区激情视频| 18禁国产床啪视频网站| 91国产中文字幕| 久久国产精品人妻蜜桃| 亚洲三区欧美一区| 国产91精品成人一区二区三区| 视频在线观看一区二区三区| 伦理电影免费视频| 国产亚洲av高清不卡| 黄色女人牲交| 久久精品91无色码中文字幕| 男女高潮啪啪啪动态图| 久久香蕉精品热| 别揉我奶头~嗯~啊~动态视频| 最新美女视频免费是黄的| 国产精品亚洲一级av第二区| 久久精品国产清高在天天线| 久久精品91无色码中文字幕| 久久热在线av| 久久婷婷成人综合色麻豆| 亚洲熟妇中文字幕五十中出 | 午夜免费鲁丝| 悠悠久久av| 男女床上黄色一级片免费看| 国产欧美日韩一区二区精品| 脱女人内裤的视频| 老熟妇仑乱视频hdxx| 亚洲av美国av| 极品教师在线免费播放| 中文字幕另类日韩欧美亚洲嫩草| 久久人妻福利社区极品人妻图片| 久久人妻熟女aⅴ| 一级片免费观看大全| 超色免费av| 老司机午夜十八禁免费视频| 色尼玛亚洲综合影院| 美女福利国产在线| 少妇 在线观看| 极品少妇高潮喷水抽搐| 亚洲第一青青草原| 精品免费久久久久久久清纯 | 久久国产精品影院| 婷婷精品国产亚洲av在线 | 国产aⅴ精品一区二区三区波| 欧美人与性动交α欧美软件| 日本五十路高清| 国产精品一区二区在线观看99| 18禁黄网站禁片午夜丰满| 亚洲精品av麻豆狂野| 国产99白浆流出| 国产精品久久久av美女十八| 午夜福利,免费看| 午夜福利,免费看| av网站在线播放免费| 91精品国产国语对白视频| 水蜜桃什么品种好| 好男人电影高清在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲性夜色夜夜综合| 国产精华一区二区三区| 首页视频小说图片口味搜索| 人妻丰满熟妇av一区二区三区 | 少妇猛男粗大的猛烈进出视频| 免费日韩欧美在线观看| 国产精品综合久久久久久久免费 | 精品久久蜜臀av无| 亚洲精品久久午夜乱码| 欧美日韩黄片免| 一进一出抽搐gif免费好疼 | 黄色片一级片一级黄色片| 亚洲欧美激情综合另类| 国产区一区二久久| 黄片播放在线免费| 国产精品乱码一区二三区的特点 | 国内久久婷婷六月综合欲色啪| 亚洲av成人一区二区三| 不卡一级毛片| 两个人免费观看高清视频| 狠狠婷婷综合久久久久久88av| 成人黄色视频免费在线看| 欧美国产精品一级二级三级| 精品久久久精品久久久| 国产一区二区三区视频了| 日韩欧美一区视频在线观看| 中文字幕人妻丝袜制服| 免费在线观看黄色视频的| 少妇的丰满在线观看| 高清av免费在线| 国产日韩一区二区三区精品不卡| 午夜91福利影院| 日韩中文字幕欧美一区二区| 丁香欧美五月| 久久国产精品大桥未久av| 免费观看a级毛片全部| 伦理电影免费视频| 岛国在线观看网站| 丝瓜视频免费看黄片| 国产精品二区激情视频| av超薄肉色丝袜交足视频| 国产精品亚洲av一区麻豆| 中文字幕av电影在线播放| 午夜两性在线视频| 国产不卡一卡二| 51午夜福利影视在线观看| 日本撒尿小便嘘嘘汇集6| av国产精品久久久久影院| 十八禁高潮呻吟视频| 亚洲av欧美aⅴ国产| 亚洲情色 制服丝袜| 亚洲 国产 在线| 午夜福利在线免费观看网站| 亚洲aⅴ乱码一区二区在线播放 | 丁香六月欧美| 亚洲中文字幕日韩| av欧美777| 狠狠婷婷综合久久久久久88av| 日韩有码中文字幕| 黄片播放在线免费| 亚洲三区欧美一区| 国产精品久久久人人做人人爽| 欧美激情久久久久久爽电影 | 日韩有码中文字幕| 亚洲欧美日韩高清在线视频| 国产野战对白在线观看| 午夜福利在线观看吧| 精品少妇久久久久久888优播| 人人妻,人人澡人人爽秒播| 久久99一区二区三区| 最近最新中文字幕大全电影3 | 国产欧美日韩一区二区精品| 亚洲自偷自拍图片 自拍| 18禁国产床啪视频网站| 免费不卡黄色视频| 久久国产精品影院| 亚洲av电影在线进入| 色综合欧美亚洲国产小说| 精品国产亚洲在线| 叶爱在线成人免费视频播放| 国产高清国产精品国产三级| 国产在线一区二区三区精| 69av精品久久久久久| 亚洲在线自拍视频| 波多野结衣av一区二区av| 欧美日韩成人在线一区二区| 亚洲av成人一区二区三| 精品国产美女av久久久久小说| bbb黄色大片| 一边摸一边抽搐一进一出视频| 日本精品一区二区三区蜜桃| 老司机午夜十八禁免费视频| 精品卡一卡二卡四卡免费| 精品福利观看| 午夜视频精品福利| 久久久久精品人妻al黑| 18禁黄网站禁片午夜丰满| 窝窝影院91人妻| 在线十欧美十亚洲十日本专区| 精品亚洲成a人片在线观看| 老司机影院毛片| 一进一出抽搐动态| 王馨瑶露胸无遮挡在线观看| 久久精品国产a三级三级三级| 国产男女超爽视频在线观看| 国产一卡二卡三卡精品| 国产淫语在线视频| 国产av精品麻豆| aaaaa片日本免费| 91在线观看av| 日本五十路高清| 自线自在国产av| 曰老女人黄片| 亚洲伊人色综图| 国产91精品成人一区二区三区| 中文欧美无线码| 色综合婷婷激情| 一a级毛片在线观看| 最近最新免费中文字幕在线| 久久人妻av系列| 亚洲国产精品合色在线| 人妻 亚洲 视频| 桃红色精品国产亚洲av| 国产1区2区3区精品| 亚洲伊人色综图| 欧美激情高清一区二区三区| 男人操女人黄网站| 国产精品免费大片| 在线国产一区二区在线| 男人舔女人的私密视频| 国产一区二区三区在线臀色熟女 | a在线观看视频网站| 国产高清videossex| 一级黄色大片毛片| 91麻豆精品激情在线观看国产 | 午夜精品久久久久久毛片777| 国产精品电影一区二区三区 | 一进一出抽搐动态| 亚洲成人手机| 一级a爱片免费观看的视频| 亚洲av欧美aⅴ国产| 精品第一国产精品| 婷婷成人精品国产| 亚洲av欧美aⅴ国产| 色播在线永久视频| 91精品三级在线观看| 免费久久久久久久精品成人欧美视频| 精品人妻1区二区| 男人的好看免费观看在线视频 | 黄色女人牲交| 少妇的丰满在线观看| 每晚都被弄得嗷嗷叫到高潮| 交换朋友夫妻互换小说| 搡老熟女国产l中国老女人| www.999成人在线观看| 国产激情久久老熟女| 亚洲少妇的诱惑av| 在线观看免费午夜福利视频| 国产亚洲一区二区精品| 男女之事视频高清在线观看| 国产1区2区3区精品| 国产精品免费一区二区三区在线 | 伊人久久大香线蕉亚洲五| 国产成人影院久久av| 欧美国产精品va在线观看不卡| 女人久久www免费人成看片| 在线播放国产精品三级| 午夜福利欧美成人| 日韩有码中文字幕| 久久久久久免费高清国产稀缺| 亚洲黑人精品在线| 国产亚洲精品一区二区www | 欧美不卡视频在线免费观看 | 99久久综合精品五月天人人| 亚洲全国av大片| 一区二区三区精品91| 激情在线观看视频在线高清 | 亚洲五月色婷婷综合| 欧美日韩乱码在线| 少妇猛男粗大的猛烈进出视频| 欧美黄色片欧美黄色片| 国产成人啪精品午夜网站| 国产男靠女视频免费网站| 一级作爱视频免费观看| 精品一区二区三区视频在线观看免费 | 亚洲久久久国产精品| 欧美国产精品一级二级三级| 日本欧美视频一区| 色婷婷av一区二区三区视频| 亚洲全国av大片| 国产av又大| 国产精品亚洲一级av第二区| 亚洲aⅴ乱码一区二区在线播放 | 大型黄色视频在线免费观看| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 国产成+人综合+亚洲专区| 一区福利在线观看| 在线观看午夜福利视频| 脱女人内裤的视频| 国产在线观看jvid| 母亲3免费完整高清在线观看| 夫妻午夜视频| 精品亚洲成a人片在线观看| 操出白浆在线播放| 亚洲精品中文字幕在线视频| 99riav亚洲国产免费| 欧美日韩精品网址| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 人人妻人人爽人人添夜夜欢视频| 国产欧美日韩一区二区精品| 在线观看免费高清a一片| 搡老岳熟女国产| 视频在线观看一区二区三区| 国产在线精品亚洲第一网站| 欧美亚洲日本最大视频资源| 亚洲熟妇熟女久久| 国产97色在线日韩免费| 女同久久另类99精品国产91| 19禁男女啪啪无遮挡网站| 亚洲成国产人片在线观看| 亚洲五月色婷婷综合| 日韩熟女老妇一区二区性免费视频| 婷婷丁香在线五月| 国产99白浆流出| 久久人妻福利社区极品人妻图片| 欧美精品av麻豆av| 大型黄色视频在线免费观看| cao死你这个sao货| 亚洲欧美日韩高清在线视频| 欧美精品高潮呻吟av久久| 大香蕉久久成人网| 好看av亚洲va欧美ⅴa在| 99精品欧美一区二区三区四区| 女人精品久久久久毛片| 亚洲成人手机| 丰满迷人的少妇在线观看| 国产视频一区二区在线看| 欧美av亚洲av综合av国产av| 久久天躁狠狠躁夜夜2o2o| 亚洲 欧美一区二区三区| 极品少妇高潮喷水抽搐| 亚洲国产欧美一区二区综合| 免费在线观看日本一区| 欧美在线黄色| 精品熟女少妇八av免费久了| 久久精品亚洲熟妇少妇任你| 91麻豆精品激情在线观看国产 | 久久久久久久精品吃奶| 正在播放国产对白刺激| 乱人伦中国视频| 最新的欧美精品一区二区| 在线国产一区二区在线| 亚洲熟妇熟女久久| 国产1区2区3区精品| 久久香蕉国产精品| 看黄色毛片网站| 视频区欧美日本亚洲| 美女福利国产在线| 国产精品自产拍在线观看55亚洲 | 亚洲五月天丁香| 午夜激情av网站| 亚洲精华国产精华精| 国产精品免费视频内射| 高清毛片免费观看视频网站 | 久久精品国产a三级三级三级| 国产成人一区二区三区免费视频网站| 精品人妻1区二区| 两人在一起打扑克的视频| 亚洲av成人不卡在线观看播放网| aaaaa片日本免费| 午夜免费成人在线视频| 女人被狂操c到高潮| 午夜影院日韩av| 女人被狂操c到高潮| 12—13女人毛片做爰片一| 欧美精品人与动牲交sv欧美| 香蕉丝袜av| 一级作爱视频免费观看| 精品福利观看| 国内毛片毛片毛片毛片毛片| 精品国内亚洲2022精品成人 | 在线视频色国产色| 久久久久久久久久久久大奶| 在线播放国产精品三级| 国产亚洲一区二区精品| 亚洲av电影在线进入| 91大片在线观看| 在线观看一区二区三区激情| 男人操女人黄网站| 老熟妇乱子伦视频在线观看| 五月开心婷婷网| 久久ye,这里只有精品| 精品久久久久久,| 99香蕉大伊视频| 亚洲专区国产一区二区| 亚洲九九香蕉| 精品国内亚洲2022精品成人 | 亚洲人成电影观看| 国产精品乱码一区二三区的特点 | 精品国产超薄肉色丝袜足j| 亚洲av熟女| 男男h啪啪无遮挡| 亚洲成人免费av在线播放| 女性被躁到高潮视频| 国产xxxxx性猛交| 如日韩欧美国产精品一区二区三区| 欧美不卡视频在线免费观看 | 亚洲av片天天在线观看| 欧美黑人欧美精品刺激| 美女扒开内裤让男人捅视频| 国产男靠女视频免费网站| 午夜精品久久久久久毛片777| 亚洲午夜理论影院| 久久精品国产亚洲av香蕉五月 | 一区二区三区激情视频| 90打野战视频偷拍视频| 久久久国产一区二区| 麻豆乱淫一区二区| 午夜亚洲福利在线播放| 欧美久久黑人一区二区| 在线观看www视频免费| 大陆偷拍与自拍| 制服人妻中文乱码| 亚洲精品中文字幕在线视频| 夜夜爽天天搞| 久久久久久免费高清国产稀缺| 久久香蕉激情| 无遮挡黄片免费观看| 成人三级做爰电影| 亚洲精品一二三| 亚洲精品久久午夜乱码| 99精国产麻豆久久婷婷| 国产精品国产高清国产av | 亚洲专区中文字幕在线| 少妇的丰满在线观看| 日韩欧美一区视频在线观看| 人妻丰满熟妇av一区二区三区 | 亚洲欧美精品综合一区二区三区| 欧美成狂野欧美在线观看| 妹子高潮喷水视频| 三级毛片av免费| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧洲精品一区二区精品久久久| 99国产精品免费福利视频| 真人做人爱边吃奶动态| 最近最新中文字幕大全免费视频| 国产成+人综合+亚洲专区| 午夜亚洲福利在线播放| 国产精品亚洲一级av第二区| 18禁黄网站禁片午夜丰满| 天天躁狠狠躁夜夜躁狠狠躁| 成人影院久久| av有码第一页| 亚洲国产精品合色在线| 欧美精品一区二区免费开放| 国产男靠女视频免费网站| 熟女少妇亚洲综合色aaa.| 宅男免费午夜| 一区二区三区国产精品乱码| 国产成人精品久久二区二区91| 午夜两性在线视频| 很黄的视频免费| 国产真人三级小视频在线观看| 国产高清国产精品国产三级| 国产免费现黄频在线看| 欧美不卡视频在线免费观看 | 久久国产亚洲av麻豆专区| 少妇裸体淫交视频免费看高清 | 日本五十路高清| 熟女少妇亚洲综合色aaa.| 一本一本久久a久久精品综合妖精| 热re99久久国产66热| 亚洲av熟女| 在线观看www视频免费| 老司机午夜十八禁免费视频| 亚洲国产毛片av蜜桃av| 国产精品亚洲av一区麻豆| 亚洲成人免费av在线播放| 母亲3免费完整高清在线观看| 久久午夜亚洲精品久久| 午夜激情av网站| 亚洲国产欧美一区二区综合| 久久精品亚洲熟妇少妇任你| 国产av精品麻豆| 国产99久久九九免费精品| 午夜免费观看网址| 十分钟在线观看高清视频www| 飞空精品影院首页| 久久久久久久久久久久大奶| 丝袜美足系列| 亚洲精品美女久久av网站| 黑人欧美特级aaaaaa片| 80岁老熟妇乱子伦牲交| 久久久国产成人精品二区 | 亚洲欧美激情在线| 男女免费视频国产| 国产精品乱码一区二三区的特点 | 丝袜在线中文字幕| 国产色视频综合| 亚洲精品在线美女| 欧美黄色淫秽网站| 免费在线观看完整版高清| netflix在线观看网站| 老鸭窝网址在线观看| 久久久久久人人人人人| 高清av免费在线| 亚洲精品美女久久av网站| 久久精品亚洲熟妇少妇任你| 亚洲熟女精品中文字幕| 日本欧美视频一区| 国产99久久九九免费精品| 男女免费视频国产| 悠悠久久av| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲成av片中文字幕在线观看| 美国免费a级毛片| a在线观看视频网站| 人人妻,人人澡人人爽秒播| 高清在线国产一区| 老司机影院毛片| 麻豆成人av在线观看| 国产精品秋霞免费鲁丝片| 黄色女人牲交| 王馨瑶露胸无遮挡在线观看| 亚洲一区二区三区欧美精品| 午夜免费成人在线视频| 最新美女视频免费是黄的| 叶爱在线成人免费视频播放| 丰满饥渴人妻一区二区三| 亚洲一区高清亚洲精品| 99精品欧美一区二区三区四区| 高清欧美精品videossex| 成年人黄色毛片网站| 亚洲欧美色中文字幕在线| 亚洲在线自拍视频| 高清黄色对白视频在线免费看| 欧美 日韩 精品 国产| 欧美日韩亚洲国产一区二区在线观看 | 欧美亚洲 丝袜 人妻 在线| 老熟女久久久| 国产日韩一区二区三区精品不卡| 香蕉丝袜av| 999久久久精品免费观看国产| 老司机午夜福利在线观看视频| 亚洲熟妇中文字幕五十中出 | 一区二区日韩欧美中文字幕| 亚洲成人手机| 国产成人影院久久av| 91成人精品电影| 淫妇啪啪啪对白视频| 男女免费视频国产| 999久久久国产精品视频| 国产精品永久免费网站| 亚洲av电影在线进入| 午夜两性在线视频| 久久草成人影院| 亚洲熟妇熟女久久| av免费在线观看网站| 性少妇av在线| 国产精品二区激情视频| 国产精品国产高清国产av | 成年动漫av网址| 国产成人精品久久二区二区免费| 18禁裸乳无遮挡免费网站照片 | 老司机午夜十八禁免费视频| 精品乱码久久久久久99久播| 国产精品乱码一区二三区的特点 | 欧美乱妇无乱码| 国产精品久久久av美女十八| 窝窝影院91人妻| 午夜福利在线观看吧| 精品亚洲成国产av| 国产精品免费大片| 一本一本久久a久久精品综合妖精| 真人做人爱边吃奶动态| 亚洲av日韩在线播放| 精品国产亚洲在线| 性少妇av在线| 欧美在线黄色| 两个人看的免费小视频| 成人免费观看视频高清| 一本大道久久a久久精品| 人成视频在线观看免费观看| 黄色成人免费大全| 久久九九热精品免费| 午夜精品国产一区二区电影| 欧美日韩福利视频一区二区| 伊人久久大香线蕉亚洲五| 久久久久国产一级毛片高清牌| 一区二区日韩欧美中文字幕| 国产高清国产精品国产三级| 捣出白浆h1v1| 9191精品国产免费久久| 叶爱在线成人免费视频播放| 成年版毛片免费区| 1024视频免费在线观看| 中文亚洲av片在线观看爽 | 三上悠亚av全集在线观看| 国产成人欧美在线观看 | 黑人猛操日本美女一级片| 国产一卡二卡三卡精品| 人人妻人人添人人爽欧美一区卜| 免费高清在线观看日韩| 免费日韩欧美在线观看| 波多野结衣一区麻豆| 天天躁狠狠躁夜夜躁狠狠躁| 欧美精品av麻豆av| 狠狠狠狠99中文字幕| 婷婷精品国产亚洲av在线 |