• <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)磁場下簡并電子氣體的磁化
    免费无遮挡裸体视频| 欧美日本视频| 国产一区二区激情短视频| 国产三级在线视频| 国产真实乱freesex| 午夜福利18| 午夜视频国产福利| 高清日韩中文字幕在线| 丝袜美腿在线中文| 97超视频在线观看视频| 天堂av国产一区二区熟女人妻| 男人狂女人下面高潮的视频| 国产精品久久久久久人妻精品电影| 亚洲美女视频黄频| 白带黄色成豆腐渣| 国产午夜精品论理片| 色哟哟·www| 在线观看一区二区三区| 18禁黄网站禁片免费观看直播| 2021天堂中文幕一二区在线观| 在线观看舔阴道视频| 亚洲国产日韩欧美精品在线观看| 麻豆成人午夜福利视频| 美女xxoo啪啪120秒动态图 | 国产色爽女视频免费观看| 99riav亚洲国产免费| 亚洲精华国产精华精| 国内毛片毛片毛片毛片毛片| 国产av在哪里看| eeuss影院久久| 亚洲一区二区三区色噜噜| 亚洲国产精品合色在线| 午夜免费男女啪啪视频观看 | 午夜免费男女啪啪视频观看 | 能在线免费观看的黄片| 亚洲激情在线av| 人人妻人人看人人澡| 最近在线观看免费完整版| 美女免费视频网站| 内射极品少妇av片p| 狂野欧美白嫩少妇大欣赏| 精华霜和精华液先用哪个| 国产精品永久免费网站| 淫秽高清视频在线观看| 欧美激情久久久久久爽电影| 非洲黑人性xxxx精品又粗又长| 男人舔奶头视频| 欧美一级a爱片免费观看看| a在线观看视频网站| 久久婷婷人人爽人人干人人爱| 亚洲精品亚洲一区二区| 精品久久久久久久久久久久久| 桃红色精品国产亚洲av| 香蕉av资源在线| 亚洲专区中文字幕在线| 女同久久另类99精品国产91| 精品久久久久久久久久免费视频| 欧美另类亚洲清纯唯美| 久久久久性生活片| 久久性视频一级片| 黄片小视频在线播放| 色哟哟哟哟哟哟| 欧美乱色亚洲激情| 国产探花极品一区二区| 757午夜福利合集在线观看| 亚洲中文字幕日韩| 亚洲国产精品999在线| 男女床上黄色一级片免费看| 欧美日韩乱码在线| 可以在线观看的亚洲视频| 狂野欧美白嫩少妇大欣赏| 午夜视频国产福利| 一进一出抽搐gif免费好疼| 久久精品国产自在天天线| 精品久久久久久久久av| 18禁黄网站禁片午夜丰满| 精品久久久久久久久久久久久| 欧美日本视频| 亚洲av成人不卡在线观看播放网| 久久久成人免费电影| 免费人成视频x8x8入口观看| 婷婷丁香在线五月| 日日摸夜夜添夜夜添av毛片 | 成人av一区二区三区在线看| 久久精品久久久久久噜噜老黄 | 精品人妻1区二区| 精品一区二区免费观看| 99久久无色码亚洲精品果冻| av女优亚洲男人天堂| 99热精品在线国产| 亚洲人成网站在线播放欧美日韩| 热99在线观看视频| 一夜夜www| 能在线免费观看的黄片| 成人av在线播放网站| 最近最新中文字幕大全电影3| 熟妇人妻久久中文字幕3abv| 欧美日本亚洲视频在线播放| 高清日韩中文字幕在线| 狂野欧美白嫩少妇大欣赏| 真实男女啪啪啪动态图| av天堂中文字幕网| 久久亚洲精品不卡| 99国产综合亚洲精品| 麻豆av噜噜一区二区三区| 亚洲在线观看片| 国语自产精品视频在线第100页| 日本黄大片高清| 亚洲一区二区三区色噜噜| 无人区码免费观看不卡| 免费看日本二区| 国产亚洲精品久久久com| 久久久久九九精品影院| 高清毛片免费观看视频网站| 俺也久久电影网| 午夜精品在线福利| 97热精品久久久久久| 亚洲五月婷婷丁香| 欧美乱妇无乱码| 99精品久久久久人妻精品| 成人永久免费在线观看视频| 亚洲中文日韩欧美视频| 国产在线男女| 国产高清视频在线播放一区| 久久中文看片网| 舔av片在线| www.www免费av| 免费看日本二区| 人人妻人人看人人澡| 999久久久精品免费观看国产| 国产精品av视频在线免费观看| 精品久久久久久久久av| 一个人免费在线观看电影| 人妻久久中文字幕网| 国产亚洲精品综合一区在线观看| 97碰自拍视频| 成人特级av手机在线观看| 色视频www国产| 天堂av国产一区二区熟女人妻| 精品一区二区三区视频在线观看免费| 黄片小视频在线播放| 亚洲狠狠婷婷综合久久图片| 亚洲不卡免费看| 亚洲无线观看免费| 国产三级在线视频| 我要看日韩黄色一级片| 国产真实乱freesex| 精品久久久久久久久久免费视频| 精品一区二区三区人妻视频| 97超视频在线观看视频| 日韩 亚洲 欧美在线| 成人美女网站在线观看视频| 在线天堂最新版资源| 亚洲电影在线观看av| 少妇人妻一区二区三区视频| 久久这里只有精品中国| 69av精品久久久久久| 很黄的视频免费| 国产综合懂色| 国产三级在线视频| 日本a在线网址| 免费在线观看亚洲国产| 国产精品人妻久久久久久| 欧美在线黄色| 女人十人毛片免费观看3o分钟| 一级av片app| 好看av亚洲va欧美ⅴa在| 国产午夜福利久久久久久| 国产精品一及| 国产亚洲av嫩草精品影院| 直男gayav资源| 91九色精品人成在线观看| 久久亚洲精品不卡| 欧美在线一区亚洲| 中文字幕人妻熟人妻熟丝袜美| 亚洲va日本ⅴa欧美va伊人久久| 国产蜜桃级精品一区二区三区| 91久久精品国产一区二区成人| 国产精品久久久久久久久免 | 亚洲专区国产一区二区| 香蕉av资源在线| 69av精品久久久久久| 成人美女网站在线观看视频| 成人毛片a级毛片在线播放| 国产av不卡久久| 性色avwww在线观看| 18+在线观看网站| 男人舔奶头视频| 午夜福利免费观看在线| 可以在线观看的亚洲视频| 一本精品99久久精品77| 美女免费视频网站| 一区二区三区免费毛片| 国产在线精品亚洲第一网站| 欧美三级亚洲精品| 国产精品,欧美在线| 日韩高清综合在线| 人妻丰满熟妇av一区二区三区| 小说图片视频综合网站| a级毛片免费高清观看在线播放| 成人av在线播放网站| 国产欧美日韩精品一区二区| 最好的美女福利视频网| www.www免费av| 亚洲人成网站高清观看| 又粗又爽又猛毛片免费看| 亚洲精品粉嫩美女一区| 久久天躁狠狠躁夜夜2o2o| 精品人妻一区二区三区麻豆 | 一边摸一边抽搐一进一小说| 欧美日韩福利视频一区二区| 国产成人a区在线观看| 一个人看视频在线观看www免费| av福利片在线观看| 色综合亚洲欧美另类图片| 日韩高清综合在线| 日日摸夜夜添夜夜添小说| 999久久久精品免费观看国产| 波多野结衣高清作品| 99在线视频只有这里精品首页| av专区在线播放| 午夜两性在线视频| 亚洲内射少妇av| 亚洲欧美激情综合另类| 人妻制服诱惑在线中文字幕| 午夜福利欧美成人| 99riav亚洲国产免费| 日本一本二区三区精品| 婷婷亚洲欧美| 村上凉子中文字幕在线| 91麻豆精品激情在线观看国产| 国产蜜桃级精品一区二区三区| 亚洲av美国av| 最后的刺客免费高清国语| 精品人妻一区二区三区麻豆 | 最新中文字幕久久久久| 亚洲五月天丁香| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜久久久久精精品| 国产精品99久久久久久久久| 国产成人福利小说| 欧美色欧美亚洲另类二区| 在线观看舔阴道视频| 国产午夜精品久久久久久一区二区三区 | 成人毛片a级毛片在线播放| 一二三四社区在线视频社区8| 精品人妻熟女av久视频| 在线十欧美十亚洲十日本专区| 国产精品久久久久久亚洲av鲁大| 九色成人免费人妻av| 日韩 亚洲 欧美在线| 搡老熟女国产l中国老女人| 18禁黄网站禁片免费观看直播| 国产精品久久久久久人妻精品电影| 亚洲中文日韩欧美视频| 国产一区二区在线观看日韩| 亚洲成av人片免费观看| 午夜久久久久精精品| 精品久久久久久,| 欧美成人性av电影在线观看| 精品久久久久久久久久免费视频| av在线蜜桃| 国产三级黄色录像| 首页视频小说图片口味搜索| 又紧又爽又黄一区二区| 久久6这里有精品| 丰满人妻熟妇乱又伦精品不卡| 国产精品久久久久久久电影| 哪里可以看免费的av片| 18禁黄网站禁片午夜丰满| 在线观看免费视频日本深夜| 国产高清视频在线观看网站| bbb黄色大片| 亚洲成av人片在线播放无| 亚洲精品日韩av片在线观看| 免费在线观看成人毛片| 国产伦一二天堂av在线观看| 丰满的人妻完整版| 久久午夜亚洲精品久久| 免费av不卡在线播放| 久久久久久久久久黄片| 日本黄色片子视频| 亚洲一区二区三区不卡视频| 亚洲内射少妇av| 夜夜夜夜夜久久久久| 欧美精品啪啪一区二区三区| 精品乱码久久久久久99久播| 啪啪无遮挡十八禁网站| 国产免费男女视频| 国产私拍福利视频在线观看| 国产麻豆成人av免费视频| 国产老妇女一区| АⅤ资源中文在线天堂| 老司机午夜福利在线观看视频| 成人高潮视频无遮挡免费网站| 琪琪午夜伦伦电影理论片6080| 久久久国产成人免费| 免费在线观看成人毛片| 国产精品日韩av在线免费观看| 欧美黑人欧美精品刺激| www日本黄色视频网| 久久6这里有精品| 色综合欧美亚洲国产小说| 日日夜夜操网爽| 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆| 亚洲无线在线观看| 久久久久九九精品影院| 亚洲人成网站在线播| 成人永久免费在线观看视频| 大型黄色视频在线免费观看| 国产私拍福利视频在线观看| 啪啪无遮挡十八禁网站| 蜜桃亚洲精品一区二区三区| 国产伦精品一区二区三区四那| 亚洲av五月六月丁香网| 亚洲精品日韩av片在线观看| 国产精品久久视频播放| 中文字幕av在线有码专区| 搡老岳熟女国产| 成人av在线播放网站| 成人高潮视频无遮挡免费网站| 五月玫瑰六月丁香| 别揉我奶头 嗯啊视频| 亚洲人成网站高清观看| 搡老岳熟女国产| 国产视频内射| 国产亚洲精品久久久久久毛片| 变态另类丝袜制服| 日本黄色片子视频| 日本免费一区二区三区高清不卡| 69av精品久久久久久| 窝窝影院91人妻| 自拍偷自拍亚洲精品老妇| 国产麻豆成人av免费视频| 少妇高潮的动态图| 国产精品三级大全| 淫妇啪啪啪对白视频| 亚洲真实伦在线观看| 一个人观看的视频www高清免费观看| 亚洲天堂国产精品一区在线| 久久精品国产亚洲av香蕉五月| 成人av在线播放网站| 又黄又爽又免费观看的视频| 亚洲国产精品sss在线观看| 一级黄片播放器| 亚洲欧美日韩东京热| 神马国产精品三级电影在线观看| 网址你懂的国产日韩在线| 啦啦啦观看免费观看视频高清| 欧美激情国产日韩精品一区| 看十八女毛片水多多多| 搡女人真爽免费视频火全软件 | 人人妻人人澡欧美一区二区| 国产一区二区亚洲精品在线观看| 在现免费观看毛片| 欧美最黄视频在线播放免费| 伊人久久精品亚洲午夜| 日韩成人在线观看一区二区三区| 久久精品人妻少妇| 亚洲激情在线av| 成人午夜高清在线视频| 在线免费观看不下载黄p国产 | 欧美性感艳星| 国产日本99.免费观看| 在线国产一区二区在线| 又黄又爽又刺激的免费视频.| 成年人黄色毛片网站| 日本黄色片子视频| 国产精品人妻久久久久久| 亚洲av成人不卡在线观看播放网| 亚洲av日韩精品久久久久久密| 一区福利在线观看| 1024手机看黄色片| 亚洲av成人精品一区久久| 久久精品91蜜桃| 日韩欧美三级三区| 在线观看午夜福利视频| 最近最新中文字幕大全电影3| 国产精品98久久久久久宅男小说| 午夜福利免费观看在线| 久久精品国产99精品国产亚洲性色| 国产av在哪里看| 青草久久国产| 内射极品少妇av片p| 欧美精品国产亚洲| 免费一级毛片在线播放高清视频| 亚州av有码| 91麻豆av在线| 麻豆久久精品国产亚洲av| 亚洲av熟女| 欧美绝顶高潮抽搐喷水| 嫩草影院新地址| 日韩精品中文字幕看吧| 免费人成视频x8x8入口观看| 日本成人三级电影网站| 最近中文字幕高清免费大全6 | 亚洲精品成人久久久久久| 国产69精品久久久久777片| 亚洲电影在线观看av| 久久午夜亚洲精品久久| 悠悠久久av| 精品一区二区三区视频在线| 久久欧美精品欧美久久欧美| 亚洲精品456在线播放app | 久久精品影院6| 美女被艹到高潮喷水动态| 久久久久久久午夜电影| 麻豆一二三区av精品| 中文字幕熟女人妻在线| 成人永久免费在线观看视频| 亚洲黑人精品在线| 国产精品98久久久久久宅男小说| 最后的刺客免费高清国语| 久久精品91蜜桃| 搡老妇女老女人老熟妇| 欧美日韩瑟瑟在线播放| 久久久久久久午夜电影| 国产色爽女视频免费观看| 搡老岳熟女国产| 18禁黄网站禁片午夜丰满| 亚洲成人免费电影在线观看| 90打野战视频偷拍视频| 国产午夜精品论理片| 国产成人福利小说| 亚洲天堂国产精品一区在线| 国产高清视频在线观看网站| 国产av一区在线观看免费| 亚洲美女黄片视频| 午夜福利视频1000在线观看| 成人毛片a级毛片在线播放| or卡值多少钱| 中文字幕高清在线视频| 成人三级黄色视频| 免费在线观看亚洲国产| 悠悠久久av| 免费人成视频x8x8入口观看| 久久香蕉精品热| 色播亚洲综合网| 精品久久久久久久人妻蜜臀av| 少妇被粗大猛烈的视频| 女人被狂操c到高潮| 久久精品国产亚洲av涩爱 | 亚洲片人在线观看| 夜夜夜夜夜久久久久| 在线观看av片永久免费下载| www.www免费av| 中文资源天堂在线| 欧美日韩中文字幕国产精品一区二区三区| 国产色婷婷99| 高清日韩中文字幕在线| 99久国产av精品| 久久亚洲精品不卡| 亚洲久久久久久中文字幕| 亚洲成人中文字幕在线播放| 最新在线观看一区二区三区| 国内精品一区二区在线观看| 免费观看人在逋| 国产亚洲精品综合一区在线观看| 免费电影在线观看免费观看| 国产aⅴ精品一区二区三区波| 少妇丰满av| 亚洲一区高清亚洲精品| 日本五十路高清| 亚洲成av人片免费观看| 国产精品久久久久久久电影| 国产精品嫩草影院av在线观看 | 男人的好看免费观看在线视频| 国产午夜精品论理片| 亚洲av日韩精品久久久久久密| 看免费av毛片| 精品午夜福利在线看| 天堂av国产一区二区熟女人妻| 免费高清视频大片| 直男gayav资源| 日本与韩国留学比较| 欧美色欧美亚洲另类二区| 不卡一级毛片| 女生性感内裤真人,穿戴方法视频| 一进一出好大好爽视频| 天堂动漫精品| 日韩欧美精品免费久久 | 男人狂女人下面高潮的视频| 在线a可以看的网站| 亚洲熟妇中文字幕五十中出| 国模一区二区三区四区视频| 久久国产精品人妻蜜桃| 国产免费男女视频| 亚洲激情在线av| 狂野欧美白嫩少妇大欣赏| 国产精品国产高清国产av| 国产高清视频在线观看网站| 国产精品亚洲av一区麻豆| 国产精品亚洲美女久久久| 亚洲真实伦在线观看| 老鸭窝网址在线观看| av专区在线播放| 一区福利在线观看| 午夜精品久久久久久毛片777| 91狼人影院| 久久久成人免费电影| 91av网一区二区| 18禁裸乳无遮挡免费网站照片| 美女被艹到高潮喷水动态| 亚洲国产精品999在线| 免费观看的影片在线观看| 在线国产一区二区在线| 国产精品av视频在线免费观看| 亚洲国产色片| 少妇高潮的动态图| 哪里可以看免费的av片| 极品教师在线视频| 国产一级毛片七仙女欲春2| 欧美色欧美亚洲另类二区| 成人毛片a级毛片在线播放| 一区二区三区四区激情视频 | 国产av不卡久久| 久久草成人影院| 亚洲精品色激情综合| 亚洲人成网站高清观看| 国产精品一区二区三区四区久久| 91狼人影院| 搡老妇女老女人老熟妇| 757午夜福利合集在线观看| 我要搜黄色片| 亚洲国产高清在线一区二区三| 欧美成人一区二区免费高清观看| 国产精品亚洲一级av第二区| 欧美xxxx性猛交bbbb| av天堂在线播放| 三级毛片av免费| 久久这里只有精品中国| 高清在线国产一区| 久久精品国产清高在天天线| 国产激情偷乱视频一区二区| 亚洲无线在线观看| 在线观看舔阴道视频| 一区二区三区免费毛片| 欧美三级亚洲精品| 最近在线观看免费完整版| 日韩成人在线观看一区二区三区| 午夜福利成人在线免费观看| 欧美bdsm另类| 久久久久性生活片| 亚洲国产精品999在线| 极品教师在线视频| 久久人人精品亚洲av| 国产在视频线在精品| 亚洲美女视频黄频| 内地一区二区视频在线| 成人特级av手机在线观看| 9191精品国产免费久久| 乱人视频在线观看| 小说图片视频综合网站| 久久九九热精品免费| 99久久九九国产精品国产免费| 能在线免费观看的黄片| 久久久久亚洲av毛片大全| 无人区码免费观看不卡| 国产精品久久久久久精品电影| 亚洲第一电影网av| 搡老岳熟女国产| 脱女人内裤的视频| 国产三级黄色录像| 真人做人爱边吃奶动态| 在线观看午夜福利视频| 午夜福利在线观看免费完整高清在 | av中文乱码字幕在线| 在线看三级毛片| 久久中文看片网| 欧美一级a爱片免费观看看| 琪琪午夜伦伦电影理论片6080| 国产色婷婷99| 美女高潮喷水抽搐中文字幕| 欧美最黄视频在线播放免费| 欧美性猛交黑人性爽| 国产精品野战在线观看| 欧美成人一区二区免费高清观看| 亚洲国产高清在线一区二区三| 国产69精品久久久久777片| 久久精品久久久久久噜噜老黄 | 99国产精品一区二区蜜桃av| 伦理电影大哥的女人| 色精品久久人妻99蜜桃| 亚洲精品在线观看二区| 色综合欧美亚洲国产小说| 国产一区二区在线观看日韩| 成人精品一区二区免费| 在线观看美女被高潮喷水网站 | 91午夜精品亚洲一区二区三区 | 亚洲久久久久久中文字幕| 国产不卡一卡二| 观看免费一级毛片| 欧美午夜高清在线| 在线播放无遮挡| 亚洲七黄色美女视频| 网址你懂的国产日韩在线| 亚洲狠狠婷婷综合久久图片| 动漫黄色视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲精品在线观看二区| 久久国产精品人妻蜜桃| 久久人人爽人人爽人人片va | 最新在线观看一区二区三区| 日本 av在线| 午夜久久久久精精品| 久久久成人免费电影| 99精品在免费线老司机午夜| 我要看日韩黄色一级片| 99久久99久久久精品蜜桃| 狠狠狠狠99中文字幕| 在线观看av片永久免费下载| 日韩成人在线观看一区二区三区| av中文乱码字幕在线| 久久精品国产99精品国产亚洲性色|