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

    標(biāo)量擬解析近似方法模擬立方體的異常電場(chǎng)

    2011-10-17 08:40:16張金會(huì)孫建國(guó)張亞?wèn)|
    物探化探計(jì)算技術(shù) 2011年2期
    關(guān)鍵詞:孫建國(guó)標(biāo)量立方體

    張金會(huì),孫建國(guó),張亞?wèn)|

    (1.吉林大學(xué) 地球探測(cè)科學(xué)與技術(shù)學(xué)院,長(zhǎng)春130026 2.國(guó)土資源部 應(yīng)用地球物理綜合解釋理論開(kāi)放實(shí)驗(yàn)室,長(zhǎng)春130026;3.安徽省勘查技術(shù)院,蚌埠233005;4.長(zhǎng)慶油田公司勘探部,西安710021)

    標(biāo)量擬解析近似方法模擬立方體的異常電場(chǎng)

    張金會(huì)1,2,3,孫建國(guó)1,2,張亞?wèn)|4

    (1.吉林大學(xué) 地球探測(cè)科學(xué)與技術(shù)學(xué)院,長(zhǎng)春130026 2.國(guó)土資源部 應(yīng)用地球物理綜合解釋理論開(kāi)放實(shí)驗(yàn)室,長(zhǎng)春130026;3.安徽省勘查技術(shù)院,蚌埠233005;4.長(zhǎng)慶油田公司勘探部,西安710021)

    擬解析近似方法是一種求解積分方程的一種近似方法,它可以處理強(qiáng)散射或者大擾動(dòng)的電磁散射問(wèn)題,在計(jì)算過(guò)程中避免了傳統(tǒng)微分?jǐn)?shù)值方法解決問(wèn)題時(shí)所遇到的大型矩陣或大型代數(shù)方程組的求解。孫建國(guó)[5]將其引入直流電場(chǎng)的積分方程中,并給出了求解異常電場(chǎng)積分方程的標(biāo)量擬解析近似公式。在以前的研究中,已經(jīng)驗(yàn)證了均勻場(chǎng)中異常球體的擬解析近似解的精度,這里對(duì)均勻場(chǎng)中的立方體異常體進(jìn)行數(shù)值模擬,得到了直流電場(chǎng)中異常立方體模型的標(biāo)量擬解析近似解。由于復(fù)雜地電模型可以用立方體的組合進(jìn)行模擬,因此對(duì)立方體異常電場(chǎng)擬解析近似解的研究,為三維直流電場(chǎng)中復(fù)雜地電模型的快速正反演模擬打下了基礎(chǔ)。

    積分方程;擬解析近似;直流電場(chǎng);立方體

    在2005年J.Sun推導(dǎo)出直流電場(chǎng)的標(biāo)量擬解析近似的理論公式基礎(chǔ)上,作者用標(biāo)量擬解析近似對(duì)均勻場(chǎng)中的異常立方體進(jìn)行研究。不同的立方體的組合,可以對(duì)復(fù)雜異常體進(jìn)行模擬。當(dāng)有多個(gè)異常體組合模擬復(fù)雜異常體時(shí),異常立方體直接的異常電場(chǎng)會(huì)受到其它異常體所產(chǎn)生的二次場(chǎng)的影響,因此在進(jìn)行模擬時(shí),還要考慮異常體所產(chǎn)生的二次場(chǎng)之間的相互影響。作者對(duì)單獨(dú)異常體的擬解析近似方法研究,為復(fù)雜地電模型的研究打下了基礎(chǔ),為三維直流電場(chǎng)中的復(fù)雜異常的正反演模擬奠定算法基礎(chǔ),可以加快三維電法勘探的進(jìn)程,具有較高的實(shí)際意義。

    1 基本方程及公式

    J.Sun在其文章中給出標(biāo)量擬解析近似方法的理論公式,作者利用其理論公式,推導(dǎo)得到擬解析近似方法在實(shí)現(xiàn)過(guò)程中的具體實(shí)現(xiàn)細(xì)節(jié)部份。

    在進(jìn)行三維異常體數(shù)值模擬時(shí),積分采用三維高斯積分進(jìn)行計(jì)算,并對(duì)積分中出現(xiàn)的格林函數(shù)的奇異性進(jìn)行了處理,采用處理方法簡(jiǎn)單,不增加運(yùn)算量的挖去法[18],在控制挖去小體積限度的基礎(chǔ)上,得到了很好的計(jì)算效果,實(shí)現(xiàn)了擬解析近似方法在直流電場(chǎng)模擬中的應(yīng)用。

    一次場(chǎng)中異常體內(nèi)部二次直流電場(chǎng)的積分方程[13]為:

    異常體外部的積分方程[13]為:

    在上面的方程中Ω表示異常體,ja(r)為二次電流密度,由式(3)給出:

    其中EBa(r)為Ea(r)的Born近似;(r',r)為并矢格林函數(shù),由式(7)給出

    J.Sun[13]于2005年通過(guò)假設(shè)電反射張量為0階張量,令:

    其中 Ι為單位張量,得到了方程(4)的解:

    通過(guò)對(duì)式(11)進(jìn)行整理,可以得到:

    在式(12)中

    將式(12)代入到式(1)中,可以得到:

    將式(14)得到的Ea(r)代入式(3),可以得到j(luò)a(r),將ja(r)代入式(2)中,即可得到異常體外部的異常電場(chǎng):

    通過(guò)方程(4),我們可以建立擬解析近似的迭代方程(16)[13]:

    其中m是迭代數(shù)。令

    利用方程(16),當(dāng)異常體電導(dǎo)率同圍巖相比變化劇烈,即大擾動(dòng)時(shí),擬解析近似解不能得到想要的結(jié)果,我們可以通過(guò)式(16)的擬解析近似的迭代解得到滿意的結(jié)果。在一般情況下,迭代次數(shù)m在不超過(guò)5的情況下,我們就能得到滿意結(jié)果。

    2 標(biāo)量擬解析近似立方體異常體電場(chǎng)值計(jì)算

    很多礦體可以由立方體組合進(jìn)行近似,選取立方體異常進(jìn)行研究,得到的異常體為立方體時(shí)的異常電場(chǎng)。在研究中給定異常體為立方體,利用標(biāo)量擬解析近似方法進(jìn)行計(jì)算,求出異常電場(chǎng)的擬解析近似解,并與Born近似解進(jìn)行比較,選取的均勻場(chǎng)場(chǎng)源中的異常立方體進(jìn)行計(jì)算。立方體邊長(zhǎng)為100 m,中心坐標(biāo)為(0 m,0 m,0 m),通過(guò)背景介質(zhì)和異常電導(dǎo)率變化,可以得到不同的模型。異常體外部異常電場(chǎng)計(jì)算點(diǎn),關(guān)于立方體中心對(duì)稱(chēng)分布在z=400 m處沿x方向。針對(duì)六種不同背景介質(zhì)電導(dǎo)率與立方體異常電導(dǎo)率的模型,分別進(jìn)行計(jì)算均勻場(chǎng)中不同電導(dǎo)率的外部電場(chǎng)。地電模型見(jiàn)圖1。

    圖1 均勻場(chǎng)中異常立方體模型示意圖Fig.1 Model of anomalous cube in the uniform field

    針對(duì)于均勻場(chǎng)中的異常立方體,我們進(jìn)行了下面幾組不同電導(dǎo)率的模型的試算。

    (1)從下頁(yè)圖2、圖3可看出Born近似解與擬解析近似解相差很小,與異常球體表現(xiàn)出的現(xiàn)象相同。

    (2)后面的圖4、圖5都是高阻立方體。由異常體引起的電場(chǎng)的Born近似解與擬解析近似解相差很大,這是由于Born近似解沒(méi)有考慮二次電場(chǎng)對(duì)異常電場(chǎng)的影響。當(dāng)二次電場(chǎng)對(duì)異常電場(chǎng)影響較大時(shí),使得擬解析近似解的數(shù)值較Born近似解高。

    (3)圖2至圖7(見(jiàn)后面)中的曲線,給出了均勻場(chǎng)中不同背景的電導(dǎo)率,與異常電導(dǎo)率立方體模型的異常電場(chǎng)Ea的x、z分量。由均勻場(chǎng)中的球體的計(jì)算結(jié)果分析,可以得出圖中的擬解析近似解比Born近似解更接近實(shí)際結(jié)果。

    3 結(jié)論

    積分方程法具有直接對(duì)異常體區(qū)域進(jìn)行數(shù)值積分的優(yōu)點(diǎn),避免了傳統(tǒng)數(shù)值方法解決問(wèn)題時(shí)所遇到的求解大型矩陣或大型代數(shù)方程組。利用標(biāo)量擬解析近似方法可以對(duì)直流電場(chǎng)的積分方程進(jìn)行求解,作者在本文對(duì)異常體為立方體時(shí)的異常電場(chǎng)進(jìn)行了標(biāo)量擬解析近似數(shù)值模擬,該方法對(duì)于其它形狀的異常體也同樣可以適用,如橢球體、圓柱體等。作者通過(guò)對(duì)均勻場(chǎng)中的異常立方體進(jìn)行的標(biāo)量擬解析模擬,表明標(biāo)量擬解析近似方法具有很強(qiáng)的適應(yīng)性。當(dāng)異常體電導(dǎo)率變化較大時(shí),可以通過(guò)對(duì)標(biāo)量擬解析近似解進(jìn)行迭代仍可以達(dá)到很好的效果,這表明該方法可以用于大擾動(dòng)問(wèn)題的模擬。立方體異常體的空間組合可以用來(lái)模擬復(fù)雜分布的地下礦體,本文作者對(duì)立方體異常的擬解析近似研究,為解決復(fù)雜地電模型提供了借鑒,初步表明了該方法具有處理復(fù)雜地電結(jié)構(gòu)的模型的能力。作者通過(guò)對(duì)直流電場(chǎng)標(biāo)量擬解析近似方法實(shí)現(xiàn)的研究,為三維直流電場(chǎng)的快速正反演模擬,提供了快速的精確度高的數(shù)值模擬方法。

    圖2 σ0=1.0 S,σa=0.01 S時(shí)異常電場(chǎng)曲線圖Fig.2 Whileσ0=1.0 S,σa=0.01 S,the anomalous electrical field curves

    圖3 σ0=1.0 S,σa=0.01 S時(shí)異常電場(chǎng)曲線圖Fig.3 Whileσ0=1.0 S,σa=0.1 S,the anomalous electrical field curves

    圖4 σ0=1.0 S,σa=-0.9 S時(shí)異常電場(chǎng)曲線圖Fig.4σ0=1.0 S,σa=-0.9 S,the anomalous electrical field curves

    圖5 σ0=1.0 S,σa=-0.99 S時(shí)異常電場(chǎng)曲線圖Fig.5σ0=1.0 S,σa=-0.99 S,the anomalous electrical field curves

    [1]HABASHY T M,GROOM W R,SPIESR B.Beyond the Born and Rytov Approximations:A Nonlinera Approach to Electromagnetic Scattering[J].Journal of geophysical research,1993,98(B2):1759.

    [2]ZHDANOV M S,F(xiàn)ANG S.Quasi-liner approximation in 3-D eledtromagnetic modeling[J].Geophysics,1996,61(3):646.

    [3]ZHDANOV M S,F(xiàn)ANG S,HURSAN G.Electromagnetic inversion using quasi-linear approximation[J].Geophysics,2000,65(5):1501.

    [4]ZHDANOV M S,DMITRIEV V I,F(xiàn)ANG S,et al.Quasi-analytical approximations and series in electromagnetic modeling[J].Geophysics,2000,65(6):1746.

    [5]孫建國(guó).聲波散射數(shù)值模擬的兩種新方案[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2006:863.

    [6]劉寧,孫建國(guó).起伏地表?xiàng)l件下的聲波散射數(shù)值模擬的積分方程法[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2007:61.

    [7]陳麗虹,孫建國(guó),王忠仁.?dāng)M線性近似方法在2_5維電阻率正演中的應(yīng)用[C].中國(guó)地球物理學(xué)會(huì)第十七屆學(xué)術(shù)年會(huì),昆明:云南科技出版社,2001.

    [8]陳麗虹,孫建國(guó).地球物理反演的擬線性近似方法綜述[J].地球物理學(xué)進(jìn)展,2002,17(3):464.

    [9]孫建國(guó).穩(wěn)定電流場(chǎng)的電位反射函數(shù)[C].中國(guó)地球物理學(xué)會(huì)第十九屆年會(huì).南京:南京師范大學(xué)出版社,2003.

    圖6 σ0=0.1 S,σa=0.9 S時(shí)異常電場(chǎng)曲線圖Fig.6σ0=0.1 S,σa=0.9 S,the anomalous electrical field curves

    圖7 σ0=0.01 S,σa=0.99 S時(shí)異常電場(chǎng)曲線圖Fig.7σ0=0.01 S,σa=0.99 S,The anomalous electrical field curves

    [10]孫建國(guó).穩(wěn)定電流場(chǎng)數(shù)值模擬的擬線性近似[C].中國(guó)地球物理學(xué)會(huì)第二十屆年會(huì),西安:西安地圖出版社,2004.

    [11]陳麗虹,孫建國(guó).穩(wěn)定電流場(chǎng)的擬線性近似方法研究[J].石油地球物理勘探,2004,39(5):589.

    [12]孫建國(guó).直流電位場(chǎng)數(shù)值模擬的擬解析法[C].中國(guó)地球物理學(xué)會(huì)第二十屆學(xué)術(shù)年會(huì),西安:西安地圖出版社,2004.

    [13]SUN J.On the electrical reflectivity tensor in d.c electric field modelling[J].Geophysical Prospecting,2005(53):411.

    [14]孫建國(guó).穩(wěn)定電流場(chǎng)擬解析近似精度分析[C].中國(guó)地球物理學(xué)會(huì)第二十一屆年會(huì).長(zhǎng)春:吉林大學(xué)出版社,2005.

    [15]孫建國(guó).穩(wěn)定電流場(chǎng)擬解析近似的數(shù)值評(píng)價(jià)[C].中國(guó)地球物理學(xué)會(huì)第二十二屆年會(huì),成都:四川科學(xué)技術(shù)出版社,2006.

    [16]傅良魁.電法勘探教程[M].北京:地質(zhì)出版社,1983.

    [17]張金會(huì),孫建國(guó).直流電場(chǎng)擬解析近似精度分析與評(píng)價(jià)[C].中國(guó)地球物理學(xué)會(huì)第二十四屆年會(huì).北京:中國(guó)大地出版社,2008.

    [18]張金會(huì),孫建國(guó).格林函數(shù)的奇異性處理[J].地球物理學(xué)進(jìn)展,2009,24(6):2106.

    P 631.3

    A

    1001—1749(2011)02—0147—06

    0 前言

    國(guó)家自然科學(xué)基金項(xiàng)目“40574052”;國(guó)家高技術(shù)研究發(fā)展計(jì)劃(863)(2006AA06Z109);國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(973)(2007CB209603)

    2010-09-28改回日期:2010-12-27

    張金會(huì)(1980-),男,河北衡水人,博士,研究方向?yàn)橹绷麟妶?chǎng)模擬。

    近年來(lái)我國(guó)資源形勢(shì)日益嚴(yán)峻,大規(guī)模的經(jīng)濟(jì)建設(shè)對(duì)金屬礦的需求越來(lái)越大,對(duì)于資源的需求應(yīng)該立足于國(guó)內(nèi),因此要加大國(guó)內(nèi)金屬礦勘探的力度。在金屬礦勘探中,電法勘探占有重要地位,研究直流電場(chǎng)快速數(shù)值模擬,對(duì)于推進(jìn)電法勘探的發(fā)展有著重要的意義。三維直流電場(chǎng)數(shù)值模擬的方法很多,大致可以分為微分方法和積分方法。在傳統(tǒng)的數(shù)值方法中,都會(huì)遇到大型矩陣或求解大型線性代數(shù)方程組的問(wèn)題,這對(duì)計(jì)算機(jī)硬件要求高,計(jì)算周期長(zhǎng)。尤其是在復(fù)雜地質(zhì)條件下,這些問(wèn)題更為突出,所以尋求快速精確的數(shù)值模擬方法,一直是地球物理工作者的研究重點(diǎn)。積分方程法在解決三維直流電場(chǎng)問(wèn)題中有其優(yōu)點(diǎn),如積分方程法只需對(duì)異常體進(jìn)行計(jì)算,而不需要對(duì)背景區(qū)域進(jìn)行計(jì)算,微分方程法在計(jì)算過(guò)程中存在誤差累積并傳遞到下一次迭代中,而積分方程法不存在誤差的累積和傳遞?;谶@些優(yōu)點(diǎn),前人對(duì)積分方程法進(jìn)行了很多的研究。T M Habashy等[1]提出用局部非線性近似方法,對(duì)電磁場(chǎng)問(wèn)題進(jìn)行研究,但這種方法運(yùn)算復(fù)雜不適于實(shí)際應(yīng)用推廣。Zhdanov等[2]受到TM Habashy的局部非線性近似方法的啟發(fā),在遇電磁場(chǎng)散射場(chǎng)問(wèn)題時(shí)提出了擬線性近似方法求解散射場(chǎng)的積分方程。這種方法通過(guò)假設(shè)異常電場(chǎng)Ea(r)與正常電場(chǎng)E0(r)之間是由電反射張量聯(lián)系起來(lái)的擬線性關(guān)系,即:。簡(jiǎn)化了電磁散射的積分方程,從而可以提高計(jì)算速度,電反射張量由最優(yōu)化的方法求出。Zhdanov[3]在2000年將擬線性近似方法應(yīng)用于電磁場(chǎng)反演,并用來(lái)反演美國(guó)新墨西哥州的CSAMT數(shù)據(jù)。由于擬線性近似時(shí)仍無(wú)法避開(kāi)解線性代數(shù)方程組,因此Zhdanov經(jīng)過(guò)研究后,提出假設(shè)電反射張量在異常體內(nèi)為緩變函數(shù),由一種相對(duì)簡(jiǎn)單的替代方法確定,這種方法稱(chēng)為擬解析近似方法[4],并提出了擬解析近似序列。孫建國(guó)[5、6]在通過(guò)電磁場(chǎng)中的擬線性近似方法和擬解析方法研究的基礎(chǔ)上,先后提出將擬線性近似方法和擬解析近似方法的應(yīng)用,擴(kuò)展到直流電場(chǎng)模擬及聲波散射模擬。孫建國(guó)[7~9]分別對(duì)2.5維直流電阻率正演和地球物理反演中擬線性近似方法的應(yīng)用進(jìn)行了討論,接下來(lái)又研究了穩(wěn)定電流場(chǎng)的電位反射函數(shù)。孫建國(guó)[10~12]2004年對(duì)穩(wěn)定電流場(chǎng)中的擬線性近似方法進(jìn)行了研究,給出了直流電位場(chǎng)數(shù)值模擬的擬解析方法。在2005年孫建國(guó)[13、14]研究了直流電場(chǎng)中的電反射張量,給出直流電場(chǎng)中的擬線性近似方法和擬解析近似方法的理論公式,并對(duì)電反射張量進(jìn)行了詳細(xì)的分析,對(duì)穩(wěn)定電流場(chǎng)擬解析近似方法也進(jìn)行精度分析。在2006年,孫建國(guó)[15]又對(duì)穩(wěn)定電流場(chǎng)擬解析近似方法進(jìn)行了數(shù)值評(píng)價(jià)。在2008年,張金會(huì)[17]用擬解析近似方法對(duì)均勻場(chǎng)中的異常球體進(jìn)行了研究,并對(duì)擬解析近似方法的精度進(jìn)行了分析和研究。

    猜你喜歡
    孫建國(guó)標(biāo)量立方體
    疊出一個(gè)立方體
    小麥籽粒大小相關(guān)基因TaGS2克隆及功能分析
    “三法”化解表面積與體積問(wèn)題
    一種高效的橢圓曲線密碼標(biāo)量乘算法及其實(shí)現(xiàn)
    一種靈活的橢圓曲線密碼并行化方法
    孫建國(guó):不懈挑戰(zhàn),讓生命自由呼吸
    圖形前線
    立方體星交會(huì)對(duì)接和空間飛行演示
    太空探索(2016年9期)2016-07-12 09:59:53
    折紙
    單調(diào)Minkowski泛函與Henig真有效性的標(biāo)量化
    女的被弄到高潮叫床怎么办| 国产无遮挡羞羞视频在线观看| 波野结衣二区三区在线| 欧美在线黄色| 超碰成人久久| 丝袜在线中文字幕| 亚洲,欧美精品.| 亚洲欧美精品自产自拍| av又黄又爽大尺度在线免费看| 国产视频首页在线观看| 老司机影院成人| 亚洲av成人不卡在线观看播放网 | 丁香六月欧美| 亚洲欧美精品综合一区二区三区| av在线老鸭窝| 免费高清在线观看日韩| 国产成人91sexporn| 高清不卡的av网站| 妹子高潮喷水视频| 国产在线视频一区二区| 免费在线观看黄色视频的| 亚洲情色 制服丝袜| 在线观看免费高清a一片| 久久久精品区二区三区| 免费在线观看完整版高清| 久久天躁狠狠躁夜夜2o2o | 啦啦啦在线免费观看视频4| 国产日韩欧美视频二区| 欧美黄色片欧美黄色片| 如日韩欧美国产精品一区二区三区| 日本猛色少妇xxxxx猛交久久| 不卡av一区二区三区| 丝瓜视频免费看黄片| 九草在线视频观看| 一区在线观看完整版| 久久久精品94久久精品| 国产亚洲av高清不卡| 精品视频人人做人人爽| 在线看a的网站| 爱豆传媒免费全集在线观看| av女优亚洲男人天堂| 国产日韩欧美在线精品| 天堂中文最新版在线下载| 免费不卡黄色视频| 久久精品亚洲熟妇少妇任你| 精品一品国产午夜福利视频| 1024视频免费在线观看| 一级毛片电影观看| 免费看av在线观看网站| 性高湖久久久久久久久免费观看| 国产成人91sexporn| 精品国产国语对白av| www.自偷自拍.com| 久久女婷五月综合色啪小说| 最新在线观看一区二区三区 | 国语对白做爰xxxⅹ性视频网站| 中文字幕制服av| 亚洲在久久综合| 免费人妻精品一区二区三区视频| 国产片内射在线| 少妇人妻精品综合一区二区| 国产毛片在线视频| 欧美人与性动交α欧美精品济南到| 青春草亚洲视频在线观看| 中文乱码字字幕精品一区二区三区| 国产日韩欧美亚洲二区| 亚洲一码二码三码区别大吗| 丰满乱子伦码专区| 9热在线视频观看99| 成年动漫av网址| 久久久久久久久久久久大奶| 亚洲国产精品一区三区| 另类精品久久| 午夜福利免费观看在线| 啦啦啦在线观看免费高清www| 高清欧美精品videossex| 无遮挡黄片免费观看| 99re6热这里在线精品视频| 成年人免费黄色播放视频| 香蕉丝袜av| 亚洲精品自拍成人| 街头女战士在线观看网站| 免费少妇av软件| 伦理电影大哥的女人| 亚洲一区二区三区欧美精品| 成人三级做爰电影| 啦啦啦在线免费观看视频4| 国产精品国产三级国产专区5o| 在线观看人妻少妇| 夫妻性生交免费视频一级片| 国产福利在线免费观看视频| 亚洲国产av新网站| 国产99久久九九免费精品| 国产熟女欧美一区二区| 久久精品国产亚洲av高清一级| 午夜老司机福利片| 我的亚洲天堂| 亚洲精品成人av观看孕妇| 无遮挡黄片免费观看| 国产精品一区二区在线观看99| 一区二区三区精品91| 国产一区二区三区av在线| av一本久久久久| 国精品久久久久久国模美| 日韩 亚洲 欧美在线| 波多野结衣av一区二区av| 国产一卡二卡三卡精品 | 国产成人精品久久久久久| 一级,二级,三级黄色视频| 成人手机av| 亚洲精品乱久久久久久| 国产av国产精品国产| 下体分泌物呈黄色| 乱人伦中国视频| 男人添女人高潮全过程视频| 亚洲av日韩精品久久久久久密 | 飞空精品影院首页| 国产日韩欧美在线精品| 成人免费观看视频高清| 汤姆久久久久久久影院中文字幕| 狂野欧美激情性xxxx| av在线播放精品| 国产乱来视频区| 99re6热这里在线精品视频| 少妇被粗大猛烈的视频| 超碰成人久久| 国产淫语在线视频| 亚洲精品日韩在线中文字幕| 中文字幕av电影在线播放| 九九爱精品视频在线观看| 精品国产一区二区三区久久久樱花| 狂野欧美激情性xxxx| 最近2019中文字幕mv第一页| 精品酒店卫生间| 性少妇av在线| 亚洲欧美激情在线| 国产成人精品福利久久| 色吧在线观看| 另类精品久久| 在线看a的网站| 国产高清国产精品国产三级| 久久久久久人妻| 国产熟女欧美一区二区| 久久久久久久久免费视频了| 无遮挡黄片免费观看| 久久久亚洲精品成人影院| 搡老乐熟女国产| 亚洲成人一二三区av| 亚洲av电影在线进入| 日韩中文字幕视频在线看片| 国产无遮挡羞羞视频在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲精品乱久久久久久| 热re99久久精品国产66热6| 欧美日本中文国产一区发布| 一本—道久久a久久精品蜜桃钙片| av女优亚洲男人天堂| 亚洲,欧美,日韩| 99热网站在线观看| 国产亚洲av高清不卡| 亚洲男人天堂网一区| 视频区图区小说| www日本在线高清视频| 美女午夜性视频免费| 欧美亚洲 丝袜 人妻 在线| 在线观看免费高清a一片| 午夜影院在线不卡| 成人黄色视频免费在线看| 国产免费一区二区三区四区乱码| 国产又色又爽无遮挡免| e午夜精品久久久久久久| 亚洲国产精品成人久久小说| 亚洲男人天堂网一区| 菩萨蛮人人尽说江南好唐韦庄| 男人爽女人下面视频在线观看| 亚洲男人天堂网一区| 日本vs欧美在线观看视频| 国产精品一区二区精品视频观看| 午夜福利在线免费观看网站| 亚洲情色 制服丝袜| 97在线人人人人妻| 两性夫妻黄色片| 丝袜美足系列| 黄片小视频在线播放| 极品人妻少妇av视频| 曰老女人黄片| 亚洲av电影在线进入| 国产精品欧美亚洲77777| 老汉色∧v一级毛片| 最近中文字幕2019免费版| 国产成人啪精品午夜网站| 亚洲精品国产区一区二| 亚洲国产欧美一区二区综合| 啦啦啦视频在线资源免费观看| 亚洲国产毛片av蜜桃av| 哪个播放器可以免费观看大片| 午夜福利影视在线免费观看| 九草在线视频观看| 久久ye,这里只有精品| 黑人巨大精品欧美一区二区蜜桃| 超色免费av| 精品久久久精品久久久| 欧美激情 高清一区二区三区| 少妇人妻久久综合中文| 在线观看www视频免费| 在线天堂最新版资源| 久久影院123| 亚洲图色成人| 一级,二级,三级黄色视频| 国产成人a∨麻豆精品| 精品一品国产午夜福利视频| 亚洲av男天堂| 黄色 视频免费看| 黄频高清免费视频| 我的亚洲天堂| 国产精品国产av在线观看| 久久这里只有精品19| 老司机在亚洲福利影院| 叶爱在线成人免费视频播放| 日韩成人av中文字幕在线观看| av在线老鸭窝| 精品国产一区二区久久| 最近2019中文字幕mv第一页| 成年美女黄网站色视频大全免费| 91精品国产国语对白视频| 一区福利在线观看| 99热网站在线观看| 999久久久国产精品视频| 伊人亚洲综合成人网| 91精品国产国语对白视频| 欧美变态另类bdsm刘玥| 国产精品人妻久久久影院| 久久久久久免费高清国产稀缺| 精品国产乱码久久久久久男人| 欧美久久黑人一区二区| 日韩精品免费视频一区二区三区| 操美女的视频在线观看| 中文字幕制服av| 日韩av不卡免费在线播放| 色播在线永久视频| 亚洲精品美女久久av网站| 中文精品一卡2卡3卡4更新| 999精品在线视频| 亚洲精品aⅴ在线观看| 国精品久久久久久国模美| 午夜福利一区二区在线看| 精品国产一区二区久久| 欧美日韩亚洲高清精品| 日韩大码丰满熟妇| 婷婷色综合www| 欧美日本中文国产一区发布| 欧美久久黑人一区二区| 国产精品国产av在线观看| 最新在线观看一区二区三区 | 亚洲国产精品国产精品| av福利片在线| 国产日韩欧美视频二区| 亚洲国产欧美一区二区综合| 亚洲精品日本国产第一区| 亚洲国产日韩一区二区| 青春草亚洲视频在线观看| 在线观看一区二区三区激情| 少妇人妻久久综合中文| 久久久久久人人人人人| 制服诱惑二区| 建设人人有责人人尽责人人享有的| 欧美激情高清一区二区三区 | 狠狠婷婷综合久久久久久88av| 天堂8中文在线网| 在线观看免费视频网站a站| 日日撸夜夜添| 久久久久久人妻| 18禁国产床啪视频网站| 亚洲av福利一区| 久久久久久久国产电影| 亚洲av中文av极速乱| 中文字幕精品免费在线观看视频| 成人亚洲精品一区在线观看| 久久97久久精品| kizo精华| 亚洲欧洲日产国产| 天天影视国产精品| 国产成人免费观看mmmm| 国产乱人偷精品视频| 香蕉国产在线看| 下体分泌物呈黄色| 国产福利在线免费观看视频| 久久精品亚洲熟妇少妇任你| av.在线天堂| 日日爽夜夜爽网站| 欧美人与善性xxx| 免费久久久久久久精品成人欧美视频| 亚洲欧美色中文字幕在线| bbb黄色大片| 超碰97精品在线观看| 人成视频在线观看免费观看| 欧美人与性动交α欧美软件| 性高湖久久久久久久久免费观看| 精品一区二区三区av网在线观看 | 国产一区二区在线观看av| 亚洲成人一二三区av| 日韩av不卡免费在线播放| 国产精品熟女久久久久浪| 国产成人精品在线电影| 91精品国产国语对白视频| 欧美黑人欧美精品刺激| av天堂久久9| 精品国产一区二区三区四区第35| 成年女人毛片免费观看观看9 | 日韩 欧美 亚洲 中文字幕| 99热网站在线观看| 极品少妇高潮喷水抽搐| 午夜老司机福利片| 亚洲图色成人| 日韩一卡2卡3卡4卡2021年| 久久久久久人人人人人| 在线观看免费午夜福利视频| 亚洲精品国产av蜜桃| 免费人妻精品一区二区三区视频| 黄网站色视频无遮挡免费观看| 美女大奶头黄色视频| 超碰成人久久| 亚洲七黄色美女视频| 9191精品国产免费久久| 精品国产一区二区三区久久久樱花| 国产1区2区3区精品| 人人妻人人爽人人添夜夜欢视频| 青青草视频在线视频观看| 亚洲视频免费观看视频| 中国国产av一级| 在线观看免费视频网站a站| 国产精品人妻久久久影院| 国产亚洲精品第一综合不卡| 精品久久久久久电影网| 国产在线一区二区三区精| 99精品久久久久人妻精品| 日韩中文字幕欧美一区二区 | 久久ye,这里只有精品| 欧美黑人精品巨大| 日韩av免费高清视频| 高清在线视频一区二区三区| 操美女的视频在线观看| 亚洲精品一区蜜桃| 操美女的视频在线观看| 无遮挡黄片免费观看| 久久久久精品久久久久真实原创| 高清欧美精品videossex| 女人被躁到高潮嗷嗷叫费观| 飞空精品影院首页| 亚洲成av片中文字幕在线观看| 日韩大片免费观看网站| 亚洲欧美激情在线| 最近中文字幕2019免费版| 在线观看国产h片| 免费日韩欧美在线观看| 18禁国产床啪视频网站| 国产一区二区在线观看av| 少妇人妻精品综合一区二区| 女人久久www免费人成看片| 日韩中文字幕视频在线看片| 亚洲精品国产区一区二| 亚洲欧美精品自产自拍| 久久97久久精品| h视频一区二区三区| 欧美日韩一级在线毛片| 女性被躁到高潮视频| 亚洲精品美女久久av网站| 自拍欧美九色日韩亚洲蝌蚪91| 黄色 视频免费看| 人成视频在线观看免费观看| 黄色毛片三级朝国网站| 午夜免费观看性视频| 日韩伦理黄色片| 九色亚洲精品在线播放| 精品免费久久久久久久清纯 | 国产精品一区二区精品视频观看| 最新的欧美精品一区二区| 日韩一区二区三区影片| 欧美日韩视频精品一区| 免费黄网站久久成人精品| 亚洲av成人精品一二三区| 乱人伦中国视频| 精品久久久久久电影网| 欧美精品高潮呻吟av久久| 亚洲精品乱久久久久久| 久久狼人影院| av一本久久久久| 亚洲精品一二三| 久久影院123| 人人妻人人爽人人添夜夜欢视频| 国精品久久久久久国模美| 五月开心婷婷网| 精品少妇内射三级| 国产有黄有色有爽视频| 久久久久精品性色| 亚洲精品久久久久久婷婷小说| 国产日韩欧美视频二区| 一区在线观看完整版| 国产欧美日韩一区二区三区在线| 亚洲美女视频黄频| 国产精品久久久av美女十八| 国产精品免费大片| 国产伦理片在线播放av一区| 一区二区三区四区激情视频| 日韩熟女老妇一区二区性免费视频| 丝袜美足系列| 亚洲欧美精品自产自拍| 久热爱精品视频在线9| 亚洲精品日韩在线中文字幕| 老司机靠b影院| 久久精品久久久久久噜噜老黄| 久久精品国产a三级三级三级| 99久久精品国产亚洲精品| 国产精品久久久av美女十八| 丰满乱子伦码专区| 不卡视频在线观看欧美| 色播在线永久视频| 男女免费视频国产| 一二三四在线观看免费中文在| 午夜日本视频在线| 老汉色∧v一级毛片| bbb黄色大片| 涩涩av久久男人的天堂| 美女扒开内裤让男人捅视频| kizo精华| 一二三四中文在线观看免费高清| 人人妻人人澡人人爽人人夜夜| 最近手机中文字幕大全| 天天操日日干夜夜撸| 久久人人爽人人片av| 国产精品无大码| 中文字幕亚洲精品专区| 一级毛片 在线播放| 亚洲国产看品久久| 国产一区二区 视频在线| 国产在线视频一区二区| 人妻人人澡人人爽人人| 伊人亚洲综合成人网| 建设人人有责人人尽责人人享有的| 日本午夜av视频| 亚洲国产精品国产精品| 中文精品一卡2卡3卡4更新| 国产午夜精品一二区理论片| 亚洲色图综合在线观看| 欧美日韩一区二区视频在线观看视频在线| 婷婷色综合www| 一级片免费观看大全| 在线观看www视频免费| 国产国语露脸激情在线看| 丁香六月天网| 午夜福利网站1000一区二区三区| 秋霞在线观看毛片| bbb黄色大片| 97在线人人人人妻| 国产成人啪精品午夜网站| 欧美黄色片欧美黄色片| 9191精品国产免费久久| av网站免费在线观看视频| 激情视频va一区二区三区| 看非洲黑人一级黄片| 久久久国产欧美日韩av| 国产精品 欧美亚洲| 一本一本久久a久久精品综合妖精| 中文字幕av电影在线播放| 欧美激情高清一区二区三区 | av福利片在线| 国产精品久久久av美女十八| 久久午夜综合久久蜜桃| 人人妻人人澡人人爽人人夜夜| 在现免费观看毛片| 欧美亚洲 丝袜 人妻 在线| 狠狠精品人妻久久久久久综合| 日日撸夜夜添| 国产精品国产av在线观看| 一区二区三区四区激情视频| 欧美xxⅹ黑人| 精品国产超薄肉色丝袜足j| 日韩精品有码人妻一区| 美女中出高潮动态图| 成人漫画全彩无遮挡| 人人妻人人澡人人爽人人夜夜| 亚洲精品日本国产第一区| 亚洲精品,欧美精品| 日韩制服丝袜自拍偷拍| 欧美日韩亚洲国产一区二区在线观看 | av片东京热男人的天堂| 80岁老熟妇乱子伦牲交| 久久国产亚洲av麻豆专区| 你懂的网址亚洲精品在线观看| 亚洲精品国产区一区二| 欧美日韩一级在线毛片| 国产免费又黄又爽又色| 日韩成人av中文字幕在线观看| 最近中文字幕2019免费版| 一二三四在线观看免费中文在| 最黄视频免费看| 日本91视频免费播放| 久久久国产一区二区| 高清视频免费观看一区二区| 国产在视频线精品| 免费人妻精品一区二区三区视频| 亚洲少妇的诱惑av| 成人影院久久| www.av在线官网国产| 免费久久久久久久精品成人欧美视频| 日韩,欧美,国产一区二区三区| 亚洲综合精品二区| 久久99精品国语久久久| 乱人伦中国视频| svipshipincom国产片| 看免费成人av毛片| 中文乱码字字幕精品一区二区三区| 如日韩欧美国产精品一区二区三区| videos熟女内射| 搡老乐熟女国产| 亚洲av电影在线进入| 亚洲欧美精品综合一区二区三区| 在线观看免费午夜福利视频| 日韩精品免费视频一区二区三区| 高清不卡的av网站| 欧美精品人与动牲交sv欧美| 久久久久久久大尺度免费视频| 亚洲一级一片aⅴ在线观看| 欧美老熟妇乱子伦牲交| 亚洲国产最新在线播放| 又大又爽又粗| 亚洲伊人色综图| 伦理电影免费视频| 亚洲av电影在线进入| 免费黄色在线免费观看| 午夜老司机福利片| 午夜免费男女啪啪视频观看| 色播在线永久视频| 丝袜人妻中文字幕| 高清在线视频一区二区三区| 亚洲av电影在线观看一区二区三区| 久久女婷五月综合色啪小说| 国产亚洲av高清不卡| 免费不卡黄色视频| 尾随美女入室| 免费在线观看视频国产中文字幕亚洲 | 午夜福利一区二区在线看| 中文字幕高清在线视频| 天天躁夜夜躁狠狠久久av| 狠狠精品人妻久久久久久综合| 男女边吃奶边做爰视频| 午夜福利在线免费观看网站| 亚洲熟女毛片儿| 日韩大片免费观看网站| 99国产精品免费福利视频| 中文字幕精品免费在线观看视频| 亚洲色图 男人天堂 中文字幕| 国产日韩欧美在线精品| 日韩av免费高清视频| 亚洲视频免费观看视频| 黄色怎么调成土黄色| 久久精品人人爽人人爽视色| 啦啦啦视频在线资源免费观看| 欧美日韩亚洲综合一区二区三区_| 下体分泌物呈黄色| 久久精品亚洲av国产电影网| 欧美亚洲日本最大视频资源| 国产极品粉嫩免费观看在线| 天天躁狠狠躁夜夜躁狠狠躁| 多毛熟女@视频| 亚洲国产精品999| 亚洲av电影在线观看一区二区三区| 美女国产高潮福利片在线看| 成人毛片60女人毛片免费| 亚洲中文av在线| 男女边吃奶边做爰视频| 久久精品国产亚洲av涩爱| 欧美乱码精品一区二区三区| 人成视频在线观看免费观看| 啦啦啦啦在线视频资源| 波多野结衣一区麻豆| 国产日韩欧美视频二区| 少妇被粗大猛烈的视频| 日日啪夜夜爽| 免费久久久久久久精品成人欧美视频| 黄片小视频在线播放| av卡一久久| 天天操日日干夜夜撸| 国产精品免费视频内射| 操出白浆在线播放| 日韩视频在线欧美| 午夜免费鲁丝| 国产av码专区亚洲av| 久久久久国产精品人妻一区二区| 国产精品成人在线| 精品亚洲成国产av| 欧美另类一区| 亚洲国产日韩一区二区| 热99国产精品久久久久久7| 91国产中文字幕| 亚洲成人免费av在线播放| 亚洲精品国产色婷婷电影| 国产在线免费精品| 国产片特级美女逼逼视频| 91国产中文字幕| 精品亚洲乱码少妇综合久久| 亚洲综合色网址| 午夜福利,免费看| 女人爽到高潮嗷嗷叫在线视频| 久久久久久久久久久久大奶| 另类精品久久| 亚洲精品乱久久久久久| 黄片小视频在线播放| 99九九在线精品视频| 亚洲精品aⅴ在线观看| 美国免费a级毛片| 交换朋友夫妻互换小说| 麻豆精品久久久久久蜜桃| 哪个播放器可以免费观看大片| videos熟女内射|