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

    臨夏臺(tái)鉆孔系統(tǒng)性質(zhì)的論證

    2016-11-16 00:55:03劉序儼楊錦玲陳超賢關(guān)玉梅陳光趙文波洪明泉
    地球物理學(xué)報(bào) 2016年9期
    關(guān)鍵詞:系統(tǒng)

    劉序儼, 楊錦玲, 陳超賢, 關(guān)玉梅, 陳光, 趙文波, 洪明泉

    福建省地震局, 福州 350003

    ?

    臨夏臺(tái)鉆孔系統(tǒng)性質(zhì)的論證

    劉序儼, 楊錦玲, 陳超賢, 關(guān)玉梅, 陳光, 趙文波, 洪明泉

    福建省地震局, 福州 350003

    對(duì)中國臨夏臺(tái)2013年至2014年兩年的井水位和四分量鉆孔應(yīng)變資料進(jìn)行了預(yù)處理以消除趨勢(shì)與突跳.根據(jù)文獻(xiàn)(Means,1982;Young and Budynas,2005),相互正交的兩條測(cè)線的應(yīng)變觀測(cè)值之和等于面應(yīng)變.文獻(xiàn)(劉序儼等,1988)證明近地表的面應(yīng)變的2/3等于體應(yīng)變,因此,可由4分量鉆孔應(yīng)變觀測(cè)值得到鉆孔體應(yīng)變,然后根據(jù)體應(yīng)變與井水位觀測(cè)資料,從時(shí)域和頻域?qū)υ撱@孔系統(tǒng)的性質(zhì)進(jìn)行了論證.結(jié)果表明,在時(shí)域,體應(yīng)變與井水位高度負(fù)相關(guān).鉆孔系統(tǒng)的靈敏度為—0.1620 mm/10-9.把兩年中的某兩個(gè)月份的兩者時(shí)間坐標(biāo)軸和縱軸比例尺放大,發(fā)現(xiàn)井水位曲線的峰/谷與體應(yīng)變觀測(cè)曲線的谷/峰一一對(duì)應(yīng),兩者的相位滯后非常小.在頻域內(nèi),本文采用Venedikov調(diào)和分析方法分別取得了井水位與體應(yīng)變9個(gè)月的半日波與全日波數(shù)個(gè)波群的逐月潮汐因子與相位滯后,然后作簡(jiǎn)單計(jì)算,得到了鉆孔系統(tǒng)對(duì)上述波群的靈敏度與相位滯后.結(jié)果表明9個(gè)月中大多數(shù)波群的靈敏度不但十分相近,且非常接近由時(shí)域得到的周年頻率分量的靈敏度,但相位滯后誤差較大,本文對(duì)此進(jìn)行了分析,認(rèn)為由反正切得到的相位滯后受計(jì)算誤差影響較大,應(yīng)以時(shí)域經(jīng)審視所得的相位滯后接近于零為準(zhǔn).通過時(shí)域與頻域的分析,表明井水位對(duì)體應(yīng)變的響應(yīng)基本是線性時(shí)不變的,論證了臨夏臺(tái)鉆孔系統(tǒng)基本上滿足了疊加性、齊次性與時(shí)不變性,基本上為一線性時(shí)不變系統(tǒng).

    四分量鉆孔應(yīng)變; 體應(yīng)變; 井水位; 傳遞函數(shù); 線性時(shí)不變系統(tǒng)

    1 引言

    眾所周知,重力儀、傾斜儀與應(yīng)變儀都幾乎不失真地分別把相應(yīng)的固體潮觀測(cè)出來,究其原因,是因?yàn)檫@些觀測(cè)儀器皆為人工精致研制的一種近乎線性時(shí)不變系統(tǒng).試問,在自然界中是否也存在這種系統(tǒng)?回答是肯定的.文獻(xiàn)(劉序儼等,2009;汪成民等,1988)從理論上證明了承壓含水層系統(tǒng)就是這樣一種系統(tǒng).該系統(tǒng)中的井水位變化是該系統(tǒng)對(duì)其體應(yīng)變的一種天然響應(yīng),且響應(yīng)是線性時(shí)不變的.線性時(shí)不變系統(tǒng)最大的特征是該系統(tǒng)的傳遞函數(shù)的振幅譜為一常數(shù),相位譜為零,且該系統(tǒng)的輸出(井水位)對(duì)輸入(體應(yīng)變)有非常大的放大作用.這種理論正確與否,還須作出觀測(cè)的實(shí)證.至今,筆者還未找到這方面的有關(guān)文獻(xiàn).其原因,除了一口鉆孔的井水位變化能清晰地顯示出固體潮形態(tài)以外,還要求能同時(shí)觀測(cè)到該系統(tǒng)的體應(yīng)變.在該系統(tǒng)中,體應(yīng)變?yōu)橐?,作為輸入,井水位變化為果,作為輸?我們就能利用兩者的觀測(cè)值對(duì)該系統(tǒng)是否為線性時(shí)不變系統(tǒng)進(jìn)行實(shí)證,從而對(duì)上述理論作出檢驗(yàn).在進(jìn)行實(shí)證之前,有必要簡(jiǎn)要介紹一下何謂線性時(shí)不變系統(tǒng)?承壓含水層系統(tǒng)為什么是天然線性時(shí)不變系統(tǒng)?該系統(tǒng)的傳遞函數(shù)又是怎樣的?要對(duì)該系統(tǒng)進(jìn)行實(shí)證,最主要是要找到一口能同時(shí)觀測(cè)到體應(yīng)變與井水位變化的鉆孔,且兩者的觀測(cè)曲線都能清晰地記錄到相應(yīng)的固體潮.經(jīng)大量搜索,結(jié)果在全國數(shù)十多個(gè)四分量應(yīng)變觀測(cè)鉆井中,發(fā)現(xiàn)中國甘肅省臨夏地震臺(tái)鉆孔的四分量應(yīng)變與井水位觀測(cè)資料最適合用于對(duì)鉆孔系統(tǒng)性質(zhì)進(jìn)行檢驗(yàn).一口鉆孔能鉆透承壓含水層且井水位與四分量應(yīng)變觀測(cè)資料均能清晰地觀測(cè)到一天雙峰雙谷的固體潮形態(tài),這種鉆孔是可遇而不可求的.有了鉆孔四分量應(yīng)變觀測(cè)資料,如何從這些資料中提取鉆孔的體應(yīng)變信息呢?又怎樣利用臨夏臺(tái)鉆孔的體應(yīng)變與井水位觀測(cè)資料對(duì)這口鉆孔系統(tǒng)的性質(zhì)進(jìn)行論證呢?以上這些問題,正是本文要探討的主要內(nèi)容.

    2 線性時(shí)不變系統(tǒng)的構(gòu)成要素

    如果輸出與輸入成正比,且輸出不失真也不延遲的系統(tǒng)稱為線性時(shí)不變系統(tǒng).在地震系統(tǒng),用以捕捉前兆信號(hào)的所有觀測(cè)系統(tǒng)絕大多數(shù)可歸入這種系統(tǒng).實(shí)際上,任何一種觀測(cè)系統(tǒng)都不可能是一種真正的線性時(shí)不變系統(tǒng).對(duì)于一個(gè)物理可實(shí)現(xiàn)的觀測(cè)系統(tǒng)來說,僅要求在規(guī)定的通頻帶內(nèi)并在允許的觀測(cè)誤差范圍內(nèi)能滿足不失真的要求就可視為一線性時(shí)不變系統(tǒng).

    如果一個(gè)系統(tǒng)同時(shí)滿足疊加性和齊次性,則稱該系統(tǒng)為線性系統(tǒng)(林秩盛,2008;Karu,1995).所謂疊加性,是指當(dāng)若干個(gè)輸入同時(shí)激勵(lì)系統(tǒng)時(shí),該系統(tǒng)的總響應(yīng)等于各個(gè)輸入激勵(lì)單獨(dú)作用時(shí)所產(chǎn)生的響應(yīng)之和,即

    (1)

    式中,x(t)與y(t)分別代表在t時(shí)刻的輸入與輸出,xi(t)代表輸入x(t)中的不同頻率分量,yi(t)為相應(yīng)于xi(t)的輸出.

    所謂齊次性,指的是當(dāng)s為各個(gè)輸出yi(t)對(duì)相應(yīng)各個(gè)輸入xi(t)的靈敏度,則系統(tǒng)的總輸出對(duì)系統(tǒng)的總輸入亦有同樣的靈敏度s.

    如果一個(gè)線性系統(tǒng)滿足

    (2)

    則稱該系統(tǒng)為時(shí)不變系統(tǒng)(林秩盛,2008;Karu,1995),式中τ代表時(shí)間原點(diǎn).時(shí)不變系統(tǒng)表明系統(tǒng)的響應(yīng)與時(shí)間原點(diǎn)無關(guān),系統(tǒng)輸出波形僅與輸入波形有關(guān).

    簡(jiǎn)而言之,一個(gè)線性時(shí)不變觀測(cè)系統(tǒng)的輸出只可能包含輸入中存在的頻率,不含有也不可能有新的頻率出現(xiàn),且輸出的每一頻率分量都放大了同樣的倍數(shù),且都無延遲或延遲相同的時(shí)間,自然它們之和所表示的波形不變.相似性是衡量一個(gè)系統(tǒng)是否為線性時(shí)不變系統(tǒng)的一個(gè)關(guān)鍵且又直觀的指標(biāo).任何一個(gè)觀測(cè)系統(tǒng)都是一種物理可實(shí)現(xiàn)的因果系統(tǒng),無輸入則無輸出,即輸出相對(duì)于輸入不會(huì)也不可能發(fā)生超前現(xiàn)象(劉序儼等,2010).

    3 承壓含水層系統(tǒng)為天然線性時(shí)不變系統(tǒng)的理論依據(jù)

    鉆透承壓含水層的井水稱為承壓水.該水體被隔水層封閉并充滿于含水層各個(gè)部分,由承壓水的埋藏條件,決定它有如下特點(diǎn)(王吉易等,1997):承壓水的分布區(qū)域與補(bǔ)給區(qū)是不一致的;地下水面承受靜水壓力是不自由的;由于受到隔水層的限制,承壓水的水位、流量、水溫及水質(zhì)等受水文氣象因素季節(jié)變化的影響遠(yuǎn)遠(yuǎn)小于潛水;承壓含水層的厚度穩(wěn)定不變,不受降水季節(jié)變化的影響.承壓含水層觀測(cè)系統(tǒng)可看作是這樣一個(gè)由隔水層構(gòu)成的地殼大容器,該容器的頂部有一鉆孔與地下水相通,承壓含水層內(nèi)流體所傳播的靜水壓強(qiáng)P為

    (3)

    式中,靜水壓強(qiáng)P的方向與作用于承壓含水層的外部應(yīng)力方向相反,大小數(shù)值相同.h為井水位,ρ為地下水的密度,g為重力加速度,ρg為地下水的重度.由于承壓含水層的靜水壓強(qiáng)的變化會(huì)導(dǎo)致含水層的體積變化,從而引起了鉆孔中井水位的變化.因此,井水位可視作為一種水壓計(jì).

    根據(jù)彈性力學(xué)與流體靜力學(xué)原理,承壓含水層這種彈性空隙介質(zhì)在靜水壓力下井水位對(duì)體應(yīng)變響應(yīng)的動(dòng)力學(xué)方程(劉序儼等,2009)為

    (4)

    式中,θ為體應(yīng)變,k為含水層彈性空隙介質(zhì)的體積模量之倒數(shù),即含水層彈性空隙介質(zhì)的壓縮系數(shù).其表達(dá)式為

    (5)

    式中n為含水層的孔隙度,ks與kw分別為含水層的巖體及所含的流體的體積模量的倒數(shù),即壓縮系數(shù).式中“-”表示體應(yīng)變壓縮則井水位上升,反之亦然.

    (6)

    在本文,h以mm為量綱,θ以10-9為應(yīng)變量綱,則s的量綱為mm/10-9.在地殼中巖石的體積模量Bs為(0.44~1)×1011Pa,水的Bw=2.2×109Pa,Bs/Bw為20~45.4,ρg=10 Pa·mm-1.在井水位以mm為量綱的前提下,不管是巖體還是水體的體積模量B的數(shù)值都要比水的重度ρg數(shù)值大108~1010倍,其結(jié)果是承壓含水層的靈敏度s最高達(dá)—10 mm/10-9,即承壓含水層發(fā)生10-9體應(yīng)變,井水位將產(chǎn)生10 mm的變化,負(fù)號(hào)表示兩者的變化方向是相反的.這就是為什么體應(yīng)變激勵(lì)承壓含水層系統(tǒng)時(shí)井水位具有天然放大機(jī)制的原因,因此井水位可作為一種體應(yīng)變計(jì).

    設(shè)η為井水位觀測(cè)的格值,由定義有

    (7)

    η的量綱為10-9/mm.

    根據(jù)控制論理論,若輸出與輸入成正比,且輸出不失真也不延遲反映輸入的環(huán)節(jié),稱為比例(或放大)環(huán)節(jié)(曾勁,2007;陸一心,2006).(4)式即為該環(huán)節(jié)的動(dòng)力學(xué)方程,式中,1/ρgk為一常數(shù),即為井水位觀測(cè)系統(tǒng)的靈敏度.因此承壓含水層觀測(cè)系統(tǒng)作為一個(gè)零階的放大環(huán)節(jié),其傳遞函數(shù)為

    (8)

    式中,ω為角速度,量綱為°/h,傳遞函數(shù)為常數(shù),表示該系統(tǒng)的振幅譜為常數(shù);相位譜為零,表明承壓含水層系統(tǒng)對(duì)輸入中任一頻率分量在輸出中均被放大了同一倍數(shù),且無延遲.又因?qū)嶋H物理系統(tǒng)總具有慣性,使輸出不會(huì)超前于輸入,因此承壓含水層系統(tǒng)可看作一個(gè)線性時(shí)不變的可實(shí)現(xiàn)的物理因果系統(tǒng).

    4 臨夏臺(tái)四分量應(yīng)變觀測(cè)鉆孔簡(jiǎn)介

    臨夏臺(tái)位于中國甘肅省,緯度為北緯35.6°,經(jīng)度為東經(jīng)103.2°.臨夏臺(tái)是花崗巖地層,井深45 m,探頭裝在44.7 m深處.鉆孔內(nèi)水位距井口10 m,水位計(jì)的分辨率為0.03 mm.應(yīng)變觀測(cè)探頭為圓筒形,其直徑107 mm,高450 mm.安裝探頭處鉆孔直徑130 mm,鉆頭上布設(shè)了相互夾角為45°的四條應(yīng)變測(cè)線.1#測(cè)線的方位角為92°或272°,四分量鉆孔應(yīng)變儀的分辨率為10-10.四分量鉆孔應(yīng)變測(cè)量元件分布如圖1所示.臨夏臺(tái)應(yīng)變觀測(cè)鉆井,除了能提供四分量應(yīng)變觀測(cè)外,其鉆孔的井水位觀測(cè)如同應(yīng)變分量一樣能顯示固體潮變化形態(tài).臨夏臺(tái)采用YRY-4型四分量鉆孔應(yīng)變儀.YRY-4型應(yīng)變儀的觀測(cè)頻帶非常寬,除了能記錄到固體潮,還能記錄到地震波.應(yīng)變觀測(cè)量綱為10-9.除應(yīng)變觀測(cè)外,還能進(jìn)行水位與水溫觀測(cè).水位觀測(cè)量綱為mm.本文收集了該鉆孔2013年1月1日0時(shí)至2014年12月31日23時(shí)的四分量應(yīng)變與井水位觀測(cè)資料.應(yīng)變與井水位觀測(cè)如同地傾斜、重力觀測(cè)一樣為相對(duì)觀測(cè).為此,把這兩年的四分量應(yīng)變與井水位整點(diǎn)觀測(cè)值分別減去2013年1月1日零時(shí)的相應(yīng)觀測(cè)值.圖2分別為四分量應(yīng)變與井水位相對(duì)觀測(cè)值曲線圖.

    從圖2可發(fā)現(xiàn),應(yīng)變分量與井水位相對(duì)觀測(cè)數(shù)據(jù)間或出現(xiàn)突跳,且應(yīng)變數(shù)據(jù)還存在趨勢(shì)變化.為此,本文對(duì)二者進(jìn)行了預(yù)處理,排除了趨勢(shì)與突跳.

    根據(jù)彈性力學(xué)理論(Means,1982;Young and Budynas,2005)通過同一點(diǎn)的相互正交的兩條應(yīng)變測(cè)線的應(yīng)變觀測(cè)值之和為鉆孔的面應(yīng)變,根據(jù)文獻(xiàn)(劉序儼等,1988),近地表的鉆孔面應(yīng)變的2/3為體應(yīng)變,從而可由四分量鉆孔應(yīng)變值得到鉆孔的體應(yīng)變.圖3為臨夏臺(tái)由相互正交方向的兩組應(yīng)變觀測(cè)值所取得的兩條體應(yīng)變觀測(cè)曲線圖.

    從圖3可發(fā)現(xiàn),鉆孔兩條體應(yīng)變觀測(cè)值是十分吻合的.本文取兩組體應(yīng)變觀測(cè)值的平均值作為最終鉆孔體應(yīng)變值.圖4為排除線性趨勢(shì)與突跳后的鉆孔體應(yīng)變與井水位觀測(cè)曲線圖.

    圖1 臨夏臺(tái)四分量鉆孔應(yīng)變測(cè)線分布圖Fig.1 Distribution of four component of borehole strain measurement lines of Linxia Station

    圖2 臨夏臺(tái)2013年1月到2014年12月鉆孔四分量應(yīng)變與井水位相對(duì)觀測(cè)曲線(上圖為應(yīng)變,量綱為10-9,下圖為井水位,量綱為mm)Fig.2 The observation curves of four-component borehole strain and well water level at Linxia Station (the unit of strain is 10-9in the figure above, the unit of well water level is mm in the figure below)

    圖3 2013年1月到2014年12月鉆孔兩條體應(yīng)變觀測(cè)曲線(量綱為10-9)Fig.3 Two graphs of volume strain of borehole(the unit is 10-9)

    圖4 2013年1月到2014年12月鉆孔體應(yīng)變與井水位觀測(cè)曲線圖上圖應(yīng)變量綱為10-9,下圖井水位量綱為mm.Fig.4 The graphs of volume strain and well water level of boreholeThe unit of strain is 10-9 in the figure above, The unit of well water level is mm in the figure below.

    5 觀測(cè)實(shí)證

    從圖4可發(fā)現(xiàn),井水位曲線與體應(yīng)變曲線都具有周年變化特征,且互為鏡像對(duì)稱,形狀十分相似.其相關(guān)系數(shù)R=—0.9353.相關(guān)系數(shù)與靈敏度皆為負(fù)值,兩者曲線又鏡像對(duì)稱,其原因是體應(yīng)變與井水位互為因果關(guān)系.對(duì)兩者進(jìn)行回歸分析.其回歸系數(shù)即為鉆孔系統(tǒng)的靈敏度,其值s=—0.1620 mm/10-9,其倒數(shù)即為格值η=—6.17×10-9/mm.由圖4無法判別井水位相對(duì)于體應(yīng)變是否存在相位滯后,要做到這點(diǎn),可由圖中疊加在兩條觀測(cè)曲線上半日波與全日波的峰/谷和谷/峰是否在同一時(shí)間上來確定.從圖5—6可發(fā)現(xiàn),井水位無超前現(xiàn)象,且兩者的相位滯后非常小,這就從時(shí)間域論證了鉆孔系統(tǒng)基本上為線性時(shí)不變系統(tǒng).

    按照線性時(shí)不變系統(tǒng)理論,疊加在體應(yīng)變與井水位觀測(cè)曲線上的半日波與全日波諸波群的靈敏度亦應(yīng)與從時(shí)域得到的周年分量的靈敏度相同.為了驗(yàn)證這一點(diǎn),此時(shí)要從時(shí)域轉(zhuǎn)向頻域進(jìn)行論證.

    臨夏臺(tái)鉆孔系統(tǒng)的具體情況,我們是無法了解的,僅能視為一個(gè)黑箱,體應(yīng)變與井水位可分別作為該黑箱的輸入與輸出,如圖7所示.

    根據(jù)控制論與信號(hào)處理理論(林秩盛,2008;Karu,1995),輸入與輸出之間存在如下關(guān)系:

    (9)

    式中,*表示褶積,h(t)稱為單位脈沖時(shí)間響應(yīng)函數(shù),亦稱為系統(tǒng)的權(quán)系數(shù).

    圖5 2013年5月份體應(yīng)變與井水位曲線圖Fig.5 The observation graphs of volume strain of borehole and well water level at May,2013

    圖6 2014年7月份體應(yīng)變與井水位曲線圖Fig.6 The observation graphs of volume strain of borehole and well water level at July,2014

    圖7 鉆孔系統(tǒng)的輸入與輸出示意圖Fig.7 The schematic diagram of input and output of the borehole system

    對(duì)(9)式,兩邊作傅里葉變換,則有

    (10)

    對(duì)體應(yīng)變與井水位相對(duì)觀測(cè)資料分別作傅里葉變換.圖8與圖9分別為體應(yīng)變和井水位的傅里葉譜.因周期T=360°/ω,圖中把角速度ω轉(zhuǎn)換成周期T(量綱為h).譜的絕對(duì)值即為振幅譜A(T).圖8—9顯示出在T=12 h,24 h處,即在半日波與全日波頻段,兩者振幅譜A(T)值十分顯著,揭示出鉆孔系統(tǒng)在這兩個(gè)頻段處有較大的靈敏度,也就是有較小的格值.而在其他頻段,振幅譜值則相對(duì)微小.對(duì)比圖8和圖9,可以發(fā)現(xiàn)井水位中保留了體應(yīng)變中的半日波與全日波分量,因?yàn)閮烧咴谶@兩個(gè)頻段上振幅譜都十分顯著.由于受到潛水、降水等非線性干擾,兩者在非潮汐頻段上則存在相對(duì)較小的連續(xù)譜.但前者比后者受潛水、降水等非體積干擾要小的多,其振幅譜要比潮汐譜小很多.如果沒有受到上述非線性干擾,由波群組合而成的復(fù)雜周期函數(shù)的傅里葉譜則應(yīng)為線譜.

    設(shè)體應(yīng)變與井水位的傅里葉譜分別為x(ω)=Ai(ω)e-ji(ω),y(ω)=Ao(ω)e-jo(ω),式中Ai(ω)與Ao(ω)分別為兩者的振幅譜,i(ω)與ο(ω)分別為其相位譜.由(10)式,則可得鉆孔系統(tǒng)的傳遞函數(shù)H(ω)及其振幅譜A(ω)與相位譜φ(ω)分別為

    (11)

    圖8 體應(yīng)變的傅里葉譜|x(T)|(量綱為10-9)Fig.8 Fourier spectrum of volume strain of borehole (the unit is 10-9)

    圖9 井水位的傅里葉譜|y(T)|(量綱為mm)Fig.9 Fourier spectrum of well water level(the unit is mm)

    限于篇幅,本文僅選取了2013年5—7月份與2014年5—7月份的體應(yīng)變與井水位觀測(cè)資料逐月進(jìn)行調(diào)和分析.在半日波調(diào)和分析結(jié)果中僅選取了振幅較大的M2波與S2波,在全日波中選取了O1波與P1波,將井水位的某波群的潮汐因子除以體應(yīng)變相應(yīng)波群的潮汐因子,則可得鉆孔系統(tǒng)在該波群頻段上的振幅譜A(ω),其值為Ao(ω)/Ai(ω),井水位與體應(yīng)變的相位譜之差即為鉆孔系統(tǒng)在該頻段上的相位滯后o(ω)-i(ω).具體數(shù)值見表1至表2.

    表1 體應(yīng)變與井水位調(diào)和分析結(jié)果及鉆孔系統(tǒng)傳遞函數(shù)一覽表(2013-05—2013-07)Table 1 The harmonic analysis results of volume strain and well water level and transfer function of borehole system at Linxia Station for M2、S2、O1 and P1 wave groups in 2013-05—2013-07

    表2 體應(yīng)變與井水位調(diào)和分析結(jié)果及鉆孔系統(tǒng)傳遞函數(shù)一覽表(2014-05—2014-07)Table 2 The harmonic analysis results of volume strain and well water level and transfer function of borehole system at Linxia Station for M2、S2、O1 and P1 wave groups in 2014-05—2014-07

    從表2可以發(fā)現(xiàn),其4個(gè)波群共12個(gè)靈敏度沒有像表1中的12個(gè)靈敏度那樣接近于周年頻率的靈敏度.從圖4可見,表2時(shí)段對(duì)應(yīng)的觀測(cè)曲線較表1時(shí)段對(duì)應(yīng)的觀測(cè)曲線更為粗糙,究其原因可能是受到的干擾比較多.為了證實(shí)這種情況,本文選擇了圖4中觀測(cè)曲線比較光滑的2013年2-4月的兩者的觀測(cè)資料進(jìn)行了鉆孔系統(tǒng)對(duì)4個(gè)波群的靈敏度與相位滯后的計(jì)算,結(jié)果見表3.從表3可看出,該時(shí)段的4個(gè)波群的靈敏度比表1的結(jié)果更接近周年頻率的靈敏度.因此,非體積參數(shù)干擾越小越少,其靈敏度就越接近.表1—3表明鉆孔系統(tǒng)在潮汐波頻段,其靈敏度是大致與周年分量的靈敏度相接近,從高頻到低頻鉆孔系統(tǒng)的振幅譜大致是平坦的.

    表3 體應(yīng)變與井水位調(diào)和分析結(jié)果及鉆孔系統(tǒng)傳遞函數(shù)一覽表(2013-02—2013-04)Table 3 The harmonic analysis results of volume strain and well water level and transfer function of borehole system at Linxia Station for M2、S2、O1 and P1 wave groups in 2013-02—2013-04

    6 認(rèn)識(shí)與討論

    文獻(xiàn)(劉序儼等,2009)推論出在含水層的傳遞函數(shù)的振幅譜為常數(shù),相位譜為零.這是根據(jù)流體靜力學(xué)和彈性力學(xué)的因果定律得出的,實(shí)際情況并非如此.在文獻(xiàn)(劉序儼等,2009)中,對(duì)承壓含水層作了理想化的假設(shè),并采用了刪繁就簡(jiǎn)的方法把引起井水位產(chǎn)生非體應(yīng)變的參數(shù)(如潛水、降水、氣壓等)刪除掉.只把含水層的巖石固體框架與流體的體積模量及其孔隙度保留住,這顯然是與實(shí)際情況不符合的.正如哲學(xué)家維特根斯坦所說“神秘的不是世界是怎樣的,而是世界是這樣的”(張法,2013).以中國臨夏臺(tái)體應(yīng)變與井水位觀測(cè)資料從潮汐波群高頻到氣溫周年變化低頻所作的論證,表明了鉆孔系統(tǒng)的傳遞函數(shù)的振幅譜是大致接近的,其相位譜大致接近于零,表明臨夏臺(tái)鉆孔系統(tǒng)是近似于一個(gè)線性時(shí)不變系統(tǒng).這種近似性緣于臨夏臺(tái)鉆孔系統(tǒng)并不完全符合文獻(xiàn)(劉序儼等,2009)對(duì)承壓含水層系統(tǒng)所作的理想假設(shè).該文獻(xiàn)所給出傳遞函數(shù)是一種抽象的數(shù)理演繹的結(jié)果,是抽取了承壓含水層的本質(zhì)屬性,撇開非本質(zhì)屬性所得到的.如果外界非體積參數(shù)(潛水、降水等)干擾越小,則鉆孔系統(tǒng)的傳遞函數(shù)的振幅譜就越接近于一個(gè)常數(shù),其相位譜就越接近于零.而本文(8)式所示的承壓含水層的理論傳遞函數(shù)則是一個(gè)可趨近但又不能到達(dá)的彼岸.如果對(duì)臨夏臺(tái)鉆孔系統(tǒng)撇開了外界那些干擾,其實(shí)證結(jié)果可能更接近線性時(shí)不變系統(tǒng).但目前,本文沒有更多的觀測(cè)資料作進(jìn)一步實(shí)證.從目前所作的實(shí)證,也說明文獻(xiàn)(劉序儼等,2009)給出的公式是正確的.理想的承壓含水層系統(tǒng)之所以為線性系統(tǒng),是由文獻(xiàn)(劉序儼等,2009)所假定的那樣是因?yàn)榈貧な蔷€性、均質(zhì)和各向同性的彈性體.

    “物含妙理總堪尋”,人類的好奇心促使人們對(duì)自然界某些現(xiàn)象作出理論解釋.人類的懷疑精神又促使人們對(duì)這些理論解釋產(chǎn)生質(zhì)疑,要以觀測(cè)事實(shí)對(duì)此解釋作出實(shí)證,如此互動(dòng)促進(jìn)了科學(xué)的進(jìn)步與發(fā)展,對(duì)理論進(jìn)行證明與證偽,至今已成為科學(xué)共同體一種共識(shí).如今井水位、深井應(yīng)變觀測(cè)越來越受到重視,測(cè)站也越來越多,本文對(duì)井水位、應(yīng)變觀測(cè)及其研究進(jìn)行了搜索,找到了較多有關(guān)文獻(xiàn),但這些文獻(xiàn)皆未對(duì)鉆孔系統(tǒng)的性質(zhì)進(jìn)行討論,最主要的原因是由于無法取得同一鉆井的井水位與體應(yīng)變觀測(cè)資料.為此,本文未將這些文獻(xiàn)一一列出.有幸本文取得了中國臨夏臺(tái)鉆孔的井水位與四分量鉆孔應(yīng)變資料,對(duì)臨夏臺(tái)鉆孔系統(tǒng)性質(zhì)所作的實(shí)證,不僅論證了臨夏臺(tái)鉆孔系統(tǒng)接近于線性時(shí)不變系統(tǒng),同時(shí)亦是對(duì)YRY-4型鉆孔應(yīng)變儀四分量應(yīng)變與井水位觀測(cè)資料的可靠性的一種佐證.

    Fang J. 1984. Solid Earth Tides (in Chinese). Beijing: Science Press.

    Jiang J, Li S L, Zhang Y B, et al. 2000.Seismic Precursor Information Processing System (in Chinese).Beijing:Seismological Press. Karu Z Z. 1995. Signals and Systems Made Ridiculously Simple. Cambridge, MA: ZiZi Press.

    Lin Z S. 2008. The Signal and Linear System (in Chinese). Beijing: Tsinghua University Press.

    Liu X Y, Li P, Zhang Y B. 1988. Calculation of theoretical value of cubic expansion tide and its harmonic analysis.CrustalDeformationandEarthquake(in Chinese), 8(4): 354-358.

    Liu X Y, Zheng X J, Wang L, et al. 2009. Response analysis of the well-water-level system in confined aquifer.ChineseJ.Geophys. (in Chinese), 52(12): 3147-3157, doi: 10.3969/j.issn.0001-5733.2009.12.025.

    Liu X Y, Huang S M, Liang Q Q, et al. 2010.Fidelity analysis of an observation system.JournalofGeodesyandGeodynamics(in Chinese), 30(Supp. 1): 1-6.

    Lu Y X. 2006. Modern Engineering Control Theory (in Chinese). Beijing: Chemical Industry Press.

    Means W D. 1982.Stress and Strain (in Chinese). Ding Z Y, Yin Y Q Trans. Beijing: Science Press, 74-114.

    Melchior P. 1984. The tides of the Planet Earth (in Chinese). Du P R, Wu Q P, Chen Y H Trans. Beijing: Science Press.

    The Earth Physics Department of Peking University, the Geodetic Survey Department of Wuhan Survey College, the Geophysics Staff Room of University of Science and Technology of China. 1982. Gravity and Earth Tide (in Chinese). Beijing: Seismological Press. Wang C M, Che Y T, Wan D K, et al. 1988. Study on the Micro Dynamics of Groundwater (in Chinese). Beijing: Seismological Press.

    Wang J Y, Dong S Y, Chen J M. 1997. Earthquake Prediction Method of Underground Fluid (in Chinese). Beijing: Seismological Press. Young W C, Budynas R G. 2005.Roark′s Formulas for Strain and Stress (in Chinese).Yue Z F, Gao X S, Wang F H, et al., Trans. Beijing: Science Press, 9-26.

    Zeng J. 2007. Foundation of Control Engineering (in Chinese). Beijing: Publishing House of Electronics Industry.

    Zhang F. 2013. Features of western philosophy in the perspective of

    four philosophical sayings.JournalofCUPL(in Chinese), 22(4): 142-148.

    附中文參考文獻(xiàn)

    北京大學(xué)地球物理系, 武漢測(cè)繪學(xué)院大地測(cè)量系, 中國科學(xué)技術(shù)大學(xué)地球物理教研室. 1982. 重力與固體潮教程. 北京: 地震出版社.

    方俊. 1984. 固體潮. 北京: 科學(xué)出版社.

    蔣駿, 李勝樂, 張雁濱等. 2000. 地震前兆信息處理與軟件系統(tǒng). 北京: 地震出版社.

    林秩盛. 2008. 信號(hào)與線性系統(tǒng). 北京: 清華大學(xué)出版社.

    劉序儼, 李平, 張雁濱. 1988. 地表的面應(yīng)變和體應(yīng)變固體潮理論值計(jì)算及其調(diào)和分析. 地殼形變與地震, 8(4): 354-358.

    劉序儼, 鄭小菁, 王林等. 2009. 承壓井水位觀測(cè)系統(tǒng)對(duì)體應(yīng)變的響應(yīng)機(jī)制分析. 地球物理學(xué)報(bào), 52(12): 3147-3157, doi: 10.3969/j.issn.0001-5733.2009.12.025.

    劉序儼, 黃聲明, 梁全強(qiáng)等. 2010. 觀測(cè)系統(tǒng)逼真度特性分析. 大地測(cè)量與地球動(dòng)力學(xué),30(S1): 1-6.

    陸一心. 2006. 現(xiàn)代工程控制理論. 北京: 化學(xué)工業(yè)出版社.

    Means W D. 1982. 應(yīng)力和應(yīng)變.丁中一, 殷有泉譯. 王仁校.北京: 科學(xué)出版社, 74-114.

    Melchior P. 1984. 行星地球的固體潮.杜品仁, 吳慶鵬, 陳益恵譯. 北京: 科學(xué)出版社.

    汪成民, 車用太, 萬迪堃等. 1988. 地下水微動(dòng)態(tài)研究. 北京: 地震出版社.

    王吉易, 董守玉, 陳建民. 1997. 地下流體地震預(yù)報(bào)方法.北京: 地震出版社.

    Young W C, Budynas R G.2005. 羅氏應(yīng)力應(yīng)變公式手冊(cè).岳珠峰, 高行山, 王峰會(huì)等譯.北京: 科學(xué)出版社, 9-26.

    曾勁. 2007. 控制工程基礎(chǔ). 北京: 電子工業(yè)出版社.

    張法. 2013. 從四句哲學(xué)名言看西方哲學(xué)的特質(zhì). 中國政法大學(xué)學(xué)報(bào), 22(4): 142-148.

    (本文編輯 胡素芳)

    The argumentation of properties of borehole system at Linxia station, China

    LIU Xu-Yan , YANG Jin-Ling, CHEN Chao-Xian, GUAN Yu-Mei, CHEN Guang, ZHAO Wen-Bo, HONG Ming-Quan

    EarthquakeAdministrationofFujianProvince,Fuzhou350003,China

    2 years of four-component borehole strain data and well water level data are preprocessed to remove trends and outliers at Linxia Station from 2013 to 2014. According to the literatures (Means, 1982; Young and Budinas, 2005), the sum of observation values of two strain measurement lines orthogonal to each other is equivalent to surface strain is equal to volume strain. Literature (Liu et al., 1988) proved that two thirds of near-ground surface strain is equal to volume strain. Therefore, the volume strain can be obtained from four-component borehole strain observation data. And the properties of the borehole system can be demonstrated from time and frequencies domain based on volume and underground well water level observation data. Results show that negative correlation between volume strain and well water level is high in time domain. The sensitivity of borehole system is —0.1620 mm/10-9.Enlarged the proportion of two coordinates of volume strain and well water level observation curves of some two months in the two year period, we found that the peak and valley of well water level corresponds to the valley and peak of volume strain one by one, the phase lag is very small. In the frequency domain, monthly tidal factor and phase lag of diurnal and semidiurnal wave groups of volume strain and water level have been obtained with Venedikov tidal harmonic analysis method. Then the sensitivity and phase lag of the borehole strain system is calculated simply. Results show that the sensitivity of most wave groups are not only close to each other, but also close to the sensitivity of annually periodic change obtained from the time domain, while the errors of phase lag are large. Considering the phase lag obtained by arctangent is influenced by calculation error, the phase lag closed to zero which is examined from time domain should be more reliable. The analysis from time and frequency domain show that the response of well water level to volume strain is substantially linear time-invariant. The borehole system of Linxia station basically meets the superposition, homogeneity and time invariance, therefore, borehole system at Linxia station is linear time-invariant system basically.

    Four-component borehole strain; Volume strain; Well water level; Transfer function; Linear time-invariant system

    10.6038/cjg20160918.Liu X Y, Yang J L, Chen C X, et al. 2016. The argumentation of properties of borehole system at Linxia station, China.ChineseJ.Geophys. (in Chinese),59(9):3343-3353,doi:10.6038/cjg20160918.

    福建省地震局科研基金(SF201513)資助.

    劉序儼,男,1939年生,研究員,1964年畢業(yè)于武漢測(cè)繪學(xué)院(現(xiàn)武漢大學(xué))天文大地測(cè)量專業(yè),長(zhǎng)期從事固體潮與地殼形變觀測(cè)研究. E-mail:xuyanliu@126.com

    10.6038/cjg20160918

    P631

    2016-04-15,2016-09-29收修定稿

    劉序儼,楊錦玲,陳超賢等. 2016. 臨夏臺(tái)鉆孔系統(tǒng)性質(zhì)的論證. 地球物理學(xué)報(bào),59(9):3343-3353,

    猜你喜歡
    系統(tǒng)
    Smartflower POP 一體式光伏系統(tǒng)
    WJ-700無人機(jī)系統(tǒng)
    ZC系列無人機(jī)遙感系統(tǒng)
    基于PowerPC+FPGA顯示系統(tǒng)
    基于UG的發(fā)射箱自動(dòng)化虛擬裝配系統(tǒng)開發(fā)
    半沸制皂系統(tǒng)(下)
    FAO系統(tǒng)特有功能分析及互聯(lián)互通探討
    連通與提升系統(tǒng)的最后一塊拼圖 Audiolab 傲立 M-DAC mini
    一德系統(tǒng) 德行天下
    PLC在多段調(diào)速系統(tǒng)中的應(yīng)用
    成人二区视频| 黄片无遮挡物在线观看| 亚洲综合色网址| 久久午夜福利片| 久久国产精品男人的天堂亚洲 | 波多野结衣一区麻豆| 精品午夜福利在线看| av在线播放精品| 亚洲五月色婷婷综合| 精品一区二区三卡| 国产av码专区亚洲av| 午夜av观看不卡| 中文字幕最新亚洲高清| 国产亚洲一区二区精品| 高清视频免费观看一区二区| 国产成人精品无人区| 国产成人91sexporn| 高清不卡的av网站| 建设人人有责人人尽责人人享有的| 精品国产一区二区三区四区第35| 免费女性裸体啪啪无遮挡网站| 日韩不卡一区二区三区视频在线| 毛片一级片免费看久久久久| 97在线人人人人妻| 中文字幕另类日韩欧美亚洲嫩草| 国产精品人妻久久久影院| 国产1区2区3区精品| 婷婷色av中文字幕| 女性生殖器流出的白浆| 亚洲国产av影院在线观看| 亚洲av电影在线观看一区二区三区| 国产精品一区www在线观看| 久久国产精品男人的天堂亚洲 | 99热6这里只有精品| 人妻 亚洲 视频| 高清黄色对白视频在线免费看| 七月丁香在线播放| 亚洲四区av| 黑人巨大精品欧美一区二区蜜桃 | 大香蕉97超碰在线| 涩涩av久久男人的天堂| 免费观看在线日韩| 亚洲一级一片aⅴ在线观看| 99久久综合免费| 久久久久久久国产电影| 女人久久www免费人成看片| 哪个播放器可以免费观看大片| 精品第一国产精品| 国产亚洲一区二区精品| 女人久久www免费人成看片| 免费在线观看黄色视频的| 久久97久久精品| 免费看av在线观看网站| 女性生殖器流出的白浆| 少妇的丰满在线观看| 最后的刺客免费高清国语| 亚洲国产精品国产精品| 欧美日韩综合久久久久久| 美女国产视频在线观看| 春色校园在线视频观看| 国产一区有黄有色的免费视频| 人体艺术视频欧美日本| 美女国产视频在线观看| 亚洲欧美中文字幕日韩二区| 秋霞伦理黄片| 男人操女人黄网站| 性色avwww在线观看| 涩涩av久久男人的天堂| 久久久久久久久久久久大奶| 中文字幕最新亚洲高清| 妹子高潮喷水视频| 国产高清不卡午夜福利| 亚洲国产欧美在线一区| 在线免费观看不下载黄p国产| 考比视频在线观看| 亚洲精品视频女| 激情五月婷婷亚洲| 久久精品久久久久久久性| 亚洲在久久综合| 亚洲色图 男人天堂 中文字幕 | 久热这里只有精品99| 男的添女的下面高潮视频| 午夜久久久在线观看| 欧美日韩av久久| 午夜免费男女啪啪视频观看| 如日韩欧美国产精品一区二区三区| 黄片播放在线免费| 亚洲少妇的诱惑av| 天天操日日干夜夜撸| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品电影小说| 日韩电影二区| 91国产中文字幕| 黄色怎么调成土黄色| 欧美成人精品欧美一级黄| 黄色一级大片看看| 一本—道久久a久久精品蜜桃钙片| 日本欧美国产在线视频| 两个人免费观看高清视频| 亚洲av国产av综合av卡| 日日爽夜夜爽网站| 天天影视国产精品| 黄片无遮挡物在线观看| 性色avwww在线观看| 国产色婷婷99| 伊人亚洲综合成人网| 制服人妻中文乱码| 亚洲国产成人一精品久久久| 90打野战视频偷拍视频| 2018国产大陆天天弄谢| 欧美 日韩 精品 国产| 亚洲中文av在线| av免费在线看不卡| 少妇人妻精品综合一区二区| 岛国毛片在线播放| 777米奇影视久久| 黄色视频在线播放观看不卡| 国产探花极品一区二区| 亚洲精品色激情综合| a级毛片黄视频| 两个人免费观看高清视频| 在线观看免费日韩欧美大片| 亚洲欧美日韩卡通动漫| 十八禁高潮呻吟视频| 狠狠婷婷综合久久久久久88av| 在线观看三级黄色| 一本大道久久a久久精品| 亚洲精品第二区| 欧美人与性动交α欧美软件 | 男男h啪啪无遮挡| 色视频在线一区二区三区| 欧美另类一区| 国产视频首页在线观看| 国产免费一区二区三区四区乱码| 校园人妻丝袜中文字幕| 国产精品女同一区二区软件| 欧美日韩视频高清一区二区三区二| 亚洲精品久久午夜乱码| 又黄又粗又硬又大视频| 午夜精品国产一区二区电影| 国产免费又黄又爽又色| videosex国产| 亚洲天堂av无毛| 国产男人的电影天堂91| 男男h啪啪无遮挡| 老熟女久久久| www.av在线官网国产| 日本与韩国留学比较| 午夜久久久在线观看| 久久人人爽人人片av| 少妇被粗大的猛进出69影院 | 综合色丁香网| 精品一区二区三区视频在线| 午夜福利视频在线观看免费| 国产白丝娇喘喷水9色精品| 又大又黄又爽视频免费| 久久精品久久精品一区二区三区| 亚洲人与动物交配视频| 午夜精品国产一区二区电影| 秋霞伦理黄片| 夫妻性生交免费视频一级片| 波多野结衣一区麻豆| 人人妻人人添人人爽欧美一区卜| 欧美日韩国产mv在线观看视频| 中文字幕制服av| 国产xxxxx性猛交| 大陆偷拍与自拍| av片东京热男人的天堂| 一区二区三区乱码不卡18| 夜夜爽夜夜爽视频| 成人毛片60女人毛片免费| 夫妻午夜视频| 五月天丁香电影| 精品亚洲成a人片在线观看| 亚洲精品一二三| 午夜福利网站1000一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 午夜91福利影院| 欧美国产精品va在线观看不卡| 欧美成人午夜免费资源| 日韩熟女老妇一区二区性免费视频| 日韩精品有码人妻一区| 91在线精品国自产拍蜜月| 亚洲国产成人一精品久久久| 黑丝袜美女国产一区| 亚洲精品,欧美精品| 少妇人妻久久综合中文| www.熟女人妻精品国产 | 婷婷色麻豆天堂久久| 国产精品女同一区二区软件| 久久人人97超碰香蕉20202| 日韩熟女老妇一区二区性免费视频| 中文字幕人妻丝袜制服| 91午夜精品亚洲一区二区三区| 精品久久久久久电影网| 国产成人免费无遮挡视频| 老熟女久久久| 午夜福利网站1000一区二区三区| 免费av中文字幕在线| 国产亚洲一区二区精品| 国产毛片在线视频| 亚洲精品日韩在线中文字幕| 亚洲av电影在线观看一区二区三区| 中文字幕制服av| 人人妻人人爽人人添夜夜欢视频| 99国产综合亚洲精品| 亚洲精品日韩在线中文字幕| 免费不卡的大黄色大毛片视频在线观看| 黑人欧美特级aaaaaa片| 国产精品99久久99久久久不卡 | 亚洲国产精品一区三区| 国产片特级美女逼逼视频| 欧美成人午夜免费资源| 亚洲,欧美,日韩| 两性夫妻黄色片 | 亚洲五月色婷婷综合| 人妻 亚洲 视频| 制服丝袜香蕉在线| 亚洲综合色惰| 国产成人欧美| av网站免费在线观看视频| 亚洲av电影在线观看一区二区三区| 99久久人妻综合| 亚洲精品久久久久久婷婷小说| 毛片一级片免费看久久久久| 国产探花极品一区二区| 一级a做视频免费观看| 亚洲精品456在线播放app| 精品熟女少妇av免费看| 大片免费播放器 马上看| 国产不卡av网站在线观看| 99热网站在线观看| 国产黄色视频一区二区在线观看| 天堂中文最新版在线下载| 亚洲欧洲日产国产| 少妇精品久久久久久久| 亚洲国产av影院在线观看| 人妻系列 视频| 99精国产麻豆久久婷婷| 一区二区三区四区激情视频| 亚洲精华国产精华液的使用体验| 有码 亚洲区| 黑人巨大精品欧美一区二区蜜桃 | 美女视频免费永久观看网站| 汤姆久久久久久久影院中文字幕| 国产精品一区www在线观看| 国产精品一区二区在线观看99| www.熟女人妻精品国产 | 久久久久人妻精品一区果冻| 精品一品国产午夜福利视频| av免费在线看不卡| 女的被弄到高潮叫床怎么办| 国产极品天堂在线| a级毛片黄视频| 波野结衣二区三区在线| 久久99热这里只频精品6学生| 秋霞伦理黄片| 国产男人的电影天堂91| 午夜福利视频在线观看免费| 日韩制服骚丝袜av| 中文字幕免费在线视频6| 国产精品 国内视频| 亚洲婷婷狠狠爱综合网| 久久鲁丝午夜福利片| 精品久久国产蜜桃| 高清在线视频一区二区三区| 婷婷色综合大香蕉| 夫妻午夜视频| 久久国产精品男人的天堂亚洲 | 日韩视频在线欧美| 亚洲一级一片aⅴ在线观看| 晚上一个人看的免费电影| 母亲3免费完整高清在线观看 | 亚洲五月色婷婷综合| 亚洲第一区二区三区不卡| 国产1区2区3区精品| 国产伦理片在线播放av一区| www.熟女人妻精品国产 | 久久人人爽av亚洲精品天堂| 亚洲第一av免费看| 七月丁香在线播放| 精品亚洲乱码少妇综合久久| 精品国产一区二区久久| 久久久久久人妻| 91午夜精品亚洲一区二区三区| 国产精品久久久久成人av| 日本黄色日本黄色录像| 欧美另类一区| 欧美变态另类bdsm刘玥| 亚洲内射少妇av| 久久精品国产鲁丝片午夜精品| 亚洲天堂av无毛| 狂野欧美激情性bbbbbb| 边亲边吃奶的免费视频| 美女国产高潮福利片在线看| 亚洲欧美精品自产自拍| 99久国产av精品国产电影| 免费高清在线观看日韩| 国产1区2区3区精品| 欧美+日韩+精品| 激情五月婷婷亚洲| 80岁老熟妇乱子伦牲交| 少妇的逼好多水| 久久久国产欧美日韩av| 一本—道久久a久久精品蜜桃钙片| 久久人人97超碰香蕉20202| 国产精品一国产av| 999精品在线视频| 美女大奶头黄色视频| av在线播放精品| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 制服丝袜香蕉在线| 国产在视频线精品| av.在线天堂| 午夜91福利影院| 久久久久国产网址| 日本爱情动作片www.在线观看| 国产xxxxx性猛交| av免费在线看不卡| 一个人免费看片子| av免费观看日本| 久久99热6这里只有精品| 亚洲av日韩在线播放| 亚洲欧美清纯卡通| 欧美xxⅹ黑人| 中文乱码字字幕精品一区二区三区| 国产有黄有色有爽视频| 又黄又粗又硬又大视频| 欧美xxⅹ黑人| 精品一区二区三区视频在线| 成人毛片a级毛片在线播放| 一区二区三区乱码不卡18| 亚洲成色77777| 久久国产亚洲av麻豆专区| 成年人免费黄色播放视频| av一本久久久久| 亚洲国产看品久久| 欧美精品国产亚洲| 成年动漫av网址| 极品人妻少妇av视频| 午夜精品国产一区二区电影| 五月伊人婷婷丁香| 久久久久久久大尺度免费视频| 亚洲色图 男人天堂 中文字幕 | 满18在线观看网站| 9热在线视频观看99| 男女下面插进去视频免费观看 | 香蕉国产在线看| 国内精品宾馆在线| 国产有黄有色有爽视频| 国产免费一区二区三区四区乱码| 人体艺术视频欧美日本| 大香蕉97超碰在线| 精品福利永久在线观看| 一二三四在线观看免费中文在 | 欧美 亚洲 国产 日韩一| 国产熟女欧美一区二区| 免费观看无遮挡的男女| 国产精品一区二区在线观看99| 男男h啪啪无遮挡| 国产白丝娇喘喷水9色精品| 久久 成人 亚洲| 亚洲中文av在线| 卡戴珊不雅视频在线播放| 国产一区二区激情短视频 | 黑人巨大精品欧美一区二区蜜桃 | 午夜福利乱码中文字幕| 国产亚洲最大av| 久久久精品免费免费高清| 亚洲美女视频黄频| 精品亚洲成国产av| 一级片'在线观看视频| 男女无遮挡免费网站观看| 2018国产大陆天天弄谢| 欧美变态另类bdsm刘玥| 久久人人爽av亚洲精品天堂| 黄色配什么色好看| 午夜福利视频精品| 香蕉丝袜av| 日韩欧美精品免费久久| 国产高清不卡午夜福利| 99久国产av精品国产电影| 大话2 男鬼变身卡| 国产精品久久久久久精品电影小说| 成人毛片a级毛片在线播放| 国产 精品1| av网站免费在线观看视频| 97人妻天天添夜夜摸| 亚洲精品自拍成人| 久久99热6这里只有精品| 日本av免费视频播放| 麻豆精品久久久久久蜜桃| 欧美精品一区二区大全| 亚洲综合精品二区| 人人妻人人澡人人爽人人夜夜| 国产免费视频播放在线视频| 伦精品一区二区三区| 亚洲精品色激情综合| 一级爰片在线观看| 亚洲综合色网址| 精品人妻在线不人妻| 狂野欧美激情性xxxx在线观看| 久久久欧美国产精品| 两个人看的免费小视频| 国产一区二区三区综合在线观看 | 免费观看在线日韩| 国产成人aa在线观看| 2022亚洲国产成人精品| 丝袜在线中文字幕| 亚洲国产精品专区欧美| 亚洲欧美色中文字幕在线| 2018国产大陆天天弄谢| 亚洲,欧美精品.| 久久人人爽av亚洲精品天堂| 看十八女毛片水多多多| 2021少妇久久久久久久久久久| 边亲边吃奶的免费视频| 国产亚洲一区二区精品| 视频区图区小说| 丁香六月天网| 国产亚洲欧美精品永久| 亚洲国产精品999| 日本爱情动作片www.在线观看| 久久久国产精品麻豆| 一级毛片我不卡| 久久99一区二区三区| 97精品久久久久久久久久精品| 18在线观看网站| 黄色 视频免费看| 久久这里只有精品19| 亚洲国产毛片av蜜桃av| 91久久精品国产一区二区三区| 欧美+日韩+精品| 91精品国产国语对白视频| 国产有黄有色有爽视频| 亚洲欧美一区二区三区国产| 波多野结衣一区麻豆| 久久精品国产a三级三级三级| 亚洲成国产人片在线观看| 成年女人在线观看亚洲视频| 国产综合精华液| 伦精品一区二区三区| 在线观看免费高清a一片| 久久女婷五月综合色啪小说| 国产老妇伦熟女老妇高清| 天天躁夜夜躁狠狠久久av| 国产日韩欧美视频二区| 国产一区二区在线观看日韩| av在线老鸭窝| 国产高清国产精品国产三级| 欧美精品国产亚洲| 久久婷婷青草| 黄色怎么调成土黄色| 国产成人精品在线电影| 美女国产高潮福利片在线看| 成人手机av| 在线亚洲精品国产二区图片欧美| 黑人高潮一二区| 少妇被粗大的猛进出69影院 | 另类亚洲欧美激情| 色视频在线一区二区三区| 看免费成人av毛片| 精品久久久精品久久久| 男女啪啪激烈高潮av片| 亚洲精品国产av蜜桃| 精品久久久久久电影网| 午夜福利乱码中文字幕| 国产成人av激情在线播放| 五月开心婷婷网| 久久国内精品自在自线图片| 美女大奶头黄色视频| 99热这里只有是精品在线观看| 亚洲美女黄色视频免费看| 国产精品久久久久久精品古装| 26uuu在线亚洲综合色| av在线老鸭窝| 欧美性感艳星| 99国产综合亚洲精品| 亚洲欧美一区二区三区国产| 精品国产乱码久久久久久小说| tube8黄色片| 成人影院久久| 欧美最新免费一区二区三区| 国产精品一区www在线观看| 亚洲熟女精品中文字幕| 极品少妇高潮喷水抽搐| 国产亚洲最大av| 亚洲中文av在线| 老司机影院成人| 欧美精品一区二区免费开放| 国产亚洲av片在线观看秒播厂| 男女高潮啪啪啪动态图| 久久久国产精品麻豆| 美女国产高潮福利片在线看| 欧美人与性动交α欧美精品济南到 | 一区二区三区乱码不卡18| 国产精品 国内视频| 超碰97精品在线观看| 日本wwww免费看| 亚洲伊人色综图| 桃花免费在线播放| 国产精品不卡视频一区二区| 99九九在线精品视频| 国产女主播在线喷水免费视频网站| 久久久久国产精品人妻一区二区| av福利片在线| 精品视频人人做人人爽| 午夜免费鲁丝| 日本黄大片高清| 亚洲av综合色区一区| 日本wwww免费看| 人妻 亚洲 视频| 国产免费视频播放在线视频| 美女大奶头黄色视频| 亚洲国产成人一精品久久久| 国产亚洲欧美精品永久| 九草在线视频观看| 国产精品久久久久久精品电影小说| 日本色播在线视频| 日韩视频在线欧美| 啦啦啦中文免费视频观看日本| 欧美xxⅹ黑人| 国产成人欧美| 大片电影免费在线观看免费| 免费大片黄手机在线观看| 久久国内精品自在自线图片| 国产精品一国产av| 国产xxxxx性猛交| 天美传媒精品一区二区| 在线看a的网站| 国产69精品久久久久777片| 亚洲一级一片aⅴ在线观看| 国产伦理片在线播放av一区| 婷婷色综合www| 街头女战士在线观看网站| 免费黄频网站在线观看国产| 人人澡人人妻人| 亚洲第一区二区三区不卡| 成人二区视频| 大片电影免费在线观看免费| 少妇高潮的动态图| 最近2019中文字幕mv第一页| 极品人妻少妇av视频| 亚洲精品一二三| 人妻人人澡人人爽人人| av播播在线观看一区| 一级爰片在线观看| 一区二区三区四区激情视频| 男的添女的下面高潮视频| 成年av动漫网址| 有码 亚洲区| 亚洲经典国产精华液单| 最近2019中文字幕mv第一页| 亚洲综合色惰| 高清视频免费观看一区二区| 女性被躁到高潮视频| 日本色播在线视频| 成年女人在线观看亚洲视频| 毛片一级片免费看久久久久| 这个男人来自地球电影免费观看 | 美女视频免费永久观看网站| 少妇的逼好多水| 国产免费现黄频在线看| 人妻系列 视频| 各种免费的搞黄视频| 深夜精品福利| a级毛片黄视频| 精品国产一区二区久久| 欧美老熟妇乱子伦牲交| av不卡在线播放| av福利片在线| 热99久久久久精品小说推荐| 久久精品国产综合久久久 | 国产亚洲精品第一综合不卡 | 狠狠婷婷综合久久久久久88av| 国产精品久久久久久精品古装| 又粗又硬又长又爽又黄的视频| 国产精品.久久久| 最近2019中文字幕mv第一页| 亚洲av.av天堂| 久久精品国产亚洲av涩爱| 9热在线视频观看99| 极品人妻少妇av视频| 国产成人一区二区在线| 亚洲国产色片| 日韩不卡一区二区三区视频在线| kizo精华| 春色校园在线视频观看| 国产一区二区在线观看日韩| 国产又色又爽无遮挡免| 国产精品久久久久成人av| 欧美bdsm另类| 97超碰精品成人国产| 国产精品熟女久久久久浪| 又大又黄又爽视频免费| 宅男免费午夜| 国产老妇伦熟女老妇高清| 亚洲国产av新网站| 欧美激情 高清一区二区三区| 色婷婷久久久亚洲欧美| 久久这里只有精品19| 深夜精品福利| 熟女人妻精品中文字幕| 最近最新中文字幕免费大全7| 成年人免费黄色播放视频| 亚洲精品成人av观看孕妇| 精品久久蜜臀av无| 色94色欧美一区二区| 韩国av在线不卡| 亚洲综合色网址| 久久久久久久久久人人人人人人| 亚洲av综合色区一区| 久久久精品94久久精品| 天堂中文最新版在线下载| 免费黄网站久久成人精品|