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

    層狀全空間介質(zhì)中重疊回線的瞬變電磁響應(yīng)

    2023-12-25 12:35:52劉建鵬車(chē)海清
    科技和產(chǎn)業(yè) 2023年22期
    關(guān)鍵詞:全空間回線層狀

    劉建鵬, 車(chē)海清, 朱 華, 劉 飛, 孫 巖

    (核工業(yè)二九〇研究所, 廣東 韶關(guān) 512029)

    地面瞬變電磁法是發(fā)射和接收裝置均位于地面的時(shí)間域方法,接收回線中得到隨時(shí)間衰減的感應(yīng)電動(dòng)勢(shì),能有效地反映該點(diǎn)豎直深度方向的電性特征[1-3]。但隨著復(fù)雜環(huán)境的地質(zhì)勘探任務(wù)需求增加,特別是在工程隧道、煤炭巷道超前探測(cè)、全空間地質(zhì)災(zāi)害等方面,全空間瞬變電磁法具有勘探效率高、對(duì)低阻體靈敏和無(wú)須接地等優(yōu)點(diǎn),全空間瞬變電磁法已經(jīng)成為首選的地球物理方法,并逐漸成為礦井勘探領(lǐng)域研究的熱點(diǎn)[4-9]。尤其受工作條件等因素限制,重疊回線裝置應(yīng)用較為廣泛。

    伴隨著煤炭等非金屬礦的開(kāi)采和交通運(yùn)輸?shù)拈_(kāi)發(fā),國(guó)內(nèi)已經(jīng)取得了許多的研究成果。部分研究通過(guò)物理實(shí)驗(yàn)的方法,探究出重疊回線的分辨率和信號(hào)強(qiáng)度主要跟發(fā)射回線參數(shù)有關(guān),并利用等磁矩的方式將多匝回線等效為同等面積的大半徑回線,得出小回線裝置同樣具有較大的探測(cè)深度[10-12]。理論研究則根據(jù)電磁場(chǎng)原理,結(jié)合數(shù)理方面已經(jīng)能模擬多維瞬態(tài)擴(kuò)散場(chǎng),而對(duì)于一維瞬態(tài)擴(kuò)散場(chǎng)研究較緩,與傳統(tǒng)磁偶源理論不同,可通過(guò)研究層狀海洋模型的回線源解析解,來(lái)探究全空間與半空間瞬變電磁法的聯(lián)系與區(qū)別[13-22]。

    從電磁場(chǎng)擴(kuò)散原理出發(fā),研究全空間中回線源產(chǎn)生的瞬變電磁場(chǎng),分析了直接求解的響應(yīng)和柱坐標(biāo)下求解的響應(yīng)表達(dá)式,并利用貝塞爾函數(shù)的積分式將這兩種解進(jìn)行了統(tǒng)一,推導(dǎo)了層狀全空間重疊回線的響應(yīng)和計(jì)算方法,從而探究發(fā)射半徑對(duì)全空間瞬變電磁場(chǎng)的影響規(guī)律。將該計(jì)算方法應(yīng)用于實(shí)測(cè)礦井?dāng)?shù)據(jù)中,經(jīng)擬二維斷面分析,計(jì)算結(jié)果證實(shí)了方法的有效性。

    1 均勻全空間的求解

    1.1 直角坐標(biāo)系下的求解

    如圖1所示,在電導(dǎo)率為σ、磁導(dǎo)率為μ0的均勻全空間中,存在一個(gè)發(fā)射電流為I、半徑為r的重疊回線,在t時(shí)刻的感應(yīng)電壓為[3]

    圖1 全空間中的回線源

    (1)

    式中:

    (2)

    式中:R為觀測(cè)點(diǎn)到原點(diǎn)的距離,R=2rsin(φ/2)。則式(2)為

    (3)

    式(1)變?yōu)?/p>

    (4)

    考慮第一類(lèi)修正的貝塞爾的積分表達(dá)式[18],令v=1有:

    (5)

    結(jié)合式(4),則可以得到感應(yīng)電壓表達(dá)式為

    (6)

    1.2 柱坐標(biāo)系下的求解

    如圖1所示,在柱坐標(biāo)系(r′,φ,z)下,回線源的中心處為坐標(biāo)原點(diǎn)(0,0,0),參數(shù)與圖1一致,得到的水平方向電場(chǎng)為[14]

    (7)

    利用法拉第電磁感應(yīng)定律,得

    V(t)=∮Eφ

    (8)

    將(7)式代入(8)式中,則有

    (9)

    對(duì)比式(6)和式(9)可知,兩種方法得出的感應(yīng)電壓完全相同,其中u1=θr,兩者可以相互驗(yàn)證。

    2 全空間層狀介質(zhì)的響應(yīng)計(jì)算方法

    如圖2所示,建立全空間層狀模型,Ti、Tj分別為第i層、j層的厚度(i=1,2,…,N-1,N;j=1,2,…,M-1,M);σi、σj分別為第i層、j層的電導(dǎo)率;半徑為r的水平重疊回線源位于薄層中(層厚非常小,可忽略影響),拉普拉斯域下瞬變電磁感應(yīng)電壓為[13]。

    圖2 位于層狀介質(zhì)中的重疊回線

    (10)

    式中:s為拉普拉斯變量,s與F(λ)有關(guān);J1為第一類(lèi)一階貝塞爾函數(shù);磁矩m=Iπr2。其他參數(shù)與圖1一致,F(λ)的推導(dǎo)公式為

    (11)

    式中:P、Q為每個(gè)層位之間的遞推值;Ki=λ2+sμ0σi,Kj=λ2+sμ0σj;Ki、Kj為第i層、j層的波數(shù)。

    遞推的關(guān)鍵步驟是得到末端層的P、Q值,即當(dāng)i=N、j=M時(shí)有Qi=N=μ0/Ki=N,Pj=M=μ0/Kj=M。即波數(shù)為Ki=N=λ2+sμ0σi=N,Kj=M=λ2+sμ0σj=M。

    由上述可知,求解全空間重疊回線瞬變電磁響應(yīng)的算法步驟為:從上下末端層(第N層,第M層)開(kāi)始,逐步往中間遞推,得到相應(yīng)第一層的P1、Q1,代入式(10),最終經(jīng)拉普拉斯逆變換得到時(shí)間域感應(yīng)電壓為

    (12)

    式中:L-1為拉普拉斯逆變換。與常規(guī)的地面瞬態(tài)響應(yīng)計(jì)算方法相同,式(12)可聯(lián)合漢克爾變換和正余弦變換得到時(shí)間域的瞬態(tài)響應(yīng)[19-20]。

    3 瞬變電磁響應(yīng)規(guī)律

    3.1 均勻全空間瞬變電磁場(chǎng)擴(kuò)散規(guī)律

    根據(jù)全空間瞬變場(chǎng)的擴(kuò)散規(guī)律,回線源關(guān)斷后產(chǎn)生的瞬變電磁場(chǎng)為水平方向的電場(chǎng)[9](后稱(chēng)為電場(chǎng)),研究電場(chǎng)在空間上的變化規(guī)律能直觀地了解場(chǎng)的擴(kuò)散特征。由式(7)可知,考慮一個(gè)垂直于φ平面且包含回線軸線的平面,該平面有利于直觀展示瞬變電磁場(chǎng)隨延時(shí)時(shí)間的擴(kuò)散規(guī)律。另一方面,為研究電場(chǎng)在r方向、z方向的傳播規(guī)律,可以給定一個(gè)半徑r(m)或者深度z(m),得到電場(chǎng)在深度或者半徑上的變化特征。

    如圖3所示,水平回線源在0.3 ms、1 ms和3 ms時(shí)的電場(chǎng)分布,其半徑為40 m,全空間電導(dǎo)率為1/100 S/m,回線中的電流為3 A。在均勻全空間介質(zhì)中,二次場(chǎng)的擴(kuò)散規(guī)律依然滿足“煙圈效應(yīng)”,并且“煙圈”會(huì)在x軸、z軸方向上運(yùn)動(dòng),而電場(chǎng)極值中心位置保持在r軸上。另外,電場(chǎng)的擴(kuò)散形態(tài)并非規(guī)則的圓形,說(shuō)明在瞬態(tài)場(chǎng)在空間中傳播存在一定的方向性,不同方向的速度不同。在不均勻介質(zhì)中,受+z、-z、+r方向的共同作用,似橢圓的電場(chǎng)將變得不光滑、不對(duì)稱(chēng),受地下結(jié)構(gòu)復(fù)雜多變的影響,滿足擴(kuò)散方程的電場(chǎng)形成更復(fù)雜的擴(kuò)散,難以判斷引起電場(chǎng)變化的異常來(lái)自+z還是-z方向。

    圖3 全空間介質(zhì)不同時(shí)刻的電場(chǎng)分布

    3.2 全空間響應(yīng)驗(yàn)證

    回線半徑是重疊回線的重要參數(shù),決定了發(fā)射、接收回線的幾何尺寸和磁矩??紤]半徑對(duì)應(yīng)電壓的影響,建立半徑分別為2、4、10、20、40、80、200 m的均勻全空間模型,同時(shí)為檢驗(yàn)全空間層狀解的精確性,利用相同參數(shù)的層狀解來(lái)計(jì)算均勻介質(zhì),并呈現(xiàn)半徑為4 m時(shí)的相對(duì)誤差曲線。由圖4可知,對(duì)于歸一化感應(yīng)電壓而言,均勻全空間響應(yīng)曲線形態(tài)基本相同;結(jié)合式(6)和式(9)對(duì)延時(shí)時(shí)間t的兩個(gè)極限狀態(tài)(早期t→0,晚期t→∞),在早期的響應(yīng)遵循t-1衰減,經(jīng)過(guò)渡期,晚期響應(yīng)遵循t-5/2衰減;半徑對(duì)早期的幅值影響大,將Kaufman[3]中的早期響應(yīng)進(jìn)行歸一化,可得到響應(yīng)與半徑成反比關(guān)系,例如在半徑為80和200 m時(shí),前者的早期響應(yīng)大于后者,到晚期則相反。另一方面,如圖4和表1所示,層狀解的響應(yīng)相對(duì)誤差較大值分布在延時(shí)時(shí)間始末時(shí)段,且相對(duì)誤差都在4%以?xún)?nèi),說(shuō)明算法有效且較準(zhǔn)確。

    表1 不同半徑的平均相對(duì)誤差

    圖4 不同半徑尺寸的解析解和層狀解

    4 瞬變電磁響應(yīng)特征

    4.1 典型全空間響應(yīng)特征

    上文已經(jīng)討論尺寸對(duì)均勻全空間的影響,得出計(jì)算的層狀解具有較高的精確度,為研究尺寸與層狀介質(zhì)響應(yīng)的影響規(guī)律提供了基礎(chǔ)。建立全空間的四層模型,發(fā)射回線的上半空間為均勻介質(zhì),下半空間為類(lèi)似地面的H型,通過(guò)Swidinsky[21]分別得到全空間瞬態(tài)響應(yīng)以及晚期視電阻率,如圖5所示。

    圖5 四層全空間不同半徑尺寸的響應(yīng)和視電阻率曲線

    由圖5(a)可知,對(duì)比均勻全空間,回線半徑小于層厚時(shí),層狀解響應(yīng)脫離均勻解析解的時(shí)間基本一致;反之,則隨半徑的增加而延后,但其尾支形態(tài)的“回歸時(shí)間”是一致的。“早期”響應(yīng)的強(qiáng)度不隨回線源半徑的改變而改變,僅延時(shí)時(shí)間不同,即回線半徑越大其低阻引起的響應(yīng)時(shí)間越延后。另一方面,延時(shí)時(shí)間在0.1~100 ms時(shí),通過(guò)增加回線源的半徑能顯著地增大響應(yīng)值,進(jìn)而提高響應(yīng)的信噪比,這段時(shí)間正是礦井瞬變電磁數(shù)據(jù)采集的主要測(cè)量時(shí)間段,例如延時(shí)時(shí)間為1 ms時(shí),半徑從2 m增加到200 m,層狀解感應(yīng)電壓從5.02×10-9V增加到0.30 V。由圖5b可知,均勻全空間的晚期視電阻率應(yīng)用范圍隨半徑的增加而延后;當(dāng)發(fā)射半徑大于等于目標(biāo)層厚度時(shí),晚期視電阻率上移,不利于晚期各個(gè)公式的應(yīng)用。

    考慮一個(gè)較為典型的七層模型,下半空間地電參數(shù)保持不變,上半空間與下半空間關(guān)于發(fā)射層對(duì)稱(chēng),即σi=1=σj=1=1/300 S/m,σi=2=σj=2=1/10 S/m,σN=σM=1/300 S/m,僅中間夾有空氣層,其厚度為4 m,其他參數(shù)與圖5一致。同樣得到全空間瞬態(tài)響應(yīng)以及對(duì)應(yīng)的晚期視電阻率曲線,如圖6所示。

    圖6 七層全空間不同半徑尺寸的(a)響應(yīng)和(b)視電阻率曲線

    由圖6(a)可知,上半空間的變化,對(duì)于小尺寸的早期影響比較大,說(shuō)明空氣層作為一種真實(shí)存在的影響因素,不能直接忽略,這與實(shí)際應(yīng)用中的結(jié)果一致。另一方面,受電磁波傳播的對(duì)稱(chēng)性,上半空間低阻層的表征基本與下半空間重合,其響應(yīng)表現(xiàn)在同一時(shí)段且疊加后的強(qiáng)度更大,但是不改變層狀解脫離均勻解析解的延時(shí)時(shí)間,僅對(duì)表征響應(yīng)的結(jié)果分析無(wú)法判斷低阻層來(lái)自上半空間或者下半空間,與前文結(jié)果相同,這正是全空間效應(yīng)的表現(xiàn)形式之一。由圖6(b)可知,半徑越小,巷道的影響越明顯,但對(duì)于薄層的巷道不影響目標(biāo)層的表征。

    4.2 不同尺寸多層全空間響應(yīng)特征

    在考慮裝置半徑對(duì)響應(yīng)的影響時(shí),可以換一個(gè)角度:在改變層位中覆蓋層、低阻層層厚的條件下,探討相同半徑在不同層厚尺度時(shí)對(duì)響應(yīng)的影響,即半徑為40 m,厚度分別為2、4、10、20、40、80、200 m及∞。圖7、圖8分別為不同層厚下的全空間響應(yīng)和視電阻率結(jié)果。綜合這些響應(yīng)結(jié)果可以得出:隨著低阻層層厚小于、等于到大于半徑改變時(shí),全空間瞬變電磁響應(yīng)對(duì)異常持續(xù)時(shí)間越長(zhǎng),偏離解析解的程度越明顯;隨著異常體離發(fā)射位置越來(lái)越遠(yuǎn),視電阻率也越 “偏離”真實(shí)值,且值偏大。

    圖7 四層全空間不同低阻層層厚尺寸的(a)響應(yīng)和(b)視電阻率曲線

    圖8 四層全空間不同覆蓋層層厚尺寸的(a)響應(yīng)和(b)視電阻率曲線

    5 應(yīng)用實(shí)例

    湖南省內(nèi)煤礦眾多,為了解某煤礦開(kāi)采中掘進(jìn)面上頂板的礦體分布情況,利用重疊回線發(fā)射源的瞬變電磁法開(kāi)展超前探測(cè)。圖9(a)為實(shí)測(cè)數(shù)據(jù)的視電阻率斷面圖,視深度結(jié)果采用于景邨[23]提出的礦井瞬變電磁法時(shí)間-深度轉(zhuǎn)換公式;視電阻率采用楊海燕等[24]提出的全區(qū)視電阻率計(jì)算公式。由圖可知,該區(qū)域橫向變化明顯,在x軸上的35 m和50~80 m位置上存在相對(duì)的低阻體,且數(shù)值較煤礦而言相近,該區(qū)域存在煤的可能性大。在此基礎(chǔ)上,建立了全空間模型,如圖9(b)所示,在x軸上布置點(diǎn)距為4.4 m的測(cè)線,長(zhǎng)度為110 m,測(cè)點(diǎn)的z方向?yàn)槿臻g層狀結(jié)構(gòu),其他參數(shù)與采集裝置相同。圖9(c)為理論數(shù)據(jù)的視電阻率斷面圖。由圖9可知,兩個(gè)斷面圖結(jié)構(gòu)上較為吻合,理論數(shù)據(jù)的視電阻率斷面明顯反映兩個(gè)低阻異常體的位置和形態(tài)。斷面結(jié)果說(shuō)明采用研究介紹的方法得出的全空間響應(yīng)能反映瞬態(tài)場(chǎng)在地下的擴(kuò)散規(guī)律。

    圖9 實(shí)測(cè)數(shù)據(jù)斷面、模型圖斷面和理論數(shù)據(jù)斷面

    6 結(jié)論

    位于全空間中的回線源,通過(guò)直接計(jì)算和采用柱坐標(biāo)系計(jì)算的方法能得出相同的瞬變電磁響應(yīng);并以此為基礎(chǔ),得到垂直平面上的水平方向電場(chǎng),該電場(chǎng)呈橢圓狀,證實(shí)回線源裝置探測(cè)具有一定的方向性,同時(shí)這也構(gòu)成了全空間瞬變電磁場(chǎng)難以分辨異常體位置的復(fù)雜性。

    在不同半徑條件下的瞬變電磁響應(yīng),相對(duì)“早期”響應(yīng),全空間響應(yīng)與半徑成反比,且按延時(shí)時(shí)間t-1衰減。

    采用層狀解公式計(jì)算的均勻全空間響應(yīng),具有相對(duì)誤差較小和精度較高等特點(diǎn),能夠用來(lái)計(jì)算層狀空間的瞬變電磁響應(yīng);在“早期”響應(yīng)中,強(qiáng)度不隨回線半徑的變化而改變,只是延時(shí)時(shí)間不同,即回線半徑越大響應(yīng)時(shí)間越延后。在“晚期”響應(yīng)中,通過(guò)增加回線源的半徑能顯著地增大響應(yīng)值,可以提高地下瞬變電磁法的抗干擾性能。

    對(duì)于全空間的晚期視電阻率應(yīng)用范圍,當(dāng)發(fā)射半徑大于等于目標(biāo)層厚度時(shí)其視電阻率值上移,不利于晚期公式的應(yīng)用。半徑越小,巷道對(duì)早期響應(yīng)的影響越明顯,但其對(duì)目標(biāo)層的表征影響較小。

    猜你喜歡
    全空間回線層狀
    無(wú)接地極直流輸電線路金屬回線選型設(shè)計(jì)
    吉林電力(2022年1期)2022-11-10 09:20:24
    軋制復(fù)合制備TA1/AZ31B/TA1層狀復(fù)合材料組織與性能研究
    全空間瞬變電磁場(chǎng)低阻層屏蔽效應(yīng)數(shù)值模擬研究
    多層圓線圈在電磁計(jì)算中全空間磁感應(yīng)強(qiáng)度B的分布
    兩級(jí)結(jié)構(gòu)層狀Ti-TiBw/Ti復(fù)合材料擴(kuò)散焊及其拉伸行為
    焊接(2016年9期)2016-02-27 13:05:22
    ±800?kV特高壓直流金屬回線斷路器保護(hù)誤動(dòng)分析
    8字形載流方形回線的空間磁場(chǎng)分布
    高韌性抗層狀撕裂Q345FTE-Z35鋼板開(kāi)發(fā)
    新疆鋼鐵(2015年2期)2015-11-07 03:27:52
    全空間上一類(lèi)半線性雙調(diào)和方程正解的衰減
    層狀雙氫氧化物處理工業(yè)廢水的研究進(jìn)展
    成人鲁丝片一二三区免费| 高清日韩中文字幕在线| 日韩,欧美,国产一区二区三区 | 哪里可以看免费的av片| 国产一区二区在线观看日韩| 亚洲第一区二区三区不卡| 色综合色国产| 小蜜桃在线观看免费完整版高清| 亚洲人成网站高清观看| 国产精品久久电影中文字幕| a级毛色黄片| 久久久久网色| 国产成人影院久久av| 亚洲av第一区精品v没综合| 哪里可以看免费的av片| 一区福利在线观看| 日本-黄色视频高清免费观看| 三级经典国产精品| 国内精品一区二区在线观看| 久久久久久久久久成人| 国产伦精品一区二区三区四那| 老女人水多毛片| 中文字幕精品亚洲无线码一区| 亚洲av免费高清在线观看| 丝袜美腿在线中文| 亚洲av一区综合| 免费看美女性在线毛片视频| 欧美日韩乱码在线| 有码 亚洲区| 日韩欧美三级三区| 国产一区亚洲一区在线观看| eeuss影院久久| 国产真实伦视频高清在线观看| 亚洲精品久久久久久婷婷小说 | 看片在线看免费视频| 日韩国内少妇激情av| 精品不卡国产一区二区三区| 三级男女做爰猛烈吃奶摸视频| 欧美激情久久久久久爽电影| 九草在线视频观看| 亚洲最大成人中文| 国产 一区 欧美 日韩| 亚洲国产精品久久男人天堂| 99久久成人亚洲精品观看| 在线播放无遮挡| kizo精华| 亚洲成av人片在线播放无| 国产精品人妻久久久久久| 亚洲在线自拍视频| 大又大粗又爽又黄少妇毛片口| 美女被艹到高潮喷水动态| 亚洲精华国产精华液的使用体验 | 亚洲熟妇中文字幕五十中出| 亚洲精品粉嫩美女一区| 97在线视频观看| 欧美性感艳星| 一区二区三区四区激情视频 | 久久久久免费精品人妻一区二区| 一个人看视频在线观看www免费| 久久久久久久久久久免费av| 熟女电影av网| 亚洲欧美精品自产自拍| 国产乱人视频| 日本色播在线视频| 亚洲熟妇中文字幕五十中出| 丰满乱子伦码专区| 日本av手机在线免费观看| 一个人观看的视频www高清免费观看| 一本久久精品| 麻豆国产av国片精品| 欧美三级亚洲精品| 18+在线观看网站| 亚洲成av人片在线播放无| 看免费成人av毛片| 免费观看在线日韩| 亚洲欧美日韩无卡精品| 热99在线观看视频| 久久这里只有精品中国| 爱豆传媒免费全集在线观看| 有码 亚洲区| 国产黄色小视频在线观看| 国模一区二区三区四区视频| 一本一本综合久久| 成人鲁丝片一二三区免费| 亚洲美女搞黄在线观看| 成人综合一区亚洲| 国产在线男女| 人妻制服诱惑在线中文字幕| 一级黄色大片毛片| 国产熟女欧美一区二区| 国内少妇人妻偷人精品xxx网站| 在线观看av片永久免费下载| 中文欧美无线码| 久久精品夜夜夜夜夜久久蜜豆| 精品国产三级普通话版| 99热只有精品国产| 淫秽高清视频在线观看| 狠狠狠狠99中文字幕| 一本久久中文字幕| 美女 人体艺术 gogo| 国产美女午夜福利| 亚洲一区高清亚洲精品| 高清毛片免费看| 久久久久久久久久久丰满| 亚洲天堂国产精品一区在线| 啦啦啦观看免费观看视频高清| 久久久久久久久中文| 看免费成人av毛片| 免费观看在线日韩| av天堂在线播放| 日韩中字成人| 国产精品免费一区二区三区在线| 一级毛片电影观看 | 精品熟女少妇av免费看| 神马国产精品三级电影在线观看| 97人妻精品一区二区三区麻豆| 国产精品久久久久久精品电影| 床上黄色一级片| 观看美女的网站| 亚洲高清免费不卡视频| av女优亚洲男人天堂| 成年免费大片在线观看| 免费黄网站久久成人精品| 级片在线观看| 小蜜桃在线观看免费完整版高清| 在线天堂最新版资源| 搡女人真爽免费视频火全软件| 国产一区亚洲一区在线观看| 狂野欧美激情性xxxx在线观看| 99九九线精品视频在线观看视频| 夜夜爽天天搞| 一区二区三区免费毛片| 国产美女午夜福利| 99热这里只有精品一区| 成人美女网站在线观看视频| 日日摸夜夜添夜夜添av毛片| 国产精品一区二区性色av| 91av网一区二区| 国产精品免费一区二区三区在线| 激情 狠狠 欧美| 中国美女看黄片| 午夜激情欧美在线| 国产av一区在线观看免费| 国产激情偷乱视频一区二区| 日韩国内少妇激情av| 欧美三级亚洲精品| 国产免费一级a男人的天堂| 大香蕉久久网| 全区人妻精品视频| 人妻夜夜爽99麻豆av| 熟妇人妻久久中文字幕3abv| 国产又黄又爽又无遮挡在线| 亚洲国产精品成人久久小说 | 变态另类成人亚洲欧美熟女| 一区二区三区高清视频在线| 3wmmmm亚洲av在线观看| 九色成人免费人妻av| 日日摸夜夜添夜夜爱| 亚洲成人久久爱视频| 亚洲人成网站在线播| 久久久a久久爽久久v久久| 免费一级毛片在线播放高清视频| 久久久欧美国产精品| 九草在线视频观看| 寂寞人妻少妇视频99o| 2022亚洲国产成人精品| 26uuu在线亚洲综合色| 国产精品人妻久久久影院| 激情 狠狠 欧美| 插阴视频在线观看视频| 国产黄色视频一区二区在线观看 | 久久午夜亚洲精品久久| 亚洲不卡免费看| 亚洲,欧美,日韩| av.在线天堂| 精华霜和精华液先用哪个| 看十八女毛片水多多多| 亚洲色图av天堂| 午夜福利高清视频| av免费观看日本| 亚洲内射少妇av| 日本爱情动作片www.在线观看| 久久九九热精品免费| 桃色一区二区三区在线观看| 亚洲第一电影网av| 天堂av国产一区二区熟女人妻| 婷婷亚洲欧美| 国产女主播在线喷水免费视频网站 | 国产伦一二天堂av在线观看| 一个人看的www免费观看视频| 日本五十路高清| 亚州av有码| 久久6这里有精品| 国产人妻一区二区三区在| 亚洲国产精品久久男人天堂| 联通29元200g的流量卡| 日韩,欧美,国产一区二区三区 | 久久久欧美国产精品| 人体艺术视频欧美日本| 少妇人妻一区二区三区视频| 波多野结衣高清作品| 国产高清有码在线观看视频| 国产黄色视频一区二区在线观看 | 免费黄网站久久成人精品| 美女被艹到高潮喷水动态| 国产精品1区2区在线观看.| 日韩欧美精品v在线| 国产成人福利小说| 18禁在线无遮挡免费观看视频| 美女高潮的动态| 波多野结衣巨乳人妻| 1000部很黄的大片| 日韩欧美一区二区三区在线观看| 精品免费久久久久久久清纯| 亚洲欧美中文字幕日韩二区| av在线老鸭窝| 三级男女做爰猛烈吃奶摸视频| 国产伦精品一区二区三区视频9| 国产黄色小视频在线观看| 久久精品国产亚洲网站| 成人综合一区亚洲| 国产精品久久久久久久电影| 深夜a级毛片| 亚洲精品粉嫩美女一区| 日韩欧美精品免费久久| 高清毛片免费看| 亚洲欧美成人精品一区二区| 在线观看av片永久免费下载| АⅤ资源中文在线天堂| 男人的好看免费观看在线视频| 99久久人妻综合| 国产精品.久久久| 亚洲乱码一区二区免费版| a级毛片免费高清观看在线播放| 午夜老司机福利剧场| 99热这里只有精品一区| 特级一级黄色大片| 午夜福利高清视频| 亚洲国产欧洲综合997久久,| 日日干狠狠操夜夜爽| 久久久国产成人精品二区| 日韩欧美在线乱码| 成人毛片60女人毛片免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜a级毛片| 午夜亚洲福利在线播放| 熟女人妻精品中文字幕| 国产成人a区在线观看| 亚洲人成网站在线播放欧美日韩| 久久婷婷人人爽人人干人人爱| 国产成人a∨麻豆精品| 欧美3d第一页| 日日摸夜夜添夜夜添av毛片| 久久欧美精品欧美久久欧美| 精品久久久久久久久av| 国产中年淑女户外野战色| 啦啦啦观看免费观看视频高清| 日本-黄色视频高清免费观看| 国产私拍福利视频在线观看| 好男人视频免费观看在线| 又黄又爽又刺激的免费视频.| 婷婷精品国产亚洲av| www.色视频.com| 国产精品久久视频播放| av黄色大香蕉| 婷婷亚洲欧美| 最新中文字幕久久久久| 免费搜索国产男女视频| 一卡2卡三卡四卡精品乱码亚洲| 内射极品少妇av片p| 日本与韩国留学比较| 国产一区二区三区在线臀色熟女| 亚洲国产色片| 日韩欧美在线乱码| 在线观看av片永久免费下载| 九色成人免费人妻av| av女优亚洲男人天堂| 99热精品在线国产| 亚洲成人久久性| 国产成人精品婷婷| 插阴视频在线观看视频| 又爽又黄a免费视频| 美女脱内裤让男人舔精品视频 | 偷拍熟女少妇极品色| 免费看日本二区| 18+在线观看网站| 国产精品福利在线免费观看| 有码 亚洲区| 国产精品嫩草影院av在线观看| 国产成人a区在线观看| 身体一侧抽搐| 一边亲一边摸免费视频| 亚洲av熟女| 久久99热6这里只有精品| 久久久久免费精品人妻一区二区| 97人妻精品一区二区三区麻豆| 免费看美女性在线毛片视频| 秋霞在线观看毛片| 亚洲人成网站在线观看播放| 亚洲精品久久久久久婷婷小说 | 国产三级中文精品| 中国国产av一级| 国产精品蜜桃在线观看 | 菩萨蛮人人尽说江南好唐韦庄 | 床上黄色一级片| 国产黄色视频一区二区在线观看 | 少妇的逼水好多| 久久欧美精品欧美久久欧美| 国产精品精品国产色婷婷| 18禁在线无遮挡免费观看视频| 亚洲人成网站高清观看| 亚洲中文字幕日韩| 免费电影在线观看免费观看| 99在线视频只有这里精品首页| 熟妇人妻久久中文字幕3abv| 91aial.com中文字幕在线观看| 狠狠狠狠99中文字幕| 日本在线视频免费播放| 亚洲电影在线观看av| 床上黄色一级片| 久久人人爽人人片av| 免费av观看视频| 精品人妻视频免费看| 亚洲婷婷狠狠爱综合网| 啦啦啦韩国在线观看视频| 国产精品久久久久久亚洲av鲁大| 亚洲av一区综合| 国产91av在线免费观看| 黄色一级大片看看| 久久99热这里只有精品18| 插阴视频在线观看视频| 国产探花在线观看一区二区| 22中文网久久字幕| 久久综合国产亚洲精品| 午夜福利成人在线免费观看| 黄色日韩在线| 国产精品,欧美在线| 亚洲国产色片| 国产成人一区二区在线| 婷婷色综合大香蕉| 好男人视频免费观看在线| 99热精品在线国产| 男女啪啪激烈高潮av片| 亚洲七黄色美女视频| 日本三级黄在线观看| av天堂中文字幕网| 黑人高潮一二区| 桃色一区二区三区在线观看| 波野结衣二区三区在线| 色综合亚洲欧美另类图片| av女优亚洲男人天堂| 亚洲,欧美,日韩| 久久婷婷人人爽人人干人人爱| 国内精品宾馆在线| 亚洲成av人片在线播放无| 色噜噜av男人的天堂激情| 高清毛片免费观看视频网站| 日韩视频在线欧美| 自拍偷自拍亚洲精品老妇| 欧美日韩乱码在线| 中文字幕免费在线视频6| 在线观看免费视频日本深夜| 麻豆国产97在线/欧美| 中出人妻视频一区二区| 国产精品国产高清国产av| 国产av不卡久久| 干丝袜人妻中文字幕| 欧美区成人在线视频| 国产国拍精品亚洲av在线观看| 久久久国产成人精品二区| 欧美日韩综合久久久久久| 久久久精品94久久精品| 蜜桃亚洲精品一区二区三区| 小说图片视频综合网站| 九色成人免费人妻av| 一级av片app| 高清午夜精品一区二区三区 | 内地一区二区视频在线| 欧美日韩一区二区视频在线观看视频在线 | 中文字幕免费在线视频6| 国产老妇女一区| 97在线视频观看| 乱系列少妇在线播放| 国产伦理片在线播放av一区 | 69人妻影院| 大又大粗又爽又黄少妇毛片口| 亚洲欧美精品专区久久| 天天躁日日操中文字幕| 久久这里有精品视频免费| 亚洲自拍偷在线| 日韩高清综合在线| 中文字幕av在线有码专区| 哪个播放器可以免费观看大片| 久久6这里有精品| 国产单亲对白刺激| 欧美最新免费一区二区三区| 亚洲精品色激情综合| 久久精品久久久久久久性| 成人永久免费在线观看视频| 日韩欧美 国产精品| 日本av手机在线免费观看| 欧美3d第一页| 欧美+日韩+精品| 菩萨蛮人人尽说江南好唐韦庄 | 97在线视频观看| 中文字幕免费在线视频6| 12—13女人毛片做爰片一| 久久久精品欧美日韩精品| 国产精品永久免费网站| 天天躁夜夜躁狠狠久久av| 两个人视频免费观看高清| 精品国内亚洲2022精品成人| 哪个播放器可以免费观看大片| 亚洲最大成人手机在线| 国产午夜精品论理片| 99久久中文字幕三级久久日本| 久久精品国产亚洲网站| 免费av不卡在线播放| 国产一区二区亚洲精品在线观看| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区免费观看| 亚洲五月天丁香| 国产精品野战在线观看| 一区福利在线观看| 国产真实乱freesex| 网址你懂的国产日韩在线| 一级毛片电影观看 | 久久久久九九精品影院| 亚洲性久久影院| 久久精品国产清高在天天线| 人妻系列 视频| 国产亚洲精品久久久久久毛片| 欧美一区二区亚洲| 亚洲在线自拍视频| 深夜a级毛片| 午夜精品在线福利| 成人高潮视频无遮挡免费网站| 欧美一区二区国产精品久久精品| 亚洲第一电影网av| 村上凉子中文字幕在线| 久久久色成人| 少妇高潮的动态图| 噜噜噜噜噜久久久久久91| 国产精品麻豆人妻色哟哟久久 | 久久午夜亚洲精品久久| 午夜精品国产一区二区电影 | 欧美日韩综合久久久久久| 床上黄色一级片| 亚洲国产精品国产精品| 欧美+亚洲+日韩+国产| 午夜福利视频1000在线观看| 日本熟妇午夜| 九色成人免费人妻av| 国产av一区在线观看免费| 亚洲欧美日韩无卡精品| 亚洲人与动物交配视频| 天堂中文最新版在线下载 | 亚洲欧洲日产国产| 99久久人妻综合| 日产精品乱码卡一卡2卡三| 给我免费播放毛片高清在线观看| 欧美+日韩+精品| 日韩三级伦理在线观看| 国产v大片淫在线免费观看| 男的添女的下面高潮视频| 好男人视频免费观看在线| 搞女人的毛片| 国产一区二区激情短视频| 国产又黄又爽又无遮挡在线| 婷婷色综合大香蕉| 麻豆乱淫一区二区| 亚洲欧美精品自产自拍| 性插视频无遮挡在线免费观看| 91狼人影院| a级毛片免费高清观看在线播放| 99热这里只有是精品50| 久久精品国产清高在天天线| 内地一区二区视频在线| 2021天堂中文幕一二区在线观| 午夜福利视频1000在线观看| 人人妻人人澡人人爽人人夜夜 | 少妇丰满av| 日日撸夜夜添| 69av精品久久久久久| 日本黄色片子视频| 精品久久久噜噜| 国产黄片美女视频| 亚洲欧美日韩无卡精品| 美女cb高潮喷水在线观看| 啦啦啦观看免费观看视频高清| 亚洲成人精品中文字幕电影| 毛片女人毛片| 久久久久久久久久成人| 亚洲精品粉嫩美女一区| 亚洲人成网站高清观看| 男人舔女人下体高潮全视频| 免费看av在线观看网站| 亚洲精品亚洲一区二区| 久久国内精品自在自线图片| 九草在线视频观看| 在线观看美女被高潮喷水网站| 小蜜桃在线观看免费完整版高清| 国产高清视频在线观看网站| 蜜桃久久精品国产亚洲av| 亚洲成人精品中文字幕电影| 日本一二三区视频观看| 小蜜桃在线观看免费完整版高清| 99久久中文字幕三级久久日本| 亚洲人成网站在线播| 亚洲激情五月婷婷啪啪| 亚洲高清免费不卡视频| 免费看光身美女| 免费看美女性在线毛片视频| 麻豆一二三区av精品| 长腿黑丝高跟| 日韩视频在线欧美| 国产一区亚洲一区在线观看| 亚洲图色成人| 久久婷婷人人爽人人干人人爱| 内地一区二区视频在线| 久久精品国产清高在天天线| 99久久人妻综合| 亚洲欧美成人综合另类久久久 | 99热这里只有是精品在线观看| 大型黄色视频在线免费观看| 精品欧美国产一区二区三| 成人特级黄色片久久久久久久| 18禁在线无遮挡免费观看视频| 欧美不卡视频在线免费观看| 一个人免费在线观看电影| 日日摸夜夜添夜夜添av毛片| 综合色av麻豆| 亚洲三级黄色毛片| 又粗又硬又长又爽又黄的视频 | 国产熟女欧美一区二区| 亚洲国产精品成人综合色| 中国美白少妇内射xxxbb| 狂野欧美激情性xxxx在线观看| 国产av麻豆久久久久久久| 国产男人的电影天堂91| 久久久精品大字幕| 国产高清视频在线观看网站| 国产高清三级在线| 国产精品免费一区二区三区在线| 直男gayav资源| 亚洲综合色惰| 午夜久久久久精精品| 久久99热这里只有精品18| 色尼玛亚洲综合影院| 男女做爰动态图高潮gif福利片| 亚洲欧美精品综合久久99| 久久国产乱子免费精品| 国产一级毛片在线| 精华霜和精华液先用哪个| 人妻久久中文字幕网| eeuss影院久久| 蜜桃久久精品国产亚洲av| 国内揄拍国产精品人妻在线| 爱豆传媒免费全集在线观看| 国产 一区 欧美 日韩| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国模一区二区三区四区视频| 国产精品1区2区在线观看.| 好男人视频免费观看在线| 国产老妇女一区| 精品欧美国产一区二区三| 成人毛片60女人毛片免费| 国产老妇伦熟女老妇高清| 寂寞人妻少妇视频99o| 国产精品久久久久久精品电影小说 | 变态另类丝袜制服| 精品国内亚洲2022精品成人| 久久精品综合一区二区三区| 国产午夜精品一二区理论片| 欧美性猛交黑人性爽| 中文资源天堂在线| 亚洲av一区综合| 岛国毛片在线播放| 六月丁香七月| 日本-黄色视频高清免费观看| 国产亚洲精品久久久com| 人妻夜夜爽99麻豆av| 国产麻豆成人av免费视频| 亚洲精品久久久久久婷婷小说 | 亚洲精品影视一区二区三区av| 国产熟女欧美一区二区| 亚洲国产色片| 亚洲av中文av极速乱| 亚洲一区高清亚洲精品| 婷婷色综合大香蕉| 亚洲人成网站在线播| 亚洲综合色惰| 久久久精品大字幕| 97在线视频观看| 国产精品人妻久久久久久| 在线观看美女被高潮喷水网站| 九草在线视频观看| 最近手机中文字幕大全| 国产免费一级a男人的天堂| 国产老妇伦熟女老妇高清| 又粗又爽又猛毛片免费看| 亚洲一区高清亚洲精品| 久久这里只有精品中国| 午夜精品一区二区三区免费看| a级毛片免费高清观看在线播放| 欧美日韩国产亚洲二区| 成人午夜高清在线视频| 九九在线视频观看精品| 精品欧美国产一区二区三| 国产精品久久久久久亚洲av鲁大| 久久草成人影院| 99热网站在线观看| 国产色爽女视频免费观看| 少妇的逼好多水| 欧美+日韩+精品|