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

    地殼應(yīng)變場觀測中體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算*

    2012-09-15 08:15:18張凌空牛安福吳利軍
    地震學(xué)報 2012年4期
    關(guān)鍵詞:空孔應(yīng)變儀鋼筒

    張凌空牛安福吳利軍

    1)中國北京100045中國地震臺網(wǎng)中心

    2)中國北京100036中國地震局地震預(yù)測研究所

    地殼應(yīng)變場觀測中體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算*

    張凌空1),牛安福1)吳利軍2)

    1)中國北京100045中國地震臺網(wǎng)中心

    2)中國北京100036中國地震局地震預(yù)測研究所

    將鉆孔體應(yīng)變儀與分量式應(yīng)變儀的觀測資料在面應(yīng)變平臺上進行整合,對于解決目前兩種應(yīng)變儀觀測數(shù)據(jù)的評估、分析、處理和地震預(yù)報問題有實際意義.根據(jù)潘立宙-陳沅俊和Evertson理論及彈性力學(xué)知識,分別建立了平面應(yīng)力作用下體應(yīng)變和面應(yīng)變觀測的力學(xué)模型,推導(dǎo)了觀測鉆孔、空孔和無孔巖石體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算公式,發(fā)現(xiàn)它們都可以歸結(jié)為同一公式描述,差異僅在于k(體應(yīng)變儀鋼筒內(nèi)壁或無孔巖石面應(yīng)變與空孔巖石面應(yīng)變之比)的取值不同.用Evertson理論推導(dǎo)的公式與空孔巖石情形相近,當(dāng)巖石彈性模量為4 ×1010—8×1010Pa時,二者都可以看成是對潘-陳公式的一種簡化、近似計算;無孔巖石的情況則相當(dāng)于在巖石彈性模量為1×1010Pa時對該式的一種估計.文中結(jié)合實際情況對轉(zhuǎn)換系數(shù)的各種影響因素進行了詳細分析.

    體應(yīng)變 面應(yīng)變 轉(zhuǎn)換系數(shù) 彈性模量 泊松系數(shù)

    Abstract:This study aims at integrating the observational data recorded by borehole volume strainmeters and component borehole strainmeters on the area strain basis.This is helpful for current data evaluation,data processing and analysis,and earthquake prediction.Based on the Pan-Chen model(Pan;Chen and Yang),Evertson theory and theory of elasticity,both volume strain and area strain models were set up under area stress.In addition,we have derived similar formulas for calculating transform coefficient between volume strain and area strain in observation hole,empty hole and non-porous rock condition,respectively.The difference only lies in the factor k(ratio of volume strain gauge on the steel cylinder wall or non-porous rock surface strain to area strain of the empty hole rock).The theoretical formula derived from Evertson model is close to the case of empty hole rock with the rock elastic modulus ranged from 4×1010Pa to 8×1010Pa.Both results may be regarded as a simplified Pan-Chen model.For non-porous rock,the factor k may be predicted by taking rock modulus as 1×1010Pa.This paper discussed some factors related to transform coefficient based on real observation data.

    Key words:volume strain;area strain;transform coefficient;elastic modulus;Poisson coefficient

    引言

    鉆孔應(yīng)變儀具有觀測精度高、響應(yīng)頻帶寬和產(chǎn)出數(shù)據(jù)質(zhì)量好等突出優(yōu)點,現(xiàn)已成為國際上地殼應(yīng)變場觀測的主流儀器.該種儀器分為兩大類:一類是體應(yīng)變儀,另一類是分量式應(yīng)變儀.“十五”期間受美國“板塊邊界觀測計劃(PBO)”的激勵(張寶紅,2004,2010),中國地震局進行了大規(guī)模數(shù)字化地震觀測網(wǎng)絡(luò)建設(shè),鉆孔應(yīng)變儀作為重要的形變前兆觀測儀器開始了大發(fā)展時期,從2007年起陸續(xù)在全國許多臺站投入使用.據(jù)中國地震臺網(wǎng)中心前兆數(shù)據(jù)庫統(tǒng)計,TJ-2型體應(yīng)變儀建立了80個測點,YRY-4型和RZB型分量式應(yīng)變儀分別建立了40個和12個測點,主要分布在首都圈、山西地震帶、東南沿海和南北地震帶等重要地震監(jiān)測區(qū),目前臺點還在不斷增加.如果根據(jù)有關(guān)數(shù)理模型(潘立宙,1981;陳沅俊,楊修信,1990),將體應(yīng)變觀測數(shù)據(jù)轉(zhuǎn)換成面應(yīng)變,分量式應(yīng)變儀經(jīng)過兩組相互垂直的測項加組合和校正(邱澤華等,2005a)也換算成面應(yīng)變,就可以將這兩種應(yīng)變儀測得的巖石面應(yīng)變進行比較,這對解決鉆孔應(yīng)變儀觀測數(shù)據(jù)的評估、分析、處理和地震預(yù)報問題有實際意義.目前,將體應(yīng)變轉(zhuǎn)換成面應(yīng)變的工作在國內(nèi)外還基本是個空白,主要原因是以前使用的分量式應(yīng)變儀比較少,中國大陸大約有8套,美、日兩國大約分別有4套和8套,因而當(dāng)時開展這項研究工作沒有太大的實用價值.有關(guān)學(xué)者(蘇愷之等,2003;邱澤華,石耀霖,2004;邱澤華,2010;歐陽祖熙等,2009;牛安福等,2011)對國內(nèi)外鉆孔應(yīng)變儀的發(fā)展?fàn)顩r已分別進行了系統(tǒng)的總結(jié)和闡述.近些年隨著分量式應(yīng)變儀在中、美兩國的大量布設(shè),體應(yīng)變與面應(yīng)變的轉(zhuǎn)換問題隨之提到議事日程,尤其是2008年汶川地震以后,如何將體應(yīng)變觀測資料在地震監(jiān)測、預(yù)報中更加合理、有效地使用就愈發(fā)顯得重要和緊迫.

    1 觀測井和儀器概況

    目前,我國使用的鉆孔體應(yīng)變儀分為兩種:一種是20世紀80年代中期從美國引進的Sacks-Evertson型(Sacks,Evertson,1971;劉瀾波等,1986),應(yīng)變靈敏度可達1×10-11,井孔深度在250—400 m;另一種是國產(chǎn)TJ型,“九五”以后大量使用的是TJ-2型,它是在Sacks-Evertson型基礎(chǔ)上加以改進的產(chǎn)品(蘇愷之等,2003),其應(yīng)變靈敏度達0.5×10-9,鉆孔深度一般在60—100 m.兩種探頭與井下巖石耦合均采用膨脹水泥固結(jié)方式,當(dāng)巖石對儀器鋼筒產(chǎn)生擠壓作用時,筒內(nèi)硅油的壓力會發(fā)生改變,通過波紋管或液壓傳感器的微小位移可測得體應(yīng)變變化.

    2 鉆孔應(yīng)變觀測力學(xué)模型的建立

    2.1 基于潘-陳理論建立的力學(xué)模型

    設(shè)有一塊無限大的巖石平板,該平板在無窮遠處承受兩項均勻拉應(yīng)力(σ1,σ2),平板中有一半徑為r3的井孔,安置其中的體應(yīng)變儀鋼筒內(nèi)半徑為r1,外半徑為r2,中間通過膨脹水泥與巖石耦合,因而鉆孔的實際情況為雙襯套結(jié)構(gòu).設(shè)鋼筒、水泥、巖石的彈性模量和泊松比分別為E1,ν1,E2,ν2,E3,ν3.假定鉆孔所處介質(zhì)近似為各向同性彈性體,遵從虎克定律,且井口及井底對傳感器的影響可以忽略,鉆孔的軸向應(yīng)力為0.在這種條件下,陳沅俊和楊修信(1990)根據(jù)潘立宙(1981)的工作進一步推導(dǎo)出鋼筒內(nèi)壁上的面應(yīng)變?yōu)?/p>

    M1是與各層材料的彈性參數(shù)及半徑有關(guān)的一個常數(shù),稱為測量系統(tǒng)的特征系數(shù),其值為

    x4由下列四元一次線性方程組確定

    其中系數(shù)為

    解得

    將式(2)代入式(1),得

    考慮到鋼筒、水泥和巖孔耦合的邊界條件,并忽略鋼筒形狀和尺寸等因素的影響,可以認為鋼筒沿其軸向的相對變化量近似等于其周圍巖體的軸向應(yīng)變,即

    則鉆孔體應(yīng)變儀測得的體應(yīng)變?yōu)?/p>

    又知,在平面應(yīng)力作用下,空孔巖石的徑向位移為(王仁,1979)

    θ為任意方向孔徑與最大主應(yīng)力之間的夾角,求得的徑向應(yīng)變?yōu)?/p>

    因為面應(yīng)變等于任意兩組相互垂直的線應(yīng)變之和(θ,θ+90°),所以空孔巖石的面應(yīng)變?yōu)?/p>

    于是

    該式表明k的物理含義是鋼筒內(nèi)壁面應(yīng)變與空孔巖石面應(yīng)變之比.略作變換,有

    將式(11)分別代入式(6)和式(10)

    式(12)、式(13)就是在平面應(yīng)力作用下,基于潘-陳理論建立的鉆孔面應(yīng)變和體應(yīng)變觀測的力學(xué)模型.

    2.2 基于Evertson理論建立的力學(xué)模型

    美國Sacks-Evertson體應(yīng)變儀研制者之一Evertson(1977)總結(jié)出,在三向空間應(yīng)力σ1,σ2,σ3作用下,探頭鋼筒內(nèi)壁上的體應(yīng)變可以寫成

    式中,εV,εH分別是空孔垂直應(yīng)變和面應(yīng)變,0.9是平面應(yīng)力作用下鋼筒內(nèi)壁面應(yīng)變(εm)與空孔巖石面應(yīng)變(εm0)之比(精確值為0.87195),也就是前面所說的k值(k=εm/εm0).如果0.9用k表示,則式(14)可以改寫成

    在平面應(yīng)力作用下(σ3=0),式(15)變?yōu)?/p>

    其面應(yīng)變?yōu)?/p>

    式(16)、式(17)就是基于Evertson理論建立的鉆孔體應(yīng)變和面應(yīng)變觀測的力學(xué)模型.

    2.3 兩組模型的比較

    基于兩種理論建立的面應(yīng)變和體應(yīng)變的力學(xué)模型表述形式完全一樣,差異僅在于k的取值不同.一個是由潘-陳理論推導(dǎo)出的比較復(fù)雜的函數(shù)關(guān)系式k=(x4E3)/E1;另一個是Evertson根據(jù)拉梅厚壁圓筒方程得出的常數(shù)0.87195,比較的焦點在于哪個k值更符合實際情況.將Sacks體應(yīng)變井有關(guān)參數(shù)(E1=19.6×1010Pa,ν1=0.3;E2=2×1010Pa,ν2=0.35;r1=54mm,r2=57mm,r3=76mm)代入k式,則k只與巖石彈性參數(shù)有關(guān).圖1顯示k值并非常量,當(dāng)E3不變時k與ν3近似呈線性關(guān)系,隨其增加而減?。划?dāng)ν3不變時,k隨E3增加而增大;當(dāng)E3=1×1010—3×1010Pa時,k隨E3增加快速上升,幅度達到0.3525,之后速率逐漸減小,各條曲線相互靠近.Evertson理論假定,鋼的彈性模量(E1)是巖石(E3)的3倍(E1=3E3),并且兩者的泊松系數(shù)(ν1=ν3)相等.事實上,E1=19.6× 1010Pa,ν1=0.3,E3=1×1010—8×1010Pa,ν3=0.20—0.35,假設(shè)條件不是嚴格成立,尤其是E3波動范圍較大,因此模型只是通常情況下的平均估計,0.87195近似等于由圖1得到的平均k值(0.8453).當(dāng)巖石較軟,即E3=1×1010—3×1010Pa時,兩數(shù)值明顯偏離很多.潘-陳模型考慮了各種因素對k值的影響,理論公式推導(dǎo)比較嚴格,不存在此類假設(shè)問題,從這個角度分析,潘-陳模型應(yīng)更為合理.

    圖1 鋼筒內(nèi)壁面應(yīng)變和空孔巖石面應(yīng)變之比(k)與巖石彈性參數(shù)(E3,ν3)的關(guān)系Fig.1 Relation between the ratio of steel tube wall plane strain to strain outside a borehole in an infinite rock plate(k)and the rock elastic parameters(E3,ν3)

    3 鉆孔體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算

    根據(jù)式(12)、(13)或式(16)、(17),可以得到

    式中,f是鉆孔體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù),利用這一關(guān)系可以將體應(yīng)變測量值換算成面應(yīng)變,從而可與分量式應(yīng)變儀測得的面應(yīng)變進行比較.

    當(dāng)k=(x4E3)/E1時,就是基于潘-陳力學(xué)模型下的體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)

    當(dāng)k=0.9時,就是基于Evertson力學(xué)模型下的體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)

    將Sacks井有關(guān)參數(shù)代入式(20),以ν3為自變量、f1為因變量作曲線(圖2a),可知f1隨ν3增加而增大.當(dāng)E3=4×1010—8×1010Pa時,f1各條曲線近于重合;而由式(21)算得的f2曲線(虛線)也恰恰分布于這一區(qū)域,說明f2是f1的一種特殊情況.通常巖石的E3在4×1010—8×1010Pa,ν3在0.20—0.35,這時k值接近0.87195,f2≈f1,所以一般情況下f2可以取代f1,并且計算簡單,只涉及一個參數(shù).但當(dāng)巖石較軟、易于發(fā)生彈性變形時,即E3<3×1010Pa,f2與f1誤差較大.國產(chǎn)體應(yīng)變儀研制者蘇愷之等(2003)給出了TJ-2型體應(yīng)變觀測井的有關(guān)參數(shù),E1=21×1010Pa,ν1=0.3;E2=3×1010Pa,ν2=0.25;r1=42mm,r2=44.5mm,r3=65mm,將這些參數(shù)代入式(20),同樣可以做出類似圖2a的曲線(圖2b),只是k值略小、f1值略大.

    圖2 兩種體應(yīng)變儀的f1與f2(圖中4與8之間的虛線)比較(a)Sacks-Evertson型;(b)TJ-2型Fig.2 Comparison between f1and f2(dashed line in two figures)for two kinds of borehole volume strainmeters(a)Sacks-Evertson type;(b)TJ-2 type

    根據(jù)寶坻等臺巖石彈性參數(shù)實際測定結(jié)果,分別對各臺站f1和f2進行計算(表1).昌平臺k值明顯偏小,導(dǎo)致f1大于f2,E3過小是造成這種差距的主要原因;其余6個臺站k值與0.87195比較接近,因而求得的f2值與f1誤差較小.

    表1 用兩個公式計算的實際臺站體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的對比Table 1 Comparison between conversion factors from volume to area strain estimated by two formulas for several stations

    4 鉆孔體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的影響因素分析

    4.1 與巖石彈性參數(shù)的關(guān)系

    圖2b表明f1隨ν3的增加而增大,隨E3的增加而減小,說明巖石的剛度越小、橫向變形越大,體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)越大;當(dāng)E3=1×1010—2×1010Pa時,f1隨ν3的增加上升速率很快,之后速率迅速減?。划?dāng)E3=4×1010—8×1010Pa時,f1各曲線接近線性,且近于重合.一般實測鉆孔巖石E3=4×1010—8×1010Pa,ν3=0.25—0.32,則f1=1.138—1.226,f2=1.167—1.224,平均值f1=1.182,f2=1.195.

    4.2 與膨脹水泥彈性參數(shù)的關(guān)系

    假定巖石的彈性參數(shù)不變,取E3=5×1010Pa,ν3=0.25,分析f1與E2,ν2的關(guān)系.圖3表明,f1隨ν2的增加而減小,隨E2的增加而增大,且基本呈線性關(guān)系,即水泥剛度越大、橫向變形越小,體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)越大.一般E2=2×1010Pa(水泥較稀)—3.5×1010Pa(水泥較硬),ν2=0.25—0.35.水泥彈性參數(shù)發(fā)生變化時,f1=1.167—1.157,最大改變量是0.01,說明膨脹水泥對體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的影響不是很大.

    圖3 轉(zhuǎn)換系數(shù)(f1)與膨脹水泥彈性參數(shù)(E2,ν2)的關(guān)系Fig.3 Relation between conversion factor(f1)and elastic parameters(E2,ν2)of expanding cement

    4.3 與儀器鋼筒半徑和鉆孔孔徑的關(guān)系

    兩代TJ型體應(yīng)變儀鋼筒內(nèi)半徑變化范圍是r1=42—54mm,外半徑r2=44—57mm,鉆孔半徑r3=55—75mm,若保持n=r2/r1=1.060不變(蘇愷之等,2003),則由式(19)可做出圖4.f1與r1近似呈線性關(guān)系,隨其增大而增大,在1.164—1.174之間變化;f1與r3呈非線性關(guān)系,隨其增大而減小,在1.161—1.173之間變化.可見,r1和r3增減對f1最多影響0.01.

    圖4 轉(zhuǎn)換系數(shù)(f1)與儀器鋼筒半徑(r1)(a)和鉆孔孔徑(r3)(b)的關(guān)系Fig.4 Relation between conversion factor(f1)and strainmeter steel tube radius(r1)(a)and borehole radius(r3)(b)

    5 空孔巖石體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算

    根據(jù)式(9)可知,在平面應(yīng)力作用下空孔巖石的面應(yīng)變?yōu)?/p>

    又根據(jù)式(5)可知,空孔巖石體應(yīng)變?yōu)?/p>

    故體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)為

    顯然,f0也可以寫成式(18)的形式,只是這時k=1.

    6 完整的無孔巖石體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算

    在平面應(yīng)力作用下,無孔巖石的線應(yīng)變有

    其面應(yīng)變?yōu)?/p>

    體應(yīng)變?yōu)?/p>

    故體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)為

    說明fw同樣也可以寫成式(18)形式,這時k=(1-ν3)/2,其物理含義是巖石面應(yīng)變與空孔巖石面應(yīng)變之比.如果取ν3=0.25,則fw=1.5,故εT=2εm/3,2/3就是由面應(yīng)變固體潮理論值計算體應(yīng)變固體潮理論值時常用的轉(zhuǎn)換系數(shù).

    7 三組體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)之間的比較

    圖5顯示fw與ν3呈非線性關(guān)系,隨其增加而增大,并且當(dāng)ν3>0.32后速率明顯加快.與圖2b比較,fw曲線類似于E3=1×1010Pa時f1的情況.當(dāng)ν3=0.2—0.32時,fw與f1近于重合,之后兩曲線差距才逐漸拉大;f2與f0曲線很接近,數(shù)值平均比f0大0.032,兩者都可以近似看成是E3=4×1010—8×1010Pa時f1的特殊情況.

    圖5 三種體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的對比Fig.5 Comparison between three conversion factors from volume to area strain in crustal strain observation

    8 討論與結(jié)論

    無論是地殼觀測鉆孔、空孔還是無孔巖石,其體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)都可以表述成fi=1/[1-ν3/(2k)]的形式(i=1,2,0,w),k的物理意義是儀器鋼筒內(nèi)壁面應(yīng)變或無孔巖石面應(yīng)變與空孔巖石面應(yīng)變之比.

    1)當(dāng)k=(x4E3)/E1,且觀測鉆孔有關(guān)參數(shù)(r1,r2,r3,E1,E2,ν1,ν2)為定值時,k只與巖石彈性模量(E3)和泊松系數(shù)(ν3)有關(guān),為一變量,f1就是基于潘立宙-陳沅俊力學(xué)模型推導(dǎo)出的鉆孔體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算式.

    2)當(dāng)k=0.87195時,f2就是基于Evertson力學(xué)模型推導(dǎo)出的鉆孔體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)計算式.0.87195近似等于用k=(x4E3)/E1公式計算得到的平均k值,f2相當(dāng)于E3=4×1010—8×1010Pa時f1的一種特殊情況.當(dāng)E3<3×1010Pa時f2與f1誤差較大,這主要是源于“鋼的彈性模量是巖石的3倍,并且二者的泊松系數(shù)相等”的假定.因此,f2不夠準確、全面地反映實際情況,它是在一般情況下對f1的一種近似、簡化計算.

    3)當(dāng)k=1時,f0是空孔巖石體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算式.f0與f2曲線很接近,但數(shù)值略小.

    4)當(dāng)k=(1-ν3)/2時,fw是無孔巖石體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算式.在E3=1 ×1010Pa,ν3=0.2—0.32時,fw與f1曲線比較接近.

    5)fi公式是在平面應(yīng)力作用下推導(dǎo)出來的,對潮汐觀測研究完全適用,但當(dāng)巖石平板存在垂直向應(yīng)力時(如大氣壓)需要另作考慮.大氣壓對鉆孔應(yīng)變觀測是一個不間斷的長期干擾源,這是今后在進行臺站數(shù)據(jù)處理時需要深入加以研究的問題,如作潮汐分析以前要首先排除日波、半日波氣壓的干擾.已知氣壓影響系數(shù)與氣壓周期、鉆孔深度和巖石彈性參數(shù)等因素有關(guān)(張凌空等,2011),下一階段的任務(wù)是如何以現(xiàn)有的氣壓干擾模型為基礎(chǔ),進一步推導(dǎo)在三向空間大氣壓力作用下體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算式.美國實施PBO計劃后,近年來在西部地區(qū)已安裝了80多臺GTSM型分量式鉆孔應(yīng)變儀,Roeloffs(2010)對華盛頓、俄勒岡和南加州12個臺站的資料進行了分析,發(fā)現(xiàn)應(yīng)變儀4個分量元件對氣壓變化的響應(yīng)范圍約在1.5×10-9—11.7×10-9hPa,這與中國和日本體應(yīng)變儀(蘇愷之等,2003;未廣重二,1979)的觀測結(jié)果大體一致.隨著體應(yīng)變儀和分量式應(yīng)變儀觀測資料的不斷積累,中美兩國今后都應(yīng)加強兩種應(yīng)變儀面應(yīng)變觀測資料的比對研究,深入分析二者之間的共性和細節(jié)上的差異,并做出合理解釋.

    6)當(dāng)儀器鋼筒、膨脹水泥和鉆孔的有關(guān)參數(shù)確定后,f1只與巖石彈性模量E3和泊松系數(shù)ν3有關(guān),f1隨ν3的增加而增大,隨E3的增加而減小,并且在E3<3×1010Pa時對f1影響很大;當(dāng)膨脹水泥的彈性模量E2和泊松系數(shù)ν2有變化,而其它參數(shù)確定后,f1隨ν2的增加而減小,隨E2的增加而增大,并且對f1的最大改變量是0.01;若保持鋼筒的內(nèi)外半徑比不變,則f1隨內(nèi)半徑r1增大而增大,隨孔徑r3增大而減小,并且二者對f1最大影響量都是0.01.

    7)體應(yīng)變轉(zhuǎn)換成面應(yīng)變后還不能直接與其它臺站的面應(yīng)變資料進行比較(除非它們的觀測條件相近),還需要求出各自的耦合系數(shù)(理論觀測面應(yīng)變/無孔巖石面應(yīng)變=2k/(1-ν3)),再將實測面應(yīng)變除以耦合系數(shù)就得到真實的巖石面應(yīng)變,這一步又稱作面應(yīng)變觀測的實地絕對標定(邱澤華等,2005b),至此各臺站之間的面應(yīng)變才具有絕對的可比性.

    8)井下探頭附近巖石彈性參數(shù)的測定是一項非常重要的基礎(chǔ)工作,不僅鉆孔體應(yīng)變與面應(yīng)變的轉(zhuǎn)換系數(shù)計算要用到它,而且二者的耦合系數(shù)計算也要用到它.目前除少數(shù)臺站外大部分臺站還沒有開展這項工作,這是今后需要專項基金給與大力支持的科研項目.

    邱澤華研究員和審稿專家對本文進行了認真審閱,并提出了建設(shè)性的修改意見,在此謹致誠摯謝意.

    陳沅俊,楊修信.1990.雙襯套鉆孔應(yīng)變測量的計算[J].華北地震科學(xué),8(4):80--89.

    劉瀾波,何世海,寧長平,劉北順.1986.Sacks-Evertson鉆孔應(yīng)變儀觀測網(wǎng)在京津地區(qū)安裝[J].中國地震,2(4):102--104.

    牛安福,張凌空,閆偉,吉平.2011.中國鉆孔應(yīng)變觀測能力及在地震預(yù)報中的應(yīng)用[J].大地測量與地球動力學(xué),31(2):48--52.

    歐陽祖熙,張鈞,陳征.2009.地殼形變深井綜合觀測技術(shù)的新進展[J].國際地震動態(tài),(11):1--13.

    潘立宙.1981.與地應(yīng)力測量有關(guān)的幾個公式的推導(dǎo)和討論[M]∥地應(yīng)力測量原理與應(yīng)用.北京:地質(zhì)出版社:163--202.

    邱澤華,石耀霖.2004.國外鉆孔應(yīng)變觀測的發(fā)展現(xiàn)狀[J].地震學(xué)報,26(增刊):162--168.

    邱澤華,石耀霖,歐陽祖熙.2005a.四分量鉆孔應(yīng)變觀測的實地相對標定[J].大地測量與地球動力學(xué),25(1):118--122.

    邱澤華,石耀霖,歐陽祖熙.2005b.四分量鉆孔應(yīng)變觀測的實地絕對標定[J].地震,25(3):27--34.

    邱澤華.2010.中國分量鉆孔地應(yīng)力-應(yīng)變觀測發(fā)展重要事件回顧[J].大地測量與地球動力學(xué),30(5):42--47.

    蘇愷之,李海亮,張均,李秀環(huán),馬相波.2003.鉆孔地應(yīng)變觀測新進展[M].北京:地震出版社:1--56.

    王仁.1979.固體力學(xué)基礎(chǔ)[M].北京:地質(zhì)出版社:188.

    張寶紅.2004.美國的板塊邊界觀測(PBO)計劃[J].大地測量與地球動力學(xué),24(3):105--108.

    張寶紅.2010.美國板塊邊界觀測(PBO)中的鉆孔應(yīng)變觀測設(shè)備[J].大地測量與地球動力學(xué),30(增刊Ⅰ):40--44.

    張凌空,王廣才,牛安福.2011.周期氣壓波對地殼應(yīng)變場觀測影響的若干因素分析[J].地震學(xué)報,33(3):351--361.

    周龍壽,邱澤華,唐磊.2008.地殼應(yīng)變場對氣壓短周期變化的響應(yīng)[J].地球物理學(xué)進展,23(6):1717--1726.

    未廣重二.1979.地殻變動連續(xù)観測と埋入式歪計(Ⅰ)[J].測候時報,46(1--2):9--25.

    Evertson D W.1977.Borehole strainmeters for seismology[R]∥Rep.ARL-TR-77-62.Applied research Lab,University of Texas,Audtin,Texas:1--141.

    Roeloffs E.2010.Tidal calibration of plate boundary observatory borehole strainmeters:Roles of vertical and shear coupling[J].J Geophys Res,115(B06405):25.

    Sacks I S,Evertson D W.1971.Sensitive borehole strain-rate meters[M]∥Carnegie Institution,Washington Yearbook.Washington:Carnegie Institution,68:448--453.

    Conversion factors from volume to area strain in crustal strain observation

    Zhang Lingkong1),Niu Anfu1)Wu Lijun2)
    1)China Earthquake Networks Center,Beijing 100045,China
    2)Institute of Earthquake Science,China Earthquake Administration,Beijing 100036,China

    10.3969/j.issn.0253-3782.2012.04.005

    P315.72+7

    A

    地震行業(yè)科研專項(201108009)資助.

    2011-08-24收到初稿,2011-12-23決定采用修改稿.

    e-mail:zhll1023@163.com

    張凌空,牛安福,吳利軍.2012.地殼應(yīng)變場觀測中體應(yīng)變與面應(yīng)變轉(zhuǎn)換系數(shù)的計算.地震學(xué)報,34(4):476--486.

    Zhang Lingkong,Niu Anfu,Wu Lijun.2012.Conversion factors from volume to area strain in crustal strain observation.Acta Seismologica Sinica,34(4):476--486.

    猜你喜歡
    空孔應(yīng)變儀鋼筒
    體應(yīng)變儀傳感器更換前后數(shù)據(jù)關(guān)聯(lián)性分析
    不耦合裝藥下炮孔—空孔距離對預(yù)裂爆破效果的影響探究
    工程爆破(2023年5期)2023-11-02 12:59:44
    TJ-2型體積式應(yīng)變儀故障檢測裝置
    水利工程預(yù)應(yīng)力鋼筒混凝土管耐久性研究及對策
    3種鉆孔應(yīng)變儀記錄的瑪多MS7.4地震同震響應(yīng)特征對比分析
    中國地震(2021年3期)2022-01-12 09:42:00
    空孔作用下定向漲裂破巖試驗研究
    光纖應(yīng)變儀的研制與應(yīng)用前景
    等體積空孔直眼掏槽槽腔形成過程及其分析*
    爆破(2020年4期)2020-12-16 09:09:52
    定向斷裂控制爆破的空孔效應(yīng)實驗分析*
    爆炸與沖擊(2015年3期)2015-06-07 11:38:12
    內(nèi)部爆炸作用下鋼筒變形過程的電探針測量技術(shù)*
    爆炸與沖擊(2014年1期)2014-12-12 06:24:36
    最黄视频免费看| 国产成人啪精品午夜网站| 亚洲天堂av无毛| 桃花免费在线播放| 丝袜美腿诱惑在线| 脱女人内裤的视频| 狂野欧美激情性xxxx| 男人操女人黄网站| 精品一区在线观看国产| 日韩 亚洲 欧美在线| 国产一区有黄有色的免费视频| 成人亚洲精品一区在线观看| 精品久久久久久电影网| 啦啦啦视频在线资源免费观看| 一边摸一边做爽爽视频免费| 国产午夜精品一二区理论片| 嫁个100分男人电影在线观看 | 国产野战对白在线观看| 亚洲久久久国产精品| 1024视频免费在线观看| 免费av中文字幕在线| 欧美精品啪啪一区二区三区 | 久久久精品免费免费高清| 99国产综合亚洲精品| 亚洲一卡2卡3卡4卡5卡精品中文| 丁香六月欧美| 欧美日韩亚洲高清精品| 看十八女毛片水多多多| 免费不卡黄色视频| 又黄又粗又硬又大视频| 亚洲精品国产区一区二| 超碰成人久久| 国产一区二区激情短视频 | av线在线观看网站| 性高湖久久久久久久久免费观看| 日本vs欧美在线观看视频| 91精品国产国语对白视频| 免费在线观看影片大全网站 | 免费久久久久久久精品成人欧美视频| 狂野欧美激情性bbbbbb| 国产精品久久久人人做人人爽| 国产成人精品无人区| 亚洲成国产人片在线观看| av视频免费观看在线观看| 欧美激情极品国产一区二区三区| 18禁黄网站禁片午夜丰满| 看免费av毛片| 久久久久久久大尺度免费视频| www.999成人在线观看| 日本av免费视频播放| 久久国产精品人妻蜜桃| 国产精品一二三区在线看| kizo精华| 久9热在线精品视频| 精品一区二区三区av网在线观看 | 国产一区二区 视频在线| 麻豆乱淫一区二区| 亚洲,欧美精品.| 视频在线观看一区二区三区| 两个人看的免费小视频| 亚洲成人免费电影在线观看 | 亚洲国产最新在线播放| 久久国产精品男人的天堂亚洲| 成人亚洲欧美一区二区av| 菩萨蛮人人尽说江南好唐韦庄| 99久久综合免费| 黑人巨大精品欧美一区二区蜜桃| 叶爱在线成人免费视频播放| 欧美日韩亚洲国产一区二区在线观看 | 久久久久久久国产电影| 欧美日韩综合久久久久久| 极品人妻少妇av视频| 欧美乱码精品一区二区三区| 日本五十路高清| 亚洲成色77777| 日本wwww免费看| 国产免费一区二区三区四区乱码| 国产在线观看jvid| 亚洲国产欧美一区二区综合| 成人影院久久| 亚洲欧美日韩另类电影网站| 国产精品二区激情视频| 黄色片一级片一级黄色片| 99国产精品99久久久久| 欧美日韩黄片免| 人人妻人人爽人人添夜夜欢视频| 国产色视频综合| 欧美精品av麻豆av| 婷婷色av中文字幕| 成人免费观看视频高清| xxx大片免费视频| 在现免费观看毛片| 欧美黄色片欧美黄色片| 午夜老司机福利片| 丝袜脚勾引网站| 久热这里只有精品99| 国产黄色免费在线视频| 国精品久久久久久国模美| 在线亚洲精品国产二区图片欧美| 久久女婷五月综合色啪小说| 中文字幕色久视频| 亚洲欧美激情在线| 成人免费观看视频高清| 啦啦啦在线观看免费高清www| 午夜激情av网站| 国产精品久久久人人做人人爽| 亚洲激情五月婷婷啪啪| 50天的宝宝边吃奶边哭怎么回事| 人人妻,人人澡人人爽秒播 | 日日摸夜夜添夜夜爱| 久久国产亚洲av麻豆专区| 黄色视频不卡| 精品福利永久在线观看| 嫁个100分男人电影在线观看 | 亚洲人成电影免费在线| 精品国产乱码久久久久久小说| 日韩中文字幕视频在线看片| 国产成人av激情在线播放| 夜夜骑夜夜射夜夜干| 999精品在线视频| 在线观看免费视频网站a站| videosex国产| 汤姆久久久久久久影院中文字幕| 777米奇影视久久| 丁香六月欧美| 国产日韩欧美视频二区| 日韩精品免费视频一区二区三区| 久久精品人人爽人人爽视色| 青草久久国产| www.精华液| 秋霞在线观看毛片| 性高湖久久久久久久久免费观看| 欧美中文综合在线视频| 黄网站色视频无遮挡免费观看| 久久精品国产亚洲av涩爱| 黄色片一级片一级黄色片| 我要看黄色一级片免费的| 男人添女人高潮全过程视频| 在线观看免费午夜福利视频| 国产老妇伦熟女老妇高清| 欧美变态另类bdsm刘玥| 又大又黄又爽视频免费| 国产在视频线精品| kizo精华| 精品国产一区二区久久| 精品福利观看| 在线精品无人区一区二区三| 国产成人欧美在线观看 | 视频区图区小说| 免费在线观看视频国产中文字幕亚洲 | 国产91精品成人一区二区三区 | 18在线观看网站| 国产成人av激情在线播放| 欧美+亚洲+日韩+国产| 中文字幕人妻丝袜制服| 制服人妻中文乱码| 在线观看人妻少妇| 国产精品成人在线| www日本在线高清视频| 国产精品国产av在线观看| 亚洲人成电影免费在线| 狂野欧美激情性bbbbbb| 99精国产麻豆久久婷婷| 久久影院123| 老熟女久久久| 美女大奶头黄色视频| 中国国产av一级| 老熟女久久久| 成人手机av| 精品国产国语对白av| 午夜免费男女啪啪视频观看| 国产欧美日韩综合在线一区二区| 777久久人妻少妇嫩草av网站| 人妻一区二区av| 午夜老司机福利片| 午夜两性在线视频| 亚洲七黄色美女视频| 搡老乐熟女国产| 国产爽快片一区二区三区| 国产福利在线免费观看视频| 精品高清国产在线一区| 尾随美女入室| 大码成人一级视频| 中文字幕制服av| 性色av乱码一区二区三区2| www.熟女人妻精品国产| 国产日韩欧美亚洲二区| 欧美人与性动交α欧美精品济南到| 天天躁日日躁夜夜躁夜夜| 亚洲第一av免费看| 中文字幕人妻熟女乱码| 一级a爱视频在线免费观看| 成年人免费黄色播放视频| 精品久久蜜臀av无| 国产高清视频在线播放一区 | www.精华液| 中文字幕制服av| 欧美在线一区亚洲| 亚洲中文日韩欧美视频| 中文字幕人妻丝袜一区二区| 91麻豆av在线| 丝袜在线中文字幕| 亚洲精品自拍成人| 午夜福利影视在线免费观看| 国产精品九九99| 中文精品一卡2卡3卡4更新| 蜜桃国产av成人99| 久久久久久久久久久久大奶| 91九色精品人成在线观看| 亚洲人成网站在线观看播放| 少妇的丰满在线观看| 精品国产超薄肉色丝袜足j| 大陆偷拍与自拍| 国产高清videossex| 伊人久久大香线蕉亚洲五| 天天操日日干夜夜撸| 国精品久久久久久国模美| avwww免费| 最黄视频免费看| 亚洲欧美中文字幕日韩二区| 19禁男女啪啪无遮挡网站| 大码成人一级视频| 亚洲精品一区蜜桃| 亚洲天堂av无毛| a级毛片黄视频| 成人免费观看视频高清| 国产一区二区 视频在线| 波多野结衣av一区二区av| 亚洲熟女精品中文字幕| 亚洲欧美日韩另类电影网站| 男女下面插进去视频免费观看| 看免费av毛片| 亚洲,欧美精品.| 久久久久国产一级毛片高清牌| 欧美成人精品欧美一级黄| 大话2 男鬼变身卡| 亚洲av片天天在线观看| 久久久久精品人妻al黑| 激情五月婷婷亚洲| 国产精品香港三级国产av潘金莲 | 国产日韩欧美在线精品| 一区二区av电影网| 老司机靠b影院| 中文字幕人妻丝袜制服| 国产免费一区二区三区四区乱码| 成年人免费黄色播放视频| 国产精品 欧美亚洲| 成人国产一区最新在线观看 | 亚洲av成人不卡在线观看播放网 | 青春草亚洲视频在线观看| 无遮挡黄片免费观看| 欧美精品一区二区免费开放| 亚洲精品一卡2卡三卡4卡5卡 | 麻豆av在线久日| 免费一级毛片在线播放高清视频 | 1024香蕉在线观看| 精品一区二区三区四区五区乱码 | 丝袜喷水一区| 日韩,欧美,国产一区二区三区| 女性生殖器流出的白浆| 亚洲人成电影免费在线| 青春草亚洲视频在线观看| 老汉色∧v一级毛片| 久久鲁丝午夜福利片| 亚洲色图 男人天堂 中文字幕| 啦啦啦视频在线资源免费观看| 久久精品国产a三级三级三级| 亚洲九九香蕉| 婷婷丁香在线五月| 又紧又爽又黄一区二区| 亚洲av美国av| 欧美日韩视频精品一区| 免费一级毛片在线播放高清视频 | 欧美日韩视频精品一区| 王馨瑶露胸无遮挡在线观看| 国产一卡二卡三卡精品| 国产精品免费大片| 亚洲国产看品久久| 日本wwww免费看| 蜜桃国产av成人99| 久久九九热精品免费| 曰老女人黄片| 免费女性裸体啪啪无遮挡网站| 免费在线观看日本一区| 欧美老熟妇乱子伦牲交| 欧美少妇被猛烈插入视频| 在线观看一区二区三区激情| 色婷婷av一区二区三区视频| 成年av动漫网址| 亚洲美女黄色视频免费看| 老司机深夜福利视频在线观看 | 免费不卡黄色视频| 欧美日韩亚洲高清精品| 免费观看a级毛片全部| 亚洲欧美精品综合一区二区三区| 男女下面插进去视频免费观看| 99九九在线精品视频| 国产精品二区激情视频| 国产人伦9x9x在线观看| 80岁老熟妇乱子伦牲交| av又黄又爽大尺度在线免费看| 亚洲国产精品国产精品| 少妇人妻久久综合中文| 日韩制服丝袜自拍偷拍| 国产欧美日韩综合在线一区二区| 在线观看www视频免费| av网站免费在线观看视频| 我要看黄色一级片免费的| 午夜福利免费观看在线| 校园人妻丝袜中文字幕| 久久国产精品男人的天堂亚洲| 国产有黄有色有爽视频| 午夜福利一区二区在线看| 视频在线观看一区二区三区| 中文精品一卡2卡3卡4更新| 波多野结衣一区麻豆| 脱女人内裤的视频| 人人妻,人人澡人人爽秒播 | 最新的欧美精品一区二区| 天堂8中文在线网| 真人做人爱边吃奶动态| 精品少妇一区二区三区视频日本电影| 在线观看人妻少妇| 国产精品一区二区在线观看99| 久久99精品国语久久久| 欧美另类一区| 秋霞在线观看毛片| 久久久精品免费免费高清| 精品人妻熟女毛片av久久网站| 后天国语完整版免费观看| 国产高清国产精品国产三级| a 毛片基地| 国产淫语在线视频| 亚洲成色77777| 国产男女内射视频| 老司机靠b影院| 成年美女黄网站色视频大全免费| 亚洲国产中文字幕在线视频| 美女视频免费永久观看网站| 精品久久久久久久毛片微露脸 | cao死你这个sao货| xxxhd国产人妻xxx| 精品国产一区二区三区四区第35| 男人操女人黄网站| 午夜福利影视在线免费观看| 欧美人与善性xxx| 欧美日韩亚洲国产一区二区在线观看 | 18在线观看网站| 国产精品 国内视频| av天堂久久9| 成年av动漫网址| 免费女性裸体啪啪无遮挡网站| 久久久久久久久免费视频了| 乱人伦中国视频| 午夜av观看不卡| 丝袜脚勾引网站| 日本黄色日本黄色录像| 欧美国产精品一级二级三级| 精品人妻1区二区| 中文欧美无线码| 欧美国产精品va在线观看不卡| 90打野战视频偷拍视频| 激情五月婷婷亚洲| 午夜福利视频在线观看免费| 女警被强在线播放| 一级黄色大片毛片| 叶爱在线成人免费视频播放| 女警被强在线播放| 国产精品久久久久久精品电影小说| 亚洲,欧美精品.| 国产精品秋霞免费鲁丝片| 九色亚洲精品在线播放| 女人久久www免费人成看片| 伦理电影免费视频| 999久久久国产精品视频| 视频区欧美日本亚洲| 国产片内射在线| 在线观看人妻少妇| 日本一区二区免费在线视频| 少妇 在线观看| 美女高潮到喷水免费观看| 亚洲欧洲日产国产| 黑人巨大精品欧美一区二区蜜桃| 十八禁人妻一区二区| 国产免费一区二区三区四区乱码| 国产欧美日韩一区二区三 | 国产精品久久久久成人av| 久久精品国产a三级三级三级| 免费高清在线观看日韩| 成年女人毛片免费观看观看9 | 欧美成人午夜精品| 老司机影院毛片| 国产一区二区三区av在线| 欧美日韩福利视频一区二区| 丰满少妇做爰视频| 赤兔流量卡办理| 国产精品国产三级专区第一集| 两个人看的免费小视频| av国产久精品久网站免费入址| 在线观看人妻少妇| 男女边摸边吃奶| 不卡av一区二区三区| www.999成人在线观看| 久久久久视频综合| 国产老妇伦熟女老妇高清| 亚洲欧美精品综合一区二区三区| 性高湖久久久久久久久免费观看| 青青草视频在线视频观看| 18禁观看日本| 天天操日日干夜夜撸| 男女国产视频网站| 男人爽女人下面视频在线观看| 成人免费观看视频高清| 国产人伦9x9x在线观看| 国产在线观看jvid| 婷婷色av中文字幕| 黄色一级大片看看| 欧美精品av麻豆av| 国产成人系列免费观看| 麻豆国产av国片精品| 啦啦啦啦在线视频资源| 国产成人精品久久久久久| 美国免费a级毛片| 久久亚洲国产成人精品v| 又粗又硬又长又爽又黄的视频| 亚洲av综合色区一区| 永久免费av网站大全| av国产久精品久网站免费入址| 国产亚洲欧美在线一区二区| 91老司机精品| 日韩电影二区| 各种免费的搞黄视频| 丝瓜视频免费看黄片| 日本午夜av视频| 妹子高潮喷水视频| 色精品久久人妻99蜜桃| 天天躁日日躁夜夜躁夜夜| 久久av网站| 中文精品一卡2卡3卡4更新| 亚洲av在线观看美女高潮| 中文字幕另类日韩欧美亚洲嫩草| 1024香蕉在线观看| av福利片在线| 日韩 亚洲 欧美在线| 精品久久久久久久毛片微露脸 | 国产免费福利视频在线观看| 亚洲自偷自拍图片 自拍| 色网站视频免费| 欧美乱码精品一区二区三区| 人体艺术视频欧美日本| 大片免费播放器 马上看| 深夜精品福利| 亚洲av电影在线观看一区二区三区| 亚洲精品久久成人aⅴ小说| 精品久久久久久电影网| 亚洲精品日本国产第一区| 国产精品人妻久久久影院| av一本久久久久| 亚洲av在线观看美女高潮| 69精品国产乱码久久久| 亚洲国产精品成人久久小说| 天天添夜夜摸| 黑人欧美特级aaaaaa片| 国产一卡二卡三卡精品| av一本久久久久| 欧美激情 高清一区二区三区| 亚洲欧洲精品一区二区精品久久久| 国产不卡av网站在线观看| 国产精品欧美亚洲77777| 国产精品香港三级国产av潘金莲 | 热re99久久精品国产66热6| 男人操女人黄网站| 免费在线观看影片大全网站 | 午夜日韩欧美国产| 国产在线免费精品| 2018国产大陆天天弄谢| 999久久久国产精品视频| 国产熟女午夜一区二区三区| 一区二区三区激情视频| 欧美日韩亚洲综合一区二区三区_| 啦啦啦视频在线资源免费观看| 赤兔流量卡办理| 成人免费观看视频高清| 午夜福利一区二区在线看| 亚洲精品乱久久久久久| 亚洲熟女毛片儿| 日本a在线网址| av国产久精品久网站免费入址| 黄色视频不卡| 50天的宝宝边吃奶边哭怎么回事| 欧美黑人精品巨大| 大片电影免费在线观看免费| 日本欧美视频一区| 国产激情久久老熟女| h视频一区二区三区| 精品欧美一区二区三区在线| av线在线观看网站| 人人妻人人澡人人看| 国产精品久久久久久精品古装| 久久精品成人免费网站| 一边摸一边抽搐一进一出视频| 国产一区二区三区综合在线观看| 国产片内射在线| tube8黄色片| 中文字幕色久视频| 国产真人三级小视频在线观看| 免费久久久久久久精品成人欧美视频| 青春草视频在线免费观看| 又粗又硬又长又爽又黄的视频| 日本欧美视频一区| 电影成人av| 国产亚洲精品久久久久5区| 久久99热这里只频精品6学生| 免费在线观看完整版高清| 丰满人妻熟妇乱又伦精品不卡| 精品一品国产午夜福利视频| 最新在线观看一区二区三区 | 丝袜喷水一区| 最新在线观看一区二区三区 | 久久天堂一区二区三区四区| 最黄视频免费看| 又大又爽又粗| 中文乱码字字幕精品一区二区三区| 香蕉丝袜av| 超碰成人久久| 亚洲av日韩在线播放| 国产免费一区二区三区四区乱码| videosex国产| 精品少妇久久久久久888优播| 久久久国产精品麻豆| 色播在线永久视频| 女人久久www免费人成看片| 精品少妇久久久久久888优播| 久久鲁丝午夜福利片| 中文字幕高清在线视频| 亚洲成色77777| 咕卡用的链子| 大话2 男鬼变身卡| 一本大道久久a久久精品| 国产一区二区三区综合在线观看| 日本av手机在线免费观看| 国产又色又爽无遮挡免| 亚洲激情五月婷婷啪啪| 日本av免费视频播放| 亚洲国产精品成人久久小说| 交换朋友夫妻互换小说| 午夜激情av网站| 国产色视频综合| 欧美成人精品欧美一级黄| 日本色播在线视频| 飞空精品影院首页| 久久精品久久久久久噜噜老黄| 三上悠亚av全集在线观看| 中文欧美无线码| 久久精品久久久久久久性| 中文欧美无线码| 亚洲图色成人| 成年女人毛片免费观看观看9 | 另类精品久久| 又大又黄又爽视频免费| 男女午夜视频在线观看| 欧美黄色淫秽网站| 少妇精品久久久久久久| 成人国语在线视频| 欧美av亚洲av综合av国产av| 又粗又硬又长又爽又黄的视频| tube8黄色片| 国产一区亚洲一区在线观看| 99久久人妻综合| av片东京热男人的天堂| 国产欧美日韩一区二区三区在线| 亚洲美女黄色视频免费看| 一级黄色大片毛片| 丝袜人妻中文字幕| 国产日韩欧美视频二区| 国产精品一国产av| 女人被躁到高潮嗷嗷叫费观| 日韩一卡2卡3卡4卡2021年| 国产一区有黄有色的免费视频| 国产精品亚洲av一区麻豆| 免费av中文字幕在线| 亚洲五月色婷婷综合| 国产精品成人在线| 欧美成人精品欧美一级黄| 欧美日韩一级在线毛片| 97在线人人人人妻| 一级毛片 在线播放| 国产亚洲精品第一综合不卡| 久久久精品94久久精品| 日本黄色日本黄色录像| 婷婷色综合www| 美女大奶头黄色视频| 中文字幕制服av| 高清av免费在线| 99国产精品免费福利视频| 午夜91福利影院| 母亲3免费完整高清在线观看| 亚洲欧美激情在线| 午夜老司机福利片| 久久狼人影院| 亚洲精品久久午夜乱码| 性色av乱码一区二区三区2| 秋霞在线观看毛片| 精品久久久久久电影网| 赤兔流量卡办理| 日日爽夜夜爽网站| 欧美xxⅹ黑人| 建设人人有责人人尽责人人享有的| 国产色视频综合| 亚洲国产精品国产精品| 亚洲精品国产一区二区精华液| 久久av网站| 丝袜美足系列| 赤兔流量卡办理| 国产成人精品无人区|