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

    二維重力數據徑向反演及應用

    2018-04-09 05:48:36林寶澤王明常
    石油地球物理勘探 2018年2期
    關鍵詞:場源正則多邊形

    林寶澤 肖 鋒 王明常

    (吉林大學地球探測科學與技術學院,吉林長春 130000)

    1 引言

    地球物理的反演問題主要包含兩個方面:一是地質體的幾何參數反演[1,2]; 二是地質體的物性參數反演[3,4]。針對地質體的幾何參數反演,建模一般分為三個方案:一是將場源表示為一系列垂直并列僅厚度未知的棱柱體[5-8]; 二是將場源表示為頂點坐標未知的多邊形或多面體[9-12]; 三是將包含場源的地下空間剖分成尺寸已知而密度未知的基本矩形單元[13-15]。

    徑向反演是一種最早由Silva等[16]提出的場源幾何參數反演方法; 賈真[17]研究了基于二維模型的重磁梯度分量徑向反演算法; Vanderlei等[18,19]將徑向反演方法擴展至重力異常的三維反演以及重力梯度張量的三維聯合反演。

    建模對反演來說至關重要,模型過于復雜會導致解的不穩(wěn)定性,可以通過Tikhonov正則化方法引入約束條件解決這一問題。徑向反演的優(yōu)點在于將多邊形頂點位置用極坐標表示,反演參數(即多邊形各頂點對應的半徑)具有同一量綱且反演參數相對于直角坐標系減少一半。利用Tikhonov正則化引入相關約束條件(場源幾何形態(tài)等先驗信息)時能取得較好的效果。

    徑向反演方法雖然已經發(fā)展得較為完善,但在國內關于該方法的研究和討論幾乎還是空白。此外,該方法的研究多見于模型實驗,對于影響反演結果的正則化參數選取更是缺乏討論。本文介紹了二維徑向反演方法的基本原理,且在此基礎上將自適應正則化因子選取引入該方法,并將該方法應用到實際資料處理,取得了較好的效果。

    2 方法原理

    2.1 正演方法

    徑向模型的中心點O(x0,z0)必須設定在場源內部,并給定徑向多邊形的頂點數M,將圓周平均分成M份,此時在直角坐標系下,多邊形第k個頂點的坐標(xk,zk)可以表示為

    (1)

    式中:rk為第k個頂點對應的半徑;θk為第k個頂點對應的角度(弧度,沿x軸正方向順時針旋轉為正)。模型及參數如圖1所示。

    圖1 徑向模型示意圖

    從圖1不難推斷,當頂點數M足夠大時,徑向模型能取得很好的近似效果。但是從式(1)來看,應用徑向模型時存在以下限制:一是包含場源的區(qū)域必須是單連通域;二是從模型中心點出發(fā)的射線與場源輪廓線有且僅有一個交點。

    反演過程中均勻多邊形截面二維體的重力異常正演采用文獻[20]給出的公式

    (2)

    式中:G為萬有引力常數; Δρ為密度差;Rk、Rk+1分別為測點到多邊形第k和k+1個頂點的距離;αk、αk+1分別為測點與第k和k+1個頂點連線與水平方向的夾角;φk為多邊形第k條邊與水平方向的夾角,如圖2所示。

    圖2 多邊形截面二度體正演模型示意型圖[20]

    2.2 約束反演

    設重力異常實測數據構成的N×1維向量為d,徑向多邊形各頂點對應的半徑構成的M×1維向量為r,由模型r得到的正演重力異常為N×1維向量g(r)。為了保證反演問題的適定性,引入約束條件后,該問題轉化為求解最優(yōu)化問題

    minψ=‖d-g(r)‖22+μ(Φ1+Φ2)

    (3)

    式中

    (4)

    (5)

    (6)

    (7)

    k=1,2,…,M;j=1,2,…,J

    其中:ψ為總目標函數;μ為正則化因子;r0為初始模型多邊形截面各頂點對應的半徑構成的M×1維向量;Φ1為相對鄰近性約束,其作用是使模型r的各個分量彼此趨于一致,但不能單獨使用;當W為單位矩陣時,Φ2為絕對鄰近性約束,作用是穩(wěn)定解;當W為非單位矩陣時,Φ2為會聚方向約束,使場源物性沿著某一特定方向集中,本文中將場源質量會聚于某一方向,能起到引入與場源幾何形態(tài)有關的先驗信息和穩(wěn)定解的作用;βj為指定的會聚方向;J表示會聚方向的總數;ε為較小的正數,與場源最小尺寸和最大尺寸之比有關。顯然,相對鄰近性約束與會聚方向約束是相互矛盾的,絕對鄰近性約束是會聚方向約束的特殊形式。除上述約束條件外,還有不等性約束和凸性約束。不等性約束的表達式為

    (8)

    不等性約束的作用是限制場源的幾何尺寸和保證頂點z坐標為非負。凸性約束就是要在每次迭代中滿足

    (9)

    式中ck為中心點O到半徑rk與多邊形第k-1、k+1個頂點連線交點之間的距離,如圖3所示。

    圖3 凸性約束示意圖

    由式(2)可以看出,g(r)是關于r的非線性函數,直接求解式(3)較困難,故采用Gauss-Newton迭代法將該非線性反演問題線性化。因此,在r0處對g(r)進行Taylor展開并取一級近似,則有

    =g(r0)+JΔr

    (10)

    式中J為Jacobi矩陣

    (11)

    令r=r0+Δr,并將式(10)代入式(3),則

    ψ=‖[d-g(r0)]-JΔr‖22+μ‖B(r0+Δr)‖22+

    μ‖WΔr‖22

    =μ(WΔr)TWΔr+μ[B(r0+Δr)]TB(r0+Δr)+

    {[d-g(r0)]-JΔr}T{[d-g(r0)]-JΔr}

    (12)

    對式(12)兩端求導,可得

    2μ(BTBΔr+BTBr0+WTWΔr)

    (13)

    令上式等于0,則迭代修正量Δr的表達式為

    Δr=[JTJ+μ(BTB+WTW)]-1×

    {JT[d-g(r0)]-μBTBr0}

    (14)

    2.3 自適應正則化因子選取

    正則化因子是數據擬合泛函和模型約束泛函之間的折衷系數,故它的選取直接影響反演結果。在誤差水平未知的情況下,常用的參數選取方法有GCV準則和L曲線準則。但是,對于非線性問題的迭代解法來說,這兩種準則不太適用。反演時總是將數據擬合置于優(yōu)先的地位,而在迭代過程中數據擬合泛函值越來越小,故模型約束泛函在總目標泛函中的權重也必須隨之減小,即正則化因子是隨著迭代過程衰減的。為了適應復雜的反演情況,正則化因子的調整過程應當與反演計算相聯系。采用陳小斌[21]提出的自適應調整方案,其表達式為

    (15)

    式中:μi為正則化因子經驗初值;k為迭代次數;φdf、φmc分別為數據擬合泛函與模型約束泛函。上述正則化因子的選取方式雖然顯著加快了收斂速度,但在將數據擬合置于優(yōu)先地位的情況下,迭代次數過多會導致過擬合,故本文將擬合標準差不再明顯改變時對應的迭代次數作為最佳迭代次數,所以當迭代次數達到最佳迭代次數或擬合標準差達到閾值時,迭代中止。

    3 巖蓋模型實驗

    采用巖蓋模型[17]分析不同約束條件對重力異常反演的影響。在有噪反演中,所有模型數據采用的噪聲變化范圍是原始異常幅值的±2%。

    巖蓋模型大致呈菌狀穹窿狀,傘蓋部分沿剖面方向寬約1km,傘柄部分與剖面方向大致成45°向下延伸。觀測剖面垂直于模型走向(走向垂直于紙面),其中心位于傘蓋中心正上方,沿剖面方向向兩側各延伸1000m,總長度為2000m,點距為20m,密度差為0.5g/cm3,徑向模型參數個數為100,初始正則化因子為0.1,會聚方向由虛線表示,與剖面方向分別成0°、45°、180°,會聚因子ε=0.167。圖4~圖9為巖蓋模型的反演結果。

    對比圖4與圖5、圖6與圖7、圖8與圖9可見,隨機噪聲對反演結果的影響不大。從圖4~圖7可以看出,絕對鄰近約束條件下的反演結果能反映模型上邊界的大致范圍,但是整體上反映的上邊界偏深,對下邊界特征則完全沒有反映。對比圖4與圖6、圖5與7可見,絕對和相對鄰近聯合約束下得到的反演結果比絕對鄰近約束下得到的結果更光滑,但是整體形態(tài)并沒有明顯變化。從圖8與圖9可以看出,會聚約束條件下的反演結果較為準確地反映了模型的真實形態(tài),準確的先驗信息有助于得到更好的反演結果。

    圖4 絕對鄰近約束下無噪反演結果

    圖5 絕對鄰近約束下有噪反演結果

    圖6 絕對鄰近約束和相對鄰近約束下無噪反演結果

    圖7 絕對鄰近約束和相對鄰近約束下有噪反演結果

    圖8 會聚約束下無噪反演結果

    圖9 會聚約束下有噪反演結果

    以上模型實驗表明,在會聚方向約束條件下,徑向反演方法效果最好;絕對鄰近約束條件下,反演結果能大致反映模型的上邊界,而下邊界總是趨于弧狀;相對鄰近約束能讓絕對鄰近約束條件下得到的反演結果更加光滑;凸性約束在控制反演結果的凹凸性上能起到一定的效果。

    上文已經提到自適應正則化方法雖然在一定程度上提高了數據擬合度,但是可能會造成過擬合。由于正則化因子伴隨著迭代修正不斷改變,但不是一味衰減而是伴隨著振蕩,這也就導致了擬合誤差在迭代到一定次數后也會發(fā)生振蕩。這種振蕩可能達到多個數量級,這也就導致了反演結果出現很大的變化,這種變化的好壞是未知的,所以需要確定最優(yōu)迭代次數。下面的例子僅僅是反演實驗中得到的一個特例。

    圖10 重力異常擬合標準差隨迭代次數變化

    圖11 巖蓋模型迭代100次反演結果

    圖10和圖11是絕對和相對鄰近聯合約束下得到的誤差收斂曲線及無噪數據誤差收斂曲線和反演結果。在25次迭代之后,擬合誤差出現了劇烈振蕩。對于迭代100次得到的結果,已經接近于會聚約束下才能得到的結果。但這僅僅是一個特例,筆者認為迭代中止原則還是遵循上文給出的最佳迭代次數原則為宜。

    4 應用實例

    實際資料選取的是在新墨西哥州里奧格蘭德峽谷大橋(圖12)上測量得到的重力數據[22]。里奧格蘭德峽谷大橋長約400m,高約175m。峽谷周邊地層大致可分為三種,分別為玄武巖地層,沖積沉積層,夾雜極薄湖成沉積層,其中玄武巖密度約為2.70g/cm3、沉積巖密度約為2.17g/cm3,具體情況如圖12所示。

    圖12 奧格蘭德峽谷大橋截面圖

    假定峽谷的周邊地層的85%為玄武巖、 15%為沖積沉積物,則峽谷與圍巖及空氣的密度差為-2.62g/cm3。觀測剖面長度為800m,點距為20m,徑向模型參數個數為40,初始正則化因子為0.1,初始模型半徑為25m,會聚因子ε=0.32。反演結果如圖13所示。

    從圖13a中的異常擬合圖可見,反演得到的模型正演異常很好地擬合了實測數據。從反演的場源形態(tài)看,反演結果較好地反映了場源的深度范圍,頂部和左側較好地反映了場源的邊界;右側底部與實際差異較大, 可能是由于實測異常并不是嚴格對稱的,從而導致反演結果在底部整體向左偏移。借鑒界面反演的相關做法,計算反演多邊形的半徑各頂點標準差為17.1186m,說明反演結果較可靠。

    圖13 會聚及凸性約束下反演結果

    5 結論

    徑向反演只能應用于單連通域,且初始模型中心必須位于場源內部。對于不太復雜的單個場源,在利用歐拉反褶積及相關成像等方法得到場源中心近似位置后,利用徑向反演可以圈定場源的大致形態(tài);如果能進一步獲得關于場源分布的半定量信息,利用會聚約束能得到較好的反演結果。

    正則化因子的自適應調整能提高數據擬合精度和收斂速度,但是由于實際數據的復雜性和反演的多解性,片面追求數據擬合效果可能會造成過擬合,利用最佳迭代次數選擇原則雖然能解決這一問題,但不可避免地引入人為因素。

    此外,由于本文采用的是Gauss-Newton迭代法,故模型初值的選擇對反演結果有一定的影響。本方法的初始模型一般采用給定中心埋深的無限長水平圓柱體。在實際數據反演中,地質體中心埋深往往是利用其他方法(如歐拉反褶積、特征點法等)得到的估計值,估計值越接近真實值,最終反演結果精度也越高。采用更加合理、精度更高的計算中心埋深方法與本方法結合將是后續(xù)研究的方向。當然,如果有更多的地球物理資料并建立較無限長水平圓柱體更加復雜的初始模型,反演的精度會進一步得到提高。

    [1]谷文彬,陳清禮,王余泉等.饒陽凹陷虎8北潛山重力三維松約束反演.石油地球物理勘探,2016,51(6):1219-1226.

    Gu Wenbin,Chen Qingli,Wang Yuquan et al.Part-constrained 3D gravity inversion for the Hubabei buried hill in Raoyang Sag.OGP,2016,51(6):1219-1226.

    [2]李婧,崔永謙,宋景明等.重力—地震聯合建模正演技術在SH地區(qū)潛山勘探中的應用.石油地球物理勘探,2016,51(增刊1):152-160.

    Li Jing,Cui Yongqian,Song Jingming et al.Gravity forward modeling with model built based on drilling and seismic data:an example of buried hill exploration in Area SH.OGP,2016,51(S1):152-160.

    [3]耿美霞,楊慶節(jié).應用RBF神經網絡反演二維重力密度分布.石油地球物理勘探,2013,48(4):651-657.

    Geng Meixia,Yang Qingjie.2-D density inversion with the RBF neural network method.OGP,2013,48(4):651-657.

    [4]王浩然,陳超,杜勁松.重力梯度張量數據的三維反演方法與應用.石油地球物理勘探,2013,48(3):474-481.

    Wang Haoran,Chen Chao,Du Jinsong.3-D inversion of gravity gradient tensor data and its applications.OGP,2013,48(3):474-481.

    [5]Cordell L and Henderson R G.Iterative three-dimensional solution of gravity anomaly data using a digital computer.Geophysics,1968,33(4):596-601.

    [6]Dyrelius D and Vogel A.Improvement of convergence in iterative gravity interpretation.Geophysical Journal of Royal Astronomical Society,1972,27(2):195-205.

    [7]Pedersen L B.Interpretation of potential field data - A generalized inverse approach.Geophysical Prospecting,1977,25(2):199-230.

    [8]Barbosa V C F,Silva J B C and Medeiros W E.Stable inversion of gravity anomalies of sedimentary basins with non-smooth basement reliefs and arbitrary density contrast variations.Geophysics,1999,64(3):754-764.

    [9]Corbato CE.A least-squares procedure for gravity interpretation.Geophysics,1965,30(2):228-233.

    [10]Al-Chalabi M.Interpretation of gravity anomalies by nonlinear optimization.Geophysical Prospecting,1972,20(1):1-16.

    [11]Moraes R A V and Hansen R O.Constrained inver-sion of gravity fields for complex 3D structures.Geophysics,2001,66(2):501-510.

    [12]Wildman R A and Gazonas G A.Gravitational and magnetic anomaly inversion using a tree-based geometry representation.Geophysics,2009,74(3):I23-I35.

    [13]Li Y and Oldenburg D W.3-D inversion of gravity data.Geophysics,1998,63(1):109-119.

    [14]Portniaguine O and Zhdanov M S.Focusing geophysical inversion images.Geophysics,1999,64(3):874-887.

    [15]Silva J B C and Barbosa V C F.Interactive gravity inversion.Geophysics,2006,71(1):1-9.

    [16]Silva J B C and Barbosa V C F.Generalized radial inversion of 2D potential field data.Geophysics,2004,69(6):1405-1413.

    [17]賈真.均勻二度體位場正演與徑向反演方法[學位論文].吉林長春:吉林大學,2009.

    Jia Zhen.Forward and Radial Inverse Modeling of Potential Fields Produced by a 2D Homogeneous Source[D].Jilin University,Changchun,Jilin,2009.

    [18]Vanderlei C,Oliveira J,Valeria C F et al.Source geo-metry estimation using the mass excess criterion to constrain 3-D radial inversion of gravity data.Geophysical Journal International,2011,187(2):754-772.

    [19]Vanderlei C,Oliveira J,Valeria C F et al.3-D radial gravity gradient inversion.Geophysical Journal International,2013,195(2):883-902.

    [20]賈真,孟令順.均勻多邊形截面二度體重力異常計算公式的改進.地球物理學進展,2009,24(2):462-467.

    Jia Zhen,Meng Lingshun.Some improvements on the formula for calculating the gravity anomaly due to a

    2D homogeneous polygonal source.Progress in Geophysics,2009,24(2):462-467.

    [21]陳小斌.大地電磁正反演新算法研究及資料處理與解釋的可視化集成系統(tǒng)開發(fā)[學位論文].北京:中國地

    震局地質研究所,2003.

    Chen Xiaobin.New Forward and Inversion Algorithms and a Visual Integrated System for MT data [D].Institute of Geology,China Earthquake Administration,Beijing,2003.

    [22]Titus W J,Titus S J and Davis J R.A Bayesian approach to modeling 2D gravity data using polygons.Geophysics,2017,82(1):G1-G21.

    猜你喜歡
    場源正則多邊形
    例談求解疊加電場的電場強度的策略
    基于深度展開ISTA網絡的混合源定位方法
    信號處理(2022年10期)2022-11-16 00:50:56
    多邊形中的“一個角”問題
    基于矩陣差分的遠場和近場混合源定位方法
    雷達學報(2021年3期)2021-07-05 11:31:08
    多邊形的藝術
    解多邊形題的轉化思想
    多邊形的鑲嵌
    趣味(數學)(2019年11期)2019-04-13 00:26:32
    剩余有限Minimax可解群的4階正則自同構
    類似于VNL環(huán)的環(huán)
    數學雜志(2018年5期)2018-09-19 08:13:48
    一種識別位場場源的混合小波方法
    黑人巨大精品欧美一区二区mp4| 黄片大片在线免费观看| 人妻一区二区av| 国产欧美日韩一区二区三区在线| 亚洲伊人色综图| 欧美成人午夜精品| av电影中文网址| 变态另类成人亚洲欧美熟女 | 久久久国产精品麻豆| 欧美日韩乱码在线| 99在线人妻在线中文字幕 | 91九色精品人成在线观看| 久久热在线av| 在线国产一区二区在线| 精品第一国产精品| 国产成人精品久久二区二区免费| 超碰成人久久| av国产精品久久久久影院| 女人爽到高潮嗷嗷叫在线视频| 午夜两性在线视频| 丰满饥渴人妻一区二区三| 又黄又粗又硬又大视频| 成人免费观看视频高清| 99re6热这里在线精品视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲成a人片在线一区二区| 少妇的丰满在线观看| 国内久久婷婷六月综合欲色啪| 黄色怎么调成土黄色| 一进一出抽搐动态| 最新在线观看一区二区三区| 久久久国产精品麻豆| 亚洲国产欧美一区二区综合| 日本vs欧美在线观看视频| 黄色a级毛片大全视频| 亚洲欧美一区二区三区久久| 国产99白浆流出| 在线观看一区二区三区激情| 岛国毛片在线播放| 男女下面插进去视频免费观看| 精品熟女少妇八av免费久了| 国产精品永久免费网站| 成人精品一区二区免费| 看黄色毛片网站| 色综合婷婷激情| 欧美激情久久久久久爽电影 | 国产成人欧美在线观看 | 亚洲av美国av| 99国产综合亚洲精品| 丝袜人妻中文字幕| 无人区码免费观看不卡| 久久久国产一区二区| 亚洲 欧美一区二区三区| 欧美性长视频在线观看| 国产99白浆流出| 波多野结衣av一区二区av| 黄片小视频在线播放| 香蕉丝袜av| 国产色视频综合| 在线av久久热| 国产成人av激情在线播放| 欧美乱色亚洲激情| 一区二区三区激情视频| 无遮挡黄片免费观看| 男女床上黄色一级片免费看| 91av网站免费观看| 欧美日韩福利视频一区二区| 亚洲国产中文字幕在线视频| 777米奇影视久久| 国产xxxxx性猛交| 韩国精品一区二区三区| 一级毛片高清免费大全| 一级a爱视频在线免费观看| 999久久久国产精品视频| 国产日韩欧美亚洲二区| 精品国产一区二区三区四区第35| 成年人免费黄色播放视频| 免费在线观看完整版高清| 久久国产精品影院| 精品国产亚洲在线| 欧美日韩精品网址| 变态另类成人亚洲欧美熟女 | 亚洲成国产人片在线观看| 91老司机精品| 欧美激情久久久久久爽电影 | 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲av第一区精品v没综合| 亚洲精品美女久久av网站| 女人被躁到高潮嗷嗷叫费观| 欧美乱色亚洲激情| 亚洲成国产人片在线观看| 久久久国产成人免费| 日韩大码丰满熟妇| 亚洲av熟女| 精品人妻1区二区| 色老头精品视频在线观看| 窝窝影院91人妻| 成人影院久久| 久久久久久亚洲精品国产蜜桃av| 亚洲熟女精品中文字幕| 日本a在线网址| 欧美日本中文国产一区发布| 久久国产精品男人的天堂亚洲| 一级黄色大片毛片| 91在线观看av| 一级毛片女人18水好多| 国产男女内射视频| 看片在线看免费视频| 一边摸一边抽搐一进一出视频| 久久婷婷成人综合色麻豆| 在线视频色国产色| 欧美激情 高清一区二区三区| 人人妻人人澡人人爽人人夜夜| 国产精品 国内视频| 欧美日韩一级在线毛片| 欧美黄色淫秽网站| 国产在线精品亚洲第一网站| 男男h啪啪无遮挡| av免费在线观看网站| 自线自在国产av| 日本欧美视频一区| 日本黄色日本黄色录像| 国产成人精品无人区| 久久九九热精品免费| 狂野欧美激情性xxxx| 脱女人内裤的视频| 一级毛片精品| 9热在线视频观看99| 美女视频免费永久观看网站| 亚洲欧美日韩高清在线视频| 日韩有码中文字幕| 亚洲第一欧美日韩一区二区三区| 亚洲欧洲精品一区二区精品久久久| 成人精品一区二区免费| 国产精品免费视频内射| 精品亚洲成国产av| 男女午夜视频在线观看| av天堂久久9| 精品第一国产精品| 黄色怎么调成土黄色| 美女福利国产在线| 一级毛片精品| 亚洲av片天天在线观看| 国产三级黄色录像| 久久中文字幕人妻熟女| 我的亚洲天堂| 国产成人精品在线电影| 亚洲熟女毛片儿| 国产亚洲精品一区二区www | 久久人人爽av亚洲精品天堂| 少妇的丰满在线观看| 在线观看免费日韩欧美大片| 久久人妻福利社区极品人妻图片| 人人妻人人爽人人添夜夜欢视频| 村上凉子中文字幕在线| 女同久久另类99精品国产91| 91在线观看av| 黄片播放在线免费| 校园春色视频在线观看| 美女 人体艺术 gogo| 午夜精品在线福利| 日本vs欧美在线观看视频| 久久久国产成人精品二区 | 人人澡人人妻人| 欧美日韩国产mv在线观看视频| 9热在线视频观看99| 国产精品综合久久久久久久免费 | 最近最新中文字幕大全免费视频| 最新的欧美精品一区二区| 大码成人一级视频| 成熟少妇高潮喷水视频| 精品国内亚洲2022精品成人 | 巨乳人妻的诱惑在线观看| 十八禁网站免费在线| 最近最新中文字幕大全免费视频| 在线观看日韩欧美| 亚洲精品国产一区二区精华液| 久久 成人 亚洲| 亚洲欧美日韩另类电影网站| 亚洲国产毛片av蜜桃av| 露出奶头的视频| 成年人黄色毛片网站| 高清av免费在线| 热re99久久国产66热| 女性被躁到高潮视频| 免费在线观看影片大全网站| 国产免费av片在线观看野外av| 黄色女人牲交| 久久人妻熟女aⅴ| bbb黄色大片| 欧美国产精品一级二级三级| 成人精品一区二区免费| 久久久久久人人人人人| 一边摸一边做爽爽视频免费| 老汉色∧v一级毛片| 精品一品国产午夜福利视频| 久久中文字幕人妻熟女| 在线观看日韩欧美| 午夜成年电影在线免费观看| 中文字幕av电影在线播放| 国产精品1区2区在线观看. | 视频区欧美日本亚洲| 国产精品偷伦视频观看了| e午夜精品久久久久久久| 黄色怎么调成土黄色| 免费不卡黄色视频| 久久午夜综合久久蜜桃| 亚洲av片天天在线观看| 黄色 视频免费看| 啦啦啦免费观看视频1| 午夜福利免费观看在线| 母亲3免费完整高清在线观看| x7x7x7水蜜桃| 精品久久久久久久久久免费视频 | av网站在线播放免费| 亚洲av片天天在线观看| 国产亚洲av高清不卡| 亚洲精品久久成人aⅴ小说| 丝袜美腿诱惑在线| 国产精品美女特级片免费视频播放器 | 夜夜夜夜夜久久久久| 色综合婷婷激情| 飞空精品影院首页| 男女床上黄色一级片免费看| 免费观看a级毛片全部| 亚洲国产中文字幕在线视频| 国产高清国产精品国产三级| 人人澡人人妻人| 在线十欧美十亚洲十日本专区| 国产精品国产高清国产av | 欧美精品人与动牲交sv欧美| 亚洲欧美激情综合另类| 色播在线永久视频| 精品一区二区三区视频在线观看免费 | 高清av免费在线| 精品少妇一区二区三区视频日本电影| 国产精品久久久人人做人人爽| 极品教师在线免费播放| 一级黄色大片毛片| 国产不卡一卡二| 国产一区二区三区在线臀色熟女 | 亚洲av美国av| ponron亚洲| 欧美精品高潮呻吟av久久| 19禁男女啪啪无遮挡网站| 国产免费av片在线观看野外av| 国产有黄有色有爽视频| 99精品久久久久人妻精品| 久久青草综合色| 国产无遮挡羞羞视频在线观看| 在线观看一区二区三区激情| 丝袜美腿诱惑在线| 99久久国产精品久久久| 国产精品永久免费网站| 国产精品国产高清国产av | 午夜亚洲福利在线播放| 午夜老司机福利片| 久久影院123| 国产亚洲av高清不卡| 欧美黑人精品巨大| 国产亚洲精品久久久久5区| 国产成人免费观看mmmm| 日韩免费高清中文字幕av| 99久久人妻综合| 岛国毛片在线播放| bbb黄色大片| 看片在线看免费视频| 亚洲精品中文字幕在线视频| 757午夜福利合集在线观看| 午夜日韩欧美国产| 中文字幕高清在线视频| av天堂在线播放| 亚洲精品国产精品久久久不卡| 欧美激情高清一区二区三区| 又大又爽又粗| 99香蕉大伊视频| 久久久久视频综合| 色综合欧美亚洲国产小说| 多毛熟女@视频| 日日爽夜夜爽网站| 岛国在线观看网站| 久久久精品区二区三区| 精品电影一区二区在线| 波多野结衣av一区二区av| 久久精品91无色码中文字幕| 怎么达到女性高潮| 老司机靠b影院| 日韩视频一区二区在线观看| 中出人妻视频一区二区| 欧美日韩av久久| 99热只有精品国产| 99热国产这里只有精品6| 精品卡一卡二卡四卡免费| 在线视频色国产色| 精品国产乱码久久久久久男人| 亚洲avbb在线观看| 成年版毛片免费区| 国产一卡二卡三卡精品| 两性夫妻黄色片| 日日夜夜操网爽| 村上凉子中文字幕在线| a级毛片在线看网站| 久久热在线av| 三上悠亚av全集在线观看| 午夜福利在线免费观看网站| 99久久精品国产亚洲精品| 久久精品亚洲熟妇少妇任你| 亚洲精品久久成人aⅴ小说| 一边摸一边抽搐一进一小说 | 国产一卡二卡三卡精品| av福利片在线| 91九色精品人成在线观看| 国产又爽黄色视频| 夜夜爽天天搞| e午夜精品久久久久久久| 欧美日韩瑟瑟在线播放| 日韩欧美一区二区三区在线观看 | 日日夜夜操网爽| 嫩草影视91久久| 18禁美女被吸乳视频| 法律面前人人平等表现在哪些方面| 亚洲一区中文字幕在线| 91麻豆精品激情在线观看国产 | 亚洲一区高清亚洲精品| 亚洲午夜精品一区,二区,三区| 国产91精品成人一区二区三区| 欧美乱妇无乱码| 啦啦啦在线免费观看视频4| 黑人巨大精品欧美一区二区蜜桃| 日韩欧美免费精品| 国产高清videossex| 免费在线观看影片大全网站| 俄罗斯特黄特色一大片| 制服人妻中文乱码| ponron亚洲| 免费观看人在逋| 成熟少妇高潮喷水视频| 中文亚洲av片在线观看爽 | 淫妇啪啪啪对白视频| 国产高清视频在线播放一区| 欧美日韩乱码在线| 多毛熟女@视频| 在线观看www视频免费| 久久久久久久午夜电影 | 国产一区有黄有色的免费视频| 国产成人精品久久二区二区91| 国产亚洲欧美精品永久| 亚洲三区欧美一区| 日韩熟女老妇一区二区性免费视频| 中出人妻视频一区二区| 亚洲成人免费电影在线观看| 久久久国产精品麻豆| 高清av免费在线| 欧美亚洲日本最大视频资源| 老司机影院毛片| www.精华液| 亚洲精品国产区一区二| 亚洲av电影在线进入| 日韩欧美在线二视频 | 一区二区三区精品91| 香蕉久久夜色| 国产淫语在线视频| 欧美日韩视频精品一区| 好看av亚洲va欧美ⅴa在| 国产一卡二卡三卡精品| 久99久视频精品免费| 欧美日韩精品网址| 亚洲人成77777在线视频| 大码成人一级视频| 激情视频va一区二区三区| 久久国产精品大桥未久av| 国产成人啪精品午夜网站| 久久久久久人人人人人| 91成年电影在线观看| 国产精品永久免费网站| 伊人久久大香线蕉亚洲五| 久久亚洲真实| 丁香六月欧美| 亚洲国产看品久久| 91国产中文字幕| av电影中文网址| 免费在线观看完整版高清| 亚洲欧美激情综合另类| 国产精品1区2区在线观看. | 极品人妻少妇av视频| 麻豆乱淫一区二区| 午夜福利免费观看在线| 夫妻午夜视频| 欧美日韩瑟瑟在线播放| 捣出白浆h1v1| 国产精品免费大片| 欧美 日韩 精品 国产| 男女之事视频高清在线观看| 精品久久久久久久久久免费视频 | 亚洲五月色婷婷综合| 欧美激情极品国产一区二区三区| 久久久国产一区二区| 国产99白浆流出| 90打野战视频偷拍视频| 777米奇影视久久| 国产99白浆流出| а√天堂www在线а√下载 | 亚洲av成人av| 色综合欧美亚洲国产小说| 日韩免费av在线播放| 一级,二级,三级黄色视频| 精品卡一卡二卡四卡免费| 两个人免费观看高清视频| 国产成人一区二区三区免费视频网站| 变态另类成人亚洲欧美熟女 | 老汉色∧v一级毛片| 久久99一区二区三区| 国产精品电影一区二区三区 | 老司机福利观看| 亚洲精品中文字幕在线视频| 免费人成视频x8x8入口观看| 热re99久久国产66热| 国产精品免费一区二区三区在线 | 中文字幕精品免费在线观看视频| 啦啦啦免费观看视频1| 波多野结衣一区麻豆| 777久久人妻少妇嫩草av网站| 亚洲精品国产精品久久久不卡| 成人影院久久| 国产野战对白在线观看| 超色免费av| 欧美国产精品va在线观看不卡| 日本黄色日本黄色录像| 亚洲第一欧美日韩一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 中文字幕人妻丝袜一区二区| 女警被强在线播放| 久久国产精品影院| 天天添夜夜摸| 精品一区二区三区av网在线观看| 日韩中文字幕欧美一区二区| 亚洲,欧美精品.| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品一区二区www | 国内久久婷婷六月综合欲色啪| 男女午夜视频在线观看| 在线观看午夜福利视频| 老司机深夜福利视频在线观看| 亚洲精品国产精品久久久不卡| 伦理电影免费视频| 在线观看一区二区三区激情| 一区在线观看完整版| 一二三四在线观看免费中文在| 黄色视频,在线免费观看| 国产精品1区2区在线观看. | 少妇 在线观看| 久久人妻熟女aⅴ| 成人精品一区二区免费| 欧美精品高潮呻吟av久久| 美女 人体艺术 gogo| 国产亚洲一区二区精品| 亚洲精品中文字幕一二三四区| 美女福利国产在线| 男女床上黄色一级片免费看| 国产精品永久免费网站| 淫妇啪啪啪对白视频| 日韩欧美一区二区三区在线观看 | 亚洲自偷自拍图片 自拍| 国产精品九九99| 成人影院久久| 国产淫语在线视频| 女人精品久久久久毛片| 大型黄色视频在线免费观看| 超色免费av| 国产精品98久久久久久宅男小说| 国产精品自产拍在线观看55亚洲 | 欧美午夜高清在线| 国产成人免费观看mmmm| 午夜视频精品福利| 亚洲精华国产精华精| 亚洲一区二区三区不卡视频| 80岁老熟妇乱子伦牲交| 国产成人欧美在线观看 | 午夜精品久久久久久毛片777| 一级黄色大片毛片| 日本a在线网址| 国产精品一区二区在线不卡| 精品国产一区二区三区四区第35| 国产麻豆69| 大型黄色视频在线免费观看| 欧美日韩视频精品一区| 麻豆成人av在线观看| 天天操日日干夜夜撸| 精品国产国语对白av| 国产午夜精品久久久久久| 欧美久久黑人一区二区| av欧美777| 女人被狂操c到高潮| 黑人巨大精品欧美一区二区蜜桃| 欧美成狂野欧美在线观看| 又紧又爽又黄一区二区| 人成视频在线观看免费观看| 在线观看免费高清a一片| 国产精品久久久av美女十八| 18在线观看网站| 91精品国产国语对白视频| 国产一区有黄有色的免费视频| 1024视频免费在线观看| 久久人人爽av亚洲精品天堂| 久久人人97超碰香蕉20202| а√天堂www在线а√下载 | 不卡一级毛片| 久久精品人人爽人人爽视色| 大码成人一级视频| 亚洲欧美激情在线| 人人妻人人澡人人爽人人夜夜| 亚洲免费av在线视频| 亚洲专区国产一区二区| 亚洲精品粉嫩美女一区| 777米奇影视久久| 黄网站色视频无遮挡免费观看| 天堂动漫精品| x7x7x7水蜜桃| 激情在线观看视频在线高清 | 精品国产一区二区久久| 欧美日韩瑟瑟在线播放| 高清视频免费观看一区二区| 日韩一卡2卡3卡4卡2021年| av国产精品久久久久影院| 精品亚洲成a人片在线观看| 电影成人av| 真人做人爱边吃奶动态| 国产精品.久久久| 我的亚洲天堂| 一本一本久久a久久精品综合妖精| 国产日韩一区二区三区精品不卡| 亚洲精品久久成人aⅴ小说| 欧美一级毛片孕妇| 波多野结衣av一区二区av| 首页视频小说图片口味搜索| 99久久人妻综合| 国产麻豆69| 亚洲 国产 在线| 国产欧美日韩综合在线一区二区| 婷婷成人精品国产| 如日韩欧美国产精品一区二区三区| 欧美日韩黄片免| 欧美色视频一区免费| 亚洲精品国产区一区二| 欧美日韩成人在线一区二区| 国产亚洲精品久久久久久毛片 | 老司机在亚洲福利影院| 久久久久久久久久久久大奶| 成年版毛片免费区| cao死你这个sao货| 亚洲色图av天堂| 欧美精品人与动牲交sv欧美| 欧美日韩av久久| 亚洲aⅴ乱码一区二区在线播放 | 精品午夜福利视频在线观看一区| 波多野结衣av一区二区av| 女人被狂操c到高潮| 精品亚洲成a人片在线观看| 日本vs欧美在线观看视频| 国产日韩欧美亚洲二区| 欧美黄色淫秽网站| 亚洲成人免费av在线播放| 久久久国产成人精品二区 | 一区二区日韩欧美中文字幕| 国产精品永久免费网站| 丰满的人妻完整版| 久久久久国产一级毛片高清牌| 一级a爱片免费观看的视频| 久久久国产一区二区| 天天影视国产精品| 久久草成人影院| 一级a爱片免费观看的视频| 免费久久久久久久精品成人欧美视频| 亚洲国产看品久久| 亚洲综合色网址| 日韩欧美一区二区三区在线观看 | 亚洲欧美激情在线| 99riav亚洲国产免费| 国产亚洲精品久久久久久毛片 | 精品少妇一区二区三区视频日本电影| 久久精品91无色码中文字幕| 777久久人妻少妇嫩草av网站| 国产精品久久视频播放| 精品少妇一区二区三区视频日本电影| 香蕉国产在线看| 三上悠亚av全集在线观看| 99久久精品国产亚洲精品| 波多野结衣av一区二区av| av网站在线播放免费| 国产精品久久久人人做人人爽| 中文字幕最新亚洲高清| 国产精品二区激情视频| 国产精品电影一区二区三区 | 亚洲精品久久午夜乱码| 99国产精品免费福利视频| 国产在线精品亚洲第一网站| 国产精品九九99| 亚洲欧美一区二区三区久久| 国产av精品麻豆| 亚洲va日本ⅴa欧美va伊人久久| 在线看a的网站| 又紧又爽又黄一区二区| 男人操女人黄网站| 久久人人97超碰香蕉20202| 久久精品国产99精品国产亚洲性色 | 悠悠久久av| 女警被强在线播放| 久久国产乱子伦精品免费另类| 别揉我奶头~嗯~啊~动态视频| 伊人久久大香线蕉亚洲五| 久久久久久亚洲精品国产蜜桃av| 国产又爽黄色视频| 久久午夜综合久久蜜桃| 国产一区二区三区综合在线观看| 热re99久久国产66热|