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

    固定翼電磁數據雙分量聯(lián)合電導率深度成像

    2015-03-06 11:58:56朱凱光李冰冰王凌群程宇奇
    關鍵詞:總場固定翼數據表

    朱凱光,李冰冰,王凌群,謝 賓,王 琦,程宇奇

    吉林大學儀器科學與電氣工程學院,長春 130026

    ?

    固定翼電磁數據雙分量聯(lián)合電導率深度成像

    朱凱光,李冰冰,王凌群,謝 賓,王 琦,程宇奇

    吉林大學儀器科學與電氣工程學院,長春 130026

    傳統(tǒng)固定翼航空電磁探測采用總場dB/dt進行電導率深度成像,不僅損失多分量測量信息,dB/dt與電導率的非唯一性也影響數據的成像精度。筆者提出了一種基于磁場雙分量(Bx,Bz)查表的聯(lián)合電導率深度成像算法,根據固定翼電磁響應的正演計算,建立按時間道劃分的Bx-Bz-電導率-飛行高度數據表,利用磁場雙分量聯(lián)合查表與插值算法確定視電導率,避免了由于電磁數據二值性引起的視導率不確定性問題;根據擴散深度公式得到視深度,并計算成像深度,從而得到雙分量聯(lián)合電導率深度成像結果?;谝痪S大地模型正演數據與準二維大地模型正演加噪數據,分別采用磁場雙分量聯(lián)合查表法、總場查表法和單分量查表法對仿真數據進行電導率深度成像,結果表明磁場雙分量聯(lián)合查表法優(yōu)于單分量與總場查表法,較單分量電導率深成像精度提高了7%。

    固定翼航空電磁;磁場雙分量查表法;聯(lián)合電導率深度成像;擴散深度

    0 引言

    固定翼航空電磁探測系統(tǒng)以固定翼飛機為載體,采用偶極-偶極收發(fā)裝置,發(fā)射線圈安置于飛機周圍,接收線圈由電纜線拖拽懸掛于飛機下方,可測量三分量電磁響應,具有勘探效率高、探測深度大、勘探面積廣等優(yōu)勢,廣泛應用于地質勘查中[1-3]。固定翼航空電磁數據量龐大,反演解釋耗時費力,一般采用電導率深度成像 ( conductivity-depth imaging,CDI) 技術快速處理海量航空電磁數據,得到地下電導率分布圖,從而確定地下異常分布[4-5]。

    目前國內外航空電磁探測的電導率深度成像方法很多,如Liu等[6]提出了基于薄板模型的成像方法,其視電導率及視深度由最小二乘迭代法導出;Wolfgram等[7]采用解卷積算法將任意波形的電磁響應轉換為階躍響應,通過τ域分解,得到地下電導率分布;Huang等[8]提出了基于假層-半空間查表的電導率深度成像方法,引入電導率為0的假想層代表人文圈;朱凱光等[9]利用神經網絡獲得更為準確的視電導率;陳小紅等[10]采用改進后的一階、二階電導率深度成像算法,提高了對地層電導率的靈敏性;毛立峰[11]建立了σ-β數據表查詢視電導率,消除了電磁響應的二值性。

    筆者根據磁場對導電異常體探測與單調性的優(yōu)勢[16]提出了一種基于B場響應的雙分量聯(lián)合電導率深度成像算法。該算法基于固定翼時間域航空電磁正演理論,采用層狀大地模型,推導出固定翼航空電磁系統(tǒng)x與z分量的B場響應,建立水平分量磁場(即Bx)-垂直分量磁場(即Bz)-電導率-飛行高度數據表,給出雙分量Bx-Bz聯(lián)合查表算法,最后通過理論仿真驗證雙分量查表法聯(lián)合電導率深度成像算法的有效性。

    1 正演算法

    圖1 固定翼時間域航空電磁系統(tǒng)坐標系Fig.1 Coordinate of the fixed-wing time-domain electromagnetic system

    根據電磁計算理論[17],接收線圈處的二次場頻率域水平分量磁場響應Bx(w)和垂直分量磁場響應Bz(w)的表達式分別為

    (1)

    (2)

    其中:J0和J1分別為零階和一階貝塞爾函數;w為角頻率;λ為積分變量;μ0為空氣中的磁導率,μ0=4×10-7π H/m;反射系數rTE為

    其中:μj=μ0。分別將式(1)、式(2)變換到s(復頻率)域中,并除以s得到x分量和z分量的磁場瞬態(tài)響應表達式:

    (3)

    (4)

    分別對式(3)、(4)做拉氏逆變換(L-1)得到系統(tǒng)雙分量B場響應,作為電導率深度成像的數據表來源:

    (5)

    (6)

    式中,零階和一階貝塞爾函數的無限積分(漢克爾變換)利用數值濾波算法(Guptasarma120點和140點)[18-19]計算。

    2 Bx-Bz分量查表法聯(lián)合電導率深度成像

    2.1 電磁響應樣本集的建立

    設定固定翼航空電磁探測系統(tǒng)保持平穩(wěn)飛行狀態(tài),飛行高度 在80~180 m等間隔選取101個,半空間模型的電導率在0.001~20.000 s/m按等對數間隔選取55個,形成101×55組半空間探測模型。利用上述正演算法,計算14個時間道(0.01~10.00 ms)各模型對應的雙分量B場響應。按時間道劃分,建立14個Bx-Bz-電導率-飛行高度數據表,表中任一點為由對應飛行高度及半空間模型確定的x和z分量B場響應,記為(Bx,Bz,σ,h)。圖2給出了第k(1≤k≤14)道的部分Bx-Bz-電導率-飛行高度數據表示意圖,其中電導率σ1<σ2<σ3,飛行高度h1>h2>h3>h4>h5>h6>h7>h8>h9。

    xk和zk為實際飛行測量的第k道x和z分量的B場響應。圖2 第k道Bx-Bz-電導率-飛行高度數據表示意圖Fig.2 Diagram of Bx-Bz-conductivity-altitude for channel k

    2.2Bx-Bz分量聯(lián)合查表算法

    傳統(tǒng)的總場查表法與單分量查表法電導率深度成像[20-22]采用總場響應或單分量響應進行匹配搜索查詢視電導率,在搜索過程中需要提供實際的飛行高度以選擇所需的數據表,而實際飛行探測過程中很難獲取準確的飛行高度,給電導率的查詢帶來一定困難。

    Bx-Bz分量聯(lián)合查表法建立了Bx-Bz-電導率-飛行高度數據表,將飛行高度引入數據表中,在無法準確提供高度信息時不但可以快速獲得視電導率,還可以得到對應的高度信息。同時采用x和z兩個分量的電磁響應進行匹配搜索,不但增加了搜索時的約束條件,而且綜合利用了x和z兩個分量電磁響應的特性,從而提高成像的精確度??倛霾楸黼m然也綜合利用了多分量的響應,但對于兩組不同的Bx、By和Bz分量,只要滿足

    就會得到相同的電導率查詢結果。因此,總場查表雖然利用了多個分量,但由于進行了二對一的轉換,損失了多分量測量信息。

    (7)

    (8)

    其中:q1--q4為線段長度。

    圖3 插值算法示意圖Fig.3 Diagram of interpolation algorithm

    利用擴散深度公式計算得到第k個時間道對應的視深度dk:

    式中,tk為第k個時間道對應的時間。在相鄰2個時間道深度變化的基礎上定義成像深度,即第k個時間道的成像深度是由上一時間道的成像深度加上利用趨膚深度公式計算得到的相鄰兩個時間道深度變化值的絕對值,從而獲得第k個時間道對應的成像深度Dk[11]:

    (9)

    綜上,Bx-Bz分量查表法聯(lián)合電導率深度成像的步驟如下:

    1)利用正演算法,建立Bx-Bz-電導率-飛行高度數據表表格;

    2)通過磁場雙分量聯(lián)合查表與插值算法,確定視電導率;

    3)利用擴散深度公式得到視深度,并計算成像深度,最終導出雙分量聯(lián)合電導率深度成像結果。

    3 仿真結果

    3.1 一維大地模型的電導率深度成像結果

    設定固定翼航空電磁系統(tǒng)的z分量發(fā)射線圈高度為h=120 m,三分量接收線圈位于(-70,0,70) m處,系統(tǒng)參數同前。針對不同層狀大地模型對Bx-Bz分量聯(lián)合查表法、總場查表法與單分量查表法CDI進行對比研究。

    以三層大地模型為例,選取了K型、H型典型地電模型作為參考模型對算法進行對比分析。表1為K、H型大地模型參數。

    表1 K、H型大地模型參數表

    Table 1 Parameter table of K and H style of geoelectric models

    模型參數電導率/(S/m)第一層第二層第三層厚度/m第一層第二層K型0.020.200.025050H型0.200.020.202030

    圖4給出了K、H型兩種三層大地模型采用Bx-Bz分量聯(lián)合查表法、總場查表法與單分量查表法的CDI結果。可以看到,總場查表法成像結果與Bx和Bz單分量的查表結果一致,大致給出了三層大地模型的變化趨勢,而Bx-Bz分量聯(lián)合查表算法對K型大地模型(圖4a)中間高導層的電導率和深度的成像結果與理論模型接近,即使H型大地的高阻層(圖4b),Bx-Bz分量聯(lián)合查表算法也獲得了優(yōu)于總場與單分量查表法的成像結果。而由于查表算法是利用半空間大地模型的電磁響應作為數據表,而觀測響應是含有異常體的三維大地模型;因此,利用觀測值在半空間模型表中查詢時,可能會出現部分點不在數據表中,使得成像結果不太光滑。

    a.K型地電模型;b.H型地電模型。圖4 三層典型大地模型CDI結果Fig.4 Results of four CDI methods on three-layer models

    a.Bx分量查表法;b. Bz分量查表法;c. 總場響應查表法;d. Bx-Bz聯(lián)合查表法。圖6 四種方法的電導率深度成像結果對比圖Fig.6 Results of four CDI methods on Pseudo-2D model

    綜上,總場查表法雖然利用了Bx和Bz兩個分量進行電導率深度成像,但是由于進行了二對一的轉換,僅獲得了與單分量查表法近似的成像結果,損失了多分量測量獲取地下信息的能力。而Bx-Bz分量聯(lián)合查表算法充分利用Bx和Bz兩個分量的信息,進行聯(lián)合查表,獲得了與真實大地模型更為接近的地下信息,提高了成像精度。

    3.2 準二維大地模型的電導率深度成像結果

    為評價Bx-Bz分量查表法聯(lián)合電導率深度成像算法在實際中的應用,仿真構建了準二維大地模型,由500個5層大地模型相互連接構成。如圖5所示,在測點100~250處的地下埋深2 m有一個厚度為5 m、電導率為0.2 s/m的異常,同時在埋深57 m處有一個厚度40 m、電導率0.2 s/m的異常,周圍基巖的電導率為0.02 s/m。對上述大地模型進行正演計算獲得其電磁響應剖面數據,并加入均方根相對誤差為5%的白噪聲,仿真野外實測數據。

    圖5 理論大地模型Fig.5 Pseudo-2D theoretical model

    為進一步對比成像效果,分別采用Bx分量查表法、Bz分量查表法、總場響應查表法及Bx-Bz分量聯(lián)合查表法對數據進行電導率深度成像。圖6分別展示了上述四種方法對于含噪準二維數據的CDI結果。

    可以看到,Bx分量查表法(圖6a)、Bz分量查表法(圖6b)與總場查表法(圖6c)的準二維含噪數據的CDI結果都較好地反映了埋深2 m的導電異常體,屏蔽了其下埋深57 m處的導電異常體,而圖6d給出的Bx-Bz分量聯(lián)合CDI結果卻較為清晰地揭示了導電薄層下的低阻異常體,與理論模型一致。筆者計算了測點175處電導率深度成像結果的相對誤差平均值得知,Bx分量查表、Bz分量查表、總場查表法以及Bx-Bz分量聯(lián)合查表法的相對誤差分別為22.07%、21.65%、21.84%、14.72%,即本文研究算法的成像精度較單分量電導率深度成像提高了7%左右。

    4 結論

    1)本文研究的基于磁場雙分量的聯(lián)合電導率深度成像算法,充分利用了多分量電磁探測數據,取得了優(yōu)于總場與單分量查表法的成像結果。

    2)一維大地模型的仿真結果表明該算法提高了對導電薄層的識別能力;準二維大地模型加噪數據的成像結果進一步證明在低阻覆蓋層存在的情況下,聯(lián)合電導率深度成像算法能夠增強對其下導電異常體的識別能力,綜合成像精度較單分量電導率深度成像提高了7%左右,為固定翼航空瞬變電磁資料解釋提供了新思路。

    吉林大學地球信息探測儀器教育部重點實驗室航空電磁組的全體成員對此文進行了幫助和指導,在此表示衷心的感謝!

    [1] Fountain D. Airborne Electromagnetic Systems:50 Years of Development[J]. Exploration Geophysics, 1998, 29(2): 1-11.

    [2] Annan A P.Effect of Differential Transmitter/Receiver Motion on Airborne Transient EM Interpretation[J]. Geophysics, 1999, 49(1): 670-670.

    [3] 胡平, 李文杰, 李軍峰, 等.固定翼時間域航空電磁勘查系統(tǒng)研發(fā)進展[J]. 地球學報,2012, 33(1): 7-12. Hu Ping, Li Wenjie, Li Junfeng, et al. The Advances in the Development of Fixed-Wing Airborne Time-Domain Electromagnetic System[J]. Acta Geoscientica Sinica, 2012, 33(1): 7-12.

    [4] Macnae J C, Smith R, Polzer B D, et al. Conductivity-Depth Imaging of Airborne Electromagnetic Step-Response Data[J]. Geophysics, 1991, 56(1): 102-114.

    [5] 嵇艷鞠, 馮雪, 于明媚, 等.基于多元線性回歸的HTEM三維異常體電導率-深度識別[J]. 吉林大學學報:地球科學版,2014, 44(5), 1687-1694. Ji Yanju, Feng Xue, Yu Mingmei, et al. Conductivity-Depth Identification of HTEM 3D Anomalies Based on Multiple Linear Regression[J]. Journal of Jilin University: Earth Science Edition, 2014, 44(5), 1687-1694.

    [6] Liu G, Asten M. Conductance-Depth Imaging of Airborne TEM data[J]. Exploration Geophysics,1994, 24(4): 655-662.

    [7] Wolfgram P, Karlik G. Conductivity-Depth Transform of GEOTEM Data[J]. Exploration Geophysics, 1995, 26(2/3): 179-185.

    [8] Huang H, Rudd J. Conductivity-Depth Imaging of Helicopter-Borne TEM Data Based on a Pseudolayer Half-Space Model[J]. Geophysics, 2008,73(3): 115-120.

    [9] 朱凱光, 林君, 韓悅慧, 等.基于神經網絡的時間域直升機電磁數據電導率深度成像[J]. 地球物理學報,2010, 53(3): 743-750. Zhu Kaiguang, Lin Jun, Han Yuehui, et al. Research on Conductivity Depth Imaging of Time Domain Helicopter-Borne Electromagnetic Data Based on Neural Network[J]. Chinese Journal of Geophysics, 2010, 53(3): 743-750.

    [10] 陳小紅, 段奶軍.時間域航空電磁快速成像研究[J]. 地球物理學進展, 2012,27(5): 2123-2127. Chen Xiaohong, Duan Naijun. Study on Fast Imaging of Airborne Time-Domain Electromagnetic Data[J]. Progress in Geophys., 2012, 27(5): 2123-2127.

    [11] 毛立峰.中心回線式直升機TEM資料的電導率-深度成像方法[J]. CT 理論與應用研究, 2013,22(3): 429-437. Mao Lifeng. Conductivity-Depth Imaging Algorithm for Central-Loop Helicopter TEM[J]. CT Theory and Applications, 2013, 22(3): 429-437.

    [12] Smith R S, Keating P B. The Usefulness of Multicomponent, Time-Domain Airborne Electromagnetic Measurements[J]. Geophysics, 1996, 61(1): 74-81.

    [13] 席振銖, 劉劍, 龍霞, 等.瞬變電磁法三分量測量方法研究[J]. 中南大學學報:自然科學版,2010,41(1):272-276. Xi Zhenzhu, Liu Jian, Long Xia, et al.Three-Component Measurement in Transient Electromagnetic Method[J]. Journal of Central South University: Science and Technology, 2010, 41(1): 272-276.

    [14] 張瑩瑩.水平電偶源地空系統(tǒng)瞬變電磁法多分量解釋技術及全域視電阻率定義研究[D]. 西安: 長安大學,2013. Zhang Yingying. Study on Multi-Component Interpretation and Full Field Apparat Resistivity Definition of Semi-Airborne Transient Electromagnetic Method with Electric Dipole on the Surface[D]. Xi’an: Chang’an University, 2013.

    [15] 王琦, 林君, 于生寶, 等. 固定翼航空電磁系統(tǒng)的線圈姿態(tài)及吊艙擺動影響研究與校正[J]. 地球物理學報,2013,56(11):3741-3750. Wang Qi, Lin Jun, Yu Shengbao, et al. Study on Influence and Correction of Coil Attitude and Bird Swing for Fixed-Wing Time-Domain Electromagnetic System[J]. Chinese Journal of Geophysics, 2013, 56(11):3741-3750.

    [16] Huang H. Locating Good Conductors by Using the B-Field Integrated From Partial dB/dtWaveforms of Time-Domain EM Systems[C]. SEG/San Antonio Annual Meeting. San Antonio: SEG, 2007: 688-692.

    [17] Wait J R. Geo-Electromagnetism[M]. [S.l.]: Academic Press Inc,1982.

    [18] Guptasarma D, Singh B. New Digital Linear Filters for Hankel J0and J1Transforms[J]. Geophysical Prospecting, 1997, 45(5): 745-762.

    [19] Guptasarma D, Singh B. Computation of the Time-Domain Response of a Polarizable Ground[J]. Geophysics, 1982, 47(11): 953-963.

    [20] Wolfgram P, Karlik G. Conductivity-Depth Transform of GEOTEM Data[J].Exploration Geophysics, 1995, 26(3): 179-185.

    [21] Vrbancich J, Fullagar P K, Macnae J. Bathymetry and Seafloor Mapping via One Dimensional Inversion and Conductivity Depth Imaging of AEM[J]. Exploration Geophysics, 2000,31(4): 603-610.

    [22] 馬銘遙.時間域航空電磁數據全波成像與處理技術研究[D]. 長春: 吉林大學,2013. Ma Mingyao. Study on Full-Waveform Imaging and Processing of Time-Domain Airborne Electromagnetic Data[D]. Changchun: Jilin University, 2013.

    Conductivity-Depth Imaging of Fixed-Wing Airborne Electromagnetic Data

    Zhu Kaiguang,Li Bingbing,Wang Lingqun,Xie Bin,Wang Qi,Cheng Yuqi

    CollegeofInstrumentationandElectricalEngineering,JilinUniversity,Changchun130026,China

    The traditional fixed-wing airborne electromagnetic survey lost multi-component measurement information for the reason of using the total field as the data base for the conductivity depth imaging. The imaging precision is decreased due to the non-uniqueness of the dB/dt. A joint conductivity depth imaging approach based on the lookup method of two-componentBfield response is proposed. Based on the fixed-wing electromagnetic forward calculation, this paper establishes a table ofBx-Bz-conductivity-altitude, usesBfield to avoid the uncertainty of conductivity due to the binary value of electromagnetic data, uses the joint lookup method of two-componentBfield and interpolation algorithm to determine the conductivity, and then takes the advantage of the diffusion depth formula to obtain the depth, and finally to process the joint conductivity depth imaging. Compared to the total field and simple component lookup method applied in the one-dimensional and quasi-two-dimensional forward simulation, the CDI of two-componentBfield gets more ideal fitting effects; and the imaging precision is improved by 7%.

    fixed-wing airborne electromagnetic;lookup method of two-component magnetic field;joint conductivity-depth imaging;diffusion depth

    10.13278/j.cnki.jjuese.201506301.

    2015-02-08

    國家自然科學基金項目(41274076);國家高技術研究發(fā)展計劃項目(2013AA063904);國家重大科研裝備研制項目(ZDYZ2012-1-03)

    朱凱光(1970--),女,教授,主要從事地球物理電磁探測及隨機信號處理研究,E-mail:zhukaiguang@jlu.edu.cn。

    10.13278/j.cnki.jjuese.201506301

    P631.3

    A

    朱凱光,李冰冰,王凌群,等. 固定翼電磁數據雙分量聯(lián)合電導率深度成像.吉林大學學報:地球科學版,2015,45(6):1839-1845.

    Zhu Kaiguang,Li Bingbing,Wang Lingqun, et al. Conductivity-Depth Imaging of Fixed-Wing Airborne Electromagnetic Data.Journal of Jilin University:Earth Science Edition,2015,45(6):1839-1845.doi:10.13278/j.cnki.jjuese.201506301.

    猜你喜歡
    總場固定翼數據表
    垂直起降固定翼無人機串聯(lián)混電系統(tǒng)優(yōu)化設計
    湖北省新冠肺炎疫情數據表
    黨員生活(2020年2期)2020-04-17 09:56:30
    綜合施策打好棉花田管“組合拳”
    一種新型固定翼無人機的研制
    電子制作(2019年7期)2019-04-25 13:18:24
    基于列控工程數據表建立線路拓撲關系的研究
    “V”系列固定翼無人機
    航空模型(2016年10期)2017-05-09 06:25:59
    前向雷達目標回波成分與特性分析
    水下滑翔機水平固定翼設計
    石總場早播棉花出苗顯行
    圖表
    网址你懂的国产日韩在线| 久久中文看片网| 一个人看的www免费观看视频| 好男人在线观看高清免费视频| 国产免费av片在线观看野外av| 村上凉子中文字幕在线| 欧美乱色亚洲激情| 国产欧美日韩一区二区精品| av国产免费在线观看| 国产成人a区在线观看| 亚洲不卡免费看| 成人鲁丝片一二三区免费| 一个人免费在线观看电影| 久久精品国产亚洲av涩爱 | 久9热在线精品视频| 亚洲国产欧美人成| 国产精品99久久久久久久久| 能在线免费观看的黄片| 一个人观看的视频www高清免费观看| 听说在线观看完整版免费高清| 午夜福利免费观看在线| 黄色女人牲交| 国产在视频线在精品| 久久久国产成人免费| 高清毛片免费观看视频网站| 少妇人妻一区二区三区视频| 自拍偷自拍亚洲精品老妇| 久久热精品热| 人人妻人人澡欧美一区二区| 精品一区二区三区人妻视频| 欧美绝顶高潮抽搐喷水| 两性午夜刺激爽爽歪歪视频在线观看| 99久久成人亚洲精品观看| 午夜激情福利司机影院| 日韩 亚洲 欧美在线| 简卡轻食公司| 精品一区二区三区视频在线| av视频在线观看入口| 一级黄色大片毛片| 成人美女网站在线观看视频| 久久久精品大字幕| 国产精品久久电影中文字幕| 18美女黄网站色大片免费观看| 欧美zozozo另类| 黄色日韩在线| 日日夜夜操网爽| 99热这里只有精品一区| 亚洲欧美日韩高清专用| 性色av乱码一区二区三区2| 久久精品综合一区二区三区| 草草在线视频免费看| 色噜噜av男人的天堂激情| 一级a爱片免费观看的视频| 亚洲第一区二区三区不卡| 成人欧美大片| 欧美极品一区二区三区四区| 99在线视频只有这里精品首页| 好男人在线观看高清免费视频| 岛国在线免费视频观看| 国产极品精品免费视频能看的| 麻豆久久精品国产亚洲av| 91久久精品电影网| 国内揄拍国产精品人妻在线| 嫁个100分男人电影在线观看| 久久午夜福利片| 精品无人区乱码1区二区| 特大巨黑吊av在线直播| 日本免费一区二区三区高清不卡| 小蜜桃在线观看免费完整版高清| 欧美另类亚洲清纯唯美| 可以在线观看的亚洲视频| 两个人视频免费观看高清| 99在线人妻在线中文字幕| av国产免费在线观看| 一区二区三区高清视频在线| 黄色一级大片看看| 欧美潮喷喷水| 亚洲精品乱码久久久v下载方式| 精品熟女少妇八av免费久了| 久久6这里有精品| 亚洲熟妇熟女久久| 黄色女人牲交| 少妇人妻一区二区三区视频| 日本免费一区二区三区高清不卡| 网址你懂的国产日韩在线| 欧美日本视频| 亚洲成av人片在线播放无| 天天一区二区日本电影三级| 在线天堂最新版资源| 国产精品一区二区性色av| 午夜福利在线在线| 国产亚洲精品久久久久久毛片| 天堂av国产一区二区熟女人妻| 可以在线观看毛片的网站| 国内久久婷婷六月综合欲色啪| 人妻丰满熟妇av一区二区三区| 国产在线男女| 国产白丝娇喘喷水9色精品| 又紧又爽又黄一区二区| 亚洲一区二区三区不卡视频| 国产成+人综合+亚洲专区| 欧美激情久久久久久爽电影| 99视频精品全部免费 在线| 日韩人妻高清精品专区| 日韩精品青青久久久久久| 久久精品91蜜桃| 国产免费男女视频| 国产探花在线观看一区二区| 日韩大尺度精品在线看网址| 99久久99久久久精品蜜桃| 欧美日韩国产亚洲二区| 麻豆成人av在线观看| 色综合婷婷激情| 级片在线观看| 欧美国产日韩亚洲一区| 国产综合懂色| 欧美国产日韩亚洲一区| 国产色婷婷99| 757午夜福利合集在线观看| 在线观看免费视频日本深夜| 久久欧美精品欧美久久欧美| 色噜噜av男人的天堂激情| 色在线成人网| 久久国产乱子免费精品| 欧美午夜高清在线| 久久这里只有精品中国| 国产蜜桃级精品一区二区三区| 九九热线精品视视频播放| 嫩草影院新地址| 久久久久九九精品影院| 男人狂女人下面高潮的视频| 色综合欧美亚洲国产小说| 99久久无色码亚洲精品果冻| 91久久精品国产一区二区成人| 国产麻豆成人av免费视频| 高清在线国产一区| 在线观看舔阴道视频| 欧美成狂野欧美在线观看| 亚洲最大成人中文| 伊人久久精品亚洲午夜| 特级一级黄色大片| 好看av亚洲va欧美ⅴa在| 国产精品99久久久久久久久| 国产探花极品一区二区| 成人特级黄色片久久久久久久| 日韩免费av在线播放| 不卡一级毛片| 午夜福利在线观看免费完整高清在 | 国产aⅴ精品一区二区三区波| 免费人成在线观看视频色| 搞女人的毛片| 日本黄色视频三级网站网址| 日本一二三区视频观看| 最近视频中文字幕2019在线8| 欧美精品国产亚洲| 九色成人免费人妻av| 欧美日韩中文字幕国产精品一区二区三区| 黄色视频,在线免费观看| 两个人视频免费观看高清| 日本 欧美在线| 1024手机看黄色片| 最近最新免费中文字幕在线| 色综合站精品国产| 欧美xxxx黑人xx丫x性爽| 亚洲成av人片免费观看| 简卡轻食公司| 51午夜福利影视在线观看| 久久精品人妻少妇| 午夜日韩欧美国产| 亚洲中文日韩欧美视频| 日日夜夜操网爽| 亚洲经典国产精华液单 | 免费av不卡在线播放| 精品国内亚洲2022精品成人| 国产综合懂色| 国产久久久一区二区三区| 婷婷丁香在线五月| 成人鲁丝片一二三区免费| 日本黄色片子视频| 亚洲精品456在线播放app | 婷婷色综合大香蕉| 精品国产三级普通话版| 婷婷色综合大香蕉| 亚洲国产精品sss在线观看| 美女免费视频网站| 精品久久久久久久久亚洲 | 欧美三级亚洲精品| 亚洲成av人片在线播放无| 国产aⅴ精品一区二区三区波| 午夜视频国产福利| 日本成人三级电影网站| 精品国产亚洲在线| 国语自产精品视频在线第100页| 国产亚洲精品综合一区在线观看| 国产在线精品亚洲第一网站| 日韩大尺度精品在线看网址| 精品久久久久久久久久免费视频| 九色国产91popny在线| 黄色丝袜av网址大全| 亚洲av电影不卡..在线观看| 成年女人毛片免费观看观看9| 日韩欧美国产在线观看| 欧美潮喷喷水| 成人美女网站在线观看视频| 色尼玛亚洲综合影院| 又爽又黄无遮挡网站| 男人的好看免费观看在线视频| 亚洲成av人片免费观看| 亚洲黑人精品在线| 日韩欧美免费精品| 乱人视频在线观看| 国产精品三级大全| 久久久久久久久久黄片| 色综合亚洲欧美另类图片| 久久亚洲精品不卡| 国产午夜精品久久久久久一区二区三区 | 性插视频无遮挡在线免费观看| 国产av不卡久久| a级毛片a级免费在线| 亚洲精华国产精华精| 久久久色成人| 亚洲人成电影免费在线| 国产av麻豆久久久久久久| 在线观看66精品国产| 亚洲国产精品sss在线观看| 成年女人永久免费观看视频| 综合色av麻豆| 国产欧美日韩一区二区精品| 成人美女网站在线观看视频| 国产精品亚洲一级av第二区| 男女下面进入的视频免费午夜| 日韩欧美免费精品| 久久久国产成人免费| 国产单亲对白刺激| 亚洲avbb在线观看| 高清在线国产一区| 在线观看舔阴道视频| 熟妇人妻久久中文字幕3abv| 久久精品人妻少妇| 国产精品日韩av在线免费观看| 嫁个100分男人电影在线观看| 尤物成人国产欧美一区二区三区| 看十八女毛片水多多多| 欧美又色又爽又黄视频| 精品久久久久久久久亚洲 | 美女大奶头视频| 欧美一区二区亚洲| 欧美乱妇无乱码| 欧美激情久久久久久爽电影| 欧美区成人在线视频| 精品久久久久久久久亚洲 | 91字幕亚洲| 国产不卡一卡二| 亚洲精品色激情综合| 国产伦人伦偷精品视频| 日日干狠狠操夜夜爽| 日本一二三区视频观看| 高清在线国产一区| 在线国产一区二区在线| 熟女电影av网| 两个人的视频大全免费| 日韩中字成人| 欧美在线黄色| 国产91精品成人一区二区三区| 国产一区二区激情短视频| 日韩欧美免费精品| 亚洲成人久久爱视频| 最新在线观看一区二区三区| 老司机午夜十八禁免费视频| 欧美黄色淫秽网站| 亚洲av电影不卡..在线观看| 国产久久久一区二区三区| 中文字幕av在线有码专区| 天天一区二区日本电影三级| 三级男女做爰猛烈吃奶摸视频| 欧美不卡视频在线免费观看| 成人毛片a级毛片在线播放| 热99在线观看视频| 亚洲精品粉嫩美女一区| 他把我摸到了高潮在线观看| 亚洲熟妇中文字幕五十中出| 国产精品爽爽va在线观看网站| 欧美性猛交黑人性爽| 国产伦人伦偷精品视频| 欧美成人一区二区免费高清观看| 三级毛片av免费| 国产精品亚洲av一区麻豆| 极品教师在线视频| 99riav亚洲国产免费| 少妇裸体淫交视频免费看高清| 极品教师在线免费播放| 男女之事视频高清在线观看| 成人无遮挡网站| 欧美日韩福利视频一区二区| 国产成人啪精品午夜网站| 亚洲成av人片在线播放无| 成年女人看的毛片在线观看| 最新中文字幕久久久久| 欧美黄色片欧美黄色片| 国产黄a三级三级三级人| 国产精品女同一区二区软件 | 九色成人免费人妻av| 在现免费观看毛片| 午夜福利成人在线免费观看| 极品教师在线视频| 可以在线观看毛片的网站| 欧美三级亚洲精品| 男女做爰动态图高潮gif福利片| а√天堂www在线а√下载| 欧美极品一区二区三区四区| 中文在线观看免费www的网站| 亚洲久久久久久中文字幕| 日本三级黄在线观看| 日韩欧美在线二视频| 丁香六月欧美| 亚洲成人免费电影在线观看| 少妇被粗大猛烈的视频| 有码 亚洲区| a在线观看视频网站| 简卡轻食公司| 欧美另类亚洲清纯唯美| 亚洲一区二区三区不卡视频| 国产视频一区二区在线看| 1000部很黄的大片| 亚洲片人在线观看| 性插视频无遮挡在线免费观看| 很黄的视频免费| 久久人妻av系列| 国产精品久久久久久久电影| 久久精品国产自在天天线| 一进一出好大好爽视频| 中亚洲国语对白在线视频| 欧美激情在线99| 老司机午夜十八禁免费视频| 日本一本二区三区精品| 国产视频一区二区在线看| 亚洲片人在线观看| 深爱激情五月婷婷| 能在线免费观看的黄片| 成人亚洲精品av一区二区| 特级一级黄色大片| 麻豆久久精品国产亚洲av| 免费观看精品视频网站| 色视频www国产| 国产 一区 欧美 日韩| 黄色日韩在线| 一a级毛片在线观看| 男女之事视频高清在线观看| 一夜夜www| 国模一区二区三区四区视频| 99热这里只有是精品在线观看 | 久久久久国产精品人妻aⅴ院| 国产免费一级a男人的天堂| 精品国产三级普通话版| 国产欧美日韩一区二区三| 我要搜黄色片| 中文字幕精品亚洲无线码一区| 国产精品一区二区三区四区免费观看 | 亚洲av二区三区四区| 人妻制服诱惑在线中文字幕| 日本黄色视频三级网站网址| 午夜福利18| 精品午夜福利在线看| 老鸭窝网址在线观看| 欧美色欧美亚洲另类二区| 尤物成人国产欧美一区二区三区| 久久久久久久久中文| 熟女电影av网| 18禁黄网站禁片免费观看直播| 亚洲熟妇熟女久久| 国产精品久久久久久人妻精品电影| 精品熟女少妇八av免费久了| 在线播放国产精品三级| 好看av亚洲va欧美ⅴa在| 一本综合久久免费| 亚洲精品一区av在线观看| 欧美黑人欧美精品刺激| 免费搜索国产男女视频| 午夜激情福利司机影院| 日本在线视频免费播放| 日本三级黄在线观看| h日本视频在线播放| 久久精品国产亚洲av香蕉五月| 亚洲成av人片免费观看| 男人狂女人下面高潮的视频| 欧美潮喷喷水| 麻豆成人av在线观看| 亚洲国产精品sss在线观看| 色哟哟·www| 一级黄片播放器| 欧美精品啪啪一区二区三区| 国产91精品成人一区二区三区| 欧美成人性av电影在线观看| 精品人妻熟女av久视频| 成人欧美大片| 国产精品影院久久| 亚洲精品久久国产高清桃花| 欧美国产日韩亚洲一区| 精品不卡国产一区二区三区| 国产欧美日韩精品亚洲av| 高清毛片免费观看视频网站| 又爽又黄无遮挡网站| 精品国内亚洲2022精品成人| 婷婷精品国产亚洲av| 国产精品免费一区二区三区在线| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久亚洲av鲁大| 99国产综合亚洲精品| 日韩亚洲欧美综合| 黄色丝袜av网址大全| 精品一区二区三区视频在线观看免费| 精华霜和精华液先用哪个| 少妇裸体淫交视频免费看高清| 久久99热6这里只有精品| 欧美日韩综合久久久久久 | 国产精品一区二区免费欧美| ponron亚洲| 九九热线精品视视频播放| 美女高潮的动态| 中文资源天堂在线| 婷婷色综合大香蕉| 久久性视频一级片| 高清在线国产一区| 国产极品精品免费视频能看的| 精品国产三级普通话版| 亚洲精华国产精华精| 国产又黄又爽又无遮挡在线| 好男人电影高清在线观看| 超碰av人人做人人爽久久| 亚洲五月婷婷丁香| 99久久成人亚洲精品观看| 国产老妇女一区| 国内久久婷婷六月综合欲色啪| 日本黄色视频三级网站网址| 香蕉av资源在线| 国产一区二区激情短视频| 国产精品永久免费网站| 免费看a级黄色片| 最后的刺客免费高清国语| 特级一级黄色大片| 俄罗斯特黄特色一大片| 久久人妻av系列| 91字幕亚洲| 欧美乱妇无乱码| 久久久久久久久久成人| 变态另类丝袜制服| 婷婷丁香在线五月| 老熟妇仑乱视频hdxx| 欧美激情在线99| 一个人看视频在线观看www免费| 我的女老师完整版在线观看| 97超级碰碰碰精品色视频在线观看| 国产免费一级a男人的天堂| 中亚洲国语对白在线视频| 内射极品少妇av片p| 国产三级黄色录像| 尤物成人国产欧美一区二区三区| 国内揄拍国产精品人妻在线| 国产成人aa在线观看| 麻豆一二三区av精品| 男人舔女人下体高潮全视频| 两人在一起打扑克的视频| 国产探花极品一区二区| 国产欧美日韩精品亚洲av| 成人鲁丝片一二三区免费| 日韩免费av在线播放| 最近最新免费中文字幕在线| 欧美乱色亚洲激情| 久久99热这里只有精品18| 99久国产av精品| xxxwww97欧美| 免费av不卡在线播放| 九九热线精品视视频播放| 中出人妻视频一区二区| 神马国产精品三级电影在线观看| 欧美精品国产亚洲| www.熟女人妻精品国产| 最新中文字幕久久久久| 久久99热6这里只有精品| 国产野战对白在线观看| 美女cb高潮喷水在线观看| 久久精品国产亚洲av涩爱 | 成年版毛片免费区| 国产私拍福利视频在线观看| 最新在线观看一区二区三区| 每晚都被弄得嗷嗷叫到高潮| a级毛片免费高清观看在线播放| xxxwww97欧美| 国产极品精品免费视频能看的| 制服丝袜大香蕉在线| 成人欧美大片| 他把我摸到了高潮在线观看| 最后的刺客免费高清国语| 欧美+日韩+精品| 久久人人爽人人爽人人片va | 久久精品国产亚洲av天美| 国产精品久久久久久亚洲av鲁大| 狂野欧美白嫩少妇大欣赏| 国语自产精品视频在线第100页| 一夜夜www| 欧美潮喷喷水| 少妇被粗大猛烈的视频| aaaaa片日本免费| 亚洲精品在线美女| 一区二区三区免费毛片| 麻豆一二三区av精品| 五月玫瑰六月丁香| 国产精品久久久久久久久免 | 国产一区二区在线av高清观看| 老司机福利观看| 亚洲人成网站高清观看| 欧美乱色亚洲激情| 黄色视频,在线免费观看| 免费在线观看日本一区| 国产精品98久久久久久宅男小说| 露出奶头的视频| 精品久久久久久久久久久久久| 欧美不卡视频在线免费观看| 琪琪午夜伦伦电影理论片6080| 老女人水多毛片| 精品国产三级普通话版| 久久久精品欧美日韩精品| 黄色丝袜av网址大全| 欧美乱色亚洲激情| 久久精品人妻少妇| 亚洲成人久久爱视频| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久久av| 美女 人体艺术 gogo| 白带黄色成豆腐渣| 久久人人精品亚洲av| 色精品久久人妻99蜜桃| 日本一本二区三区精品| 国产精品影院久久| 免费大片18禁| 欧美最黄视频在线播放免费| 亚洲熟妇熟女久久| av黄色大香蕉| 午夜福利欧美成人| 午夜福利欧美成人| 亚洲人成电影免费在线| 黄色配什么色好看| 嫩草影院新地址| 真人做人爱边吃奶动态| АⅤ资源中文在线天堂| 日韩欧美在线二视频| 一a级毛片在线观看| 欧美日韩福利视频一区二区| 波野结衣二区三区在线| 免费搜索国产男女视频| 久久久久免费精品人妻一区二区| 国产精品自产拍在线观看55亚洲| 人妻丰满熟妇av一区二区三区| 亚洲精品色激情综合| 国产精品伦人一区二区| 国内少妇人妻偷人精品xxx网站| 国产欧美日韩一区二区精品| 欧美日韩综合久久久久久 | 久久人妻av系列| 99热只有精品国产| 天堂影院成人在线观看| 97超视频在线观看视频| 亚洲男人的天堂狠狠| 麻豆久久精品国产亚洲av| 97热精品久久久久久| 黄色一级大片看看| 嫩草影院新地址| 啦啦啦韩国在线观看视频| 国产精品一区二区三区四区久久| 欧美极品一区二区三区四区| 国产男靠女视频免费网站| 欧美色欧美亚洲另类二区| 在线天堂最新版资源| 一本久久中文字幕| 国产精品98久久久久久宅男小说| 午夜老司机福利剧场| 高潮久久久久久久久久久不卡| 午夜福利欧美成人| 亚洲美女搞黄在线观看 | 毛片一级片免费看久久久久 | 亚洲国产精品999在线| 国产69精品久久久久777片| 国产精品一区二区性色av| 最后的刺客免费高清国语| 99久久九九国产精品国产免费| 欧美日韩乱码在线| 少妇的逼水好多| 黄片小视频在线播放| 淫妇啪啪啪对白视频| 亚洲av免费在线观看| 亚洲av中文字字幕乱码综合| 哪里可以看免费的av片| 日日摸夜夜添夜夜添小说| 男女床上黄色一级片免费看| 人妻丰满熟妇av一区二区三区| 精品一区二区免费观看| 国产成人a区在线观看| 亚洲精品日韩av片在线观看| 亚洲男人的天堂狠狠| 精华霜和精华液先用哪个| 精品99又大又爽又粗少妇毛片 | 午夜免费男女啪啪视频观看 | 久久久精品大字幕| 亚洲熟妇中文字幕五十中出| 天堂av国产一区二区熟女人妻| 亚洲国产欧美人成| 欧美+日韩+精品| 久久久久久久午夜电影| 搡老熟女国产l中国老女人| or卡值多少钱| 国产成人a区在线观看| 久久久久久久久久黄片| 国产在视频线在精品| 亚洲自偷自拍三级|