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

    中國“采水型”地裂縫特征和成因分析

    2017-01-10 08:06:13王藝偉葉淑君龔緒龍
    高校地質(zhì)學報 2016年4期
    關鍵詞:基巖含水層成因

    王藝偉,葉淑君*,于 軍,龔緒龍

    1.南京大學地球科學與工程學院水科學系,南京210023;2.國土資源部地裂縫地質(zhì)災害重點實驗室,江蘇省地質(zhì)調(diào)查研究院,南京210000

    中國“采水型”地裂縫特征和成因分析

    王藝偉1,葉淑君1*,于 軍2,龔緒龍2

    1.南京大學地球科學與工程學院水科學系,南京210023;2.國土資源部地裂縫地質(zhì)災害重點實驗室,江蘇省地質(zhì)調(diào)查研究院,南京210000

    過量開采地下水誘發(fā)的地裂縫(“采水型”地裂縫)災害對社會造成嚴重損失,受到人們普遍關注。近幾十年來,國內(nèi)外學者對其成因機制進行了大量研究,取得了豐碩成果。本文系統(tǒng)分析了前人的研究成果,歸納出“采水型”地裂縫的差異沉降成因、拉張性水平應變成因和四種成因模式;然后針對中國三個典型的“采水型”地裂縫發(fā)育區(qū),即汾渭盆地、華北平原、蘇錫常地區(qū),闡述了各地地裂縫分布、發(fā)育特征和成因機理,并選擇典型地段地裂縫進行成因分析,認為汾渭盆地和華北平原地區(qū)部分地裂縫屬于“先期斷層模式”,而蘇錫常地區(qū)地裂縫總體屬于“基巖起伏模式”,并闡述了該模式下導致地裂縫發(fā)育的差異性變形、水平應變作用機理。

    地裂縫;地下水開采;成因機制;汾渭盆地;華北平原;蘇錫常

    1 引言

    地裂縫是在內(nèi)外動力作用下,巖土層發(fā)生破裂的一種表生的地質(zhì)現(xiàn)象,是地質(zhì)災害中地面變形災害的一種(王景明,2000)。自20世紀20年代在美國德克薩斯州的Goose Creek油田發(fā)現(xiàn)第一條剪切性質(zhì)的地裂縫之后(Prattand Johnson,1926),至今已在世界多個地區(qū)發(fā)現(xiàn)了多條規(guī)模、形態(tài)、啟裂性質(zhì)各異的地裂縫(武強和陳佩佩,2003),如美國西南部的亞利桑那州(Leonard,1929;Feth,1951;Carpenter,1993)、加利福尼亞州、內(nèi)瓦達州、德克薩斯州、新墨西哥州(Holzer,1976,1984;Holzer et al.,1979;Holzer and Pampeyan, 1981)、猶他州艾恩縣(Inkenbrandtet al.,2014)等、沙特阿拉伯(Bankher and Al-Harthi,1999;Youssefetal.,2014)、墨西哥(Pacheco-Martínez et al.,2013;Chaussard etal.,2014)、土耳其的安納托利亞和利比亞薩里爾(Rothenburgetal.,彭建兵等,2007)、中國的汾渭盆地(彭建兵,1992,2012;張家明,1990;王景明,1989)、華北平原(李昌存,2003;徐繼山,2012;安海波,2011)、蘇錫常地區(qū)(于軍等,2004;劉聰?shù)龋?004;Wangetal.,2009,2010)等都為地裂縫嚴重發(fā)育區(qū)。

    地裂縫作為一種地質(zhì)災害,嚴重影響人類生命、財產(chǎn)安全和社會穩(wěn)定與發(fā)展。常見的危害包括房屋坍塌、道路變形、地下管道破裂、對地鐵、高鐵等大型基礎設施造成嚴重破壞、地面開裂導致地下水污染等。在我國地裂縫發(fā)育的汾渭盆地、華北地區(qū)和蘇南地區(qū),上述破壞現(xiàn)象普遍存在,嚴重制約了當?shù)亟?jīng)濟和社會的可持續(xù)發(fā)展。

    防治地裂縫災害,必須首先查明其成因機制。美國學者關于地裂縫的研究開展較早,對地裂縫成因機制的研究也較為成熟。最早關于地裂縫成因研究的學說是Leonard(1929)提出的構(gòu)造成因?qū)W說,其認為亞利桑納州Picacho城附近的地裂縫的產(chǎn)生為附近地震活動所致。隨后,一些學者提出地下水開采成因觀點,如Flectcher等(1954)的滲透變形機理、Feth(1951)的不均勻沉降和Schumann與Poland(1970)的差異壓密變形機理、Neal等(1968)的土層失水收縮變形機理、Lofgren(1971)的滲透應力拖曳作用機理、Bouwer(1977)的剛性折裂機理等。隨著對地裂縫研究的不斷深入,又形成了構(gòu)造與地下水開采復合成因觀點。Sheng等(2003)將地質(zhì)條件和地下水開采聯(lián)系起來,提出不同地質(zhì)條件下可能產(chǎn)生地裂縫的四種成因模式。國內(nèi)地裂縫的成因研究最早始于西安地裂縫,目前該地區(qū)也是國內(nèi)地裂縫研究最為深入和系統(tǒng)的地區(qū)。張家明(1990)、王景明(1989)、彭建兵(1992)等先后就西安地裂縫成因提出多種學說,如:斷塊掀斜成因論、隱伏斷層蠕動成因論、斷裂伸展活動成因論等。他們的主要觀點為:構(gòu)造斷裂帶為地裂縫的產(chǎn)生提供地質(zhì)背景,過量抽取地下水是加速地裂縫活動并使其顯現(xiàn)于地表的直接誘發(fā)因素。2006-2011年,長安大學彭建兵教授對西安地裂縫進行了更為深入的研究,總結(jié)西安地裂縫成因機理為區(qū)域引張環(huán)境下的斷層伸展活動以及疊加抽水作用的內(nèi)外動力耦合成因模式(彭建兵,2012)。針對華北平原地裂縫成因也形成多種學說,如:古河道成因說(吳忱等,1991)、微裂隙閉鎖-開啟成因說(王景明,2000)、差異沉降說(邢忠信等,2004)、斷層蠕滑成因(李世雄等,2006)、優(yōu)勢面理論(羅國煜,2006)等。在對蘇錫常地裂縫的研究中,于軍等(2004)、伍洲云等(2003)認為蘇錫常地裂縫是特定構(gòu)造地質(zhì)背景下地下水強烈開采引發(fā)的不均勻地面沉降所致,基巖面的起伏及第四紀沉積物結(jié)構(gòu)差異是主要控制因素,長期過量開采地下水導致地面沉降與地裂縫相伴生。

    本文以過量開采地下水誘發(fā)的地裂縫(“采水型”地裂縫)為研究對象,在總結(jié)歸納國內(nèi)外學者對地裂縫成因機制研究成果的基礎上,重點分析了我國三個主要地裂縫發(fā)育區(qū)地裂縫的特征和成因機制,為科學布設地裂縫監(jiān)測系統(tǒng)、構(gòu)建模擬地裂縫的數(shù)學模型提供必要的理論依據(jù)。

    2 “采水型”地裂縫的成因機制

    19世紀50年代開始,美國學者對地下水開采產(chǎn)生的地裂縫成因機制進行一系列的研究(Feth,1951;Wolff,1970;Schumann and Poland,1970;Holzer,1976,1981,1984;Bouwer,1977;Helm,1994a,b;Budhu,2008,2011,Budhu and Adiya?man,2012)。由于發(fā)生地裂縫的地區(qū)一般也伴隨著地下水位下降及地面沉降災害,因此很早就有學者提出地下水開采成因,但對其影響地裂縫產(chǎn)生的機理解釋卻不盡相同。綜合前人研究,“采水型”地裂縫成因機制主要有差異沉降機制和水平拉張應變機制。

    2.1 差異沉降機制

    “采水型”地裂縫成因源自于差異沉降學說,其主要觀點是:地下水過量開采導致土層的壓縮變形,而在基巖起伏、壓縮層厚度不均勻等情況下,就會產(chǎn)生差異沉降,使上覆土體受彎或受剪導致破壞,并在地表產(chǎn)生地裂縫。

    差異沉降學說被許多學者的研究成果所證實。Feth(1951)對亞利桑那州Picacho盆地地裂縫進行了研究,認為由于含水層厚度變化引起不均勻沉降產(chǎn)生拉張應力,導致地裂縫產(chǎn)生,為拉張性破裂;Schumann與Poland(1970)對該地區(qū)進行了更為深入的研究,認為該地區(qū)下部可能存在隱伏斷層或是基巖面形態(tài)的突變,導致上覆土層的差異性沉降,從而在差異沉降最大部位產(chǎn)生最大的拉張應力,導致土層開裂變形產(chǎn)生地裂縫。Bouwer(1977)提出沉降盆地中的地裂縫產(chǎn)生的剛性折裂機理,由于沉降中心壓縮量大于沉降盆地邊緣,導致上覆土層在沉降盆地邊緣的斷層陡坎處可能發(fā)生整體的剛性反轉(zhuǎn),導致盆地邊緣發(fā)生折裂而形成地裂縫。

    從19世紀70年代開始,Holzer等分別對亞利桑那州東南部Picacho盆地和加里福利亞州的Fremont谷地分布的幾條主要地裂縫進行了一系列的調(diào)查研究,發(fā)現(xiàn)大多數(shù)地裂縫發(fā)生在沉降盆地中,且地裂縫最活動的地區(qū)分布在基巖起伏曲率最大處,此處由于差異變形易形成拉張區(qū)(Holzer,1976,1981,1984)。此外,為了進一步研究地下水位、地面沉降與地裂縫三者之間的關系,Holzer和Jachens(1982)對亞利桑那州中南部的卡薩格蘭卡進行了較為詳細的勘察和研究,了解該地區(qū)1940-1980地下水位和地面沉降歷史變化,并利用重力測量數(shù)據(jù)反演地下空間基巖的分布情況,研究結(jié)果顯示:重力異常反映出地下基巖面的傾斜或凸起的部位與地下水位變化及地裂縫產(chǎn)生位置具有很好的匹配性。Budhu(2008)提出一種單井抽水情況下潛水含水層中地裂縫的產(chǎn)生機制,認為潛水含水層中地裂縫的形成是彎曲和剪切共存,水位漏斗形成后,將有限單元體所承受力分解為壓應力和指向沉降中心的力矩,該力矩使相鄰土體單元間形成差異沉降,利用摩爾-庫倫準則作為判定準則,在力矩存在情況下土體所受剪應力更易達到抗剪強度而發(fā)生破壞。由這種理論得出,這種類型的地裂縫既可產(chǎn)生于地表并向下發(fā)展,也可由剪切力產(chǎn)生于地層深部而向上發(fā)展,但該理論僅適用于潛水含水層。Budhu和Adiyaman(2012)通過計算得出若地下無不連續(xù)結(jié)構(gòu)的存在(如先期斷層、基巖起伏等),當沉降漏斗的沉降梯度大于8×10-5,且與沉降中心的距離大于含水層厚度的1.4倍時,地裂縫才可能形成;若存在不連續(xù)結(jié)構(gòu),則在不連續(xù)界面處易發(fā)生地裂縫。

    以上研究者從不同方向證實抽取地下水引起的差異沉降是引發(fā)地裂縫活動和發(fā)展的重要機制。

    2.2 水平拉張應變機制

    隨著對地裂縫研究的進一步深入,人們發(fā)現(xiàn)有些地區(qū)的地裂縫發(fā)生在地表以下,并向上擴展,而根據(jù)差異沉降機理,地裂縫一般在地表發(fā)生張性破裂并自上而下發(fā)展。由此人們便得出差異沉降不可能是地裂縫的唯一誘發(fā)因素。隨著地面沉降由垂向一維變形向三維變形的研究進展,人們開始關注含水層在水平方向上的壓縮變形,由此提出水平拉張應變機制。據(jù)研究,這種水平方向上的拉張應變機制可能來自于兩個方面:一種是含水層中的水平滲流力對固體顆粒的拖拽作用(Lofgren,1971;Wolff,1970;Holzer&Jachens, 1982),另一種是含水層在水平方向上的壓縮釋水,引發(fā)水平應變所致(Burbey,1999)

    Lofgren指出:含水層中的水平滲流力可以提供足夠的水平應變以產(chǎn)生地裂縫,即抽水過程中沿著滲流方向產(chǎn)生的滲流力對固體顆粒產(chǎn)生牽引拖拽作用,該作用在含水層中逐漸積累并導致上覆土層發(fā)生張裂變形(Lofgren,1971)。Wolff(1970)在未固結(jié)含水層單井定流量抽水條件下,在抽水井附近利用測量裝置測量地面的徑向應變和切向應變,指出抽水所產(chǎn)生的滲流力可在近距離范圍內(nèi)對應變具有重要影響,且在抽水井附近徑向應變?yōu)閴嚎s應變,隨著距離抽水井位置的逐漸變遠,徑向應變由壓縮變?yōu)槔瓘?,此觀點在后來的研究中得到證實。

    Holzer與Jachens(1982)在對卡薩格蘭卡地裂縫進行研究時,估算出發(fā)生地裂縫的拉張應變的范圍為0.02%~0.2%。Rothenburg(1995)通過對利比亞薩里爾農(nóng)業(yè)區(qū)的水位下降歷史、地面沉降狀況及地裂縫特征進行分析,發(fā)現(xiàn)該地區(qū)地下水位下降和地面沉降量并不大,然而水平滲流力的存在使含水層產(chǎn)生較大的拉張應變,據(jù)他估算,每1×10-5數(shù)量級的拉張應變足以產(chǎn)生地裂縫。Helm也提出含水層徑向水平拉張應變可能是地裂縫產(chǎn)生的一種機理,并以Darcy-Gersevanov理論為基礎,假定水和固體顆粒均不可壓縮,結(jié)合含水層質(zhì)量守恒方程,推導出承壓含水層中徑向位移和徑向應變的解析解表達式,解釋了Wolff(1970)提出的未固結(jié)含水層中抽水井附近出現(xiàn)壓縮應變而遠離抽水井區(qū)域出現(xiàn)拉張應變的現(xiàn)象(Helm,1994a,b)。王慶良和劉玉海(2002)估算抽水所引起的含水層水平拉張應變可以達到10-6量級,推測抽水活動引起的含水層拉張應變可能是地裂縫活動的誘發(fā)因素之一,并通過大同機車工廠現(xiàn)場抽水實驗證實這一觀點。

    Burbey(1999,2001,2008)指出含水層壓縮產(chǎn)生的水平應變在地裂縫的形成中發(fā)揮著重要作用。含水層的釋水壓縮在垂直方向和水平方向均有體現(xiàn),將僅考慮垂向應變的一維數(shù)值模型和考慮水平應變的二維模型作比較發(fā)現(xiàn),隨著時間推移,水平方向上的釋水比例高達70%,由此產(chǎn)生的水平應力應變在盆地邊緣、斷層或基巖凸起附近易發(fā)生累積,容易誘發(fā)地裂縫形成。

    2.3 “采水型地裂縫”的成因模式

    通過對不同地區(qū)采水型地裂縫形成機制的研究發(fā)現(xiàn),抽取地下水引發(fā)地面沉降災害的地區(qū),如果沒有構(gòu)造作用形成的隱伏斷層,或下伏基巖面起伏等地質(zhì)背景,單純的地面沉降很難產(chǎn)生地裂縫。據(jù)此,Sheng等(2003)綜合分析國內(nèi)外各處地裂縫產(chǎn)生的各種地質(zhì)條件,提出不同地質(zhì)條件下產(chǎn)生地裂縫的四種成因模式:

    “薄弱水平面模式”,含水層存在薄弱水平層面。在含水層中,層與層之間的接觸面相對比較薄弱,當抽取地下水時,由于水平滲流力在接觸面處產(chǎn)生剪切應力而發(fā)生剪切位移,當達到抗剪強度時便產(chǎn)生剪切破壞并使其差異變形,導致深部地裂縫產(chǎn)生。

    “先期斷層模式”,含水層中存在具有壓性結(jié)構(gòu)面的斷層(逆斷層)。當不透水性斷層通過含水層時,由抽水引起含水層的水平壓縮和垂直壓縮,斷盤兩側(cè)產(chǎn)生差異性位移,并導致先期隱伏斷層復活或在斷層附近形成新的垂向裂隙構(gòu)造,若斷層位于抽水形成的拉張區(qū)域內(nèi),則更容易形成裂縫例如:美國拉斯維加斯谷地、Picacho盆地和中國的西安地區(qū),這些地區(qū)地裂縫明顯受到斷裂帶的制約。

    “基巖起伏模式”,含水層中存在基巖面起伏的情況。這種模型下有兩種可能產(chǎn)生地裂縫的機理:一是基巖起伏導致上覆粘土層發(fā)生彎曲而在隆起兩側(cè)產(chǎn)生差異性垂向壓縮,引起含水層出現(xiàn)水平拉伸并在含水層的最薄處形成拉伸帶,當拉伸強度達到一定極限強度時便發(fā)生張性破裂。二是基巖起伏導致含水層水平位移受到阻礙,從而發(fā)生深部土體破壞。這種破裂可能出現(xiàn)在毗鄰抽水井一方的基巖起伏附近,如果有兩個抽水中心,破裂也可能出現(xiàn)在基巖起伏的上部,當這些破裂向上擴展到地表就產(chǎn)生地裂縫。例如美國亞利桑那州的卡薩格蘭卡地區(qū)和中國的蘇錫常地區(qū),地裂縫的發(fā)生地區(qū)都位于基巖起伏曲率較大部位的上方。

    “結(jié)構(gòu)突變模式”,含水層結(jié)構(gòu)突變。該模型中,含水層結(jié)構(gòu)突變指的是含水層介質(zhì)厚度在某處發(fā)生突然變化,該變化可能會引起含水層深部產(chǎn)生不同程度的垂直剪切或水平拉伸,并產(chǎn)生裂縫。一般來說,先期斷裂和基巖起伏的存在也會引起含水層結(jié)構(gòu)的非均一性。

    Sheng等(2003)通過數(shù)值模擬對上述模型中影響地裂縫產(chǎn)生的因素進行敏感性分析得出,在不同的地質(zhì)環(huán)境中影響地裂縫的形成因素是不同的。在模式一中,土層周圍壓力和含水層的埋深是影響地裂縫形成的兩個重要因素。當含水層應力狀態(tài)為拉張時,易產(chǎn)生張性斷裂構(gòu)造;在相同的抽水條件下,淺部含水層比深部含水層更容易發(fā)生拉張破壞。在模式二中,土層周圍壓力、斷層的摩擦角和斷層傾角成為引發(fā)地裂縫的重要因素。模擬計算證實,斷層傾角為60°,斷層的摩擦臨界角為18°時,有利于地裂縫形成;若摩擦角增大可能會阻礙上覆地層產(chǎn)生拉張破壞,因此,低摩擦角有利于地裂縫的形成。在模式三中,含水層深度是影響地裂縫形成的重要因素。在其假定的模式條件下,當含水層深度大于200米時,地裂縫很難產(chǎn)生。在模式四中,土層周圍壓力和含水層的深度發(fā)揮重要作用,由于含水層的臨界深度在不同的地質(zhì)環(huán)境和水力條件下是不同的,因此深部含水層相對淺部含水層來說,不易產(chǎn)生地裂縫。

    圖1 中國地面沉降狀況和地裂縫發(fā)育分布(徐繼山,2012)Fig.1 Land subsidence and earth fissure distribution in China

    3 中國“采水型”地裂縫的特征與成因分析

    地裂縫在中國自古就有且分布范圍甚廣,早在幾千年前就有關于地裂縫的記載,這些裂縫記載均為自然因素所致(陳立偉,2007)。1958年,中國第一條采水型地裂縫在西北大學校園內(nèi)被發(fā)現(xiàn)(李永善,1986),其產(chǎn)生與當?shù)卮罅砍槿〉叵滤嘘P。由于發(fā)育的地裂縫數(shù)量少,不具有區(qū)域性特征,因而未引起足夠重視。隨著社會經(jīng)濟不斷發(fā)展,一些地區(qū)開始大量開采地下水來滿足經(jīng)濟發(fā)展需求,導致地面沉降災害的發(fā)生,并在地質(zhì)結(jié)構(gòu)中存在先期斷裂、基巖凸起、含水層介質(zhì)巖性突變的地區(qū)引發(fā)地裂縫災害。目前在我國的山西、陜西、河北、河南、山東、安徽、江蘇和北京等十余省市均有發(fā)育(圖1),分布廣泛。汾渭盆地、華北平原、蘇錫常三個地區(qū)是地裂縫災害較為嚴重的地區(qū),同時也是地面沉降現(xiàn)象最為嚴重的地區(qū)。

    3.1 汾渭盆地

    3.1.1 地裂縫分布

    汾渭盆地位于青藏、鄂爾多斯、華南、華北四大塊體之間,總體呈NE-SW向的“S”型,由一系列新生代斷陷盆地組成,是中國地裂縫發(fā)育最為強烈地區(qū)。自20世紀中、后期以來,汾渭盆地已有50余個縣市,180處發(fā)現(xiàn)地裂縫,總計約500余條(彭建兵等,2007),其中尤以西安、大同、榆次、臨汾、運城、三原、清徐等城市最為典型,同時也是由過量開采地下水引起地面沉降的嚴重地區(qū)。

    西安市是汾渭盆地乃至全國范圍內(nèi)地裂縫災害最為嚴重的地區(qū)之一。自1958年西安市發(fā)現(xiàn)第一條地裂縫以來,長安大學及西安市地調(diào)局等單位在地裂縫方面做了大量的地質(zhì)研究工作。目前在西安市發(fā)現(xiàn)的地裂縫達14條,地表出露總長度約70 km,面積約155 km2,均近東西向橫穿西安市區(qū)和郊區(qū)(圖2)(彭建兵,2012),且地裂縫帶在沉降區(qū)內(nèi)幾個較大沉降槽邊緣的陡變地帶顯現(xiàn),其發(fā)育過程與地下水開采密切相關。

    3.1.2 地裂縫成因研究

    在西安地裂縫研究的基礎上,提出了多種成因?qū)W說,主要有隱伏斷層蠕動成因(王景明,1989)、基底斷裂活動成因(李永善,1986)、斷塊掀斜成因(張家明,1990)、主斷裂伸展活動成因(彭建兵,1992)、內(nèi)外動力耦合成因論(彭建兵,2012)和構(gòu)造重力擴展成因(趙其華和王蘭生,1995)等。上述學說中大多以構(gòu)造作用作為背景來探索地裂縫的成因,其基本論點為:地裂縫形成于基底伸展構(gòu)造的橫向拉張應力場作用下,過量開采地下水為其重要的誘發(fā)因素。

    圖2 西安地裂縫分布圖Fig.2 Earth fissure distribution in Xi’an

    通過對西安地裂縫進行的大量野外調(diào)查、槽探以及淺層鉆探等一系列地質(zhì)工作,發(fā)現(xiàn)西安地裂縫與基底隱伏斷層緊密相連(陳立偉,2007,王景明,1989),因此西安地裂縫的成因符合Sheng等(2003)提出的“先期斷層模式”。但是控制西安地裂縫的基底隱伏斷層具有一定的活動性,這與第四紀地層中發(fā)現(xiàn)的斷層,其產(chǎn)狀(走向NE-NEE,傾角70°)一致得到印證。王啟耀等(2013)認為抽水作用下,先期斷層傾角在50°~60°區(qū)間時,上盤巖土體穩(wěn)定系數(shù)最低,最容易發(fā)生滑移,從而加劇了地裂縫的活動。在大量抽取深層地下水之前,西安基底隱伏斷層的活動速率很小,不足以引起地裂縫。但過量開采地下水之后,水位快速下降,逐步出現(xiàn)沉降漏斗區(qū),基底隱伏斷層重新復活,導致地裂縫開始出現(xiàn)并呈現(xiàn)持續(xù)加速活動狀態(tài)。

    位于西安交警總隊的F11地裂縫走向NE60°,傾角為80°,是西安地區(qū)活動最為強烈的地裂縫之一(彭建兵,2012)。根據(jù)GPS監(jiān)測數(shù)據(jù),2005-2008年,該地裂縫南北兩側(cè)的平均沉降速率分別為63mm/a和32mm/a,地裂縫兩側(cè)的差異沉降達到31mm/a,并分別以17.0 mm/a和8.4mm/a的年平均速率向西南、西北方向擴展。F11地裂縫的快速發(fā)展與地下水開采具有密切的關系,是地裂縫擴展的外動力因素。根據(jù)有效應力原理,抽水導致的有效應力的增加使含水層產(chǎn)生垂向和橫向變形,由于先期斷層的存在使斷層兩側(cè)地層產(chǎn)生差異性位移,這些位移會導致先期隱伏斷層重新復活,并在斷層處形成新的垂向裂隙,并且向地表發(fā)展;若斷層位于抽水形成的拉張區(qū)域內(nèi),則土體更容易產(chǎn)生破壞。

    3.2 華北平原

    3.2.1 地裂縫分布

    華北平原由北部的河北平原與南部的黃淮平原組成,是中國地裂縫分布數(shù)量最多的地區(qū)(王景明,2000)。華北平原地裂縫分布具有區(qū)域上的成帶性、局部的相對叢集性、擴展的方向性等特征,規(guī)模較大的地裂縫主要分布在近EW向(NWW向或NEE向)斷裂帶附近,地裂縫的發(fā)育受到構(gòu)造斷裂帶的控制(徐繼山,2012)。據(jù)調(diào)查顯示,截止到2013年,僅北部的河北平原就有9市共出現(xiàn)裂縫839條(圖3)(呂鳳蘭等,2014),主要分布在北京、保定、滄州、衡水、廊坊、邯鄲、邢臺、石家莊、河間等市,這些地區(qū)都是地面沉降的重災區(qū),具有多個大型水位降落漏斗,也是采水型地裂縫主要分布區(qū)。

    圖3 華北平原地裂縫分布圖(據(jù)中國地調(diào)局內(nèi)部資料,2013)Fig.3 Earth fissure distribution in North China Plain

    3.2.2 地裂縫成因研究

    針對華北平原“采水型”地裂縫成因問題,曾提出兩種假說:古河道成因說(吳忱等,1991)與差異沉降成因說(邢忠信等,2004)。近幾年,一些研究者對河北平原幾處典型地裂縫進行詳細的勘察,研究發(fā)現(xiàn),由于斷層長期緩慢活動形成與諸多斷層相連的隱伏裂縫,過量地下水開采等人類活動使隱伏裂縫活動性增強并逐漸擴展至地表(李志明等,2010;徐繼山,2012a,b)。李俊等(2003)、呂鳳蘭等(2014)、李世雄等(2006)等認為河北平原地裂縫形成與地震活動、斷層蠕滑、地下水開采、巖溶塌陷等幾種主要因素密切相關。邵長慶(2015)則從坳陷伸展構(gòu)造環(huán)境出發(fā),認為其為地裂縫發(fā)育提供了有利的拉張應力場條件,長期大量抽取地下水等工程活動為外部因素。綜合來說,華北平原“采水型”地裂縫的成因?qū)儆跀嗔褬?gòu)造與地下水開采共同作用,即斷裂構(gòu)造為地裂縫的發(fā)育提供孕裂背景,一些地區(qū)大規(guī)模的抽取地下水改變土層中的應力應變狀態(tài),在拉張應力或剪切應力超過最大極限后沿斷層面發(fā)生開裂。

    華北平原隆堯地裂縫長度約36 km,是該地區(qū)規(guī)模較大的地裂縫之一,地裂縫走向主要為EW和NEE向,這與隆堯斷裂構(gòu)造走向基本吻合(圖4)。宋偉等(2011)認為邢臺斷裂的長期蠕滑活動,以及人類對地下水的開采利用是形成隆堯地裂縫的主要原因。在構(gòu)造上,隆堯地裂縫受控于下覆斷裂帶,其緩慢的長期的構(gòu)造活動導致在一定深度范圍內(nèi)產(chǎn)生附加的伸展變形,這是一定的物質(zhì)基礎,隨著地下水不斷被開采,先期斷層兩側(cè)產(chǎn)生不均勻沉降,造成深部斷層附近拉張應力集中,當拉張應力超過一定極限便形成隱伏裂縫,隨著抽水活動持續(xù)進行,不均勻沉降進一步發(fā)育,隱伏地裂縫向上延伸,最終出露于地表。據(jù)資料顯示,隆堯地裂縫具有分段性特征,每段發(fā)生的時間點不同。東段出現(xiàn)在1966年邢臺7.2級地震后,部分地裂縫曾經(jīng)閉合,至今有些地段裂縫正逐漸趨于閉合狀態(tài),其影響并未持續(xù)到近期。中段、柏舍裂縫帶20世紀80年代后出現(xiàn),西段裂縫2000年顯現(xiàn),這三段地裂縫的的產(chǎn)生時間處于水位快速下降期(曹建玲等,2015)。據(jù)資料,20世紀70年代后至2010年,降落漏斗的地下水水位埋深已經(jīng)從3~5m下降到60m左右(馬學軍等,2011),且地裂縫發(fā)生地位于地下水降落漏斗的邊緣。駱祖江等(2013)通過全耦合模型模擬土體水平位移發(fā)現(xiàn),在地面沉降漏斗的邊緣區(qū)水平位移較大,屬于地裂縫的高發(fā)區(qū)。綜上所述,隆堯中段、柏舍裂縫帶和西段地裂縫的發(fā)育主要由地下水的過量開采誘發(fā),可以歸為“先期斷層模式”。

    圖5 蘇錫常地裂縫分布圖(劉聰?shù)龋?004)Fig.5 Earth fissure distribution in Su-Xi-Chang

    3.3 蘇錫常平原

    3.3.1 地裂縫分布及成因研究

    20世紀70年代后,蘇錫常地區(qū)開始過量開采地下水,從而引發(fā)了嚴重的地面沉降災害,造成嚴重經(jīng)濟損失。隨著地面沉降的發(fā)展,于1989年在常州武進區(qū)發(fā)生第一條地裂縫,90年代地裂縫進入高發(fā)期,共發(fā)育20處,2000年以后又發(fā)生4處。截止到目前,蘇錫常共發(fā)育地裂縫25處(圖5),主要分布在常州武進區(qū)、無錫錫山區(qū)、江陰市及蘇州張家港等地區(qū)。區(qū)內(nèi)的地裂縫在地表的展布方向不一致,主要有NE(NNE)向,近EW、NW向三組方向,并呈密集的帶狀展布,同一組地裂縫帶表現(xiàn)出平行斷續(xù)延伸。17處NE向地裂縫與下覆基巖凸起走向基本一致,約占總數(shù)的68%。地裂縫帶的寬度一般在30~100m,長度可達200~600 m,個別達1000 m以上(劉聰?shù)龋?004)。

    蘇錫常地區(qū)發(fā)育較厚的松散、高壓縮性粘土層,由于受長期超量開采地下水的影響,第II承壓含水層已形成較大范圍的降落漏斗,漏斗中心水位埋深由20m(二十世紀六十年代)降至80m左右(1994年),最大沉降量已接近2000 mm(2000年)(Zhang etal.,2010)。1995年,政府開始實施地下水限采措施,之后水位下降速率變緩,地面沉降現(xiàn)象得到明顯控制,最大沉降速率由1994年的110 mm/a逐漸降至2006年的10~20 mm/a,局部地裂縫活動速率減緩,地裂縫數(shù)目不再增多,但部分地裂縫仍處于活動之中。曾經(jīng)有過強烈活動的江陰河塘地裂縫及石塘灣因果岸地裂縫現(xiàn)已基本得到控制,目前活動最強且持續(xù)擴張的地裂縫位于無錫市錫山區(qū)光明村楊墅里一帶。

    劉聰?shù)龋?004)將優(yōu)勢面理論結(jié)合地質(zhì)統(tǒng)計分析法應用到蘇錫常地裂縫成因機理研究中,查明了該地區(qū)地裂縫發(fā)生的影響因素后,采用聚類分析方法進行成因分類,將該地區(qū)分為基巖潛山、埋藏階地、巖溶塌陷、土層結(jié)構(gòu)差異、地下水綜合開采五種類型。Wang等(2010)認為江陰河塘地裂縫的形成是基巖面起伏不平與長期過量開采地下水共同作用所致。Yu等(2006)利用三維地震勘探技術對因果岸地裂縫發(fā)生區(qū)進行地質(zhì)勘探,運用地震數(shù)據(jù)反演地層結(jié)構(gòu)及基巖面形態(tài),查明粘彈性地層、基巖面起伏形態(tài)、含水層空間分布差異是地裂縫產(chǎn)生的地質(zhì)條件,而人為過量開采地下水則是外部因素。第四紀海侵海退影響使該地區(qū)海陸過渡相沉積的硬塑-可塑狀粘土層分布較厚,構(gòu)造背景控制著地裂縫的發(fā)育。

    3.3.3 光明村楊墅里地裂縫研究

    (1)地裂縫分布

    光明村楊墅里地裂縫最初發(fā)現(xiàn)于1998年,2000年后在水渠巷、湯家里、徐家住基、基家壩4個村相繼發(fā)現(xiàn)地裂縫。2007年前后,地裂縫災害處于最嚴重時期,原有的地裂縫進一步變寬。楊墅里地裂縫呈北東向貫穿于楊墅里村的南部和東部,在村南的41號住宅和村東的25號住宅附近,地裂縫破壞最為嚴重(圖6,7)。位于村南的楊墅里村41~44號是主地裂縫區(qū),該處地裂縫走向48°,長度100m,影響寬度達20m,地面呈凹槽狀,裂縫帶西側(cè)比東側(cè)下沉近0.3m,造成多處房屋墻體出現(xiàn)裂縫,部分房屋傾斜,門窗變形,造成了大量的危房。

    圖6 光明村楊墅里地裂縫分布及水準線示意圖Fig.6 The distribution ofearth fissureand leveling survey line in Yang Shuli,Guangming Village

    圖7 光明村楊墅里地裂縫形跡及房屋破壞現(xiàn)象(左:41號和右:43號)Fig.7 Earth fissuresand damaged houses in Yang Shuli,Guangm ing Village

    (2)地裂縫動態(tài)

    為獲取地裂縫活動的定量數(shù)據(jù),江蘇省地質(zhì)調(diào)查研究院自2009年6月開始對楊墅里地裂縫進行定期監(jiān)測,并在橫跨地裂縫方向上布置了3條水準測量剖面(圖6)。由楊墅里A-A’剖面水準測量結(jié)果顯示(圖8),楊墅里43號西側(cè)(圖6,DL4點)地面沉降速率最大,可達5mm/a,2009年6月至2013年12月累積沉降量約20mm,而相鄰的DL2點僅有微量沉降,兩點間差異沉降達16mm(圖8)。此外,DL4和DL6點之間也存在約7mm的差異沉降。水準測量成果充分反映地面沉降由東向西具有明顯的差異分布規(guī)律。由于本地沒有絕對穩(wěn)定的沉降參考點,上述觀測到的沉降量只是相對值,實際沉降量應更大。

    圖8 光明村楊墅里地裂縫形變A-A’剖面Fig.8 The deformation of A-A'crosssection in Yang Shuli, Guangming Village

    4 地裂縫成因分析

    通過對無錫光明村歷史地下水位,基巖面起伏及地裂縫的發(fā)生及發(fā)展進行研究,發(fā)現(xiàn)三者具有較好的時空相關性(Wang et al.,2009,2010,武健強等,2013,劉聰?shù)龋?004)。根據(jù)地層剖面資料(圖9),光明村基底巖性為石英砂巖,基巖面埋深在70~160m間,起伏較大。楊墅里所在位置是基巖隆起區(qū),山脊最淺埋深70m左右。光明村基底分別向NWW和SEE方向傾斜,地裂縫發(fā)育的走向NE,與下伏基巖山脊走向一致。受基底起伏影響,含水層厚度分布不均,楊墅里以西一帶第Ⅱ承壓含水層相對發(fā)育,楊墅里淺部地下含水層系統(tǒng)(微承壓層)相對發(fā)育,第II承壓含水層(埋深70~160m)在基巖凸起附近發(fā)生尖滅乃至消失,含水層空間結(jié)構(gòu)表現(xiàn)出較大的差異性,這為差異性沉降的形成奠定了基礎。

    學術期刊是發(fā)布和傳播科研成果的重要平臺和媒介,以推動社會進步為己任,對傳播科研成果、促進科學技術發(fā)展發(fā)揮了重要作用。

    光明村地下水位下降受到周圍地下水開采的影響,光明村以北的河塘鎮(zhèn)一帶第Ⅱ承壓含水層超采導致水位降落漏斗不斷擴大,逐漸影響到光明村境內(nèi)。1998年,光明村地表發(fā)生第一條地裂縫,此時正處地下水位快速下降期(下降速率為10m/a)(圖10),同時還受含水層厚度分布的差異性以及基底起伏影響,三重因素共同作用下導致了差異性沉降,從而誘發(fā)地裂縫產(chǎn)生,之后水位下降速率變緩,地裂縫也處于緩慢活動中。地裂縫的產(chǎn)生受到基巖起伏等地質(zhì)背景條件的控制,再加上土層的沉積條件復雜,含水層本身結(jié)構(gòu)存在很大的空間差異性,含水層的厚度和富水性的不同,導致人類開采活動特征及強度時空差異性。在含水層富水性好的區(qū)域,相應的開采強度也較大,水位下降速度快,進而出現(xiàn)應力應變分布的空間差異性,導致垂向或水平方向的差異性變形,進而誘發(fā)地裂縫形成。

    圖10 光明村附近觀測孔地下水位變化Fig.10 Groundwater level changeswith time in Guangm ing Village

    關于光明村地裂縫成因機理主要有兩種:其一,基巖起伏導致上覆土層結(jié)構(gòu)和厚度突變,過量開采地下水導致基巖凸起兩側(cè)產(chǎn)生明顯的差異性垂向壓縮,而差異性垂向壓縮引起土層產(chǎn)生水平拉伸并在土層最薄處形成拉伸帶,當拉伸達到破壞強度時便會發(fā)生土層破裂。其二,當開采地下水后,含水層產(chǎn)生壓縮釋水,包括水平方向上和垂直方向上的壓縮釋水,隨著開采時間的推移,含水層的水平方向的釋水量逐漸增大,由此產(chǎn)生的水平應變力傳播到基巖隆起處而受阻,為地裂縫的產(chǎn)生創(chuàng)造了有利條件。從光明村地裂縫的形成和發(fā)展過程,可以充分說明該區(qū)地裂縫受控于基底巖層的起伏并與水位動態(tài)變化密切相關,屬于典型的采水型地裂縫。

    5 結(jié)論

    中國“采水型”地裂縫的成因可歸納為由隱伏斷層、或基底基巖起伏與過量開采地下水共同作用的結(jié)果,兩種因素缺一不可。其中先期斷層或基底基巖起伏對地裂縫的發(fā)育和分布起著控制作用,而過量開采地下水是地裂縫產(chǎn)生的直接誘發(fā)因素。但不同地區(qū)地裂縫特征和成因又存在明顯的差異性,總體來說遵循兩種成因模式,即“先期斷層模式”和“基巖起伏模式”。

    汾渭盆地和華北地區(qū)采水型地裂縫分布與先期斷層的分布在空間上具有高度一致性,抽水引起地裂縫可能由兩種機理分別或共同作用導致:一是斷層的不透水性使含水層系統(tǒng)的水平壓縮量和垂直壓縮量在斷層兩側(cè)產(chǎn)生差異,差異位移會導致先期隱伏斷層重新復活或在斷層處形成新的裂隙構(gòu)造,隨著抽水繼續(xù),這些裂隙向上發(fā)展,一般在暴雨過后裸露于地表形成地裂縫。二是含水層的釋水壓縮更趨向于水平方向的壓縮,過量開采地下水產(chǎn)生的水平應力應變在遇到阻礙其傳播的斷層邊界條件時,便會發(fā)生應變累積導致地裂縫形成。

    蘇錫常地區(qū)采水型地裂縫成因總體屬于“基巖起伏模式”,可能存在兩種機理分別或共同作用導致地裂縫的發(fā)生:一是基巖起伏導致上覆土層結(jié)構(gòu)和厚度突變,過量開采地下水導致基巖凸起兩側(cè)產(chǎn)生明顯的差異性垂向壓縮所致;二是過量開采地下水產(chǎn)生的水平應力應變在遇到阻礙其傳播的基巖凸起時,便會在此發(fā)生累積導致形成地裂縫萌發(fā)的有利條件所致。此外,水平滲流力土體顆粒的拖拽作用也可能在土層深部形成拉張性應力區(qū),在上述兩種情況下,會進一步加劇地裂縫的形成。

    “采水型”地裂縫的成因機制是非常復雜的,不論受控于先期斷層、還是基底起伏,都還包含了受控于含水層系統(tǒng)結(jié)構(gòu)的空間差異性,進而控制了人類采水活動特征及強度的時空差異性,地裂縫的形成受到多重因素的制約。本文綜合前人研究結(jié)果對中國采水型地裂縫的特征和成因進行了定性分析,這將為科學布設地裂縫監(jiān)測系統(tǒng)、構(gòu)建模擬地裂縫的數(shù)學模型提供必要的理論依據(jù)。

    陳立偉.2007.地裂縫擴展機理研究[D].西安:長安大學.

    曹建玲,白淑芳,姜正義.2015.隆堯地裂縫形變特征及成因分析[J].華北地震科學,33(2):13-19.

    李昌存.2003.河北平原地裂縫研究[D].北京:中國地質(zhì)大學.

    李俊,劉金峰,莫多聞.2003.河北平原地裂縫的分布規(guī)律及成因初探[J].水土保持研究,10(3):62-65.

    李世雄,李守定,郜洪強.2006.河北平原地裂縫分布特征及成因機制研究[J].工程地質(zhì)學報,14(2):178-183.

    李永善.1986.西安地裂縫[M].北京:地震出版社.

    李志明,楊旭東,蘭劍梅,等.2010.河北邢臺柏鄉(xiāng)地裂縫成因分析[J].水文地質(zhì)工程地質(zhì),37(2):135-138.

    劉聰,袁曉軍,朱錦旗.2004.蘇錫常地裂縫[M].武漢:中國地質(zhì)大學出版社.

    羅國煜.2006.關于地裂縫地質(zhì)災害的優(yōu)勢面分析[J].地質(zhì)學刊,30(1):1-5.

    駱祖江,王琰,田小偉,等.2013.滄州市地下水開采與地面沉降地裂縫模擬預測[J].水利學報,44(2):198-204.

    呂鳳蘭,馬學軍,邵長慶,等.2014.河北平原地裂縫分布特征及成因分析[J].上海國土資源,35(4):49-52.

    馬學軍,宋偉,王洪磊,等.2011.隆堯地裂的成因[J].地質(zhì)災害與環(huán)境保護,22(3):43-46.

    彭建兵,范文,李喜安,等.2007.汾渭盆地地裂縫成因研究中的若干關鍵問題[J].工程地質(zhì)學報,15(4):433-440.

    彭建兵.1992.渭河斷裂帶的構(gòu)造演化與地震活動[J].地震地質(zhì),14(2): 113-120.

    彭建兵.2012.西安地裂縫災害[M].北京:科學出版社.

    宋偉,馬學軍,賈國欣.2011.邢臺隆堯典型斷層蠕滑型地裂縫成因初析[J].地質(zhì)學刊,35(4):405-412.

    邵長慶.2015.地裂縫構(gòu)造背景與坳陷伸展構(gòu)造環(huán)境的關系[J].中國水土保持,(4):53-56.

    王景明.2000.地裂縫及其災害的理論與應用[M].西安:陜西科學技術出版社.

    王景明.1989.論西安地裂縫[J].地震地質(zhì),11(3):85-93.

    王啟耀,蔣臻蔚,彭建兵.2013.抽水作用下先存斷裂活化滑移機制研究[J].水文地質(zhì)工程地質(zhì),40(2):108-112.

    王慶良,劉玉海.2002.抽水引起的含水層應變—地裂縫活動新機理[J].工程地質(zhì)學報,10(1):46-50.

    吳忱,朱宣清,何乃華,等.1991.華北平原古河道的形成研究[J].中國科學:化學生命科學地學,21(2):188-197.

    伍洲云,余勤,張云.2003.蘇錫常地區(qū)地裂縫形成過程[J].水文地質(zhì)工程地質(zhì),30(1):67-72.

    武健強,周鋼軍,吳曙亮,等.2013.無錫光明村地裂縫研究[J].地質(zhì)學刊,37(2):203-207.

    武強,陳佩佩.2003.地裂縫災害研究現(xiàn)狀與展望[J].中國地質(zhì)災害與防治學報,14(1):22-27.

    邢忠信,李和學,張熟,等.2004.滄州市地面沉降研究及防治對策[J].地質(zhì)調(diào)查與研究,27(3):157-163.

    徐繼山,馬潤勇,彭建兵,等.2012a.河間地區(qū)地裂縫分形特征研究[J].自然災害學報,21(3):177-183.

    徐繼山,彭建兵,馬學軍,等.2012b.邢臺市隆堯地裂縫發(fā)育特征及成因分析[J].工程地質(zhì)學報,20(2):160-169.

    徐繼山.2012.華北陸緣盆地地裂縫成因機理研究[D].西安:長安大學.

    于軍,王曉梅,蘇小四,等.2004.蘇錫常地區(qū)地裂縫地質(zhì)災害形成機理分析[J].吉林大學學報(地球科學版),34(2):236-241.

    張家明.1990.西安地裂縫研究[M].西安:西北大學出版社.

    趙其華,王蘭生.1995.構(gòu)造重力擴展機制的地質(zhì)力學模擬研究[J].地質(zhì)災害與環(huán)境保護,3(1):21-27.

    Bankher K A and Al-Harthi A A.1999.Earth fissuring and land subsidence inwestern SaudiArabia[J].NaturalHazards,20(1):21-42.

    Bouwer H.1977.Land subsidence and cracking due to ground-water depletiona[J].Groundwater.,15(5):358-364.

    Budhu M and Adiyaman I.2012.Earth fissure formation from groundwater pumping and the influence of a stiff upper cemented layer[J].Quarterly JournalofEngineeringGeology and Hydrogeology,45(2):197-205.

    Budhu M.2010.Earth fissure formation from themechanics of groundwater pumping[J].International JournalofGeomechanics,11(1):1-11.

    Budhu M.2008.Mechanicsof earth fissures using the Mohr-Coulomb failure criterion[J].Environmental&EngineeringGeoscience,14(4):281-295. Burbey T J.1999.Effects of horizontal strain in estimating specific storage and compaction in confined and leaky aquifer systems[J].Hydrogeology Journal,7(6):521-532.

    Burbey T J.2001.Stress-strain analyses for aquifer-system characterization [J].Groundwater,39(1):128-136.

    Burbey T J.2008.The influence of geologic structures on deformation due to groundwaterwithdrawal[J].Groundwater,46(2):202-11.

    Carpenter M C.1993.Earth-fissuremovements associated with fluctuations in ground-water levels near the Picacho Mountains,south-central Arizona,1980-84[J].ProfessionalPaper.

    Chaussard E,Wdowinski S,Cabral-Cano E,et al.2014.Land subsidence in central Mexico detected by ALOS InSAR time-series[J].Remote Sensing of Environment,140(1):94-106.

    Feth J H.1951.Structural reconnaissance of the Red Rock quadrangle, Arizona[J].Geological Society of America Bulletin Part Pages: 2331-1258.

    Fletcher JE,Harris K,Peterson H B,etal.1954.Piping[J].Eos Transactions American GeophysicalUnion,35(2):258–263.

    Helm D C.1994a.Horizontal aquifermovement in a theis‐thiem confined system[J].Water Resources Research,30(4):953-964.

    Helm DC.1994b.Hydraulic forces thatplay a role in generating fissures at depth[J].Environmental&Engineering Geoscience,31:3.

    Hernandez-Marin M and Burbey T J.2010.Controls on initiation and propagation of pumping-induced earth fissures:insights from numerical simulations[J].Hydrogeology Journal,18(8):1773-1785.

    Holzer TL,DavisSN and Lofgren BE.1979.Faulting caused by groundwater extraction in southcentral Arizona[J].Journal of Geophysical Research Atmospheres,84(B2):603-612.

    Holzer T L.1976.Ground failure in areas of subsidence due to ground-water decline in the United States[J].Proceedings of the Anaheim Symposium;

    Holzer TL and Pampeyan EH.1981.Earth fissuresand localized differential subsidence[J].WaterResourcesResearch,17(1):223-227.

    Holzer T L.1984.Ground failure induced by ground-water withdrawal from unconsolidated sediment[J].Reviews in Engineering Geology,67-105.

    Inkenbrandt P,Lund W,Lowe M,et al.2014.Investigation of land subsidence and earth fissures in Cedar Valley,Iron County,Utah[J]. Utah GeologicalSurvey.

    Jachens R C and Holzer TL.1982.Differential compaction mechanism for earth fissures near Casa Grande,Arizona[J].Geological Society of America Bulletin,93(10):998-1012.

    Leonard R.1929.An earth fissure in southern Arizona[J].The Journal of Geology,37(8):765-774.

    Lofgren B.1971.Significant role of seepage stresses in compressible aquifer systems[J].Transactions American GeophysicalUnion,52(11).

    Neal J T,Langer A M and Kerr P E.Giant Desiccation Polygons Of Great Basin Playas[J].Geological Society of America Bulletin,79(1):69.

    Pacheco-Martínez J,Hernandez-Marín M,Burbey T J,et al.2013.Land subsidence and ground failure associated to groundwater exploitation in the Aguascalientes Valley,México[J].Engineering Geology,164(18): 172-186.

    PrattW E and Johnson DW.1926.Local subsidence of the Goose Creek oil field[J].The JournalofGeology,34(7):577-590.

    Rothenburg L,Obah A and El B S.1995.Horizontal ground movements due to water abstraction and formation of earth fissures[J].International Association ofHydro logical Sciences,234:239-249.

    Schumann H and Poland J.1970 Land subsidence,earth fissures and groundwater withdrawal in south-central arizona,U.S.A[J].USENIX Systems Administration Conference.

    Sheng Z,Helm D C,Li J.2003.Mechanisms of earth fissuring caused by groundwaterwithdrawal[J].Environmental&Engineering Geoscience., 9(9):313-324.

    Wang G Y,You G,Shi B,et al.2010.Earth fissures in Jiangsu Province, China and geological investigation of Hetang earth fissure[J]. Environmental Earth Sciences,60(1):35-43.

    Wang G Y,You G,Shi B,et al.2009.Earth fissures triggered by groundwaterwithdrawal and coupled by geological structures in Jiangsu Province,China[J].EnvironmentalGeology,57(5):1047-1054.

    Wolff R G.1970.Relationship between horizontal strain near a well and reverse water level fluctuation[J].Water Resources Research,6(6): 1721-1728.

    Youssef A M,Sabtan A A,Maerz N H,et al.2014.Earth fissures in Wadi Najran,Kingdom of Saudi Arabia[J].Natural Hazards,71(3): 2013-2027.

    Yu J,Chen G M,Wang X M,et al.2006.Formation mechanism analysis on ground fissure based on 3D seismic exploration data in Yinguo'an area, WuxiCity[J].Hydrogeology&EngineeringGeology,33(4):117-120.

    Zhang Y,Xue Y Q,Wu JC,et al.2010.Excessive groundwater withdrawal and resultant land subsidence in the Su-Xi-Chang area,China[J]. EnvironmentalEarth Sciences,61(6):1135-1143.

    Featuresand M echanism sof Earth Fissures Induced GroundwaterW ithdrawalin China

    WANG Yiwei1,YE Shujun1*,YU Jun2,GONG Xulong2
    1.Departmentof Hydrosciences,Schoolof Earth Science and Engineering,Nanjing University,Nanjing 210023,Chian; 2.Key Laboratory of Earth FissureGeologicalGazard of Land and Resources,instituteofgeologicalsurvey in Jiangsu,Nanjing 210000,China

    Excessive exploitation of groundwater in China has induced serious earth fissure hazards,which can cause significant damage and heavy economic losses,and has received significant attention.In recent decades,mechanisms of earth fissures induced by groundwaterwithdrawalhave been studied by researchersworldwide andmany remarkableachievementshave been accomplished.The results of previous studies were analyzed in this article firstly and then the mechanisms of earth fissure induced by excessive exploitation of groundwater were summarized and categorized as differential subsidence,horizontal tensile strain,and Four Earth Fissure Formation ConceptualModels.Then the distribution ofearth fissures,development featuresandmechanismsof earth fissures in three regions(i.e.,the Fenwei Basin,North China Plain and Su-Xi-Chang in China)where typical“groundwater-withdraw-induced”fissures developed,were elaborated,and a typical earth fissurewas chosen to analyze its cause.Results indicate that earth fissures in the FenweiBasin and North China Plain weremainly caused by pre-existing faultsand excessively ground water exploitation,while the fissures in the Su-Xi-Chang region were generally induced by a“bedrock ridge”below the aquifer and excessively ground water exploitation.Differentialcompactionand horizontalstrainmechanismscausingthedevelopmentofsuchearth fissureswere thenexplained.

    earth fissure;groundwaterwithdrawal;mechanism;FenweiBasin;North China Plain;Su-Xi-Changarea

    YEShujun,Professor;E-mail:sjye@nju.edu.cn

    P641

    A文獻標識碼:1006-7493(2016)04-0741-12

    10.16108/j.issn1006-7493.2015259

    2016-01-08;

    2016-07-27

    本文為國土資源部地裂縫地質(zhì)災害重點實驗室開放課題(2014003);國土資源部公益項目(201411096);江蘇省自然科學基金(批BK2012730);國家自然科學基金項目(41272259)聯(lián)合資助

    王藝偉,女,1991年生,碩士研究生,主要從事地下水和地面沉降數(shù)值模擬的研究;E-mail:wangyiwei1005@126.com

    *通訊作者:葉淑君,教授;E-mail:sjye@nju.edu.cn

    猜你喜歡
    基巖含水層成因
    全球多個含水層里的水正快速流失
    說說兩種『白氣』的成因
    輸水渠防滲墻及基巖滲透系數(shù)敏感性分析
    暈紋石成因解讀(上)
    寶藏(2017年7期)2017-08-09 08:15:19
    基于改進物元的大壩基巖安全評價
    美國西部奧加拉拉含水層水位下降原因初探
    翻譯實踐問題及成因
    唐山文學(2016年11期)2016-03-20 15:25:54
    河北省基巖熱儲開發(fā)利用前景
    全球地下含水層下降驚人:要被抽干了
    地理教學(2015年14期)2015-03-31 20:04:53
    巖溶含水層水流模型研究進展
    国产精品免费一区二区三区在线| 日本在线视频免费播放| 国产欧美日韩精品亚洲av| 欧美又色又爽又黄视频| 日本黄色视频三级网站网址| 搞女人的毛片| 丁香六月欧美| 一级毛片精品| 日本一本二区三区精品| 天堂影院成人在线观看| 久久国产精品人妻蜜桃| 亚洲第一电影网av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲成a人片在线一区二区| 亚洲欧美日韩高清专用| 成人无遮挡网站| 亚洲欧美日韩卡通动漫| 成人午夜高清在线视频| 欧美日韩中文字幕国产精品一区二区三区| 色吧在线观看| 国产高清激情床上av| 欧美av亚洲av综合av国产av| 我的老师免费观看完整版| 精品福利观看| 色老头精品视频在线观看| 久久精品国产综合久久久| 久久国产精品人妻蜜桃| 亚洲在线观看片| 国产高清videossex| xxxwww97欧美| 久久亚洲精品不卡| 女人高潮潮喷娇喘18禁视频| 欧美高清成人免费视频www| 国产高清videossex| 亚洲精品乱码久久久v下载方式 | 久久亚洲精品不卡| 国产精品香港三级国产av潘金莲| 国产精品久久久av美女十八| 女人高潮潮喷娇喘18禁视频| 日韩精品青青久久久久久| 日韩精品青青久久久久久| 午夜影院日韩av| 最好的美女福利视频网| 国产乱人视频| 三级国产精品欧美在线观看 | 亚洲国产欧美人成| 亚洲av成人av| 日本与韩国留学比较| av天堂在线播放| 国产精品久久久久久亚洲av鲁大| 亚洲成人久久性| 婷婷六月久久综合丁香| 99久久精品一区二区三区| 免费一级毛片在线播放高清视频| 国产综合懂色| 亚洲在线观看片| 九色国产91popny在线| 国产aⅴ精品一区二区三区波| 久久人妻av系列| 午夜久久久久精精品| 国产人伦9x9x在线观看| 亚洲欧洲精品一区二区精品久久久| 非洲黑人性xxxx精品又粗又长| tocl精华| 又大又爽又粗| 老鸭窝网址在线观看| 天天躁日日操中文字幕| a级毛片在线看网站| 欧美xxxx黑人xx丫x性爽| 中亚洲国语对白在线视频| 久久久久国内视频| 最近在线观看免费完整版| tocl精华| 亚洲成a人片在线一区二区| 在线看三级毛片| 级片在线观看| 亚洲 欧美一区二区三区| 好看av亚洲va欧美ⅴa在| 亚洲一区二区三区不卡视频| 中文字幕高清在线视频| 免费在线观看视频国产中文字幕亚洲| 国语自产精品视频在线第100页| 国产精品1区2区在线观看.| 国产黄片美女视频| 国产视频一区二区在线看| 欧美日韩黄片免| 大型黄色视频在线免费观看| 日本五十路高清| 午夜免费成人在线视频| 日韩成人在线观看一区二区三区| 亚洲成人中文字幕在线播放| 国产视频一区二区在线看| 一个人看视频在线观看www免费 | 中文字幕久久专区| 亚洲av中文字字幕乱码综合| 天堂av国产一区二区熟女人妻| 特级一级黄色大片| 成人国产综合亚洲| 国产精品久久久久久久电影 | 久久久久久久久中文| 日日摸夜夜添夜夜添小说| 少妇的逼水好多| 波多野结衣巨乳人妻| 成人高潮视频无遮挡免费网站| 亚洲电影在线观看av| 国内精品美女久久久久久| 此物有八面人人有两片| 人人妻人人澡欧美一区二区| 亚洲片人在线观看| 88av欧美| 国产伦人伦偷精品视频| 99riav亚洲国产免费| 国产精品美女特级片免费视频播放器 | 免费看美女性在线毛片视频| av视频在线观看入口| 国产精品久久视频播放| 一进一出抽搐动态| 久久久久性生活片| 亚洲精品乱码久久久v下载方式 | 国产精品永久免费网站| www.熟女人妻精品国产| 伦理电影免费视频| 黄色丝袜av网址大全| 国产伦在线观看视频一区| 国产99白浆流出| 午夜免费成人在线视频| www日本黄色视频网| 天天躁狠狠躁夜夜躁狠狠躁| 99精品欧美一区二区三区四区| 成人特级av手机在线观看| 亚洲精品中文字幕一二三四区| 久久久色成人| 欧美一级毛片孕妇| 国产精品一区二区免费欧美| 国产熟女xx| av福利片在线观看| 精品电影一区二区在线| 亚洲片人在线观看| 淫秽高清视频在线观看| 亚洲av熟女| 久久性视频一级片| 成人鲁丝片一二三区免费| 真实男女啪啪啪动态图| 黑人操中国人逼视频| av福利片在线观看| 国产三级中文精品| 草草在线视频免费看| 九色成人免费人妻av| 一夜夜www| 午夜福利视频1000在线观看| 在线免费观看不下载黄p国产 | 午夜日韩欧美国产| 国产视频一区二区在线看| 午夜福利18| 男女床上黄色一级片免费看| 一夜夜www| 国产真实乱freesex| 国产精品 欧美亚洲| 可以在线观看毛片的网站| 日韩精品中文字幕看吧| 夜夜看夜夜爽夜夜摸| 久久久久亚洲av毛片大全| 精品国产美女av久久久久小说| 舔av片在线| 麻豆av在线久日| 女人被狂操c到高潮| 真人做人爱边吃奶动态| 国产精品国产高清国产av| 午夜两性在线视频| 好男人电影高清在线观看| 黄片小视频在线播放| 制服丝袜大香蕉在线| 欧美乱妇无乱码| 欧美成人免费av一区二区三区| 中文字幕av在线有码专区| 亚洲天堂国产精品一区在线| 黄片大片在线免费观看| 国产单亲对白刺激| 夜夜躁狠狠躁天天躁| 欧美绝顶高潮抽搐喷水| 一个人看的www免费观看视频| 午夜视频精品福利| 免费av毛片视频| 黄色视频,在线免费观看| 美女cb高潮喷水在线观看 | 在线a可以看的网站| 午夜免费成人在线视频| 美女高潮喷水抽搐中文字幕| 国产成人精品无人区| 青草久久国产| 国产精品永久免费网站| 又黄又粗又硬又大视频| 成年女人看的毛片在线观看| 成人三级黄色视频| 国产欧美日韩精品亚洲av| 99久久无色码亚洲精品果冻| 最新在线观看一区二区三区| 人妻丰满熟妇av一区二区三区| 免费看光身美女| 国产极品精品免费视频能看的| 在线观看一区二区三区| 国产精品爽爽va在线观看网站| 青草久久国产| 天堂影院成人在线观看| 亚洲午夜精品一区,二区,三区| 久久欧美精品欧美久久欧美| 国产91精品成人一区二区三区| 国产成人精品久久二区二区91| 成年免费大片在线观看| 天堂网av新在线| 一个人看的www免费观看视频| 最好的美女福利视频网| 成年女人永久免费观看视频| 午夜福利在线观看吧| 三级毛片av免费| 中文在线观看免费www的网站| 天天一区二区日本电影三级| 男人舔奶头视频| 亚洲精品色激情综合| 亚洲av第一区精品v没综合| 国产成人av教育| 丰满人妻熟妇乱又伦精品不卡| 欧美精品啪啪一区二区三区| 国产三级在线视频| 国产午夜精品久久久久久| 午夜福利在线观看吧| 精品午夜福利视频在线观看一区| 国产精品野战在线观看| 色综合站精品国产| 亚洲国产欧美人成| 在线观看日韩欧美| 级片在线观看| 国产精品1区2区在线观看.| 日韩 欧美 亚洲 中文字幕| 少妇裸体淫交视频免费看高清| 国产亚洲精品久久久com| 国产高清videossex| 一个人看的www免费观看视频| 最新中文字幕久久久久 | 国产视频内射| 亚洲欧美日韩高清在线视频| 亚洲精品在线美女| 免费av毛片视频| 夜夜爽天天搞| 亚洲欧美激情综合另类| 成人一区二区视频在线观看| av在线蜜桃| 久久中文看片网| www.www免费av| 国产成人啪精品午夜网站| 不卡av一区二区三区| 中文字幕人成人乱码亚洲影| 久久久精品欧美日韩精品| 中文字幕精品亚洲无线码一区| 久久久久久久久免费视频了| 不卡一级毛片| 一区福利在线观看| 亚洲aⅴ乱码一区二区在线播放| 女人被狂操c到高潮| 日韩国内少妇激情av| 日本五十路高清| 黄色视频,在线免费观看| 动漫黄色视频在线观看| 大型黄色视频在线免费观看| 一个人免费在线观看电影 | 国产欧美日韩一区二区三| 欧美色视频一区免费| 日日摸夜夜添夜夜添小说| 国产三级在线视频| 69av精品久久久久久| 久久热在线av| 99久久久亚洲精品蜜臀av| 欧美日韩中文字幕国产精品一区二区三区| 女人被狂操c到高潮| 日本成人三级电影网站| 国产精品九九99| 少妇的丰满在线观看| 亚洲专区国产一区二区| www.999成人在线观看| aaaaa片日本免费| 日日摸夜夜添夜夜添小说| av天堂在线播放| 叶爱在线成人免费视频播放| 精品日产1卡2卡| 白带黄色成豆腐渣| 丰满的人妻完整版| 毛片女人毛片| 两个人看的免费小视频| 禁无遮挡网站| 亚洲欧美日韩高清在线视频| 三级国产精品欧美在线观看 | 久久久久久久久免费视频了| 午夜福利在线观看免费完整高清在 | 美女午夜性视频免费| 性色av乱码一区二区三区2| 亚洲avbb在线观看| 757午夜福利合集在线观看| 国产又色又爽无遮挡免费看| 国产成人精品久久二区二区91| 欧美成人免费av一区二区三区| 免费无遮挡裸体视频| 少妇丰满av| 亚洲国产欧美一区二区综合| 亚洲欧洲精品一区二区精品久久久| 国产午夜福利久久久久久| 久久久水蜜桃国产精品网| 18禁黄网站禁片午夜丰满| 久久精品国产综合久久久| 亚洲欧美日韩无卡精品| 久久久久久久精品吃奶| 特大巨黑吊av在线直播| 曰老女人黄片| 男女床上黄色一级片免费看| 国产成人av教育| 中文字幕av在线有码专区| 男人舔奶头视频| 日本a在线网址| 精品久久久久久成人av| 免费av不卡在线播放| 1024香蕉在线观看| 欧美黑人巨大hd| 伦理电影免费视频| 亚洲国产高清在线一区二区三| 亚洲精品色激情综合| 麻豆成人午夜福利视频| 国产激情欧美一区二区| 神马国产精品三级电影在线观看| 观看美女的网站| 欧美极品一区二区三区四区| 亚洲色图av天堂| 丰满的人妻完整版| 久久九九热精品免费| 精品久久久久久久久久免费视频| 老汉色∧v一级毛片| 免费大片18禁| 久久久久精品国产欧美久久久| 国产视频一区二区在线看| 琪琪午夜伦伦电影理论片6080| 两人在一起打扑克的视频| 岛国在线观看网站| 久久久成人免费电影| 欧美日韩福利视频一区二区| 久久久久免费精品人妻一区二区| 国产一区在线观看成人免费| 香蕉国产在线看| 国产伦精品一区二区三区四那| 国产野战对白在线观看| 这个男人来自地球电影免费观看| 美女被艹到高潮喷水动态| 精品国产三级普通话版| 久久久国产成人免费| 欧美日本视频| 久久久成人免费电影| 别揉我奶头~嗯~啊~动态视频| 久久精品亚洲精品国产色婷小说| 午夜福利高清视频| 亚洲熟妇熟女久久| 欧美色视频一区免费| 精品午夜福利视频在线观看一区| 亚洲 欧美 日韩 在线 免费| 午夜福利免费观看在线| 美女扒开内裤让男人捅视频| 不卡一级毛片| 最近视频中文字幕2019在线8| www.熟女人妻精品国产| 亚洲中文字幕一区二区三区有码在线看 | 欧美日韩乱码在线| 久久人妻av系列| 老熟妇仑乱视频hdxx| 搡老妇女老女人老熟妇| 中文在线观看免费www的网站| 18禁观看日本| 男女视频在线观看网站免费| 99久久精品一区二区三区| 一a级毛片在线观看| 欧美乱妇无乱码| 国产精品乱码一区二三区的特点| 久久久久久久久久黄片| xxxwww97欧美| 国产精品亚洲一级av第二区| 国产午夜精品久久久久久| 亚洲国产中文字幕在线视频| 12—13女人毛片做爰片一| 国产高清有码在线观看视频| 夜夜躁狠狠躁天天躁| 欧美av亚洲av综合av国产av| 桃红色精品国产亚洲av| 国产亚洲精品一区二区www| 久久久久久久久久黄片| 99久国产av精品| 久久香蕉精品热| 偷拍熟女少妇极品色| 国产在线精品亚洲第一网站| 人人妻,人人澡人人爽秒播| 十八禁人妻一区二区| 欧美日韩国产亚洲二区| 黄片小视频在线播放| bbb黄色大片| 特级一级黄色大片| 久久精品国产亚洲av香蕉五月| 亚洲人与动物交配视频| 又黄又爽又免费观看的视频| 男人舔女人下体高潮全视频| 成人性生交大片免费视频hd| 国产精品乱码一区二三区的特点| 午夜免费观看网址| 精品无人区乱码1区二区| 白带黄色成豆腐渣| 亚洲无线在线观看| 中文在线观看免费www的网站| av天堂在线播放| 18禁裸乳无遮挡免费网站照片| 国产精品野战在线观看| 黄色女人牲交| 亚洲人成伊人成综合网2020| 精品不卡国产一区二区三区| 午夜福利免费观看在线| 九色国产91popny在线| 婷婷精品国产亚洲av在线| 免费观看的影片在线观看| 国内精品一区二区在线观看| 国产精品九九99| 国产精品 欧美亚洲| 国产野战对白在线观看| 禁无遮挡网站| 岛国在线观看网站| 国产三级黄色录像| 国内久久婷婷六月综合欲色啪| 久久久久亚洲av毛片大全| 最近最新免费中文字幕在线| 黑人操中国人逼视频| 长腿黑丝高跟| 亚洲国产精品久久男人天堂| 亚洲真实伦在线观看| 久久久国产成人免费| 国产精品98久久久久久宅男小说| 天天添夜夜摸| 淫妇啪啪啪对白视频| 日韩成人在线观看一区二区三区| 一个人观看的视频www高清免费观看 | 亚洲av电影不卡..在线观看| 动漫黄色视频在线观看| 最近视频中文字幕2019在线8| 高清在线国产一区| 午夜福利18| 久9热在线精品视频| 欧美精品啪啪一区二区三区| 琪琪午夜伦伦电影理论片6080| 观看免费一级毛片| 日本与韩国留学比较| 久久九九热精品免费| 免费av不卡在线播放| 亚洲精品中文字幕一二三四区| 小蜜桃在线观看免费完整版高清| 日本 av在线| xxxwww97欧美| 国产乱人视频| 欧美成狂野欧美在线观看| 亚洲中文字幕日韩| 天天躁日日操中文字幕| 91老司机精品| 国产精品99久久久久久久久| 亚洲国产欧美人成| 国产在线精品亚洲第一网站| 国产高清激情床上av| 亚洲欧美精品综合久久99| 国产精品久久久久久亚洲av鲁大| 青草久久国产| 欧美日韩综合久久久久久 | 变态另类丝袜制服| 色视频www国产| 不卡一级毛片| 久久精品91无色码中文字幕| 老司机深夜福利视频在线观看| 亚洲狠狠婷婷综合久久图片| 美女被艹到高潮喷水动态| 欧美在线黄色| 很黄的视频免费| 亚洲成人久久性| 国产精华一区二区三区| 老熟妇乱子伦视频在线观看| 国产精品av久久久久免费| 国产成人精品久久二区二区免费| 亚洲18禁久久av| 久久午夜综合久久蜜桃| 啦啦啦免费观看视频1| 国产亚洲欧美98| 国产三级中文精品| 真人一进一出gif抽搐免费| 国产视频一区二区在线看| cao死你这个sao货| 久久精品亚洲精品国产色婷小说| 视频区欧美日本亚洲| 久久久久久久久免费视频了| 在线永久观看黄色视频| 色综合欧美亚洲国产小说| 亚洲乱码一区二区免费版| 欧美成狂野欧美在线观看| 日韩中文字幕欧美一区二区| 嫩草影院入口| 日本熟妇午夜| 制服丝袜大香蕉在线| 国产精品野战在线观看| 男女床上黄色一级片免费看| 亚洲欧洲精品一区二区精品久久久| 国产精品久久久久久精品电影| 成人亚洲精品av一区二区| 美女大奶头视频| 国内精品久久久久久久电影| 大型黄色视频在线免费观看| 久久久色成人| 国产午夜福利久久久久久| 亚洲精品久久国产高清桃花| 国产精品国产高清国产av| 岛国在线免费视频观看| 又黄又粗又硬又大视频| 精品久久久久久久毛片微露脸| 免费在线观看日本一区| 亚洲激情在线av| 成人18禁在线播放| 男女做爰动态图高潮gif福利片| 成熟少妇高潮喷水视频| 亚洲精品一区av在线观看| 很黄的视频免费| 午夜免费成人在线视频| av片东京热男人的天堂| 婷婷精品国产亚洲av在线| 欧美色欧美亚洲另类二区| 好男人电影高清在线观看| 青草久久国产| 亚洲一区二区三区色噜噜| 又大又爽又粗| 亚洲人成网站在线播放欧美日韩| 欧美xxxx黑人xx丫x性爽| 日韩中文字幕欧美一区二区| 最近最新中文字幕大全免费视频| 窝窝影院91人妻| 国内揄拍国产精品人妻在线| 很黄的视频免费| 嫁个100分男人电影在线观看| 中文字幕精品亚洲无线码一区| 伦理电影免费视频| 久久久久久九九精品二区国产| 日韩三级视频一区二区三区| www.熟女人妻精品国产| 国产免费av片在线观看野外av| 午夜久久久久精精品| 香蕉av资源在线| 99久久成人亚洲精品观看| 在线永久观看黄色视频| 99久久成人亚洲精品观看| 国产精品,欧美在线| 观看免费一级毛片| 国产高清视频在线观看网站| 国产成人av激情在线播放| 欧美黑人欧美精品刺激| 亚洲在线自拍视频| АⅤ资源中文在线天堂| 一级毛片精品| 亚洲av成人av| 99热精品在线国产| 日韩三级视频一区二区三区| 国产1区2区3区精品| 国产精品九九99| 色综合婷婷激情| 亚洲精品美女久久久久99蜜臀| 高潮久久久久久久久久久不卡| 免费电影在线观看免费观看| 级片在线观看| 亚洲国产欧洲综合997久久,| 两个人看的免费小视频| 亚洲精品美女久久久久99蜜臀| 午夜精品一区二区三区免费看| 少妇丰满av| 99久久成人亚洲精品观看| 亚洲国产看品久久| 亚洲第一欧美日韩一区二区三区| 美女免费视频网站| 久久伊人香网站| 亚洲激情在线av| 久久热在线av| 亚洲精品国产精品久久久不卡| 日韩欧美在线乱码| 欧美一区二区国产精品久久精品| 亚洲av片天天在线观看| 国产精品日韩av在线免费观看| 国产精品亚洲美女久久久| 国产伦人伦偷精品视频| 老熟妇乱子伦视频在线观看| 性色avwww在线观看| 成年女人看的毛片在线观看| 欧美成狂野欧美在线观看| 我的老师免费观看完整版| 桃红色精品国产亚洲av| 久久久久久九九精品二区国产| x7x7x7水蜜桃| 日韩 欧美 亚洲 中文字幕| 丰满人妻熟妇乱又伦精品不卡| 色吧在线观看| 成人欧美大片| 在线观看美女被高潮喷水网站 | 99热只有精品国产| 黑人欧美特级aaaaaa片| 日韩三级视频一区二区三区| 午夜日韩欧美国产| 亚洲无线在线观看| 黑人操中国人逼视频| 国产高清有码在线观看视频| 三级男女做爰猛烈吃奶摸视频| 久久久色成人| 国产三级在线视频| 成人亚洲精品av一区二区| 黄片小视频在线播放| 一级毛片精品| 亚洲av电影在线进入|