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

    地下水對氣壓和固體潮響應研究進展

    2021-09-23 16:18:26何冠儒史浙明
    地震研究 2021年4期
    關鍵詞:固體潮大氣壓力潮汐

    何冠儒 史浙明

    摘要:分析地下水對氣壓和固體潮的響應可以獲得多項水文地質參數(shù)和含水系統(tǒng)特性指標。通過收集國內外相關研究,系統(tǒng)論述了基于地下水對大氣壓和固體潮響應來獲取含水層系統(tǒng)參數(shù)的方法。介紹了傳遞函數(shù)計算、諧波分析及回歸反卷積等計算氣壓效率的方法;介紹了利用氣壓效率與固體潮響應模型計算含水層水力參數(shù)的方法。結果表明:相比于傳統(tǒng)水文地質調查手段,地下水氣壓和固體潮響應方法更簡便高效,能夠深入了解水力特性的空間和時間變化規(guī)律,可用于水位校正、含水層性質評估、水力參數(shù)計算以及評估地震對含水層造成的影響等。

    關鍵詞:地下水位;氣壓響應;固體潮;水力特性

    中圖分類號:P315.72? 文獻標識碼:A文章編號:1000-0666(2021)04-0541-09

    0引言

    研究地下水位的動態(tài)變化對于深入了解地下水系統(tǒng)特性和資源的開發(fā)利用具有重要意義。地下水位變化受到多種因素影響,其中大氣壓力和固體潮對地下水的作用是最為常見的現(xiàn)象(車用太,1984;張邵棟等,1989)。固體潮是太陽或月球對地球引力產生周期性變化,令固體地球發(fā)生周期形變的現(xiàn)象,從而引起觀測井中水位產生周期性波動(Doan et al ,2006)。同樣,大氣壓力的波動也會使井水位發(fā)生變化。由于晝夜溫度變化,氣體產生周期性的膨脹,此時大氣壓力載荷從地表傳遞到隔水層,部分荷載由含水層孔隙水承受,部分由含水層基質承受,造成含水層與井筒之間存在壓力不平衡,從而引起井孔與含水層之間的地下水流動,導致井水位出現(xiàn)相應的變化(Tur-nadge et al ,2019)。通過評估地下水對固體潮和氣壓的響應,結合多孔彈性介質模型,可以確定地下水的承壓性質,量化水力特性參數(shù),如滲透系數(shù)、導水率、儲水率、孔隙度、地層壓縮系數(shù)等(MeMillan et al, 2019)

    由于傳統(tǒng)的水文地質參數(shù)獲取常常借助于抽水試驗、注水試驗等現(xiàn)場測試手段,這些方法往往耗時耗力且無法獲取連續(xù)的參數(shù)變化情況。而基于地下水位對氣壓和固體潮響應的分析則可以更便捷、高效地獲取水文地質參數(shù)信息及其動態(tài)變化,與傳統(tǒng)水文地質方法相比優(yōu)勢巨大。針對地下水對大氣壓力和固體潮的響應,本文總結了具體的分析方法、流程以及相應的應用,并對未來該方法的推廣和完善提出了建議。

    1氣壓響應分析方法

    隨著觀測手段的發(fā)展,水位、氣壓等已經實現(xiàn)高精度、自動化觀測。水位與大氣壓力信號之間的關系可用氣壓效率(BE)和荷載效率反映。氣壓效率指得是地下水位壓力變化與大氣壓力變化之比,當前其計算方法主要有5種時域的方法、3種頻域的方法(Turnadge et al ,2019)。時域方法是將氣壓的變化與地下水壓力的變化在時間域上聯(lián)系起來,較為容易實現(xiàn),但結果也容易受到其他因素干擾(如降水、蒸散等)的影響。而頻域分析方法則將不同信號從時間域轉換到頻率域,可以有效分離不同信號,不僅能得到響應的振幅,還可以得到不同信號響應間的相位滯后。

    Bredehoeft(1967)指出,井水位所有潮汐諧波分量中M2、S2、N2、K1和O1占潮汐勢的95%,為主要分析的頻率。Cutillo和Bredehoeft(2011)指出K1波易受蒸散干擾,N2(1.895982cpd)波振幅很小可以忽略,主要解析O1(0.929536cpd)和M2(1.932274cpd)的響應振幅。然而Acworth (2016)認為固體潮引起的S2波對氣壓引起的S2波會產生干擾,由此提出利用M2波校正來消除固體潮相應波段影響,從而獲得更為準確的氣壓效率這一方法。

    目前計算氣壓效率的主流方法有:傳遞函數(shù)法Quilty, Roeloffs, 1991; Rojstaczer, 1988a, b 諧波調和法(Acworth et al ,2016;rau et al ,2020)和回歸反卷積法(Butler et al ,2011;Rasmussen ,Crawford, 1997: Toll, Rasmussen, 2007)

    1.1傳遞函數(shù)法

    水位對氣壓的傳遞函數(shù)可用于擬合氣壓響應函數(shù),其振幅可作為氣壓效率。計算傳遞函數(shù)的步驟為:①確定水位、氣壓和固體潮信號的自功率譜和交叉譜;②求解下列復線性方程組(Ro-jstaczer, 1988a)

    式中:BB和TT分別表示大氣壓力和固體潮的自功率譜;BT和TB分別表示大氣載荷和固體潮之間的交叉譜和它的復共軛;BW表示大氣載荷和水位的交叉譜;TW表示固體潮和水位之間的交叉譜;HB和HT分別表示水位與大氣載荷和水位與固體潮之間的傳遞函數(shù)。

    去除固體潮效應,計算傳遞函數(shù)的方法可以簡化為:

    振幅和相位可以表示為(Hussein et al ,2013):

    在使用傳遞函數(shù)時,往往需要使用濾波器對數(shù)據(jù)進行預處理,常用的有帶通、低通、高通等濾波器,濾波器的截止頻率與采樣率有關,可以過濾掉3cpd以上的高頻噪音,或使用逐漸變小的漢寧窗,分成50%重疊率的頻帶進行平均等手段(Hussein et al ,2013)。

    1.2諧波調和法

    Acworth 和Brain (2008)通過地下水和氣壓的振幅比估計氣壓效率:

    式中:A(W)為水位波動的振幅;A(P)為氣壓波動振幅。

    但是在進行諧波調和分析計算氣壓效率時也需要識別并剔除固體潮響應,為此Acworth 等(2015)提出利用離散傅立葉變換可以提取顯著頻率的振幅,并指示地下水對固體潮的響應。水位信號中S2分量頻率(2cpd)以下的響應包括對大氣潮的響應信號Z和對固體潮響應信號z-,即水頭響應信號為:

    由于固體潮的M2和S2分量的振幅比是恒定的,水位對M2分量的響應不會受大氣壓力的干擾,于是Z可由Z計算得出:

    氣壓效率可表示為:

    在地下水存在補給、排泄或氣壓效率較低的系統(tǒng)中,氣壓波動與水位波動沒有顯著的相關性。為了更好地分解氣壓信號和固體潮信號,Acworth 等(2016)采用諧波加和定理實現(xiàn)了突破,提出了一種利用消除固體潮響應來計算氣壓效率的新方法,指出響應振幅取決于兩個驅動因素之間的相位差,只有當相位差等于0時,振幅才可以做簡單的加減。

    為了準確量化固體潮分量對地下水頭的影響,可以通過計算固體潮在S2分量下的振幅,用S2分量下大氣潮和固體潮之間的相位差來校正,氣壓效率的計算公式為:

    式中:S2和M2分量的振幅,可采用快速傅立葉變換獲得(Acworth et al ,2016)。

    Rau等(2020)指出,Acworth 等(2016)的方法是基于鉆孔水位代表含水層孔隙壓力的假設,即存在瞬時和無阻尼響應。在這種情況下,僅需要理論固體潮和氣壓潮驅動之間的相位差來校正地下水響應振幅。然而Hsieh 等(1987)指出,固體潮M2分量和井水位響應之間存在相位差,在利用地下水位對氣壓潮的響應中量化BE時,還必須考慮相位延遲(Rau et al ,2020)。

    考慮地下水位與孔隙壓力的關系,阻尼大小可以用它們的振幅比A計算,方法采用Hsieh 等(1987)和Rojstaczer(1988b)的理論公式,即:

    式中:Z為孔隙壓力對S2波響應信號;f為S2波的頻率;T為導水率;S為儲水率;rn為井半徑;rn為套管半徑。

    振幅和相位變化都受到水力傳導系數(shù)的影響,對于承壓層A>0.99,即滲透系數(shù)K>1×10-m/s,阻尼造成的誤差可以忽略,然而在滲透系數(shù)較小的條件下A顯著降低,計算氣壓效率時便需要考慮地下水位與孔隙壓力的關系,所以更客觀的氣壓效率的計算方法如下:

    值得注意的是,Rau等(2020)在提出諧波信號時采用了諧波最小二乘法,計算公式為:

    式中:N是離散樣本數(shù);yn是時間tn的樣本值;C是頻率為f的潮汐分量C的總數(shù),通過最小二乘法擬合出參數(shù)an和b作為復數(shù)信號的實部和虛部。

    Schweizer 等(2021)比較了離散傅立葉變換和諧波最小二乘法提取諧波分量特性的能力,通過驗證兩個具有不同特征的真實數(shù)據(jù)集的模型,發(fā)現(xiàn)Rau等(2020)使用的諧波最小二乘法在各方面都優(yōu)于前人使用的離散傅立葉變換。諧波最小二乘法在未來也會得到更多應用。

    1.3回歸反卷積法

    由于承壓層之上存在包氣帶,氣壓傳遞會受到氣體擴散系數(shù)的影響,從而造成氣壓效應的滯后;在較厚或低滲透性的區(qū)域,包氣帶會導致周期性信號產生較大的衰減和相移。針對這種響應的時間滯后,Rasmussen 和Crawford (1997)根據(jù)Furbish 早期的研究,提出了利用回歸反卷積方法來計算地下水對大氣壓力的響應,氣壓效率可以表示為延遲響應的函數(shù),可以用如下的卷積計算方法來計算氣壓響應函數(shù)(BRF)(Butler et al ,2011):

    式中:△w(t)為無趨勢水位;a1和B1為滯后i時的大氣壓力和固體潮的單位脈沖,通過最小二乘回歸可以計算a1、β;j為觀測時間。

    Butler 等(2011)提出某一口井的氣壓響應函數(shù)并不是一成不變的,其形式可能取決于穿過滲流區(qū)的壓力傳播的性質。對于承壓含水層,短期內(1h)大氣壓力荷載由孔隙水承擔,含水層對大氣壓力的變化不敏感,BRF(j)數(shù)值較低,隨著壓力波通過不飽和帶,井水頭與大氣壓力變化達到平衡狀態(tài),單位脈沖響應a1逐漸減弱,響應達到穩(wěn)定。

    例如云南會澤井的氣壓響應函數(shù)(圖1a),隨著滯后時間的增大逐步穩(wěn)定,反映所處的地層可能為承壓含水層。當在潛水或半承壓含水層時,井水位和孔隙水同時受到大氣壓中短周期信號產生的響應,隨后大氣壓力被孔隙水分擔一部分,導致井水位的延遲響應越來越小,在壓力平衡后達到穩(wěn)定。這類氣壓響應函數(shù)的特征表明含水層可能屬于半承壓或未承壓的狀態(tài),如新10井(圖1b)。

    2地下水潮汐響應模型

    眾所周知,固體潮可以根據(jù)已知的天文參數(shù)進行精確計算(Rau et al ,2020)。通過分析地下水位與固體潮、氣壓的響應信號,結合多孔彈性介質理論模型(McMillan et al ,2019),可以挖掘出地下水位記錄的豐富地質與地球物理信息。例如水力擴散率、導水率和儲存系數(shù),降低含水層抽水試驗分析中固有的較大不確定性(MeMillan et al,2019)。常用的多孔介質模型可分為水平流模型、垂直流模型和混合模型。

    2.1水平流模型

    Jacob (1940)通過建立一個理想均勻邊界的彈性承壓含水層,說明了氣壓和潮汐波動可以作為量化含水層彈性的指標,潮汐效率(荷載效率)或氣壓效率與儲水率之間有直接關系。

    氣壓效率為:

    潮汐效率為:

    式中:θ為孔隙度;K1為水的彈性模量;K為地層的彈性模量;ρ為水的密度,根據(jù)泰斯的儲水率計算公式,儲水率可由氣壓效率計算:

    Bredehoeft(1967)提出水位對固體潮響應的振幅可以用于量化水文地質參數(shù)。假設忽略固體的可壓縮性,近地表含水層的垂直應變取決于含水層的泊松比,含水層中鉆孔水位變化和潮汐擾動的振幅變化之比是一個與儲水率有關的定值:

    式中:dH為水位變化;dW2為太陽或月球的干擾潮汐;S為儲水率;v為泊松比;h、l為洛夫數(shù);a為地球半徑;g為重力加速度。

    當觀測到的地球潮汐振幅有超過1mm左右的水位變化時,那么地層一定是相對不可壓縮的,固體的可壓縮性是不可忽略的。所以Kamp和Gale(1983)考慮了固體可壓縮性的影響,導出了通過固體潮和大氣壓力響應計算儲水率的公式。如果固體是可壓縮的,儲水率可以寫成:

    式中:y為荷載效率;G為剛性模量;a為Biot-Willis 常數(shù)。

    Hsieh 等(1987)通過理想承壓含水層模型和開爾文函數(shù)解出了井水位對固體潮響應的解析解,發(fā)現(xiàn)振幅和相位是井半徑、儲水率和導水率的函數(shù),為水文地質學家通過固體潮波動下水位記錄計算含水層導水率提供了可能。

    潮汐引起的水位波動振幅為:

    相位為:

    式中:x為水位波動振幅;s為體應變波動的振幅;E、F由開爾文函數(shù)計算:

    式中:φ、φ由an的一階開爾文函數(shù)計算;Ker、Kei為零階開爾文函數(shù)的實部和虛部;an由儲水率、導水率、頻率、井徑計算得出:

    在實際應用中可以通過原位應變儀測得潮汐作用引起的體應變,也可以通過理論固體潮計算得到,通過對水位和體應變進行傅立葉變換得到振幅和相位,再結合上述公式,可獲得含水層的導水系數(shù)和儲水率。

    2.2垂直流模型

    Roeloffs(1996)提出,在不排水條件下,假設忽略水流橫向流動,只考慮垂向擴散,當氣體能夠擴散到整個潛水含水層時,氣壓效率為井水深和水力擴散系數(shù)的函數(shù),可表示為:

    式中:z為地表到含水層的距離;zn為地表到井內水面的距離;c為水力擴散系數(shù);vn為未排水地層的泊松比。

    孔隙壓與固體潮體應變的振幅比可表示為(Roeloffs, 1996)

    式中:po(z)為深度z時的孔隙壓力;e為固體潮體應變振幅。

    Roeloffs(1996)指出,對于特定的潮汐分量主頻率,ω(z-z)2/c可以用于反映含水層的承壓性,當ω(z-z)2/c<0.2時,如果潮汐應變產生的擾動小于1/e,那么可以指示該含水層具有良好的承壓性。

    2.3混合模型

    Rojstaczer(1988a)提出簡化井-含水層模型,考慮氣壓變化的3個流動因素:①從地表到地下水位的非飽和帶中的垂直氣流;②承壓層內的垂直地下水水流;③含水層和井孔之間的水平地下水水流。井水位對大氣壓力的響應大小取決于井的半徑、含水層的水平向水力擴散率、包氣帶厚度和垂直氣動擴散率以及潛水層(或部分承壓層)垂直水力擴散率5個參數(shù)。這些水文地質參數(shù)與響應頻率的關系如下。

    在地表到地下水位的氣流交換中,定義反映非飽和帶厚度和氣體擴散關系的頻率參數(shù)R為:

    式中:L為地表到地下水位的距離;Dn為氣體擴散系數(shù)。

    在潛水面到含水層間的垂直水流交換中,定義反映潛水層厚度和水力擴散系數(shù)關系的頻率參數(shù)Q為:

    式中:b為潛水面到承壓層上頂板的距離;D為垂直水力擴散系數(shù)。

    在鉆孔與含水層之間的水流交換中,定義反映井徑與滲透系數(shù)關系的頻率參數(shù)W為:

    式中:r為井半徑;K為含水層滲透系數(shù);b為含水層厚度。

    在一般情況下,井水位對大氣荷載的響應可以表示為:

    式中:x為水位信號的振幅;A為大氣壓力信號的振幅;p為含水層中孔隙壓力;s為井內水位降深。式(32)描述了井水位在頻率域中的響應,井筒內水位的變化加上大氣負荷(等效水位變化),等于遠場孔隙壓力(等效水位變化)減去井水位降深。

    氣壓效率為響應頻率的函數(shù),氣壓效率可計算表示為響應函數(shù)的增益:

    相位計算如下:

    Rojstaczer和Agnew (1989)隨后的研究中利用井水位對固體潮的響應來獲取地殼應變,發(fā)現(xiàn)井水位對大氣載荷和面應變響應的理論模型可以與交叉譜估計結合使用,進而估算含水層的彈性參數(shù),如骨架壓縮系數(shù)、孔隙率、儲水率。水井對荷載和固體骨架應變的靈敏度和地層壓縮性可表示為:

    式中:w為水位變化;E,為固體骨架體應變;B為Skempton系數(shù);B地層飽水時的壓縮系數(shù)。

    孔隙度可表示為:

    式中:β1為液體壓縮系數(shù);Bn為固體壓縮系數(shù);B由面應變靈敏度和荷載效率計算得出。

    由于含水層往往介于完全承壓與完全不承壓之間,Wang等(2018)提出了更具普適性的滲流模型,即考慮地下水水平流動和垂向的越流,在這種情況下,潮汐響應的振幅比為:

    相位滯后表示為:

    式中:hn表示地下水的振幅,可通過修正Bessel 函數(shù)和Laplce變化得到:

    式中:S為儲水率;k為弱透水層滲透系數(shù);b為弱透水層厚度;ξ為導水率,并徑的Bessel 函數(shù)。

    通過上述公式可以得到潮汐響應振幅、相移與導水系數(shù)、儲水率和越流系數(shù)的關系。

    3應用

    3.1水位校正

    地下水位動態(tài)經常受氣壓和固體潮變化的影響。這些干擾會對后續(xù)地下水位的動態(tài)分析造成誤差(王麗亞等,2012)。即使在低水力梯度的包氣帶或非承壓含水層中,氣壓的波動都會使地下水流動產生瞬時變化,從而影響地下水位測量的準確性。因此,在地下水流的研究中應考慮氣壓和固體潮的影響。利用回歸反卷積法從水位中去除氣壓和固體潮效應在水位校正中已經進行了廣泛應用,并取得了較好的效果(Rasmussen ,Crawford ,1997;Spane,2002; Toll, Rasmussen, 2007)

    3.2識別含水層特性

    Rahi和Halihan(2013)提出利用頻譜分析的方法可以判斷含水層的類型:①非承壓含水層的地下水信號存在S2和K1分量,但缺失M2分量;②半承壓含水層地下水信號中S2分量占主導,同時存在M2分量;③承壓含水層中則有明顯的M2分量。盡管Rahi和Halihan(2013)沒有給出定量的標準,但后續(xù)研究表明:當水位M2分量的振幅大于1mm時就可以指示含水層具備較好的承壓性(Turnadge et al ,2019)。前人的多個研究結果表明,諧波分析法是確定含水層類型的一種準確、低成本的方法(Hussein et al ,2013;rahi,Hali-han, 2013: Acworth et al, 2015)

    Acworth 等(2017)指出,當諧波分量S2的地下水位與大氣壓力的相位差達到-180°時可以說明含水層具備良好的承壓性。該研究指出地下水位的氣壓響應可以作為一種天然的示蹤劑,有助于理解地下水系統(tǒng)的時空特征和概念,量化地質單元中多個含水層的承壓狀態(tài)和壓縮儲存特性。

    3.3計算含水層水力參數(shù)

    前人通過地下水位對潮汐和氣壓的響應模型推導了許多水力參數(shù)的計算方法(Bredehoeft,1967; Hsieh et al, 1987; Jacob, 1940; Kamp Gale, 1983; Roeloffs, 1996; Rojstaczer, Agnew 1989)。在后續(xù)的改良與應用中,Hussein 等(2013)通過傳遞函數(shù)方法計算了含水層多個水文地質參數(shù),并且根據(jù)結果提出了含水層脆弱性的評價方法,使地下水的氣壓響應方法成為半承壓含水層脆弱性評估的有效工具。Xue等(2013)通過水位對氣壓潮響應來研究斷裂帶的水文地質結構,估算了斷裂帶的儲水量和滲透率。Sun和Xiang (2019)等運用頻譜分析,氣壓響應函數(shù)等手段計算了會澤井井水位對周期性荷載的響應靈敏度,說明響應靈敏度取決于其井結構和含水層滲透性。Sun和Xiang (2020)等通過比較不同模型,得出在氣壓荷載低頻區(qū),忽略非飽和帶影響的模型擬合結果更優(yōu),特別是識別淺層地殼應變或應力的井。Qu等(2020a)研究了煤礦開采區(qū)地下水位對大氣壓力的響應,計算了煤礦上覆含水層的水力特性,證明了潮汐效應計算的水力特性有效半徑與抽水試驗相似。

    基于地下水對氣壓響應的分析可以求解礦區(qū)、斷裂帶等特殊地層的水力特性參數(shù),并且隨著模型復雜性的提升,反演結果也更加客觀。未來研究中可以對水力特性的有效范圍開展進一步研究,探索其尺度效應的變化。

    3.4分析地震對含水層介質的影響

    Rau等(2018)提出了一種基于耦合跨孔地震探測與地下水位氣壓效應來估算含水層儲水率的新方法。該方法考慮了固體顆粒與骨架壓縮性對儲水率計算結果的影響,相比傳統(tǒng)方法提高了計算的精度。Sun和Xiang (2019)利用井水位對大氣壓力的響應計算了地震對含水層系統(tǒng)垂向和水平向滲透率的影響,發(fā)現(xiàn)地震造成了地層滲透性的增大。地下水對潮汐和氣壓響應方法為研究淺層地殼介質滲透率變化及了解其各向異性分析提供一種新的途徑。Zhang 等(2019)利用地下水位對大氣壓力和固體潮的響應,研究了地震對同一含水層系統(tǒng)中含水層和弱透水層滲透率和導水率的影響。說明地下水對氣壓響應分析同樣適用于弱透水層的水力參數(shù)計算。

    利用地下水對潮汐和氣壓的響應,結合地震分析可以改進水文地質概念模型(孫小龍等,2020a,b),加深對含水層斷裂結構和形變的認識。通過分析地震前后地下水位相移和振幅的變化也為監(jiān)測地震地下水前兆提供新的思路。

    4研究展望

    地下水對潮汐及大氣壓力響應的研究可以從以下兩個角度進一步開展:

    (1)分析方法研究。如結合其它領域理論對已有的計算水位地質參數(shù)的方法進行提升(Rau et al,2018)。由于實際地質條件的復雜性,現(xiàn)有的方法均進行了不同程度的簡化,因此現(xiàn)存的不同方法獲得的結果可能存在一定差異,可以探索結果差異的原因以及改良的方法(Rau et al ,2020);對現(xiàn)有模型進行比較和細分,尋找更適合野外實際的概念模型(Sun,Xiang,2020)。

    (2)模型拓展應用研究。如地震、采礦、長期地下水超采等對地下水系統(tǒng)的影響(Qu et al ,2020b; Sun, Xiang, 2019; Zhang et al, 2019)在地下水污染評價中,對半承壓含水層脆弱性進行評估(Hussein et al ,2013;Odling et al ,2015)。隨著分析方法和理論模型的不斷成熟,地下水的潮汐和氣壓響應可以和更多領域進行交叉研究,應用于地面沉降、資源開采評價等方面。

    參考文獻:

    車用太。1984。井孔水位的微動態(tài)特征綜述[J].水文地質工程地質,(4):18-22,42.

    孫小龍,劉耀煒,付虹,等。2020a.我國地震地下流體學科分析預報研究進展回顧[J].地震研究,43(2):216-231。

    孫小龍,向陽,李源。2020b.深井水位對地震波、固體潮和氣壓的水力響應——以范縣井為例[J].地震學報,42(6):719-731。

    王麗亞,郭海朋,李文鵬,等。2012。氣壓對觀測井水位的影響及校正方法[J].水文地質工程地質,39(6):29-34。

    張昭棟,鄭金涵,張廣城,等。1989。承壓井水位對氣壓動態(tài)過程的響應[J].地球物理學報,(5):539-549。

    Acworth rI. Brain T.2008. Calculation of barometric efficiency in shal-low piezometers using water levels, atmospherc and earth tide data LJ]. Hydrogeology Journal, 16:1469-1481.

    Acworth RI, Halloran L, Rau GC, et al.2016. An objective frequency do-main method for quantifying confined aquifer compressible storage u-sing Earth and atmospheric tides[] . Geophysieal Research Letters 43(11):671-611,678.

    Acworth RI, Rau GC, Halloran L et al.2017. Vertical groundwater stor-age properties and changes in confinement determined using hydrau-lie head response to atmospherietides.[ J]. Water Resources Re-search ,53(4):2983-2997.

    Acworth RI. Rau GC, MeCallum A M, et al.2015. Understanding con-nected surface-water/ groundwater systems using Fourier analysis of daily and sub-daily head fluctuations.[ J]. Hydrogeology Journal 23:143-159.

    BredehoeftJD.1967. Response of well-aquifer systems to Earth tides.[]Journal of Geophysical Research, 72(12): 3075-3087.

    Butler JJ, Jin W, Mohammed GA, et al.2011. New insights from well re-sponses to fluctuations in barometric pressure[ J]. Ground Water, 49(4):525-533.

    Cutillo PA, Bredehoeft J D.2011. Estimating Aquifer Properties from the Water Level Response to Earth Tides.[ ]Ground Water,49(4):600-610.

    Doan ML, Brodsky EE ,Prioul R, et al.2006. Tidal analysis of borehole pressure -A tutoria .[EB/OL].(2006-12-20)[2021-05-17]Santa Cruz Romain Prioul: University of Califomia. websites.pme. ussc. edw/_ motoan/ articles/ tidcu _ tutorial_ SDR. pdf.

    Hsieh PA, BredehoeftJD, FarrJM.1987. Determination of aquifer trans-missivity from Earth tide analysis[ J]. Water Resources Research, 23(10):1824-1832.

    Hussein M, Odling NE, Clark RA.2013. Borehole water level response to barometric pressure as an indicator of aquifer vulnerability [ J]. Wa-ter Resources Research, 49(10): 7102-7119.

    Jacob CE.1940. On the fow of waler in an elastie artesian aquifer[]Eos Trans Agu, 21(2): 574-586

    Kamp G, Gale JE.1983. Theory of earth tide and barometrie effects in porous formations with compressible grains [ J]. Water Resources Re-search ,19(2):538-544.

    MeMillan TC, Rau G C, Timms WA, et al.2019. Utilizing the impact of earth and atmospheric tides on groundwater systems: A review reveals the future potential[ J]. Reviews of Geophysics, 57(2): 281-315.

    Odling NE, Serrano RP, Hussein MEA, et al.2015. Detecting the vul-nerability of groundwater in semi-confined aquifers using baromet-rie response functions[ J]. Journal of Hydrology, 520:143-156

    Qu S, Shi Z, Wang G, et al 2020a. Using water-level fluctuations in re-sponse to Earth-tide and barometrie - pressure changes to measure the in - situ hydrogeological properties of an overburden aquifer in a coalfield[ J]. Hydrogeology Journal, 28:1-15

    Qu S, Wang G, Shi Z, et al 2020b. Temporal changes of hydraulic proper-ties of overburden aquifer induced by longwall mining in Ningtiaota coalfield, northwest China[ J]. Journal of Hydrology, 582:124525

    Quilty E G, Roeloffs E A.1991. Removal of barometrie pressure re-sponse from water level data[ J]. Joumnal of Geophysical Research Solid Earth and Planets 96( B6): 10209-10218.

    Rasmussen TC, Crawford LA.1997. Identifying and removing barometric pressure effects in confined and unconfined aquifers [ Ground Water ,35:1109-1109.

    Rahi KA, Halihan T.2013. Identifying aquifer type in fractured rock aq-uifers using harmonic analysis[ J]. Ground water, 51(1): 76-82.

    Rau GC, Acworth R I, Halloran L JS, et al.2018. Quantifying compres-sible groundwater storage by combining cross-hole seismic surveys and head response to atmospherie tides [ j]. Journal of Geophysical Research-Earth Surface, 123(8): 1910-1930

    Rau g c. Cuthbert m o. Acworth rlet a.2020. Technical note. Disen-tangling the groundwater response to Earth and atmospherie tides to improve subsurface characterisation[] . Hydrology and Earth System Sciences ,24(12):6033-6046.

    Roeloffs E A.1996. Poroelastic techniques in the study of earthquake related hydrologic phenomena?[ J. Advances in Geophysics, 37:135-195

    Rojstaezer S, Agnew D C.1989. The influence of formation material prop-erties on the response of waler levels in wells to earth tides and at-mospherie loading [ J]. Journal of Geophysical Research - Solid Earth and Planets, 94( B9): 12403-12411

    Rojstaezer S.1988a. Determination of fluid flow properties from the re-sponse of water levels in wells to atmospheric loadingL]. Water Re-sources Research, 24(11): 1927-1938.

    Rojstaczer S.1988b. Intermediate Period Response of Water Levels in Wells to Crustal Strain - Sensitivity And Noise-Level[J].Journal of Geophysical Research - Solid Earth and Planets, 931(B11):13619-13634

    Rahi KA, Halihan T.2013. Identifying aquifer type in fractured rock ag-uifers using harmonic analysis[ J]. Ground water, 5(1): 76-82.

    Schweizer D, Ried V, Rau GC, et al.2021. Comparing methods and defi-ning practical requirements for extracting harmonic tidal components from groundwater level measurements [ J]. Mathematical Geosciences ,doi:10.1007/s11004-020-09915-9.

    Spane FA.2002. Considering barometric pressure in groundwater flow in-vestigations[ J]. Water Resources Research, 38(6): 141-148

    Sun x, Xiang Y.2019. Heterogeneous permeability changes along a fault zone caused by the Xingwen M5.7 Earthquake in SW China[ J]. Ge-ophysical Research Letters, 46(24): 14404-14411.

    Sun X, Xiang Y.2020. Comparison of transfer function models for well aquifer system response to atmospheric loading [ J]. Journal of Hy-drology ,590:125494.

    Toll NJ, Rasmussen TC.2007. Removal of barometrie pressure effects and earth tides from observed water levels.[ J. Ground Water, 45(1):101-105.

    Turnadge C, Crosbie RS, Barron O.2019. Comparing methods of baromel-rieeffieiency characterization for specific storage estimation [ J]Groundwater, 57(6):844-859

    Wang CY, Doan ML, Xue L, et a.2018. Tidal response of groundwater in a leaky aquiferApplication to Oklahoma [ J]. Water Resour Res,54(10):8019-8033.

    Xue L, Li H, Brodsky E E, et al.2013, Continuous permeability measure-ments record healing inside the Wenchuan Earthquake fault zone [J].Science ,340(6140):1555-1559.

    Zhang H, Shi Z, Wang G, et al.2019. Large Earthquake reshapes the groundwater flow System: Insight from the water-level response to earth tides and atmospheric pressure in a deep wellJ]. Water Re-sources Research 55(7): 4207-4219.

    Advance in the Groundwater Level Response to Barometric Pressure and earth tide

    HE Guanru. SHI Zheming

    School of Water Resources and Environment, China University of Geosciences, Beijing 100083, China)

    Abstract

    Numerous observations and studies have shown that groundwater can record clear signals from atmospheric pressure and Earth tide. Analysis of groundwater response to atmospheric pressure and Earth tide will help infer hydrogeological parameters and characterize the aquifer system. In this paper we reviewed the methods of obtaining aquifer system parameters based on groundwater response to atmospheric pressure and Earth tide, introduced meth ods of calculating barometric efficiency such as transfer function calculation, harmonic analysis and regression de-convolution. We further reviewed the method of inferring aquifer hydraulic parameters by groundwater level in re-sponse to barometric efficiency and tidal response. Compared with the traditional hydrogeological investigation meth-ods, the new method is much simpler and more efficient and it can provide insight on the spatial and temporal variation patterns of hydraulic properties, and can be used to correct water level measurement, infer aquifer prop-erty, and evaluate the impact of earthquakes on aquifers

    Keywords: water level; barometric response; Earth tide; hydraulic properties

    猜你喜歡
    固體潮大氣壓力潮汐
    固體潮模型在GPS/BDS PPP 中的影響分析
    制作硅膠吸盤
    小學科學(2022年16期)2022-11-19 08:54:20
    潮汐與戰(zhàn)爭(上)
    精密水準測量中固體潮改正的應用
    測繪工程(2022年4期)2022-08-08 13:39:56
    我的一次科學實驗
    絕美海灘
    潮汐式灌溉控制系統(tǒng)的設計及應用
    電子制作(2017年9期)2017-04-17 03:00:56
    干法紙的潮汐
    生活用紙(2016年6期)2017-01-19 07:36:25
    本溪井水位固體潮畸變的研究
    控者無歡
    意林(2014年9期)2014-06-18 14:37:25
    91九色精品人成在线观看| 99riav亚洲国产免费| 女生性感内裤真人,穿戴方法视频| 天天躁狠狠躁夜夜躁狠狠躁| 中文字幕精品免费在线观看视频| 久久人妻av系列| 自拍欧美九色日韩亚洲蝌蚪91| 午夜精品在线福利| 国产成人精品久久二区二区91| 国产精品九九99| 一区福利在线观看| 久久人人精品亚洲av| 91字幕亚洲| 国产成人欧美| 黄色 视频免费看| 亚洲欧美日韩另类电影网站| 久久久久国产精品人妻aⅴ院| 色播在线永久视频| 欧美一级a爱片免费观看看 | 高清在线国产一区| 亚洲成a人片在线一区二区| 久久性视频一级片| 免费人成视频x8x8入口观看| 国产成人精品在线电影| 午夜a级毛片| 啦啦啦免费观看视频1| 国产亚洲精品第一综合不卡| 大型黄色视频在线免费观看| 亚洲欧美激情在线| 97人妻天天添夜夜摸| 一本综合久久免费| 91在线观看av| 美国免费a级毛片| 久久青草综合色| 久久午夜综合久久蜜桃| 亚洲国产中文字幕在线视频| 日本vs欧美在线观看视频| 制服诱惑二区| 这个男人来自地球电影免费观看| 啪啪无遮挡十八禁网站| 69av精品久久久久久| 香蕉丝袜av| 日韩视频一区二区在线观看| 亚洲人成77777在线视频| 久久久久久大精品| 99在线人妻在线中文字幕| 黄片小视频在线播放| 国产欧美日韩一区二区三| 亚洲片人在线观看| 在线观看免费视频网站a站| 久久国产精品男人的天堂亚洲| 国产在线观看jvid| 村上凉子中文字幕在线| av中文乱码字幕在线| 人人妻,人人澡人人爽秒播| 伊人久久大香线蕉亚洲五| 亚洲精品中文字幕在线视频| 亚洲 欧美 日韩 在线 免费| 亚洲成人免费电影在线观看| 亚洲电影在线观看av| 色婷婷久久久亚洲欧美| 少妇熟女aⅴ在线视频| 超碰成人久久| 亚洲一区二区三区不卡视频| 中文字幕高清在线视频| 国产激情欧美一区二区| 久久人人爽av亚洲精品天堂| 国产精品一区二区在线不卡| 婷婷丁香在线五月| 亚洲欧美精品综合久久99| 国产精品亚洲av一区麻豆| 成人国产综合亚洲| 激情在线观看视频在线高清| 又黄又粗又硬又大视频| 搞女人的毛片| 女性被躁到高潮视频| 在线十欧美十亚洲十日本专区| 亚洲专区字幕在线| 亚洲最大成人中文| 国产精品99久久99久久久不卡| 一二三四在线观看免费中文在| 精品福利观看| 国产精品一区二区精品视频观看| 九色亚洲精品在线播放| 免费无遮挡裸体视频| 岛国在线观看网站| 国产精品免费视频内射| 国产亚洲av高清不卡| 黄网站色视频无遮挡免费观看| 精品国产超薄肉色丝袜足j| 夜夜看夜夜爽夜夜摸| 精品不卡国产一区二区三区| 中文亚洲av片在线观看爽| 久久精品91蜜桃| 麻豆久久精品国产亚洲av| 麻豆久久精品国产亚洲av| 成人av一区二区三区在线看| 亚洲全国av大片| 亚洲精品一区av在线观看| 91成人精品电影| 男人舔女人下体高潮全视频| 男女做爰动态图高潮gif福利片 | 国产精品爽爽va在线观看网站 | 亚洲三区欧美一区| 动漫黄色视频在线观看| 日本免费一区二区三区高清不卡 | 久久精品91无色码中文字幕| 国产成人影院久久av| 一本久久中文字幕| 在线观看免费午夜福利视频| 免费不卡黄色视频| 欧洲精品卡2卡3卡4卡5卡区| 国产aⅴ精品一区二区三区波| 国产精品一区二区免费欧美| 黄色毛片三级朝国网站| 精品久久久精品久久久| 亚洲自偷自拍图片 自拍| 亚洲精品一卡2卡三卡4卡5卡| 在线观看www视频免费| 国产精品综合久久久久久久免费 | 欧美日本中文国产一区发布| 啦啦啦 在线观看视频| 日韩视频一区二区在线观看| 日日干狠狠操夜夜爽| 免费看十八禁软件| 成年人黄色毛片网站| 久久这里只有精品19| 精品久久蜜臀av无| 91大片在线观看| 午夜日韩欧美国产| 欧美激情久久久久久爽电影 | or卡值多少钱| 欧美成人性av电影在线观看| 久久中文字幕人妻熟女| 国产三级黄色录像| 婷婷六月久久综合丁香| 欧美黑人精品巨大| 亚洲成av片中文字幕在线观看| 亚洲中文av在线| 国产精品自产拍在线观看55亚洲| 欧美在线一区亚洲| 高清毛片免费观看视频网站| 亚洲色图 男人天堂 中文字幕| av福利片在线| 国产区一区二久久| 欧美另类亚洲清纯唯美| 黄色丝袜av网址大全| 精品国产乱子伦一区二区三区| 亚洲 欧美一区二区三区| 咕卡用的链子| 色婷婷久久久亚洲欧美| 精品不卡国产一区二区三区| 久久精品91蜜桃| 亚洲 国产 在线| 美女午夜性视频免费| 国产高清激情床上av| 免费在线观看完整版高清| 欧美日本中文国产一区发布| or卡值多少钱| 亚洲av熟女| 18禁裸乳无遮挡免费网站照片 | av欧美777| 国产成年人精品一区二区| 成人亚洲精品一区在线观看| 国产精品久久久久久人妻精品电影| 国产精品久久视频播放| 91国产中文字幕| 国产成人欧美在线观看| 欧美中文综合在线视频| 精品国产乱子伦一区二区三区| 久久伊人香网站| 看免费av毛片| 女人爽到高潮嗷嗷叫在线视频| 757午夜福利合集在线观看| 午夜激情av网站| 久久婷婷成人综合色麻豆| 亚洲av成人不卡在线观看播放网| 无遮挡黄片免费观看| 国产主播在线观看一区二区| 亚洲男人天堂网一区| 亚洲人成77777在线视频| 日本 欧美在线| 一区在线观看完整版| 亚洲欧美激情综合另类| 欧洲精品卡2卡3卡4卡5卡区| 成人欧美大片| 757午夜福利合集在线观看| 一边摸一边做爽爽视频免费| 每晚都被弄得嗷嗷叫到高潮| 搡老岳熟女国产| 99国产精品免费福利视频| 制服丝袜大香蕉在线| 免费看十八禁软件| 久久这里只有精品19| 国产成人啪精品午夜网站| 精品久久久久久久久久免费视频| 搡老熟女国产l中国老女人| 国产精品一区二区精品视频观看| 悠悠久久av| 国产xxxxx性猛交| 美女高潮喷水抽搐中文字幕| 91成年电影在线观看| 国产精品久久久久久人妻精品电影| 欧美人与性动交α欧美精品济南到| 在线av久久热| 级片在线观看| 久久久久久久久中文| 97人妻天天添夜夜摸| 日本免费一区二区三区高清不卡 | 日本一区二区免费在线视频| 亚洲精品国产一区二区精华液| 日日夜夜操网爽| 又黄又爽又免费观看的视频| 欧美亚洲日本最大视频资源| 最好的美女福利视频网| 久久婷婷人人爽人人干人人爱 | 亚洲天堂国产精品一区在线| 99国产精品99久久久久| 亚洲午夜精品一区,二区,三区| 亚洲成人久久性| 欧美在线黄色| АⅤ资源中文在线天堂| 亚洲一码二码三码区别大吗| 午夜老司机福利片| √禁漫天堂资源中文www| 亚洲国产欧美网| 中文亚洲av片在线观看爽| 欧美日本亚洲视频在线播放| 欧美黄色片欧美黄色片| 一夜夜www| 国产亚洲av嫩草精品影院| 老司机午夜福利在线观看视频| 我的亚洲天堂| 美女国产高潮福利片在线看| 精品国产亚洲在线| 村上凉子中文字幕在线| www.999成人在线观看| 欧美日韩乱码在线| 精品欧美国产一区二区三| aaaaa片日本免费| 午夜久久久久精精品| 久久精品亚洲精品国产色婷小说| 亚洲五月色婷婷综合| 午夜免费成人在线视频| 亚洲激情在线av| or卡值多少钱| 免费看十八禁软件| 在线av久久热| 黄色a级毛片大全视频| 成人永久免费在线观看视频| 999久久久精品免费观看国产| 亚洲国产精品sss在线观看| 国产精品,欧美在线| 久久中文字幕一级| 18美女黄网站色大片免费观看| 精品一区二区三区四区五区乱码| 久久国产精品人妻蜜桃| av电影中文网址| 999精品在线视频| 无遮挡黄片免费观看| 国产精品一区二区三区四区久久 | 亚洲三区欧美一区| 免费少妇av软件| 成在线人永久免费视频| 久久久久九九精品影院| 两人在一起打扑克的视频| 成人18禁高潮啪啪吃奶动态图| 热re99久久国产66热| 久久国产精品影院| 99re在线观看精品视频| 国产一区二区激情短视频| 免费人成视频x8x8入口观看| 久99久视频精品免费| 久久伊人香网站| 久久久久久久久中文| 国产午夜福利久久久久久| av欧美777| 免费观看人在逋| 国产亚洲av嫩草精品影院| 成人亚洲精品av一区二区| 一级毛片精品| 真人做人爱边吃奶动态| 色婷婷久久久亚洲欧美| 国产精品影院久久| 长腿黑丝高跟| 中国美女看黄片| 一卡2卡三卡四卡精品乱码亚洲| 日日摸夜夜添夜夜添小说| 动漫黄色视频在线观看| 国产主播在线观看一区二区| 男女之事视频高清在线观看| 可以在线观看的亚洲视频| 50天的宝宝边吃奶边哭怎么回事| 9热在线视频观看99| 亚洲国产精品sss在线观看| 国产av精品麻豆| 免费在线观看黄色视频的| 精品久久久久久久久久免费视频| 又紧又爽又黄一区二区| 男女做爰动态图高潮gif福利片 | 欧美日韩福利视频一区二区| 国产精品永久免费网站| 亚洲午夜精品一区,二区,三区| 成人三级黄色视频| 国产精品爽爽va在线观看网站 | 亚洲第一青青草原| 久久久久久久精品吃奶| 久久久久国产一级毛片高清牌| 亚洲中文日韩欧美视频| 999久久久精品免费观看国产| 久久天堂一区二区三区四区| а√天堂www在线а√下载| 欧美黄色淫秽网站| 精品国产国语对白av| 欧美+亚洲+日韩+国产| 亚洲精品美女久久久久99蜜臀| 丝袜人妻中文字幕| 成人三级做爰电影| 亚洲最大成人中文| 成人三级黄色视频| 久久人妻福利社区极品人妻图片| 国产精品美女特级片免费视频播放器 | 国产主播在线观看一区二区| 高清黄色对白视频在线免费看| www国产在线视频色| 久久精品国产综合久久久| 伊人久久大香线蕉亚洲五| 亚洲无线在线观看| 真人一进一出gif抽搐免费| 久久久久久久久免费视频了| 久久香蕉精品热| 久热这里只有精品99| cao死你这个sao货| 国产精品久久电影中文字幕| 老司机在亚洲福利影院| 在线观看一区二区三区| 亚洲国产欧美一区二区综合| 多毛熟女@视频| av免费在线观看网站| avwww免费| 久久精品国产亚洲av高清一级| www国产在线视频色| 美国免费a级毛片| 村上凉子中文字幕在线| 中文字幕另类日韩欧美亚洲嫩草| 欧美丝袜亚洲另类 | 欧美激情 高清一区二区三区| 这个男人来自地球电影免费观看| 日韩欧美在线二视频| 一级a爱视频在线免费观看| 宅男免费午夜| 又大又爽又粗| 久久影院123| 亚洲 欧美 日韩 在线 免费| 在线免费观看的www视频| 青草久久国产| 免费观看精品视频网站| 乱人伦中国视频| 久久人妻福利社区极品人妻图片| 欧美一级毛片孕妇| 女性被躁到高潮视频| 国产又色又爽无遮挡免费看| 男女下面插进去视频免费观看| 中文字幕色久视频| av电影中文网址| 天天躁狠狠躁夜夜躁狠狠躁| 香蕉丝袜av| 国产单亲对白刺激| 成人永久免费在线观看视频| 女人高潮潮喷娇喘18禁视频| 久久人妻福利社区极品人妻图片| 成人亚洲精品av一区二区| 手机成人av网站| 男女做爰动态图高潮gif福利片 | 99香蕉大伊视频| 一级毛片高清免费大全| 少妇粗大呻吟视频| 亚洲中文av在线| 国产成人av教育| 国产亚洲精品综合一区在线观看 | 色播亚洲综合网| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩精品网址| 日韩 欧美 亚洲 中文字幕| 亚洲精品中文字幕一二三四区| 欧美日韩亚洲国产一区二区在线观看| 搞女人的毛片| 欧美在线黄色| av网站免费在线观看视频| www.精华液| 亚洲三区欧美一区| 中文字幕人妻丝袜一区二区| 好男人在线观看高清免费视频 | √禁漫天堂资源中文www| 久热爱精品视频在线9| 91麻豆精品激情在线观看国产| 久久久久九九精品影院| 看免费av毛片| 亚洲专区中文字幕在线| 成人欧美大片| 99re在线观看精品视频| 母亲3免费完整高清在线观看| 亚洲男人的天堂狠狠| 又紧又爽又黄一区二区| 天堂动漫精品| 美女 人体艺术 gogo| 一级毛片高清免费大全| 在线视频色国产色| 久久人妻熟女aⅴ| 多毛熟女@视频| 亚洲欧美激情在线| 欧美日本中文国产一区发布| 在线观看免费午夜福利视频| 99国产精品99久久久久| 无遮挡黄片免费观看| 日韩av在线大香蕉| 正在播放国产对白刺激| 麻豆久久精品国产亚洲av| 美女高潮到喷水免费观看| 夜夜夜夜夜久久久久| 色婷婷久久久亚洲欧美| 久久天堂一区二区三区四区| 好男人在线观看高清免费视频 | 制服人妻中文乱码| 成人av一区二区三区在线看| or卡值多少钱| 亚洲第一欧美日韩一区二区三区| 九色亚洲精品在线播放| 岛国视频午夜一区免费看| 久久久国产成人免费| 中文字幕高清在线视频| av在线播放免费不卡| 日韩高清综合在线| 脱女人内裤的视频| 91成人精品电影| 日韩大码丰满熟妇| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久亚洲av鲁大| 欧美日韩亚洲综合一区二区三区_| 久久草成人影院| 国产精品一区二区精品视频观看| 美女高潮喷水抽搐中文字幕| x7x7x7水蜜桃| 老司机在亚洲福利影院| 欧美黄色片欧美黄色片| 97超级碰碰碰精品色视频在线观看| 成人永久免费在线观看视频| 色综合亚洲欧美另类图片| 精品熟女少妇八av免费久了| 精品久久久久久成人av| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲国产看品久久| 精品欧美国产一区二区三| 韩国av一区二区三区四区| 亚洲片人在线观看| 99在线视频只有这里精品首页| 国产亚洲精品综合一区在线观看 | 制服诱惑二区| 精品国产国语对白av| 亚洲国产毛片av蜜桃av| 一级,二级,三级黄色视频| 中文字幕av电影在线播放| 国产精品久久久av美女十八| 一边摸一边抽搐一进一小说| 麻豆一二三区av精品| 色在线成人网| 日韩大码丰满熟妇| 在线国产一区二区在线| 国产一区在线观看成人免费| 欧美乱妇无乱码| 欧美乱码精品一区二区三区| 国产欧美日韩一区二区三| 欧美日韩一级在线毛片| 亚洲一区中文字幕在线| 国产精品,欧美在线| 亚洲成a人片在线一区二区| 九色亚洲精品在线播放| 国产成人精品在线电影| 多毛熟女@视频| 黑人巨大精品欧美一区二区蜜桃| 久久久久久久精品吃奶| 啦啦啦观看免费观看视频高清 | 亚洲人成电影免费在线| 午夜激情av网站| 91av网站免费观看| 亚洲天堂国产精品一区在线| 神马国产精品三级电影在线观看 | 日韩大尺度精品在线看网址 | 一级作爱视频免费观看| 亚洲专区字幕在线| 日韩三级视频一区二区三区| 人妻久久中文字幕网| 女人被躁到高潮嗷嗷叫费观| 51午夜福利影视在线观看| 正在播放国产对白刺激| 女人高潮潮喷娇喘18禁视频| 午夜免费激情av| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 人妻久久中文字幕网| 午夜福利高清视频| 精品免费久久久久久久清纯| 精品国产乱子伦一区二区三区| 精品国产一区二区三区四区第35| 亚洲成av人片免费观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久久久99蜜臀| 中文字幕色久视频| 久久久久久人人人人人| 男女下面插进去视频免费观看| 久久国产亚洲av麻豆专区| 1024香蕉在线观看| 中出人妻视频一区二区| 久久久久久人人人人人| 母亲3免费完整高清在线观看| 丝袜美足系列| 免费av毛片视频| 欧美丝袜亚洲另类 | 亚洲精品粉嫩美女一区| 国产成人啪精品午夜网站| 国产av在哪里看| 男人的好看免费观看在线视频 | 日本精品一区二区三区蜜桃| 人妻久久中文字幕网| av在线天堂中文字幕| 午夜免费鲁丝| 中文字幕人成人乱码亚洲影| 久久香蕉激情| 国产亚洲av高清不卡| 久久中文字幕人妻熟女| 啦啦啦观看免费观看视频高清 | 国产欧美日韩一区二区精品| 老司机在亚洲福利影院| 久久久精品欧美日韩精品| 国产在线精品亚洲第一网站| 村上凉子中文字幕在线| 91国产中文字幕| 亚洲av五月六月丁香网| 久久人妻福利社区极品人妻图片| 一级毛片精品| 亚洲av成人av| 国产片内射在线| 国产精品香港三级国产av潘金莲| 亚洲成国产人片在线观看| 亚洲伊人色综图| 中文字幕人妻熟女乱码| 国产一区二区三区视频了| 这个男人来自地球电影免费观看| 久久久久国产精品人妻aⅴ院| 免费女性裸体啪啪无遮挡网站| 看黄色毛片网站| 19禁男女啪啪无遮挡网站| 亚洲情色 制服丝袜| 黄片小视频在线播放| 成人18禁在线播放| 在线观看一区二区三区| 精品日产1卡2卡| 久久国产乱子伦精品免费另类| 99久久精品国产亚洲精品| 成人国产一区最新在线观看| 99在线视频只有这里精品首页| 日本黄色视频三级网站网址| 十八禁人妻一区二区| 黄色成人免费大全| cao死你这个sao货| 欧美精品啪啪一区二区三区| 亚洲av第一区精品v没综合| 无限看片的www在线观看| or卡值多少钱| 亚洲中文日韩欧美视频| 久久精品影院6| 欧美最黄视频在线播放免费| 99国产精品免费福利视频| 黑人欧美特级aaaaaa片| 精品人妻在线不人妻| 色尼玛亚洲综合影院| 国产免费av片在线观看野外av| 色综合欧美亚洲国产小说| 亚洲成人精品中文字幕电影| 最新在线观看一区二区三区| av中文乱码字幕在线| 在线观看免费午夜福利视频| 妹子高潮喷水视频| 午夜福利在线观看吧| netflix在线观看网站| 亚洲成人精品中文字幕电影| 精品熟女少妇八av免费久了| 午夜福利成人在线免费观看| 12—13女人毛片做爰片一| 母亲3免费完整高清在线观看| 757午夜福利合集在线观看| 美女高潮喷水抽搐中文字幕| 日本在线视频免费播放| 在线观看免费视频网站a站| 日本一区二区免费在线视频| 不卡一级毛片| a在线观看视频网站| 天天一区二区日本电影三级 | 国产精品 国内视频| 超碰成人久久| 日韩三级视频一区二区三区| 亚洲第一电影网av| 在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 亚洲国产高清在线一区二区三 | 午夜日韩欧美国产| 国产精品,欧美在线| 久久午夜综合久久蜜桃| 国产真人三级小视频在线观看| 国产国语露脸激情在线看| 18禁观看日本| 在线观看午夜福利视频| 99国产精品一区二区三区| 中文字幕另类日韩欧美亚洲嫩草| 99国产精品99久久久久|