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

    基于磁偶極子構(gòu)造法幾何建模的鐵管磁異常正演

    2017-03-27 06:23:48朱慧慧劉得軍閆景富
    計(jì)算機(jī)測量與控制 2017年3期
    關(guān)鍵詞:磁偶極子磁化率磁化

    朱慧慧,劉得軍,馮 碩,潘 琦,閆景富

    (中國石油大學(xué)(北京) 地球物理與信息工程學(xué)院,北京 102249)

    基于磁偶極子構(gòu)造法幾何建模的鐵管磁異常正演

    朱慧慧,劉得軍,馮 碩,潘 琦,閆景富

    (中國石油大學(xué)(北京) 地球物理與信息工程學(xué)院,北京 102249)

    地下鐵磁性金屬管線受地球原磁場磁化產(chǎn)生磁場會(huì)改變原地磁場分布,而形成地磁異常;通過管道幾何建模和磁偶極子構(gòu)造法對管道磁異常進(jìn)行正演分析;首先利用近似的立方體單元對管道進(jìn)行幾何建模并輸出相應(yīng)的立方體單元體積和中心坐標(biāo);然后將立方體單元近似為磁偶極子,使用磁偶極子構(gòu)造法進(jìn)行管道磁異常正演;最后研究了距徑比,磁傾角,磁偏角,地球場強(qiáng),管道材料磁化率對管道磁異常正演的影響;研究結(jié)果表明:在6倍距徑比外,同一球面上測點(diǎn)磁異常值接近一致, 且管道磁異常正演結(jié)果不受磁傾角,磁偏角和材料磁化率的影響;隨著地磁場強(qiáng)度增大,管道磁異常呈線性增大;研究成果將為自主研發(fā)高精度磁異常地下管線探測設(shè)備奠定理論基礎(chǔ)。

    鐵質(zhì)管線,幾何建模,磁偶極子構(gòu)造法,磁異常正演

    0 引言

    自然界的磁體都可以認(rèn)為是大量形狀和磁矩大小不同的基本磁體單元的組合[1],當(dāng)磁體單元最大外形尺寸遠(yuǎn)小于磁體到磁場測點(diǎn)之間的距離時(shí),基本磁體單元可視作磁偶極子來計(jì)算磁場[2-3]。磁偶極子構(gòu)造法基于該原理,將被地磁場磁化的鐵質(zhì)管道分割成有限數(shù)量的單元塊,由單元塊在幾何結(jié)構(gòu)上重構(gòu)管道[4-7]。將每個(gè)單元塊視作磁偶極子,使用磁偶極子公式計(jì)算所有單元在測點(diǎn)的磁場,依據(jù)磁場矢量疊加原理得到磁化鐵質(zhì)管道在測點(diǎn)的磁場[4-6]。

    對于形狀和分布規(guī)則的球體、無限長水平圓柱體及板狀體等磁源,其磁異常計(jì)算存在解析解[8-9],而管道、鐵皮桶、鋼鐵構(gòu)件及未爆炸武器(UXO)鐵殼等形狀及分布不規(guī)則的鐵磁質(zhì)磁源,磁異常解析解求解較為困難[9-10],使用數(shù)值計(jì)算的方式進(jìn)行磁異常仿真較為合適[9]。Marchetti等人和Furness分別使用表面積分方程法求解了地下鐵皮桶所致探測面磁異常[10-11]。Aydin使用棱柱單元計(jì)算了磁基底的磁異常[12]。Churchill等人使用COMSOL有限單元分析軟件計(jì)算了不同形態(tài)下UXO的磁異常,并將求解結(jié)果與解析法橢球體模型計(jì)算結(jié)果進(jìn)行了對比[13]。Sanchez等人及其所在研究團(tuán)隊(duì)使用立方磁體單元構(gòu)造了UXO及Non-UXO構(gòu)件磁異常數(shù)值分析模型[9]。由于地下管道結(jié)構(gòu)的特殊性,棱柱單元不能有效的重建管道幾何結(jié)構(gòu)[9,12];而有限元軟件在偽二度體磁異常探測建模時(shí),計(jì)算內(nèi)存需求大且求解易發(fā)散[13]。為此,本文利用近似立方體單元構(gòu)造管道幾何模型,使用磁偶極子構(gòu)造法建立管道磁異常數(shù)值分析模型,分析了不同因素對管道正演磁異常的影響,為管道探測磁異常數(shù)據(jù)處理及反演提供理論基礎(chǔ)。

    1 原理方法

    1.1 磁偶極子構(gòu)造法原理

    當(dāng)磁體單元的外形尺寸遠(yuǎn)小于磁體中心到磁場測量點(diǎn)的距離時(shí),磁體在測點(diǎn)處的磁場可使用磁偶極子公式(1)近似計(jì)算[2,14]。定義測點(diǎn)到磁體中心的距離與磁體最大外形尺寸之間的比為距徑比,Zhang等[14]人在一定距徑比范圍內(nèi)進(jìn)行了磁體磁偶極子磁場計(jì)算測點(diǎn)磁場誤差的實(shí)驗(yàn)分析,從實(shí)驗(yàn)上證明了磁偶極子磁場計(jì)算公式的適用性。Yang等[3]人在論文中給出了磁偶極子公式計(jì)算磁體磁場的誤差計(jì)算公式,但未說明其具體來源。Guo等[2]人基于圓環(huán)電流磁體單元定量分析了磁偶極子計(jì)算磁體的誤差,并與前人研究成果進(jìn)行了對比分析說明了其研究結(jié)果的準(zhǔn)確性。在此,我們將管道劃分所得單元視作磁偶極子,其測點(diǎn)處的磁場計(jì)算公式為:

    (1)

    (2)

    位于地磁場中的鐵質(zhì)管線分割單元在地磁場的磁化下產(chǎn)生了磁化磁矩,每一個(gè)磁化單元能夠視作一個(gè)磁偶極子[2,14]。所有磁化單元的在測點(diǎn)處的磁場疊加,即為地下磁化管道在測點(diǎn)處產(chǎn)生的磁化磁場。管道磁異常探測過程中,Pj點(diǎn)的探測磁場Bdj為磁化磁場與地球背景場之和,即:

    Bdj=Bj+Be

    (3)

    其中:Bj為管線被地磁場磁化后產(chǎn)生的磁化磁場,Be為地球背景場。

    1.2 管道幾何建模及磁異常正演計(jì)算

    如圖1所示,管道外直徑為,管道厚度為δ,在地磁觀測坐標(biāo)系中管道軸線距離地面為0.4m。為了使得磁化管道單元在距離單元中心各個(gè)方向全球面上測點(diǎn)磁場盡量一致,將管道單元分割為盡量接近于立方體的單元。沿管道軸向?qū)⒐艿赖确譃棣拈L的管道環(huán),將每一小段管道圓環(huán)等分為N等份,N為的舍入整數(shù),則長度為Lp的管線劃分所得單元數(shù)為,這個(gè)單元在幾何上還原構(gòu)造了原管道結(jié)構(gòu)。

    圖1 磁偶極子構(gòu)造法管道幾何建模模型

    位于地磁場中的鐵質(zhì)管線分割單元在地磁場的磁化下產(chǎn)生了磁化磁矩,每一個(gè)磁化單元能夠視作一個(gè)磁偶極子。所有磁化單元的在測點(diǎn)處的磁場疊加,即為地下磁化管道在測點(diǎn)處產(chǎn)生的磁化磁場。在管線外徑與厚度確定時(shí),管道單元大小即被確定,圖1中單元塊表示由管道分割所得單元。管道空間位置確定時(shí),每個(gè)管道單元的幾何中心即被確定,圖1中點(diǎn)即為單元的幾何中心。此時(shí),在地磁坐標(biāo)系中,影響點(diǎn)磁化磁場的因素只有管道材料磁化率與背景磁場。結(jié)合式(2)與式(3),最終可表示為:

    式(4)就是磁偶極子構(gòu)造法計(jì)算地磁坐標(biāo)系中,受地磁場磁化的鐵質(zhì)管道(磁化率為)在點(diǎn)的探測磁異常。

    要想確定地面上一點(diǎn)的地磁場的強(qiáng)度與方向,至少要測出任意3個(gè)彼此獨(dú)立的地磁要素,稱為地磁三要素[15]。在本文研究過程中,以地球場強(qiáng)、磁傾角及磁偏角描述磁場測點(diǎn)背景場。由公式(4)可知,距徑比,地球場強(qiáng),磁傾角,磁偏角,材料磁化率對測點(diǎn)磁異常均有影響。本文重點(diǎn)研究這些因素對測點(diǎn)磁異常的影響。為了研究分析地下鐵管正演磁異常,建立如圖2所示的分析模型。假設(shè)地球場強(qiáng),磁傾角,磁偏角,材料磁化率為固定值,研究距徑比對球面上磁異常的影響。此外,進(jìn)一步研究不同距徑比條件下,在一定范圍內(nèi)改變地球場強(qiáng),磁傾角,磁偏角,材料磁化率,通過分析模型正演計(jì)算測點(diǎn)磁異常,并分析這些因素對測點(diǎn)磁異常的影響。

    圖2 管道磁異常正演計(jì)算模型

    2 數(shù)值分析

    2.1 距徑比對球面上磁異常的影響

    圖3 距徑比λ對測點(diǎn)磁異常B的影響

    圖4 地磁場強(qiáng)度對測點(diǎn)磁異常的影響

    圖5 磁傾角對測點(diǎn)磁異常的影響

    2.2 地磁場強(qiáng)度,磁傾角,磁偏角,材料磁化率對測點(diǎn)磁異常的影響

    設(shè)直線管道單元外徑60 cm,厚度6 mm,材料磁化率為200 SI,磁傾角為59.061°,磁偏角為-6.629°,背景磁場從22 000 nT到68 000 nT以1 000 nT為間隔變化,計(jì)算距徑比分別為2,4,6,8,10時(shí),取球面上第22、28、72、161點(diǎn)四組有代表性的磁異常數(shù)據(jù)如圖4所示。由圖4可知:同一測點(diǎn)上,隨著地磁場強(qiáng)度增大,磁異常呈線性增大;隨著距徑比的增大,同一測點(diǎn)在相同地磁場強(qiáng)度時(shí)磁異常之間的差異逐漸減小;當(dāng)λ≥6時(shí),同一測點(diǎn)磁異常隨地磁場強(qiáng)度增大同比例線性增大。

    設(shè)直線管道單元外徑60cm,厚度6mm,材料磁化率為200SI,背景磁場大小為54583.6nT,磁偏角為-6.629°,磁傾角從-90°到90°變化,計(jì)算距徑比分別為2,4,6,8,10時(shí),取球面上第22、28、72、161點(diǎn)四組有代表性的磁異常數(shù)據(jù)如圖5所示。由圖5可知:在固定距徑比球面上,同一測點(diǎn)磁異常隨磁傾角的改變成正弦波動(dòng);在同一測點(diǎn)上,不同距徑比條件下,磁異常隨磁傾角變化趨勢相同;隨著距徑比的增大,同一測點(diǎn)磁異常受磁傾角的影響逐漸減??;當(dāng)λ≥6時(shí),同一測點(diǎn)磁異常不受磁傾角的影響。

    設(shè)直線管道單元外徑60cm,厚度6mm,材料磁化率為200SI,背景磁場大小為54583.6nT,磁傾角為59.061°,磁偏角從-30°到30°變化,計(jì)算距徑比分別為2,4,6,8,10時(shí),取球面上第22、28、72、161點(diǎn)四組有代表性的磁異常數(shù)據(jù)如圖6所示。由圖6可知:同一測點(diǎn)上,不同距徑比條件下,磁異常隨磁偏角的變化趨勢相同;隨著距徑比的增大,同一測點(diǎn)磁異常受磁偏角的影響逐漸減??;當(dāng)λ≥6時(shí),同一測點(diǎn)磁異常不受磁偏角的影響。

    圖6 磁偏角對測點(diǎn)磁異常的影響

    設(shè)直線管道單元外徑60cm,厚度6mm,背景磁場大小為54583.6nT,磁傾角為59.061°,磁偏角為-6.629°,材料磁化率從100SI到1000SI以50SI為間隔變化,計(jì)算距徑比分別為2,4,6,8,10時(shí),取球面上第22、28、72、161點(diǎn)四組有代表性的磁異常數(shù)據(jù)如圖7所示。由圖7可知:同一測點(diǎn)上,不同距徑比條件下測點(diǎn)磁異常隨材料磁化率的變化趨勢相同;隨著距徑比的增大,同一測點(diǎn)磁異常受材料磁化率影響逐漸變?。划?dāng)λ≥6時(shí),同一測點(diǎn)磁異常不受材料磁化率的影響。

    圖7 材料磁化率對測點(diǎn)磁異常的影響

    3 結(jié)論

    利用近似立方體單元重構(gòu)地下鐵質(zhì)管線的幾何結(jié)構(gòu),在此基礎(chǔ)上使用磁偶極子構(gòu)造法建立管道磁異常正演模型。研究了測點(diǎn)空間位置對磁異常的影響,同時(shí)討論了磁傾角,磁偏角,地球場強(qiáng),材料磁化率對磁異常的影響。由研究分析可知,磁傾角,磁偏角,地球場強(qiáng),材料磁化率一定時(shí),在不同的距徑比球面上,管道磁異常分布趨勢相似;隨著距徑比的增大,同一球面上測點(diǎn)的磁異常變化逐漸變??;當(dāng)距徑比大于6時(shí),不同測點(diǎn)處磁異常近似相同。隨著地球場強(qiáng)的增大,測點(diǎn)磁異常線性增大。當(dāng)距徑比大于6時(shí),測點(diǎn)磁異常不受磁傾角,磁偏角,材料磁化率的影響。文章提供了一種有效的地下鐵管磁異常正演方法,為自主研發(fā)高精度磁異常地下管線探測設(shè)備奠定理論基礎(chǔ),同時(shí)為實(shí)際儀器的標(biāo)定、儀器響應(yīng)模擬及探測特性評(píng)價(jià)等提供理論依據(jù)。

    [1] 張恒磊,胡祥云,劉天佑. 基于二階導(dǎo)數(shù)的磁源邊界與頂部深度快速反演[J]. 地球物理學(xué)報(bào). 2012, 55(11):3839-3847.

    [2] 底青云,方廣有,張一鳴. 地面電磁探測系統(tǒng)(SEP)研究[J]. 地球物理學(xué)報(bào). 2013, 56(11): 3629-3639.

    [3] 張恒磊,MarangoniYR,左仁廣,等. 改進(jìn)的各向異性標(biāo)準(zhǔn)化方差探測斜磁化磁異常源邊界[J]. 地球物理學(xué)報(bào),2014, 57(8): 2724-2731.

    [4]YangWA,HuC,LiM,etal.Anewtrackingsystemforthreemagneticobjectives[J].IEEETransactionsonMagnetics, 2010, 46(12): 4023-4029.

    [5] 張朝陽, 肖昌漢, 高俊吉,等. 磁性物體磁偶極子模型適用性的試驗(yàn)研究[J]. 應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào), 2010, 18(5): 862-868.

    [6]HuC,MengMQandMandalM.Efficientmagneticlocalizationandorientationtechniqueforcapsuleendoscopy[J].InternationalJournalofInformationAcquisition, 2005, 2(01): 23-36.

    [7] 郭文波,李琳,薛國強(qiáng)等. 瞬變電磁快速成像解釋系統(tǒng)研究[J]. 地球物理學(xué)報(bào),2005, 48(6): 1400-1405.

    [8]ShawRK,AgarwalBNP,NandiBK.UseofWalshtransformsinestimationofdepthsofidealizedsourcesfromtotal-fieldmagneticanomalies[J].Computers&geosciences, 2007, 33(7): 966-975.

    [9]SanchezV,LiYG,NabighianMNetal.NumericalModelingofHigherOrderMagneticMomentsinUXODiscrimination[J].IEEETransactionsonGeoscienceandRemoteSensing, 2008, 46(9): 2568-2583.

    [10]MarchettiM,SapiaVandSettimiA.Magneticanomaliesofsteeldrums:areviewoftheliteratureandresearchresultsoftheINGV[J].AnnalsofGeophysics, 2013, 56(1):R0108.

    [11]FurnessP.Modellingmagneticfieldsduetosteeldrumaccumulations[J].GeophysicalProspecting, 2007, 55: 1-12.

    [12]AydinI,OksumE.MATLABcodeforestimatingmagneticbasementdepthusingprisms[J].Computers&Geosciences, 2012, 46: 183-188.

    [13]ChurchillKM,LinkC,YoumansCC.Acomparisonofthefinite-elementmethodandanalyticalmethodformodelingunexplodedordnanceusingmagnetometry[J].IEEETransactionsonGeoscienceandRemoteSensing, 2012, 50(7): 2720-2732.

    [14]JessellM.Three-dimensionalgeologicalmodellingofpotential-fielddata[J].Computers&Geosciences, 2001, 27(4): 455-465.

    [15] 徐文耀. 地球電磁形象物理學(xué)[M]. 合肥:中國科學(xué)技術(shù)大學(xué)出版社, 2009.

    [16] 洪立波. 城市地下管線面臨的挑戰(zhàn)與機(jī)遇[J]. 城市運(yùn)行管理, 2009(2): 32-34.

    [17] 雷林源. 城市地下管線探測與測漏[M]. 北京:冶金工業(yè)出版社, 2003.

    Magnetic Dipole Reconstruction Geometry Modeling for Underground Ferromagnetic Pipe Magnetic Abnormal Detection

    Zhu Huihui, Liu Dejun, Feng Shuo, Pan Qi, Yan Jingfu

    (College of Geophysics and Information Engineering, China University of Petroleum, Beijing 102249, China)

    Underground ferromagnetic pipe can be magneticed by earth magnetic field, superposition of magneticed field and earth magnetic field produce magnetic abnormal. By pipe geometry modeling and Magnetic dipole reconstruction we can do Magnetic abnormal forward analysis. Initially, approximate cube elements are used to reconstruct pipe geometric model and output volume and center coordinates of elements. Then cube elements are seen as magnetic dipole and magnetic dipole reconstruction is used to forward pipe magnetic abnormal. At last research the influence of Distance-to-Diameter ratio, magnetic inclination, magnetic declination, geomagnetic field and pipeline’s susceptibility on pipe magnetic abnormal forward calculation. The results show that: when Distance-to-Diameter ratio is bigger than 6, magnetic abnormal of measurement points on the same sphere are very close and pipe magnetic abnormal is not affected by magnetic inclination, magnetic declination, and pipeline’s susceptibility; pipe magnetic abnormal increases linearly with the increase of geomagnetic field. This result can lay a theory foundation for high precision underground pipeline detection equipment.

    underground ferromagnetic pipeline; geometry modeling; magnetic dipole reconstruction; magnetic abnormal forward

    2016-10-17;

    2016-11-21。

    國家自然科學(xué)基金(41374151)。

    朱慧慧(1989-),女,湖北鐘祥人,碩士研究生,主要從事地下鐵質(zhì)管線磁異常正演方法方向的研究。

    劉得軍(1965-),男,河北遷安人,教授、博士生導(dǎo)師,主要從事磁異常探測與儀器響應(yīng)數(shù)值模擬方向的研究。

    1671-4598(2017)03-0201-04

    10.16526/j.cnki.11-4762/tp.2017.03.055

    TP631

    A

    猜你喜歡
    磁偶極子磁化率磁化
    東北豐磁化炭基復(fù)合肥
    雙色球磁化炭基復(fù)合肥
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    磁偶極子跟蹤的漸進(jìn)貝葉斯濾波方法
    基于超拉普拉斯分布的磁化率重建算法
    基于磁化能量的鋰電池串模塊化均衡方法
    基于磁偶極子的磁場梯度張量縮并的試驗(yàn)驗(yàn)證及相關(guān)參數(shù)確定
    巖(礦)石標(biāo)本磁化率測定方法試驗(yàn)及認(rèn)識(shí)
    水平磁偶極子電磁場特征研究
    超強(qiáng)磁場下簡并電子氣體的磁化
    精品国内亚洲2022精品成人| 高清日韩中文字幕在线| 国产精品麻豆人妻色哟哟久久 | 婷婷色综合大香蕉| 九九爱精品视频在线观看| 免费电影在线观看免费观看| 国产精品一二三区在线看| 国产熟女欧美一区二区| 日本黄色片子视频| 毛片女人毛片| 淫秽高清视频在线观看| 日产精品乱码卡一卡2卡三| 亚洲精品aⅴ在线观看| 可以在线观看毛片的网站| 亚洲国产精品久久男人天堂| 欧美zozozo另类| 91久久精品国产一区二区三区| 日本wwww免费看| 两个人的视频大全免费| 欧美变态另类bdsm刘玥| 成人鲁丝片一二三区免费| 乱人视频在线观看| 国内揄拍国产精品人妻在线| 精品国产一区二区三区久久久樱花 | 国产亚洲精品久久久com| 精品国产三级普通话版| 永久免费av网站大全| 日本熟妇午夜| 老司机福利观看| 日韩 亚洲 欧美在线| av免费观看日本| 亚洲人成网站高清观看| 有码 亚洲区| 啦啦啦韩国在线观看视频| 亚洲欧洲日产国产| 在线观看66精品国产| 午夜精品在线福利| 26uuu在线亚洲综合色| 国产又色又爽无遮挡免| 午夜福利成人在线免费观看| 深爱激情五月婷婷| 寂寞人妻少妇视频99o| 亚洲av中文字字幕乱码综合| 欧美成人精品欧美一级黄| 日韩一区二区三区影片| 日本一二三区视频观看| 亚洲五月天丁香| 国产私拍福利视频在线观看| 中文天堂在线官网| 亚洲精华国产精华液的使用体验| 性色avwww在线观看| 亚洲怡红院男人天堂| 欧美日韩在线观看h| 三级男女做爰猛烈吃奶摸视频| 直男gayav资源| 一个人免费在线观看电影| 热99re8久久精品国产| h日本视频在线播放| 毛片女人毛片| www.色视频.com| 精品熟女少妇av免费看| 亚洲不卡免费看| 国产av不卡久久| 国产极品天堂在线| 欧美人与善性xxx| 国产伦在线观看视频一区| 中文在线观看免费www的网站| 亚洲一级一片aⅴ在线观看| 久久鲁丝午夜福利片| 内地一区二区视频在线| 色哟哟·www| 人人妻人人澡欧美一区二区| 中国美白少妇内射xxxbb| 国产又黄又爽又无遮挡在线| 色视频www国产| 性色avwww在线观看| 国产亚洲午夜精品一区二区久久 | 中国国产av一级| 亚洲av成人精品一二三区| 大话2 男鬼变身卡| 欧美日本视频| 国产精品野战在线观看| 六月丁香七月| 菩萨蛮人人尽说江南好唐韦庄 | 美女被艹到高潮喷水动态| 免费看日本二区| 三级男女做爰猛烈吃奶摸视频| 国产精品永久免费网站| 久久精品影院6| 日本黄大片高清| av福利片在线观看| 亚洲国产精品国产精品| 能在线免费看毛片的网站| 色吧在线观看| 午夜日本视频在线| 久久久久免费精品人妻一区二区| 亚洲精品乱久久久久久| 在线观看66精品国产| 韩国av在线不卡| 一边摸一边抽搐一进一小说| 爱豆传媒免费全集在线观看| 一区二区三区高清视频在线| 国产免费又黄又爽又色| 亚洲欧美中文字幕日韩二区| 日本wwww免费看| 尤物成人国产欧美一区二区三区| 桃色一区二区三区在线观看| 国产一区二区在线av高清观看| 伊人久久精品亚洲午夜| 国国产精品蜜臀av免费| 国产成人精品一,二区| 乱人视频在线观看| 三级国产精品欧美在线观看| 伦精品一区二区三区| 国产精品不卡视频一区二区| 少妇的逼水好多| 久久久a久久爽久久v久久| 在线播放国产精品三级| 国产在视频线在精品| 欧美另类亚洲清纯唯美| 一级黄片播放器| 99在线视频只有这里精品首页| 性插视频无遮挡在线免费观看| 嫩草影院新地址| 国产一区二区亚洲精品在线观看| 七月丁香在线播放| 亚洲综合色惰| 91午夜精品亚洲一区二区三区| 久久久久久久久久成人| 亚洲精品国产av成人精品| 一级毛片电影观看 | 尾随美女入室| 少妇人妻精品综合一区二区| 我的女老师完整版在线观看| av播播在线观看一区| 亚洲高清免费不卡视频| 精品免费久久久久久久清纯| 我要搜黄色片| 91精品伊人久久大香线蕉| 女人被狂操c到高潮| 美女被艹到高潮喷水动态| 99久久九九国产精品国产免费| 午夜日本视频在线| 高清日韩中文字幕在线| av又黄又爽大尺度在线免费看 | 最近最新中文字幕免费大全7| 欧美日韩精品成人综合77777| 亚洲伊人久久精品综合 | 亚洲综合精品二区| 男人的好看免费观看在线视频| 18禁裸乳无遮挡免费网站照片| 午夜免费激情av| 丰满少妇做爰视频| 国语对白做爰xxxⅹ性视频网站| 久久国内精品自在自线图片| 男女视频在线观看网站免费| 亚洲欧洲国产日韩| 国产不卡一卡二| 日本五十路高清| 亚洲欧美精品自产自拍| 黄色配什么色好看| av在线亚洲专区| 黄色配什么色好看| 国产成人精品久久久久久| 午夜福利在线观看免费完整高清在| 免费无遮挡裸体视频| 一级黄色大片毛片| 亚洲精品自拍成人| 国产三级中文精品| 亚洲成av人片在线播放无| 久久久久久久国产电影| 美女cb高潮喷水在线观看| 欧美性感艳星| av免费在线看不卡| 人体艺术视频欧美日本| www.av在线官网国产| 国产v大片淫在线免费观看| 91久久精品国产一区二区成人| 97超视频在线观看视频| 亚洲精品日韩av片在线观看| 国内精品美女久久久久久| 成人亚洲欧美一区二区av| 熟女电影av网| 日日摸夜夜添夜夜添av毛片| 少妇的逼水好多| 欧美人与善性xxx| 亚洲色图av天堂| 国产免费福利视频在线观看| 亚洲激情五月婷婷啪啪| 亚洲精品久久久久久婷婷小说 | 国产精品久久久久久精品电影| 免费观看a级毛片全部| 精品人妻偷拍中文字幕| 国产极品精品免费视频能看的| 村上凉子中文字幕在线| 麻豆国产97在线/欧美| 国产成人精品婷婷| 精华霜和精华液先用哪个| 国产私拍福利视频在线观看| 乱人视频在线观看| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩高清专用| 日本欧美国产在线视频| 亚洲伊人久久精品综合 | 一个人看视频在线观看www免费| 亚洲图色成人| 国产69精品久久久久777片| 久久久久网色| 一本一本综合久久| 国产单亲对白刺激| 中文字幕久久专区| 成人美女网站在线观看视频| 国产又色又爽无遮挡免| 日本熟妇午夜| 欧美日韩一区二区视频在线观看视频在线 | 最后的刺客免费高清国语| 七月丁香在线播放| 久久久久久伊人网av| 搡女人真爽免费视频火全软件| 五月伊人婷婷丁香| 国产一区二区三区av在线| 亚洲乱码一区二区免费版| 日韩av不卡免费在线播放| 三级国产精品欧美在线观看| 亚洲激情五月婷婷啪啪| 亚洲国产精品sss在线观看| 波多野结衣高清无吗| 国产不卡一卡二| 最近最新中文字幕免费大全7| 国产高清国产精品国产三级 | 亚洲美女视频黄频| 亚洲怡红院男人天堂| 国产国拍精品亚洲av在线观看| 亚洲精品国产成人久久av| 欧美性猛交黑人性爽| 亚洲精品影视一区二区三区av| 国内精品宾馆在线| 欧美一区二区国产精品久久精品| www日本黄色视频网| 久久99蜜桃精品久久| 深夜a级毛片| 熟女电影av网| 亚洲电影在线观看av| 联通29元200g的流量卡| 亚洲人成网站在线播| 禁无遮挡网站| 亚洲欧美清纯卡通| 亚洲精品国产成人久久av| 中文精品一卡2卡3卡4更新| 麻豆国产97在线/欧美| 美女大奶头视频| 久久精品国产自在天天线| 七月丁香在线播放| 久久婷婷人人爽人人干人人爱| 日韩欧美精品v在线| 亚洲中文字幕日韩| 欧美三级亚洲精品| 最近中文字幕2019免费版| 一个人看的www免费观看视频| 亚洲性久久影院| 麻豆一二三区av精品| 男人和女人高潮做爰伦理| 久久精品综合一区二区三区| 日本-黄色视频高清免费观看| 久久久欧美国产精品| 亚洲在线观看片| 天堂影院成人在线观看| 中文在线观看免费www的网站| 久久久久久久亚洲中文字幕| 直男gayav资源| videos熟女内射| 老女人水多毛片| 国产精品电影一区二区三区| 免费看光身美女| 日韩av不卡免费在线播放| 日本免费一区二区三区高清不卡| 尾随美女入室| 插逼视频在线观看| 日日啪夜夜撸| 爱豆传媒免费全集在线观看| 综合色av麻豆| 国产免费视频播放在线视频 | 91av网一区二区| 国产男人的电影天堂91| 欧美性猛交╳xxx乱大交人| 亚洲国产色片| 九九爱精品视频在线观看| 免费av观看视频| 一级黄色大片毛片| 精品不卡国产一区二区三区| av福利片在线观看| 免费看a级黄色片| 国产精品一及| 最近视频中文字幕2019在线8| 欧美成人a在线观看| 三级毛片av免费| 蜜桃久久精品国产亚洲av| 日韩成人伦理影院| 一级黄色大片毛片| 亚洲自拍偷在线| 欧美精品国产亚洲| 嫩草影院入口| 久久综合国产亚洲精品| 日本爱情动作片www.在线观看| 毛片女人毛片| 久热久热在线精品观看| 久久久久久久久久久丰满| 日韩一区二区视频免费看| 99久久九九国产精品国产免费| 夫妻性生交免费视频一级片| 日本-黄色视频高清免费观看| 日本与韩国留学比较| 中文字幕熟女人妻在线| 狂野欧美白嫩少妇大欣赏| 大话2 男鬼变身卡| 日本色播在线视频| 在线观看一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲18禁久久av| 狠狠狠狠99中文字幕| 国产真实乱freesex| 日韩一区二区三区影片| 欧美成人精品欧美一级黄| 内地一区二区视频在线| 欧美极品一区二区三区四区| 在线观看人妻少妇| 午夜免费男女啪啪视频观看| 久久久久精品人妻al黑| 国产毛片在线视频| 成人亚洲精品一区在线观看| 中文字幕人妻丝袜制服| 一区二区三区精品91| 久久这里只有精品19| 国产av国产精品国产| 国产在视频线精品| 99热全是精品| 中文字幕免费在线视频6| 亚洲精品一二三| 熟女电影av网| 亚洲av.av天堂| 最近中文字幕高清免费大全6| 人妻系列 视频| 亚洲,欧美精品.| 国产精品 国内视频| 夫妻午夜视频| 51国产日韩欧美| 午夜免费观看性视频| 色94色欧美一区二区| 男女下面插进去视频免费观看 | 国产精品一区二区在线不卡| 春色校园在线视频观看| 久久精品国产亚洲av天美| 日本欧美视频一区| 久久精品国产亚洲av涩爱| 久久久亚洲精品成人影院| 国产男女内射视频| 免费观看性生交大片5| 人妻人人澡人人爽人人| 国产精品久久久久久精品电影小说| 欧美丝袜亚洲另类| 在线亚洲精品国产二区图片欧美| 丝袜喷水一区| 久久精品久久久久久久性| 下体分泌物呈黄色| av不卡在线播放| 韩国高清视频一区二区三区| 高清毛片免费看| 宅男免费午夜| 久久久久久久久久人人人人人人| 天堂中文最新版在线下载| 久久99蜜桃精品久久| 久久热在线av| 黄色视频在线播放观看不卡| 午夜免费男女啪啪视频观看| 大话2 男鬼变身卡| 国产国拍精品亚洲av在线观看| 国产精品久久久久久av不卡| 激情视频va一区二区三区| 精品一品国产午夜福利视频| 伊人亚洲综合成人网| 欧美精品av麻豆av| 天堂俺去俺来也www色官网| 中国美白少妇内射xxxbb| 免费人妻精品一区二区三区视频| 看免费av毛片| 欧美成人午夜免费资源| 久久女婷五月综合色啪小说| 亚洲av男天堂| www日本在线高清视频| 国产精品人妻久久久久久| 成人黄色视频免费在线看| 一本大道久久a久久精品| 国产毛片在线视频| 国产精品一国产av| 国精品久久久久久国模美| 国产亚洲av片在线观看秒播厂| 国产有黄有色有爽视频| 成人国产av品久久久| 99视频精品全部免费 在线| 黄色一级大片看看| 久久99蜜桃精品久久| 久久热在线av| 成年人午夜在线观看视频| 亚洲精品第二区| 国产日韩欧美在线精品| 大陆偷拍与自拍| 制服诱惑二区| 国产精品免费大片| 91在线精品国自产拍蜜月| 建设人人有责人人尽责人人享有的| 九草在线视频观看| 日韩中文字幕视频在线看片| 青春草国产在线视频| 国产探花极品一区二区| 亚洲激情五月婷婷啪啪| 色哟哟·www| 国产黄色视频一区二区在线观看| 亚洲精华国产精华液的使用体验| av播播在线观看一区| 亚洲 欧美一区二区三区| 欧美成人午夜免费资源| 最近中文字幕2019免费版| 一本色道久久久久久精品综合| 又黄又粗又硬又大视频| 国产一区有黄有色的免费视频| 欧美精品一区二区大全| 中文字幕最新亚洲高清| 欧美成人午夜精品| 免费在线观看完整版高清| 国产精品无大码| 在线观看免费高清a一片| 亚洲精品aⅴ在线观看| 18禁裸乳无遮挡动漫免费视频| 9191精品国产免费久久| 妹子高潮喷水视频| 国产免费现黄频在线看| 国产一区二区在线观看av| 春色校园在线视频观看| 午夜影院在线不卡| 一级片免费观看大全| 91久久精品国产一区二区三区| 久久人人爽人人片av| 日韩精品免费视频一区二区三区 | 男的添女的下面高潮视频| 自拍欧美九色日韩亚洲蝌蚪91| 精品国产国语对白av| 欧美成人午夜免费资源| 你懂的网址亚洲精品在线观看| 大码成人一级视频| 9热在线视频观看99| 久热这里只有精品99| 国产永久视频网站| 国产日韩一区二区三区精品不卡| 天堂俺去俺来也www色官网| 亚洲,欧美,日韩| 国产精品无大码| 国产1区2区3区精品| 人人妻人人添人人爽欧美一区卜| av不卡在线播放| 国产不卡av网站在线观看| www.av在线官网国产| 最近最新中文字幕免费大全7| 高清欧美精品videossex| 日韩制服丝袜自拍偷拍| 女人久久www免费人成看片| 王馨瑶露胸无遮挡在线观看| 国产又爽黄色视频| 久久久久久久久久久久大奶| 波多野结衣一区麻豆| 亚洲三级黄色毛片| 国产欧美日韩一区二区三区在线| 亚洲成色77777| 国产熟女午夜一区二区三区| 国国产精品蜜臀av免费| 春色校园在线视频观看| 久久精品国产鲁丝片午夜精品| 一二三四在线观看免费中文在 | 日韩制服丝袜自拍偷拍| 免费在线观看完整版高清| 久久99蜜桃精品久久| 丰满迷人的少妇在线观看| 国产精品一区二区在线不卡| 男人爽女人下面视频在线观看| 18禁观看日本| 色网站视频免费| 久久久精品94久久精品| 亚洲欧美色中文字幕在线| 欧美日韩视频精品一区| 精品一区在线观看国产| 亚洲精品成人av观看孕妇| 亚洲成av片中文字幕在线观看 | 插逼视频在线观看| 久久人人爽人人爽人人片va| 成人二区视频| 熟女电影av网| 精品少妇久久久久久888优播| 在线观看免费视频网站a站| 最新中文字幕久久久久| 99热全是精品| 久久久久久久久久人人人人人人| 一本—道久久a久久精品蜜桃钙片| 一级黄片播放器| 国产熟女欧美一区二区| 美女xxoo啪啪120秒动态图| 九九爱精品视频在线观看| 一级毛片黄色毛片免费观看视频| 搡老乐熟女国产| 2021少妇久久久久久久久久久| 精品人妻在线不人妻| 亚洲美女黄色视频免费看| 又黄又粗又硬又大视频| 伦精品一区二区三区| 日韩精品有码人妻一区| 久久久久精品性色| 亚洲国产日韩一区二区| 日韩中文字幕视频在线看片| 侵犯人妻中文字幕一二三四区| tube8黄色片| 日本与韩国留学比较| 亚洲高清免费不卡视频| 国产又爽黄色视频| 美女主播在线视频| av女优亚洲男人天堂| 欧美人与善性xxx| 99久久中文字幕三级久久日本| 国产亚洲一区二区精品| 色94色欧美一区二区| 色哟哟·www| 新久久久久国产一级毛片| 国产女主播在线喷水免费视频网站| 亚洲国产av影院在线观看| 97精品久久久久久久久久精品| 国产白丝娇喘喷水9色精品| 日韩三级伦理在线观看| 国产精品一区www在线观看| 久久精品aⅴ一区二区三区四区 | 男人添女人高潮全过程视频| 日本午夜av视频| 97精品久久久久久久久久精品| 精品国产一区二区三区久久久樱花| 国产成人av激情在线播放| 美女中出高潮动态图| 久久久久精品久久久久真实原创| 日韩成人伦理影院| 777米奇影视久久| 久久久a久久爽久久v久久| 日韩中字成人| 亚洲av综合色区一区| 18在线观看网站| 99热这里只有是精品在线观看| 亚洲精品,欧美精品| 亚洲精品中文字幕在线视频| 一区在线观看完整版| 免费黄色在线免费观看| 高清黄色对白视频在线免费看| 美国免费a级毛片| 久久女婷五月综合色啪小说| 精品久久国产蜜桃| 亚洲av综合色区一区| 22中文网久久字幕| av国产精品久久久久影院| 99热全是精品| 中文天堂在线官网| 一本大道久久a久久精品| 国产在线视频一区二区| 久久ye,这里只有精品| 亚洲少妇的诱惑av| 中文欧美无线码| 男人爽女人下面视频在线观看| 日韩,欧美,国产一区二区三区| 久久久国产一区二区| 午夜福利视频在线观看免费| 少妇人妻 视频| 在线观看www视频免费| av网站免费在线观看视频| 久久午夜福利片| 国产一区二区三区综合在线观看 | 丝袜脚勾引网站| 黑丝袜美女国产一区| 久久人人爽av亚洲精品天堂| h视频一区二区三区| 女人精品久久久久毛片| 女人被躁到高潮嗷嗷叫费观| 汤姆久久久久久久影院中文字幕| av不卡在线播放| 国产免费现黄频在线看| 啦啦啦中文免费视频观看日本| 18禁国产床啪视频网站| 一本色道久久久久久精品综合| 欧美 日韩 精品 国产| 视频在线观看一区二区三区| 精品国产国语对白av| 免费看光身美女| 亚洲欧洲日产国产| 乱码一卡2卡4卡精品| 又大又黄又爽视频免费| 亚洲精品中文字幕在线视频| 国产永久视频网站| 婷婷色综合www| 黑丝袜美女国产一区| 久久97久久精品| 国产毛片在线视频| 黑人猛操日本美女一级片| 91成人精品电影| 亚洲美女视频黄频| 久久久久久久国产电影| 欧美日韩视频精品一区| 性色av一级| 亚洲精品视频女| 水蜜桃什么品种好| 一级片免费观看大全| 99九九在线精品视频| 日韩熟女老妇一区二区性免费视频| 黄片播放在线免费| 亚洲精华国产精华液的使用体验|