• <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
    前向雷達目標回波成分與特性分析
    水下滑翔機水平固定翼設計
    石總場早播棉花出苗顯行
    圖表
    亚洲国产日韩一区二区| 国产精品亚洲av一区麻豆 | 中文字幕av电影在线播放| 亚洲欧美中文字幕日韩二区| 国产亚洲av片在线观看秒播厂| 国产精品二区激情视频| av网站免费在线观看视频| 亚洲av.av天堂| 最近中文字幕高清免费大全6| 在线免费观看不下载黄p国产| 国产精品麻豆人妻色哟哟久久| 国产精品.久久久| 久久99蜜桃精品久久| 色哟哟·www| 国产精品久久久久成人av| 男人舔女人的私密视频| 91精品三级在线观看| 亚洲一区二区三区欧美精品| 国产亚洲午夜精品一区二区久久| 美女国产高潮福利片在线看| 精品一区二区三卡| 亚洲av.av天堂| 国产有黄有色有爽视频| 久久精品国产亚洲av涩爱| 老汉色av国产亚洲站长工具| 久久毛片免费看一区二区三区| 久久精品熟女亚洲av麻豆精品| 伊人久久大香线蕉亚洲五| 久久精品国产亚洲av天美| 国产精品蜜桃在线观看| 久久久久久人人人人人| 丝瓜视频免费看黄片| 欧美老熟妇乱子伦牲交| 99re6热这里在线精品视频| 黄色配什么色好看| 久久久国产欧美日韩av| 亚洲图色成人| 韩国av在线不卡| 97精品久久久久久久久久精品| 五月天丁香电影| 少妇人妻久久综合中文| 亚洲美女视频黄频| av免费观看日本| 1024视频免费在线观看| 黄频高清免费视频| 色94色欧美一区二区| 女人久久www免费人成看片| 精品酒店卫生间| 亚洲,欧美,日韩| 欧美变态另类bdsm刘玥| 少妇人妻精品综合一区二区| 精品卡一卡二卡四卡免费| 丝袜美足系列| 日韩欧美精品免费久久| 亚洲,一卡二卡三卡| 在现免费观看毛片| 精品少妇一区二区三区视频日本电影 | 久久国产精品男人的天堂亚洲| 天天躁日日躁夜夜躁夜夜| 精品久久久精品久久久| 91久久精品国产一区二区三区| 成人18禁高潮啪啪吃奶动态图| 久久热在线av| 日本午夜av视频| av免费观看日本| 人人妻人人澡人人看| 这个男人来自地球电影免费观看 | 九九爱精品视频在线观看| 国产色婷婷99| 国产麻豆69| av.在线天堂| 在线天堂最新版资源| 老汉色av国产亚洲站长工具| 五月开心婷婷网| 老汉色av国产亚洲站长工具| 国产一区二区激情短视频 | 久久精品国产综合久久久| 看免费av毛片| 看免费av毛片| 亚洲婷婷狠狠爱综合网| 性色avwww在线观看| 日韩免费高清中文字幕av| 久久久久久久国产电影| 成人免费观看视频高清| 观看av在线不卡| 制服人妻中文乱码| 亚洲,欧美精品.| 男的添女的下面高潮视频| 可以免费在线观看a视频的电影网站 | 国产精品免费视频内射| 曰老女人黄片| 精品视频人人做人人爽| 亚洲av日韩在线播放| 中国三级夫妇交换| 老司机亚洲免费影院| 老熟女久久久| 岛国毛片在线播放| 国产一区二区在线观看av| 一本久久精品| 亚洲精品久久成人aⅴ小说| 高清黄色对白视频在线免费看| 国产精品 欧美亚洲| 亚洲三区欧美一区| 久久精品人人爽人人爽视色| 国产成人精品婷婷| 婷婷色av中文字幕| 看十八女毛片水多多多| 一边亲一边摸免费视频| 妹子高潮喷水视频| 不卡视频在线观看欧美| 日韩一卡2卡3卡4卡2021年| 王馨瑶露胸无遮挡在线观看| 美女中出高潮动态图| 母亲3免费完整高清在线观看 | 精品99又大又爽又粗少妇毛片| 亚洲经典国产精华液单| 2018国产大陆天天弄谢| 欧美成人精品欧美一级黄| 亚洲国产日韩一区二区| 精品一区二区免费观看| 久久久久网色| 久久久欧美国产精品| 日产精品乱码卡一卡2卡三| 国产日韩一区二区三区精品不卡| 飞空精品影院首页| 99精国产麻豆久久婷婷| 大码成人一级视频| 国产亚洲一区二区精品| 精品国产露脸久久av麻豆| 一本—道久久a久久精品蜜桃钙片| 99久久综合免费| 国产野战对白在线观看| 电影成人av| 制服人妻中文乱码| 秋霞在线观看毛片| 九草在线视频观看| 国产精品人妻久久久影院| 成年人午夜在线观看视频| 赤兔流量卡办理| 91国产中文字幕| 亚洲,一卡二卡三卡| 伊人久久国产一区二区| 满18在线观看网站| 一区二区三区四区激情视频| 精品亚洲乱码少妇综合久久| 欧美精品人与动牲交sv欧美| 久久综合国产亚洲精品| 日韩制服丝袜自拍偷拍| 日韩av在线免费看完整版不卡| xxxhd国产人妻xxx| 日韩视频在线欧美| 成人手机av| 国产免费视频播放在线视频| 涩涩av久久男人的天堂| 久久99精品国语久久久| 天天躁日日躁夜夜躁夜夜| 久久精品熟女亚洲av麻豆精品| 视频在线观看一区二区三区| 黄色毛片三级朝国网站| 国产精品香港三级国产av潘金莲 | 大话2 男鬼变身卡| 久久精品国产a三级三级三级| 中文字幕另类日韩欧美亚洲嫩草| 亚洲综合色惰| 人妻系列 视频| 老汉色av国产亚洲站长工具| 可以免费在线观看a视频的电影网站 | 如何舔出高潮| 亚洲欧美精品自产自拍| 黄色毛片三级朝国网站| 久久久欧美国产精品| 视频区图区小说| 亚洲av免费高清在线观看| 亚洲精品一二三| 成人国产麻豆网| 美女国产高潮福利片在线看| 捣出白浆h1v1| 男女边吃奶边做爰视频| 久久精品国产鲁丝片午夜精品| 中文天堂在线官网| 亚洲国产最新在线播放| 人妻一区二区av| 亚洲av成人精品一二三区| 丝袜脚勾引网站| 欧美+日韩+精品| 国产成人91sexporn| 亚洲欧美一区二区三区国产| 亚洲av电影在线进入| 日日撸夜夜添| 亚洲欧洲精品一区二区精品久久久 | 有码 亚洲区| 久久韩国三级中文字幕| 亚洲欧美中文字幕日韩二区| 久久久久久久国产电影| 久久亚洲国产成人精品v| 寂寞人妻少妇视频99o| 街头女战士在线观看网站| 有码 亚洲区| 欧美 亚洲 国产 日韩一| 久久精品国产鲁丝片午夜精品| 欧美日韩综合久久久久久| 亚洲av男天堂| 午夜老司机福利剧场| 女人高潮潮喷娇喘18禁视频| 大陆偷拍与自拍| 爱豆传媒免费全集在线观看| 侵犯人妻中文字幕一二三四区| 黄片小视频在线播放| av在线观看视频网站免费| 看免费av毛片| 熟妇人妻不卡中文字幕| av.在线天堂| www.熟女人妻精品国产| 欧美亚洲日本最大视频资源| 99国产综合亚洲精品| 国产激情久久老熟女| 我要看黄色一级片免费的| 久久国产精品男人的天堂亚洲| 最新的欧美精品一区二区| 久久青草综合色| 精品亚洲成a人片在线观看| 国产精品av久久久久免费| 爱豆传媒免费全集在线观看| av国产久精品久网站免费入址| 亚洲国产精品国产精品| 亚洲成人av在线免费| 18禁裸乳无遮挡动漫免费视频| 老司机影院成人| 精品一区二区三卡| 天堂中文最新版在线下载| 日本猛色少妇xxxxx猛交久久| 日本欧美国产在线视频| 久久婷婷青草| 一级,二级,三级黄色视频| 精品久久蜜臀av无| 国产老妇伦熟女老妇高清| 老女人水多毛片| 日本-黄色视频高清免费观看| 熟女少妇亚洲综合色aaa.| 亚洲欧美精品综合一区二区三区 | www.av在线官网国产| 人成视频在线观看免费观看| 精品午夜福利在线看| 王馨瑶露胸无遮挡在线观看| 在线亚洲精品国产二区图片欧美| 亚洲欧美中文字幕日韩二区| 波多野结衣av一区二区av| videosex国产| 欧美精品一区二区免费开放| 中文字幕最新亚洲高清| 97人妻天天添夜夜摸| 国产精品国产av在线观看| 日韩av不卡免费在线播放| 国产精品三级大全| 久久精品夜色国产| 久久精品熟女亚洲av麻豆精品| 永久免费av网站大全| 一本色道久久久久久精品综合| 男男h啪啪无遮挡| 伊人久久国产一区二区| 91午夜精品亚洲一区二区三区| 国产成人免费无遮挡视频| 人妻 亚洲 视频| 久久久久久伊人网av| 久久久久久免费高清国产稀缺| 日本免费在线观看一区| 毛片一级片免费看久久久久| 久久婷婷青草| 亚洲色图 男人天堂 中文字幕| 99九九在线精品视频| 赤兔流量卡办理| 波多野结衣一区麻豆| 丝袜美足系列| 日韩欧美一区视频在线观看| 91aial.com中文字幕在线观看| 久久青草综合色| 国产野战对白在线观看| 一级片免费观看大全| av天堂久久9| 青春草国产在线视频| 深夜精品福利| 午夜av观看不卡| 精品第一国产精品| 欧美亚洲日本最大视频资源| 久久女婷五月综合色啪小说| 免费播放大片免费观看视频在线观看| 亚洲精品av麻豆狂野| 精品99又大又爽又粗少妇毛片| 国产成人精品在线电影| 伦理电影大哥的女人| 少妇的丰满在线观看| 建设人人有责人人尽责人人享有的| 亚洲av日韩在线播放| 亚洲欧美一区二区三区久久| 久久精品国产亚洲av涩爱| 男女高潮啪啪啪动态图| 成人午夜精彩视频在线观看| 男女边吃奶边做爰视频| 午夜福利乱码中文字幕| 又粗又硬又长又爽又黄的视频| 999久久久国产精品视频| 久久99蜜桃精品久久| 一二三四在线观看免费中文在| 制服丝袜香蕉在线| 老司机亚洲免费影院| 成人亚洲欧美一区二区av| av在线播放精品| 人妻一区二区av| 激情五月婷婷亚洲| 国产精品国产av在线观看| 国产亚洲欧美精品永久| 亚洲少妇的诱惑av| 亚洲伊人久久精品综合| 亚洲图色成人| 日韩一卡2卡3卡4卡2021年| 午夜精品国产一区二区电影| 免费人妻精品一区二区三区视频| 国产精品免费视频内射| 1024视频免费在线观看| 欧美成人午夜精品| 国产精品 国内视频| 亚洲精品日本国产第一区| 在线天堂最新版资源| 嫩草影院入口| 黑丝袜美女国产一区| 在线观看免费高清a一片| 精品少妇黑人巨大在线播放| 人人澡人人妻人| 18+在线观看网站| 国产精品免费大片| 久久这里有精品视频免费| 可以免费在线观看a视频的电影网站 | 国产乱人偷精品视频| 99热国产这里只有精品6| 亚洲成人一二三区av| av在线观看视频网站免费| 侵犯人妻中文字幕一二三四区| 伊人亚洲综合成人网| 在现免费观看毛片| 成人毛片a级毛片在线播放| 免费在线观看视频国产中文字幕亚洲 | 日韩成人av中文字幕在线观看| 黄片小视频在线播放| 亚洲欧洲精品一区二区精品久久久 | 亚洲av在线观看美女高潮| 26uuu在线亚洲综合色| 美女视频免费永久观看网站| 香蕉国产在线看| 激情视频va一区二区三区| 超色免费av| 三级国产精品片| freevideosex欧美| 在线观看www视频免费| 国产黄色免费在线视频| 久久久精品区二区三区| 街头女战士在线观看网站| 激情视频va一区二区三区| 欧美亚洲 丝袜 人妻 在线| 亚洲精品成人av观看孕妇| 1024视频免费在线观看| 丰满饥渴人妻一区二区三| 在线观看免费日韩欧美大片| 啦啦啦啦在线视频资源| 最黄视频免费看| 老汉色∧v一级毛片| 制服诱惑二区| 晚上一个人看的免费电影| 婷婷色综合大香蕉| 久久精品国产a三级三级三级| 亚洲欧美成人综合另类久久久| 熟妇人妻不卡中文字幕| 亚洲av男天堂| 日韩视频在线欧美| 黑人猛操日本美女一级片| 成年动漫av网址| 色哟哟·www| 国产av一区二区精品久久| 国产精品 欧美亚洲| 国产午夜精品一二区理论片| 在线观看国产h片| 亚洲视频免费观看视频| 人人妻人人添人人爽欧美一区卜| 老汉色av国产亚洲站长工具| 91精品伊人久久大香线蕉| 亚洲经典国产精华液单| 亚洲欧美清纯卡通| 大话2 男鬼变身卡| 97精品久久久久久久久久精品| 亚洲精品视频女| 一级黄片播放器| 超色免费av| 免费看av在线观看网站| 免费播放大片免费观看视频在线观看| 新久久久久国产一级毛片| 99热网站在线观看| 丁香六月天网| 国产成人精品婷婷| 国产不卡av网站在线观看| 精品第一国产精品| 国产一区二区激情短视频 | 最近最新中文字幕大全免费视频 | 国产精品 国内视频| 亚洲欧美清纯卡通| 人成视频在线观看免费观看| 人妻系列 视频| 亚洲美女黄色视频免费看| 久久精品久久久久久久性| 如日韩欧美国产精品一区二区三区| 超碰97精品在线观看| 免费观看a级毛片全部| 丝袜在线中文字幕| 久久这里只有精品19| 99久久人妻综合| av片东京热男人的天堂| 国产精品久久久av美女十八| 男女边吃奶边做爰视频| 亚洲成av片中文字幕在线观看 | 一个人免费看片子| 亚洲国产欧美在线一区| 日韩欧美一区视频在线观看| 中文字幕最新亚洲高清| 久久鲁丝午夜福利片| 精品少妇内射三级| 一级毛片黄色毛片免费观看视频| 国产精品久久久久久精品电影小说| av网站免费在线观看视频| 一级黄片播放器| 哪个播放器可以免费观看大片| 纯流量卡能插随身wifi吗| 色婷婷av一区二区三区视频| 久久精品熟女亚洲av麻豆精品| 成人影院久久| 成年人免费黄色播放视频| 王馨瑶露胸无遮挡在线观看| 男女下面插进去视频免费观看| 少妇 在线观看| 国产黄频视频在线观看| 亚洲美女搞黄在线观看| 国产男人的电影天堂91| 亚洲国产精品999| 一区二区三区四区激情视频| 欧美激情高清一区二区三区 | 两个人免费观看高清视频| 国产一区二区 视频在线| 黄色一级大片看看| 欧美 亚洲 国产 日韩一| 一级片免费观看大全| 精品国产超薄肉色丝袜足j| 久久久久精品久久久久真实原创| 一个人免费看片子| 国产视频首页在线观看| 五月开心婷婷网| 春色校园在线视频观看| 久久免费观看电影| 国产精品久久久久久久久免| 日韩一卡2卡3卡4卡2021年| 日韩视频在线欧美| 国产一区二区 视频在线| 女人久久www免费人成看片| 亚洲欧美清纯卡通| 这个男人来自地球电影免费观看 | 看免费成人av毛片| 免费在线观看黄色视频的| 亚洲成人一二三区av| 亚洲美女视频黄频| 男女无遮挡免费网站观看| 在线观看国产h片| 国产成人精品福利久久| 99久国产av精品国产电影| 亚洲精品自拍成人| 丝袜在线中文字幕| 亚洲一级一片aⅴ在线观看| 秋霞伦理黄片| 欧美日韩视频高清一区二区三区二| 伦理电影大哥的女人| 一本久久精品| 天天操日日干夜夜撸| a 毛片基地| 国产福利在线免费观看视频| 春色校园在线视频观看| 免费观看av网站的网址| 国产在视频线精品| 一级毛片黄色毛片免费观看视频| 天天躁日日躁夜夜躁夜夜| 老司机影院毛片| 亚洲成人手机| 欧美日韩综合久久久久久| 一级毛片我不卡| 精品国产一区二区三区久久久樱花| 啦啦啦视频在线资源免费观看| 人妻一区二区av| 久久人人爽av亚洲精品天堂| 成年女人在线观看亚洲视频| 午夜激情久久久久久久| 亚洲国产毛片av蜜桃av| 18禁裸乳无遮挡动漫免费视频| 精品亚洲成a人片在线观看| 亚洲av.av天堂| 午夜福利视频在线观看免费| 中国三级夫妇交换| 十八禁高潮呻吟视频| 日韩一卡2卡3卡4卡2021年| 亚洲美女搞黄在线观看| 久久久精品国产亚洲av高清涩受| 国产精品蜜桃在线观看| 中文字幕人妻熟女乱码| 黄色 视频免费看| 欧美日韩精品成人综合77777| 另类精品久久| 国产黄频视频在线观看| 在线观看www视频免费| 欧美精品一区二区免费开放| 啦啦啦啦在线视频资源| 国产一区二区三区av在线| 超碰97精品在线观看| 亚洲天堂av无毛| 成年av动漫网址| 精品国产乱码久久久久久小说| 国产精品久久久av美女十八| 丰满乱子伦码专区| www.av在线官网国产| 欧美国产精品va在线观看不卡| 永久免费av网站大全| 久久久久久免费高清国产稀缺| 电影成人av| 国产精品一二三区在线看| 性少妇av在线| 宅男免费午夜| 亚洲综合精品二区| 中文天堂在线官网| 亚洲欧美日韩另类电影网站| 国产黄色免费在线视频| 韩国av在线不卡| 两性夫妻黄色片| 国产一级毛片在线| 青草久久国产| 一区二区三区激情视频| 高清黄色对白视频在线免费看| 亚洲欧美清纯卡通| 久久精品国产亚洲av高清一级| 中文字幕精品免费在线观看视频| 免费女性裸体啪啪无遮挡网站| 欧美 亚洲 国产 日韩一| 边亲边吃奶的免费视频| 少妇 在线观看| 国产黄频视频在线观看| 一边摸一边做爽爽视频免费| 免费高清在线观看视频在线观看| 久久精品久久久久久久性| 少妇被粗大猛烈的视频| 国产爽快片一区二区三区| 女人精品久久久久毛片| 午夜影院在线不卡| 免费av中文字幕在线| 黄片小视频在线播放| 国产av国产精品国产| 色网站视频免费| 老司机影院成人| 18在线观看网站| 免费看av在线观看网站| 美女主播在线视频| 九九爱精品视频在线观看| 女人精品久久久久毛片| 女人高潮潮喷娇喘18禁视频| 午夜激情久久久久久久| 热99久久久久精品小说推荐| 好男人视频免费观看在线| 男女下面插进去视频免费观看| 国产精品蜜桃在线观看| 黑人猛操日本美女一级片| 国产av精品麻豆| 成人二区视频| 国产精品三级大全| 成人亚洲欧美一区二区av| 热re99久久精品国产66热6| 爱豆传媒免费全集在线观看| 国产在线视频一区二区| 亚洲欧美一区二区三区国产| 美女国产视频在线观看| 国产精品国产av在线观看| 国产在视频线精品| 日韩中字成人| 热99久久久久精品小说推荐| 夜夜骑夜夜射夜夜干| 午夜日本视频在线| 国产亚洲精品第一综合不卡| 我的亚洲天堂| 精品人妻偷拍中文字幕| 国产野战对白在线观看| 久久精品久久精品一区二区三区| 日韩人妻精品一区2区三区| 精品人妻一区二区三区麻豆| 午夜免费鲁丝| 天天躁夜夜躁狠狠久久av| 国产亚洲精品第一综合不卡| 久久99热这里只频精品6学生| 天天躁夜夜躁狠狠久久av| 国产野战对白在线观看| 亚洲 欧美一区二区三区| 男女下面插进去视频免费观看| av线在线观看网站| 欧美日韩视频高清一区二区三区二| 精品福利永久在线观看| 亚洲av在线观看美女高潮| 我的亚洲天堂| 日韩av不卡免费在线播放| 亚洲国产毛片av蜜桃av| 午夜福利视频精品| 一区二区三区精品91| 精品视频人人做人人爽| 一级毛片我不卡| av福利片在线| 黄片小视频在线播放| 亚洲国产色片| av免费在线看不卡|