• <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)展
    91精品国产九色| 国产精品一区二区三区四区免费观看| 国产爱豆传媒在线观看| 日韩欧美三级三区| 亚洲高清免费不卡视频| 亚洲精品乱码久久久v下载方式| 麻豆成人av视频| 大型黄色视频在线免费观看| a级毛片a级免费在线| 亚洲精品久久国产高清桃花| 99久久人妻综合| 亚洲欧美日韩卡通动漫| 亚洲成a人片在线一区二区| 又黄又爽又刺激的免费视频.| 免费av观看视频| 边亲边吃奶的免费视频| 国产高清视频在线观看网站| 人妻少妇偷人精品九色| 小蜜桃在线观看免费完整版高清| 欧美区成人在线视频| 精华霜和精华液先用哪个| 日韩三级伦理在线观看| 97在线视频观看| 69人妻影院| 久久九九热精品免费| 老司机福利观看| av在线亚洲专区| 日韩av不卡免费在线播放| 亚州av有码| 午夜老司机福利剧场| 亚洲久久久久久中文字幕| 欧美人与善性xxx| 蜜臀久久99精品久久宅男| 国产精品日韩av在线免费观看| 久久午夜亚洲精品久久| 久久人妻av系列| 成年版毛片免费区| 久久精品久久久久久噜噜老黄 | 99riav亚洲国产免费| 99热这里只有精品一区| 国产精品.久久久| 五月伊人婷婷丁香| 六月丁香七月| 成人美女网站在线观看视频| 欧美一区二区精品小视频在线| 国产精品一区二区三区四区久久| 久久久久网色| 青春草视频在线免费观看| 插逼视频在线观看| 美女黄网站色视频| 麻豆国产av国片精品| 国产在线精品亚洲第一网站| 夜夜看夜夜爽夜夜摸| 99久久中文字幕三级久久日本| 精品久久久久久久末码| 老女人水多毛片| 夫妻性生交免费视频一级片| 国产一区亚洲一区在线观看| 国产精品久久久久久精品电影小说 | 人妻系列 视频| 亚洲av中文字字幕乱码综合| 在线观看av片永久免费下载| 丰满乱子伦码专区| 国产精品伦人一区二区| .国产精品久久| 九九久久精品国产亚洲av麻豆| 国产精品久久久久久av不卡| 尤物成人国产欧美一区二区三区| 精品一区二区三区人妻视频| 综合色av麻豆| 又爽又黄无遮挡网站| 欧美人与善性xxx| 亚洲成人中文字幕在线播放| 一区二区三区免费毛片| 级片在线观看| 男的添女的下面高潮视频| 国产av不卡久久| 老司机影院成人| 啦啦啦观看免费观看视频高清| 三级国产精品欧美在线观看| 国产黄a三级三级三级人| 久99久视频精品免费| 国产av麻豆久久久久久久| 国内少妇人妻偷人精品xxx网站| 国产一区亚洲一区在线观看| 热99re8久久精品国产| 欧美高清成人免费视频www| 久久欧美精品欧美久久欧美| 日韩av在线大香蕉| 国产成人a区在线观看| 三级国产精品欧美在线观看| 床上黄色一级片| 国产片特级美女逼逼视频| 欧美最黄视频在线播放免费| 国产精品无大码| 最近中文字幕高清免费大全6| 性欧美人与动物交配| 中文在线观看免费www的网站| 日韩欧美在线乱码| 免费无遮挡裸体视频| av在线蜜桃| 人人妻人人澡欧美一区二区| 3wmmmm亚洲av在线观看| 国产一区亚洲一区在线观看| 国产国拍精品亚洲av在线观看| 综合色av麻豆| 成人午夜高清在线视频| 午夜福利在线观看吧| 日本一二三区视频观看| 97热精品久久久久久| 国产精品嫩草影院av在线观看| 女人被狂操c到高潮| 亚洲欧美日韩东京热| 午夜a级毛片| 亚洲婷婷狠狠爱综合网| 男女那种视频在线观看| av天堂在线播放| 国产黄色小视频在线观看| 性欧美人与动物交配| 国产白丝娇喘喷水9色精品| 国产成人91sexporn| 久久国内精品自在自线图片| 国产在线男女| 国产一区二区在线观看日韩| 99riav亚洲国产免费| 大又大粗又爽又黄少妇毛片口| 天天躁夜夜躁狠狠久久av| 青春草亚洲视频在线观看| 久久久国产成人精品二区| 99久久九九国产精品国产免费| 丰满的人妻完整版| 久久午夜亚洲精品久久| 天天一区二区日本电影三级| 亚洲国产欧美在线一区| 国产v大片淫在线免费观看| h日本视频在线播放| 国产精品av视频在线免费观看| 嫩草影院新地址| 99久久无色码亚洲精品果冻| 十八禁国产超污无遮挡网站| 国产精品福利在线免费观看| 国产成人影院久久av| 99热精品在线国产| 少妇被粗大猛烈的视频| 国产精品久久久久久精品电影小说 | 成人美女网站在线观看视频| 日韩,欧美,国产一区二区三区 | 日韩欧美精品v在线| 精品久久国产蜜桃| 国产亚洲av嫩草精品影院| 国产91av在线免费观看| 日本熟妇午夜| 国产成人精品婷婷| 非洲黑人性xxxx精品又粗又长| 久久久久久久久久成人| 五月玫瑰六月丁香| 婷婷色av中文字幕| 国内精品美女久久久久久| 综合色av麻豆| 老师上课跳d突然被开到最大视频| 边亲边吃奶的免费视频| 色播亚洲综合网| 国内久久婷婷六月综合欲色啪| 成年版毛片免费区| 国内揄拍国产精品人妻在线| av在线天堂中文字幕| 热99re8久久精品国产| a级一级毛片免费在线观看| 天堂中文最新版在线下载 | av黄色大香蕉| 日韩精品青青久久久久久| 深爱激情五月婷婷| 亚洲在线观看片| 欧美又色又爽又黄视频| 最近手机中文字幕大全| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品av视频在线免费观看| 老司机影院成人| 亚洲,欧美,日韩| 久久99热这里只有精品18| 精品一区二区三区视频在线| 三级经典国产精品| 久久久精品欧美日韩精品| 亚洲不卡免费看| 亚洲成人久久性| 嫩草影院精品99| 成年av动漫网址| 成人亚洲欧美一区二区av| 久久国产乱子免费精品| 久久久久网色| 成年版毛片免费区| 久久久精品欧美日韩精品| 别揉我奶头 嗯啊视频| 久久久久久九九精品二区国产| 日本爱情动作片www.在线观看| 国产高清三级在线| 18禁黄网站禁片免费观看直播| 久久韩国三级中文字幕| 如何舔出高潮| 国内精品一区二区在线观看| 久久人人爽人人爽人人片va| 免费黄网站久久成人精品| 色哟哟哟哟哟哟| 69人妻影院| 免费观看精品视频网站| 成人国产麻豆网| 成人性生交大片免费视频hd| 插逼视频在线观看| 欧美在线一区亚洲| 亚洲国产日韩欧美精品在线观看| 久久精品国产亚洲av香蕉五月| 成人av在线播放网站| 国产精品一区二区三区四区免费观看| 青春草亚洲视频在线观看| 天堂影院成人在线观看| 18禁黄网站禁片免费观看直播| 精品欧美国产一区二区三| 欧美激情国产日韩精品一区| 国产精品日韩av在线免费观看| 九九爱精品视频在线观看| 午夜福利高清视频| 毛片女人毛片| 国产黄片美女视频| 日本免费a在线| 十八禁国产超污无遮挡网站| 一区福利在线观看| 国产单亲对白刺激| 亚洲av中文av极速乱| 欧美xxxx黑人xx丫x性爽| 日韩欧美在线乱码| 国产黄a三级三级三级人| 欧洲精品卡2卡3卡4卡5卡区| 蜜臀久久99精品久久宅男| 亚洲欧美清纯卡通| 国语自产精品视频在线第100页| 热99在线观看视频| 一个人观看的视频www高清免费观看| 成人性生交大片免费视频hd| 一级毛片aaaaaa免费看小| 在线观看av片永久免费下载| 亚洲欧美清纯卡通| 成年版毛片免费区| 日本av手机在线免费观看| 联通29元200g的流量卡| 国产私拍福利视频在线观看| 精品久久久久久成人av| 深夜a级毛片| 亚洲欧美中文字幕日韩二区| 免费观看在线日韩| 欧美xxxx黑人xx丫x性爽| 色5月婷婷丁香| 亚洲18禁久久av| 校园春色视频在线观看| 欧美日本亚洲视频在线播放| 日韩 亚洲 欧美在线| 夫妻性生交免费视频一级片| 99国产精品一区二区蜜桃av| kizo精华| 色哟哟·www| 亚洲图色成人| 在线观看一区二区三区| 婷婷色av中文字幕| 国产三级在线视频| 99热6这里只有精品| 插逼视频在线观看| 午夜福利在线观看吧| 国产精品99久久久久久久久| 青春草亚洲视频在线观看| 国产色婷婷99| 日韩三级伦理在线观看| 国产大屁股一区二区在线视频| 中文欧美无线码| 成人性生交大片免费视频hd| 身体一侧抽搐| 99久久精品国产国产毛片| 岛国毛片在线播放| 日本免费a在线| 亚洲av中文字字幕乱码综合| 国产黄a三级三级三级人| 日韩欧美国产在线观看| 国产伦精品一区二区三区视频9| 秋霞在线观看毛片| 日韩一本色道免费dvd| 午夜福利视频1000在线观看| 女同久久另类99精品国产91| 亚洲精品影视一区二区三区av| 亚洲无线观看免费| 中文字幕熟女人妻在线| 国产在线男女| 久久精品综合一区二区三区| 青春草视频在线免费观看| 久久久久久久久大av| 国产av不卡久久| 蜜桃久久精品国产亚洲av| av在线观看视频网站免费| 国产真实乱freesex| 欧美区成人在线视频| 日本av手机在线免费观看| 少妇的逼水好多| 日韩欧美三级三区| 九九爱精品视频在线观看| 深夜a级毛片| 非洲黑人性xxxx精品又粗又长| 国产精品一区二区性色av| 欧美成人一区二区免费高清观看| 看十八女毛片水多多多| 最近中文字幕高清免费大全6| 国产精品人妻久久久久久| 亚洲熟妇中文字幕五十中出| 女的被弄到高潮叫床怎么办| 蜜臀久久99精品久久宅男| 成人av在线播放网站| 成年女人看的毛片在线观看| 久久国内精品自在自线图片| 国产极品精品免费视频能看的| 校园人妻丝袜中文字幕| 中国美白少妇内射xxxbb| 国产一区二区激情短视频| 99热这里只有是精品50| 久久久久久久久久黄片| 精品人妻一区二区三区麻豆| 国产私拍福利视频在线观看| 国产亚洲av片在线观看秒播厂 | 日韩国内少妇激情av| 乱人视频在线观看| 久久久久久久久中文| 久久久久久伊人网av| 高清毛片免费看| 黄片wwwwww| 免费看日本二区| 自拍偷自拍亚洲精品老妇| 热99re8久久精品国产| 边亲边吃奶的免费视频| 国产成人aa在线观看| 12—13女人毛片做爰片一| 给我免费播放毛片高清在线观看| 国产亚洲精品av在线| 欧美性猛交黑人性爽| 免费人成在线观看视频色| 波多野结衣高清作品| 久久人妻av系列| 高清日韩中文字幕在线| 中文字幕人妻熟人妻熟丝袜美| 热99在线观看视频| 精品久久久久久久久久久久久| 日韩精品有码人妻一区| 国产久久久一区二区三区| 久久精品国产自在天天线| 国产一级毛片在线| 男人舔女人下体高潮全视频| 别揉我奶头 嗯啊视频| 久久精品夜夜夜夜夜久久蜜豆| 欧美成人精品欧美一级黄| 国产一区二区三区av在线 | 丰满人妻一区二区三区视频av| 99久久人妻综合| 草草在线视频免费看| 人妻久久中文字幕网| 美女大奶头视频| 干丝袜人妻中文字幕| 国产精品福利在线免费观看| 国产69精品久久久久777片| 亚洲真实伦在线观看| 精品无人区乱码1区二区| 男女啪啪激烈高潮av片| 国产视频内射| 成人三级黄色视频| 欧美另类亚洲清纯唯美| 亚洲美女视频黄频| 日韩欧美国产在线观看| 亚洲18禁久久av| 精品久久久噜噜| 国产男人的电影天堂91| 免费人成在线观看视频色| 哪里可以看免费的av片| 亚洲丝袜综合中文字幕| 亚洲最大成人av| 国产成人午夜福利电影在线观看| 国产一区二区在线av高清观看| 国产探花在线观看一区二区| 亚洲无线在线观看| 欧美+日韩+精品| 观看美女的网站| 亚洲最大成人中文| 美女内射精品一级片tv| 亚洲第一区二区三区不卡| 久久久色成人| 亚洲av中文av极速乱| 日本撒尿小便嘘嘘汇集6| 国产三级在线视频| 成人亚洲精品av一区二区| 在线观看一区二区三区| 久久精品国产99精品国产亚洲性色| 人妻夜夜爽99麻豆av| videossex国产| 精品少妇黑人巨大在线播放 | av在线天堂中文字幕| 99热这里只有是精品50| 欧美bdsm另类| 午夜激情福利司机影院| 你懂的网址亚洲精品在线观看 | 卡戴珊不雅视频在线播放| 国产精品野战在线观看| 日日撸夜夜添| 亚洲欧美精品综合久久99| 又粗又爽又猛毛片免费看| 高清毛片免费观看视频网站| 免费看光身美女| a级毛片a级免费在线| 亚洲av成人精品一区久久| 国产伦精品一区二区三区视频9| 性欧美人与动物交配| 精品人妻偷拍中文字幕| 亚洲精品影视一区二区三区av| 久久久久久久久中文| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 又爽又黄无遮挡网站| 91av网一区二区| 九九爱精品视频在线观看| 91麻豆精品激情在线观看国产| 国产亚洲av片在线观看秒播厂 | 天天躁夜夜躁狠狠久久av| 精品无人区乱码1区二区| 黄片无遮挡物在线观看| 午夜精品国产一区二区电影 | 中文字幕免费在线视频6| www.色视频.com| 国产精品麻豆人妻色哟哟久久 | 高清日韩中文字幕在线| 亚洲丝袜综合中文字幕| 国产精品女同一区二区软件| 成年版毛片免费区| 99久久精品国产国产毛片| 欧美一区二区精品小视频在线| 韩国av在线不卡| 97人妻精品一区二区三区麻豆| 午夜精品国产一区二区电影 | 国产91av在线免费观看| 97超视频在线观看视频| 中文字幕av在线有码专区| 深爱激情五月婷婷| 一个人看的www免费观看视频| 18禁裸乳无遮挡免费网站照片| 身体一侧抽搐| 成年免费大片在线观看| 久久精品国产亚洲av香蕉五月| 国产成人a区在线观看| 国产一区二区亚洲精品在线观看| 午夜爱爱视频在线播放| 午夜福利在线在线| 亚洲精品456在线播放app| 国内揄拍国产精品人妻在线| 亚洲在线观看片| 狂野欧美激情性xxxx在线观看| 精品少妇黑人巨大在线播放 | 国产亚洲av片在线观看秒播厂 | 日韩成人伦理影院| 一个人观看的视频www高清免费观看| 全区人妻精品视频| 三级毛片av免费| 日本色播在线视频| 国产大屁股一区二区在线视频| 亚洲内射少妇av| 国产av麻豆久久久久久久| 免费大片18禁| 一级黄色大片毛片| 国产黄色视频一区二区在线观看 | 日本av手机在线免费观看| 波野结衣二区三区在线| 亚洲国产精品成人久久小说 | 午夜久久久久精精品| 有码 亚洲区| 日韩在线高清观看一区二区三区| 亚洲av熟女| 在线天堂最新版资源| 51国产日韩欧美| 精华霜和精华液先用哪个| 国产又黄又爽又无遮挡在线| 国产精品爽爽va在线观看网站| 黄色配什么色好看| 精品欧美国产一区二区三| 国产亚洲5aaaaa淫片| 中文字幕免费在线视频6| 亚洲自偷自拍三级| 99热这里只有是精品50| av又黄又爽大尺度在线免费看 | 夜夜夜夜夜久久久久| 国产高潮美女av| 麻豆成人午夜福利视频| 国产极品精品免费视频能看的| 亚洲精品久久国产高清桃花| www.av在线官网国产| 一个人免费在线观看电影| 尾随美女入室| 中文字幕熟女人妻在线| 欧美变态另类bdsm刘玥| av免费在线看不卡| 久久精品国产鲁丝片午夜精品| 午夜a级毛片| 中出人妻视频一区二区| 欧美日韩在线观看h| 欧美丝袜亚洲另类| av在线蜜桃| 日韩亚洲欧美综合| 欧美日韩一区二区视频在线观看视频在线 | 男人舔女人下体高潮全视频| 国产大屁股一区二区在线视频| 91久久精品国产一区二区三区| 白带黄色成豆腐渣| 精品人妻熟女av久视频| 久久亚洲精品不卡| 精品久久久久久久人妻蜜臀av| 69av精品久久久久久| 成人毛片60女人毛片免费| 91精品国产九色| 亚洲自偷自拍三级| 久久精品国产亚洲av涩爱 | 亚洲av中文av极速乱| 亚洲美女搞黄在线观看| 精品久久国产蜜桃| 国产黄片美女视频| 久久精品夜色国产| 色综合亚洲欧美另类图片| 赤兔流量卡办理| 亚洲精品影视一区二区三区av| 天天躁夜夜躁狠狠久久av| 免费黄网站久久成人精品| 色哟哟·www| 国产白丝娇喘喷水9色精品| 亚洲成人久久爱视频| 精品国产三级普通话版| 天堂av国产一区二区熟女人妻| 91在线精品国自产拍蜜月| av在线蜜桃| 国产精品蜜桃在线观看 | 午夜久久久久精精品| 激情 狠狠 欧美| 美女脱内裤让男人舔精品视频 | 午夜福利高清视频| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久大精品| 九九在线视频观看精品| 亚洲成a人片在线一区二区| 99久久中文字幕三级久久日本| 色5月婷婷丁香| 国产亚洲精品久久久com| 成人三级黄色视频| 欧美最黄视频在线播放免费| 丝袜美腿在线中文| 中文字幕av成人在线电影| 99久久人妻综合| 舔av片在线| 高清午夜精品一区二区三区 | 国产精品一及| 熟妇人妻久久中文字幕3abv| 亚洲欧美日韩高清专用| 男女下面进入的视频免费午夜| 国产伦精品一区二区三区四那| 性插视频无遮挡在线免费观看| 亚洲精品日韩在线中文字幕 | 噜噜噜噜噜久久久久久91| 人体艺术视频欧美日本| 色视频www国产| 深夜精品福利| 黄色一级大片看看| 午夜激情欧美在线| 日日撸夜夜添| 亚洲在久久综合| 国产精品久久久久久精品电影小说 | 深夜精品福利| www.av在线官网国产| 亚洲人与动物交配视频| 91麻豆精品激情在线观看国产| 伦精品一区二区三区| 国产淫片久久久久久久久| 亚洲欧美成人综合另类久久久 | 国产精品福利在线免费观看| 欧美一级a爱片免费观看看| 欧美激情在线99| 亚洲av中文av极速乱| 国产色爽女视频免费观看| 亚洲七黄色美女视频| 中文字幕av成人在线电影| 一本精品99久久精品77| 国产黄片视频在线免费观看| 在线观看av片永久免费下载| 天堂影院成人在线观看| 精品免费久久久久久久清纯| 久久热精品热| 国产又黄又爽又无遮挡在线| 亚洲七黄色美女视频| 欧美+亚洲+日韩+国产| 真实男女啪啪啪动态图| 久久6这里有精品| 亚洲av中文av极速乱| 男女那种视频在线观看| 欧美在线一区亚洲| 日本五十路高清| 啦啦啦观看免费观看视频高清| 欧美成人a在线观看| 村上凉子中文字幕在线| 成人欧美大片| 亚洲va在线va天堂va国产| 最好的美女福利视频网| 高清毛片免费看| 久久人人爽人人爽人人片va| 欧美高清成人免费视频www| 久久综合国产亚洲精品| 亚洲欧美日韩东京热| 91久久精品国产一区二区成人| 婷婷精品国产亚洲av| 日本一二三区视频观看| 午夜福利在线在线| 午夜精品一区二区三区免费看| 六月丁香七月|