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

    多面函數(shù)擬合區(qū)域垂直形變場(chǎng)中核函數(shù)的優(yōu)化選擇

    2017-11-29 05:14:24鄭洪艷王世進(jìn)
    華北地震科學(xué) 2017年4期

    田 曉,鄭洪艷,張 超,王世進(jìn)

    (中國(guó)地震局第一監(jiān)測(cè)中心,天津 300180)

    多面函數(shù)擬合區(qū)域垂直形變場(chǎng)中核函數(shù)的優(yōu)化選擇

    田 曉,鄭洪艷,張 超,王世進(jìn)

    (中國(guó)地震局第一監(jiān)測(cè)中心,天津 300180)

    應(yīng)用4種核函數(shù)對(duì)3個(gè)不同地區(qū)的垂直形變場(chǎng)進(jìn)行多面函數(shù)擬合比較,并通過(guò)外符合精度以及得到的實(shí)際等值線圖來(lái)綜合選擇較優(yōu)核函數(shù)。通過(guò)對(duì)比發(fā)現(xiàn),錐面函數(shù)是一種很好的核函數(shù),適應(yīng)性較強(qiáng),能較好地?cái)M合出區(qū)域垂直形變場(chǎng)的升降變化,具有一定的實(shí)際應(yīng)用價(jià)值。

    多面函數(shù);垂直形變場(chǎng);核函數(shù);優(yōu)化選擇

    0 引言

    目前,利用重復(fù)高精度水準(zhǔn)測(cè)量仍然是研究地殼垂直運(yùn)動(dòng)的主要方法。但由于水準(zhǔn)重合點(diǎn)是離散的,且水準(zhǔn)網(wǎng)內(nèi)部沒(méi)有數(shù)據(jù),導(dǎo)致動(dòng)態(tài)平差得到的垂直運(yùn)動(dòng)速率在時(shí)間域和空間域都是不連續(xù)的,實(shí)際上地殼垂直運(yùn)動(dòng)是一種連續(xù)的形態(tài)變化。為了實(shí)現(xiàn)在時(shí)空域內(nèi)客觀準(zhǔn)確地描述形變分布特征,往往需要基于離散的測(cè)點(diǎn)數(shù)據(jù)構(gòu)建合適的數(shù)學(xué)模型擬合整個(gè)地區(qū)的地殼垂直形變速率面。

    在地殼垂直形變場(chǎng)的擬合中,多面函數(shù)法是一種傳統(tǒng)的、得到廣泛應(yīng)用的方法。多面函數(shù)法中核函數(shù)的選擇是一個(gè)至關(guān)重要的方面,會(huì)嚴(yán)重影響整個(gè)地區(qū)的擬合效果。在我國(guó),20世紀(jì)80年代末期開始,多位學(xué)者研究了多面函數(shù)法在擬合地殼垂直形變場(chǎng)中的應(yīng)用,這其中也對(duì)核函數(shù)的選擇做了相關(guān)研究,并得出了不同的觀點(diǎn)。楊國(guó)華等[1]對(duì)多面函數(shù)的特性進(jìn)行了初步數(shù)值分析,指出就地殼垂直形變分析來(lái)講,選用冪次不高的“缽”形函數(shù)作為核函數(shù)為佳;趙承坤等[2]研究了速率面擬合法中核函數(shù)中心點(diǎn)的選擇,提出用逐步回歸理論來(lái)篩選核函數(shù)中心點(diǎn),收到了較好的效果,并指出雙曲型核函數(shù)一般優(yōu)于三次曲面核函數(shù),其中又以正雙曲型核函數(shù)為佳;陶本藻等[3-4]對(duì)多面函數(shù)擬合法做了相關(guān)試驗(yàn)研究,建議在地殼垂直形變分析中,使用正雙曲面型核函數(shù),并使平滑因子為0,應(yīng)用此核函數(shù)研究了鮮水河地區(qū)地殼垂直運(yùn)動(dòng),獲得了該地區(qū)地殼垂直形變速率面的預(yù)測(cè)結(jié)果,較好地反映了地殼垂直形變及應(yīng)變積累的特征;黃立人等[5]在多面函數(shù)用于擬合地殼垂直運(yùn)動(dòng)的研究中,指出雙曲型函數(shù)優(yōu)于三次函數(shù),且在大多數(shù)情況下,倒雙曲型更好一些;董鴻聞等[6-7]在中國(guó)大陸現(xiàn)今地殼垂直運(yùn)動(dòng)的研究中,認(rèn)為以全國(guó)范圍采用倒雙曲型核函數(shù)為宜,并且核函數(shù)中心點(diǎn)取所有重合點(diǎn)的1/3較好。另外,郝明在研究青藏高原東緣地區(qū)的垂直運(yùn)動(dòng)時(shí)采用了正雙曲面核函數(shù),平滑因子取了0.01,核函數(shù)中心點(diǎn)取了1/3水準(zhǔn)重合點(diǎn)[8];秦姍蘭和劉同文在研究西秦嶺地區(qū)的垂直形變時(shí)分別采用了倒雙曲面核函數(shù)和正雙曲面核函數(shù)[9-10]。

    總之,學(xué)者們對(duì)多面函數(shù)擬合地殼垂直形變場(chǎng)中核函數(shù)的選擇主要集中在選擇正雙曲面還是倒雙曲面核函數(shù),并且大多學(xué)者是采用回代點(diǎn)(擬合點(diǎn))或者插值點(diǎn)(檢核點(diǎn))的擬合精度這一指標(biāo)來(lái)判別選擇哪種核函數(shù)。實(shí)際應(yīng)用中這樣選擇核函數(shù)存在不妥之處,因?yàn)橐粋€(gè)好的核函數(shù)應(yīng)該對(duì)于已測(cè)點(diǎn)上的運(yùn)動(dòng)速率有良好的擬合效果,對(duì)于未測(cè)點(diǎn)的預(yù)測(cè)值較為平穩(wěn),不至于在預(yù)測(cè)點(diǎn)偏離數(shù)據(jù)點(diǎn)時(shí),內(nèi)插值有急劇的變化。但是對(duì)于實(shí)測(cè)水準(zhǔn)網(wǎng)不管采用擬合點(diǎn)(內(nèi)符合)精度還是檢核點(diǎn)(外符合)精度都只是檢核水準(zhǔn)網(wǎng)上是否有良好的擬合效果,對(duì)于水準(zhǔn)網(wǎng)內(nèi)部的擬合效果得不到驗(yàn)證。因此,本文采用檢核點(diǎn)(外符合)精度以及利用該核函數(shù)得到的實(shí)際等值線圖(驗(yàn)證水準(zhǔn)網(wǎng)內(nèi)部是否有急劇的變化)來(lái)選擇較優(yōu)的核函數(shù)。

    1 多面函數(shù)擬合法

    多面函數(shù)擬合法是美國(guó)Hardy教授于1977年提出的,并建議用于地殼垂直形變分析。該方法的基本思想是:任何數(shù)學(xué)表面和任何不規(guī)則的圓滑表面,總可以用一系列有規(guī)則的數(shù)學(xué)表面的和以任意精度逼近[10]。根據(jù)這一思想,假設(shè)地球表面上有一點(diǎn)(x,y),它的垂直速率值ξ(x,y)可以用下式來(lái)表示:

    (1)

    式中,aj為待定系數(shù),Q(x,y,xj,yj)是多面函數(shù)的核函數(shù)。核函數(shù)有多種類型,常用的是距離型核函數(shù)和錐面函數(shù),其中距離型核函數(shù)表達(dá)式為:

    Q(x,y,xj,yj)=[(x-xj)2+(y-yj)2+δ]μ

    (2)

    式中:δ為平滑因子,通常可取一小正數(shù)或0;μ一般取1/2、-1/2、3/2。當(dāng)μ=1/2時(shí),核函數(shù)稱為正雙曲面函數(shù);當(dāng)μ=-1/2時(shí),核函數(shù)稱為倒雙曲面函數(shù);當(dāng)μ=3/2時(shí),核函數(shù)稱為三次曲面函數(shù)。

    錐面函數(shù)表達(dá)式為:

    (3)

    式中:c為待定參數(shù)。

    設(shè)有m個(gè)已知水準(zhǔn)點(diǎn)(xi,yi),選取其中n(n≤m)個(gè)點(diǎn)作為核函數(shù)的中心點(diǎn)(xj,yj),令Qij=Q(xi,yi,xj,yj),則各數(shù)據(jù)點(diǎn)應(yīng)滿足:

    (4)

    由此可列出誤差方程:

    V=QX-ξ

    (5)

    根據(jù)最小二乘原理可求得待定系數(shù)X,即

    X=(QTQ)-1QTξ

    (6)

    待定系數(shù)求出后,根據(jù)式(4)可計(jì)算測(cè)區(qū)各未知點(diǎn)的垂直形變速率值。

    2 算例與分析

    2.1 算例介紹

    本文選取3個(gè)不同地區(qū)的各2期一等水準(zhǔn)網(wǎng)數(shù)據(jù),分別采用自由網(wǎng)動(dòng)態(tài)平差計(jì)算得到各個(gè)地區(qū)的垂直形變速率值,并剔除平差后速率中個(gè)別量值過(guò)大的突變點(diǎn),然后均勻選取1/5點(diǎn)作為外部檢核點(diǎn),其余點(diǎn)用于擬合多面函數(shù)模型參數(shù)。應(yīng)用多面函數(shù)中4種核函數(shù)對(duì)3個(gè)地區(qū)進(jìn)行擬合比較,其中距離型核函數(shù)的平滑因子是從0到5(約100D,D為水準(zhǔn)點(diǎn)平均間距)間隔0.01共501次計(jì)算中根據(jù)外符合精度選取的最優(yōu)值,錐面函數(shù)參數(shù)是從0到100共101次計(jì)算中根據(jù)外符合精度選取的最優(yōu)值,核函數(shù)中心點(diǎn)是均勻選取1/3點(diǎn)。平差情況如下:

    算例1:晉冀蒙地區(qū)2002年和2006年2期數(shù)據(jù),以大同基巖點(diǎn)為參考點(diǎn),平差后的單位權(quán)中誤差為1.15 mm/a,共有627個(gè)具有平差后速率的數(shù)據(jù)點(diǎn)(圖1)。

    圖1 晉冀蒙地區(qū)垂直形變速率矢量圖

    算例2:鄂爾多斯地區(qū)1990年和2013年2期數(shù)據(jù),以棋東49基為參考點(diǎn),平差后的單位權(quán)中誤差為0.84 mm/a,共有343個(gè)具有平差后速率的數(shù)據(jù)點(diǎn)(圖2)。

    算例3:蘇魯皖地區(qū)1985年和1995年2期數(shù)據(jù),以青島原點(diǎn)為參考點(diǎn),平差后的單位權(quán)中誤差為1.10 mm/a,共有796個(gè)具有平差后速率的數(shù)據(jù)點(diǎn)(圖3)。

    圖2 鄂爾多斯地區(qū)垂直形變速率矢量圖

    圖3 蘇魯皖地區(qū)垂直形變速率矢量圖

    為了選出較優(yōu)的核函數(shù),本文采用外符合精度檢核水準(zhǔn)網(wǎng)上的擬合效果,以及繪制等值線圖檢核水準(zhǔn)網(wǎng)內(nèi)部的擬合效果,進(jìn)行綜合選取。外符合精度σ的表達(dá)式為:

    (7)

    式中:Δi表示檢核點(diǎn)速率擬合值與實(shí)際觀測(cè)值之差,n為檢核點(diǎn)個(gè)數(shù)。

    2.2 不同核函數(shù)擬合比較

    分別應(yīng)用多面函數(shù)法中的4種核函數(shù)對(duì)3個(gè)不同地區(qū)的垂直形變場(chǎng)進(jìn)行擬合,得到的最優(yōu)外符合精度以及對(duì)應(yīng)的參數(shù)(表1~3)。

    表1 晉冀蒙地區(qū)擬合比較表

    表2 鄂爾多斯地區(qū)擬合比較表

    表3 蘇魯皖地區(qū)擬合比較表

    對(duì)3個(gè)地區(qū)應(yīng)用各自得到的4種核函數(shù)進(jìn)行格網(wǎng)化插值,插值間隔設(shè)為1′×1′,并繪制3個(gè)地區(qū)的垂直形變速率等值線圖(圖4~6),其中4種核函數(shù)分別對(duì)應(yīng)每幅圖中的a、b、c、d四幅圖。

    a 正雙曲面;b 倒雙曲面; c 三次曲面;d 錐面圖4 晉冀蒙地區(qū)垂直形變速率等值線圖

    a 正雙曲面;b 倒雙曲面;c 三次曲面;d 錐面圖5 鄂爾多斯地區(qū)垂直形變速率等值線圖

    a 正雙曲面;b 倒雙曲面;c 三次曲面;d 錐面圖6 蘇魯皖地區(qū)垂直形變速率等值線圖

    2.3 算例結(jié)果分析

    1)從表1和圖4可以看出,當(dāng)錐面函數(shù)作為核函數(shù)時(shí)擬合精度優(yōu)于其他3種核函數(shù),并且用錐面函數(shù)繪制的等值線圖與圖1反映的升降情況一致,不僅水準(zhǔn)網(wǎng)上具有一致的速率值,而且水準(zhǔn)網(wǎng)內(nèi)部也變化平穩(wěn),較好地反映了整個(gè)地區(qū)的升降情況。前3種核函數(shù)雖然對(duì)水準(zhǔn)網(wǎng)上具有較好的擬合效果,但是水準(zhǔn)網(wǎng)內(nèi)部偏離數(shù)據(jù)點(diǎn)的位置產(chǎn)生急劇的升降變化,不符合地殼垂直形變場(chǎng)擬合的原則,失真度較高,不應(yīng)采用。綜合來(lái)看,晉冀蒙地區(qū)應(yīng)該選擇錐面函數(shù)作為核函數(shù)。

    (2)從表2和圖5可以看出,正雙曲面函數(shù)和錐面函數(shù)的擬合精度一樣,都是最高的,而且參數(shù)也相同,實(shí)際上此情況下,2種核函數(shù)完全一致,具有相同的表達(dá)式,說(shuō)明當(dāng)平滑因子為0時(shí)的正雙曲面函數(shù)是錐面函數(shù)的一個(gè)特例。從等值線圖看,用這2個(gè)核函數(shù)繪制的等值線與圖2反映的升降情況一致,在水準(zhǔn)網(wǎng)上和內(nèi)部都表現(xiàn)出較好的擬合效果。但是倒雙曲面函數(shù)同樣能反映圖2的升降情況,并且水準(zhǔn)網(wǎng)內(nèi)部也變化比較平穩(wěn),同樣適用該地區(qū),不過(guò)從擬合精度來(lái)看,比正雙曲面和錐面函數(shù)略低。三次曲面為核函數(shù)時(shí)插值效果不好,失真度較高,不應(yīng)采用。綜合來(lái)看,鄂爾多斯地區(qū)應(yīng)該選擇正雙曲面和錐面函數(shù)為核函數(shù)。

    (3)從表3和圖6可以看出,三次曲面為核函數(shù)時(shí)的擬合精度最高,說(shuō)明此例中以三次曲面為核函數(shù)時(shí)對(duì)路線上有較好的擬合效果,由于擬合精度不能對(duì)環(huán)內(nèi)部給出真實(shí)的反映,在此例中對(duì)比等值線圖可以發(fā)現(xiàn)水準(zhǔn)網(wǎng)內(nèi)部的插值效果很差(繪圖時(shí)對(duì)于小于-10 mm/a的區(qū)域一律用-10 mm/a顏色代替,對(duì)于大于10 mm/a的區(qū)域一律用10 mm/a顏色代替),偏離數(shù)據(jù)點(diǎn)時(shí)產(chǎn)生了急劇的升降變化,失真度很高,不是好的核函數(shù)。正雙曲面和倒雙曲面存在同樣的問(wèn)題。雖然錐面函數(shù)為核函數(shù)時(shí)擬合精度最低,但是從等值線圖可以發(fā)現(xiàn)其能反映圖3的升降情況,并且水準(zhǔn)網(wǎng)內(nèi)部變化比較平穩(wěn),能較好地反映該地區(qū)垂直形變場(chǎng)的升降變化。綜合來(lái)看,蘇魯皖地區(qū)應(yīng)該選擇錐面函數(shù)作為核函數(shù)。

    3 結(jié)語(yǔ)

    1)多面函數(shù)擬合區(qū)域垂直形變場(chǎng)時(shí),核函數(shù)選擇的好壞能產(chǎn)生至關(guān)重要的作用。并且不能單靠擬合精度來(lái)判別選擇哪種核函數(shù)最好,應(yīng)該結(jié)合等值線圖中水準(zhǔn)網(wǎng)內(nèi)部的擬合情況來(lái)選擇最優(yōu)核函數(shù)。

    2)正雙曲面函數(shù)和倒雙曲面函數(shù)作為核函數(shù)有時(shí)也能有較好的擬合效果,但是這兩種核函數(shù)具有反映特別靈敏的平滑因子需要選取,如果選取不當(dāng),擬合效果會(huì)很差,實(shí)際選取比較麻煩。但是錐面函數(shù)此時(shí)同樣具有較好的擬合效果,而且其參數(shù)的選取對(duì)擬合效果靈敏度較低。

    3)從本文3個(gè)具體算例來(lái)看,錐面函數(shù)是一種很好的核函數(shù),適應(yīng)性較強(qiáng),能較好地?cái)M合出區(qū)域垂直形變場(chǎng)的升降變化。綜合來(lái)看,建議多面函數(shù)擬合垂直形變場(chǎng)時(shí)選用錐面函數(shù)作為核函數(shù)。

    [1] 楊國(guó)華, 黃立人. 速率面擬合法中多面函數(shù)幾個(gè)特性的初步數(shù)值研究[J]. 地殼形變與地震, 1990,10(4):70-82.

    [2] 趙承坤, 黃立人. 速率面擬合法中核函數(shù)中心點(diǎn)的選擇[J]. 地殼形變與地震, 1991,11(2):48-54.

    [3] 陶本藻, 王新洲, 于正林,等. 用于垂直形變模型的多面函數(shù)擬合法的試驗(yàn)研究[J]. 地殼形變與地震, 1992,12(1):1-13.

    [4] 陶本藻, 杜方. 多面函數(shù)法研究鮮水河地區(qū)現(xiàn)今地殼垂直運(yùn)動(dòng)[J]. 武漢測(cè)繪科技大學(xué)學(xué)報(bào), 1991, 16(3):12-21.

    [5] 黃立人, 陶本藻, 趙承坤. 多面函數(shù)擬合在地殼垂直運(yùn)動(dòng)研究中的應(yīng)用[J]. 測(cè)繪學(xué)報(bào), 1993,22(1):25-32.

    [6] 董鴻聞, 顧旦生, 李國(guó)智,等. 中國(guó)大陸現(xiàn)今地殼垂直運(yùn)動(dòng)研究[J]. 測(cè)繪學(xué)報(bào), 2002, 31(2):100-103.

    [7] 王文利, 陳士銀, 董鴻聞,等. 利用多面函數(shù)擬合中國(guó)陸地垂直運(yùn)動(dòng)速率圖[J]. 測(cè)繪通報(bào), 2002(8):6-8,11.

    [8] 郝明. 基于精密水準(zhǔn)數(shù)據(jù)的青藏高原東緣現(xiàn)今地殼垂直運(yùn)動(dòng)與典型地震同震及震后垂直形變研究[D]. 北京:中國(guó)地震局地質(zhì)研究所, 2012.

    [9] 秦姍蘭, 王慶良, 季靈運(yùn),等. 利用水準(zhǔn)資料研究西秦嶺地區(qū)的垂直形變[J]. 大地測(cè)量與地球動(dòng)力學(xué), 2012, 32(2):16-19.

    [10] 劉同文, 楊志強(qiáng), 王慧敏. 西秦嶺地區(qū)垂直形變的分析研究[J]. 測(cè)繪科學(xué), 2014, 39(4):61-63.

    OptimizationSelectionofKernelFunctioninFittingRegionVerticalDeformationFieldUsingMulti-surfaceFunction

    TIAN Xiao, ZHENG Hong-yan, ZHANG Chao, WANG Shi-jin

    (The First Monitoring and Application Center, CEA, Tianjin 300180, China)

    Selecting kernel function in fitting region vertical deformation field using multi-surface function is an important aspect. Four kernel functions are compared in fitting vertical deformation field of three different regions in this paper, and the best kernel function is selected by the external accuracy and the contour map. Tapered surface function is a good kernel function with a strong adaptability by comparison, and it can well fit the rise and fall of region vertical deformation field. Therefore it has some practical values.

    multi-surface function; vertical deformation field; kernel function; optimization selection

    田曉,鄭洪艷,張超,等. 多面函數(shù)擬合區(qū)域垂直形變場(chǎng)中核函數(shù)的優(yōu)化選擇[J].華北地震科學(xué),2017,35(4):11-15.

    2017-04-27

    中國(guó)綜合地球物理場(chǎng)觀測(cè)——大華北地區(qū)(201508009);中國(guó)地震局監(jiān)測(cè)、預(yù)測(cè)、科研三結(jié)合課題(CEA-JC/3JH-163711)

    田 曉(1988—),男,山東人,助理工程師,主要從事精密水準(zhǔn)和跨斷層數(shù)據(jù)處理.E-mail:tx_yice@163.com

    P315.7

    A

    1003-1375(2017)04-0011-05

    10.3969/j.issn.1003-1375.2017.04.002

    91字幕亚洲| 两个人看的免费小视频| 国产精品成人在线| 欧美激情高清一区二区三区| 国产亚洲欧美在线一区二区| 啦啦啦 在线观看视频| 色94色欧美一区二区| 丁香欧美五月| 婷婷精品国产亚洲av在线 | av网站免费在线观看视频| 19禁男女啪啪无遮挡网站| 日韩视频一区二区在线观看| 精品国产一区二区三区久久久樱花| 免费女性裸体啪啪无遮挡网站| 俄罗斯特黄特色一大片| 在线天堂中文资源库| 精品国产美女av久久久久小说| 精品视频人人做人人爽| 一a级毛片在线观看| 国产真人三级小视频在线观看| 亚洲成人免费av在线播放| 亚洲国产毛片av蜜桃av| 精品国产一区二区三区四区第35| 无遮挡黄片免费观看| 久久国产精品男人的天堂亚洲| 精品人妻1区二区| 免费黄频网站在线观看国产| 热re99久久精品国产66热6| 色精品久久人妻99蜜桃| 国产精品久久久av美女十八| 欧美日韩av久久| 在线国产一区二区在线| 久久性视频一级片| 视频在线观看一区二区三区| 免费观看a级毛片全部| 波多野结衣一区麻豆| 国产极品粉嫩免费观看在线| 日本一区二区免费在线视频| 麻豆国产av国片精品| 波多野结衣一区麻豆| 热99re8久久精品国产| 国产亚洲欧美98| 久久精品国产99精品国产亚洲性色 | 岛国毛片在线播放| 美女高潮到喷水免费观看| 欧美日本中文国产一区发布| 视频区欧美日本亚洲| 亚洲九九香蕉| 夜夜爽天天搞| 老鸭窝网址在线观看| 国产精品久久视频播放| 久久精品91无色码中文字幕| 成人黄色视频免费在线看| aaaaa片日本免费| 国产亚洲欧美精品永久| www.精华液| 久久午夜综合久久蜜桃| 精品少妇久久久久久888优播| 一区二区三区国产精品乱码| 亚洲专区字幕在线| www日本在线高清视频| av免费在线观看网站| 一级黄色大片毛片| 国产无遮挡羞羞视频在线观看| 亚洲成人免费电影在线观看| 老汉色∧v一级毛片| 亚洲精品在线观看二区| 国产xxxxx性猛交| 丰满的人妻完整版| 亚洲精品美女久久av网站| 俄罗斯特黄特色一大片| 亚洲黑人精品在线| 亚洲人成电影观看| 热99re8久久精品国产| 老司机午夜福利在线观看视频| 国产精品久久电影中文字幕 | 久久影院123| 大香蕉久久成人网| 国产精品偷伦视频观看了| 中文字幕人妻熟女乱码| 国产97色在线日韩免费| 人妻久久中文字幕网| 欧美国产精品一级二级三级| 久久青草综合色| 黑人操中国人逼视频| 又黄又粗又硬又大视频| 水蜜桃什么品种好| 日本黄色视频三级网站网址 | 男人操女人黄网站| bbb黄色大片| 成年女人毛片免费观看观看9 | 欧美日韩精品网址| 欧美 日韩 精品 国产| av中文乱码字幕在线| 人妻 亚洲 视频| 在线天堂中文资源库| 国产精品av久久久久免费| 午夜福利在线观看吧| 国产成人一区二区三区免费视频网站| 自拍欧美九色日韩亚洲蝌蚪91| 人人妻人人添人人爽欧美一区卜| 中文字幕精品免费在线观看视频| 亚洲中文字幕日韩| 日韩欧美国产一区二区入口| 欧美在线黄色| 19禁男女啪啪无遮挡网站| 男女免费视频国产| 亚洲精华国产精华精| 波多野结衣一区麻豆| 91麻豆av在线| 久久影院123| 亚洲一区二区三区不卡视频| 天堂俺去俺来也www色官网| 国产深夜福利视频在线观看| 国产免费av片在线观看野外av| 亚洲性夜色夜夜综合| 99riav亚洲国产免费| 欧美日韩福利视频一区二区| 美女扒开内裤让男人捅视频| 国产精品一区二区免费欧美| 国产在视频线精品| 激情视频va一区二区三区| 国产野战对白在线观看| 久久精品国产清高在天天线| 久久精品人人爽人人爽视色| 免费人成视频x8x8入口观看| 夜夜爽天天搞| 久久久久久久久久久久大奶| 99国产精品一区二区蜜桃av | 欧美人与性动交α欧美软件| 国产精品久久久久久精品古装| 国产深夜福利视频在线观看| 在线观看免费视频网站a站| 日本精品一区二区三区蜜桃| 香蕉丝袜av| 久久精品91无色码中文字幕| 老司机午夜十八禁免费视频| 久久人妻av系列| 久久精品国产综合久久久| 国产精品亚洲一级av第二区| 搡老岳熟女国产| 99re6热这里在线精品视频| 午夜亚洲福利在线播放| 国产精华一区二区三区| 最新在线观看一区二区三区| 人人澡人人妻人| 不卡av一区二区三区| 国产精品综合久久久久久久免费 | 十八禁人妻一区二区| 可以免费在线观看a视频的电影网站| 精品少妇久久久久久888优播| 亚洲午夜精品一区,二区,三区| 久久人人97超碰香蕉20202| 国产精品香港三级国产av潘金莲| 视频区欧美日本亚洲| 麻豆av在线久日| 欧美黄色淫秽网站| 人妻丰满熟妇av一区二区三区 | 免费黄频网站在线观看国产| 宅男免费午夜| 亚洲欧美一区二区三区黑人| 一级作爱视频免费观看| 正在播放国产对白刺激| 欧美乱色亚洲激情| 好男人电影高清在线观看| 午夜福利乱码中文字幕| 亚洲熟妇中文字幕五十中出 | 交换朋友夫妻互换小说| 视频区图区小说| 国产在线精品亚洲第一网站| 高清欧美精品videossex| 天天躁夜夜躁狠狠躁躁| 午夜免费观看网址| 少妇裸体淫交视频免费看高清 | 在线国产一区二区在线| 精品福利观看| 亚洲成人免费av在线播放| 女性生殖器流出的白浆| 一级,二级,三级黄色视频| 国产成人免费观看mmmm| 一级毛片精品| 色老头精品视频在线观看| 成人特级黄色片久久久久久久| 超碰97精品在线观看| 亚洲精品美女久久久久99蜜臀| 国产精品久久视频播放| 精品午夜福利视频在线观看一区| 美女福利国产在线| 亚洲三区欧美一区| 精品人妻在线不人妻| 叶爱在线成人免费视频播放| 波多野结衣一区麻豆| 两个人看的免费小视频| 欧美日韩av久久| 亚洲欧美色中文字幕在线| 国产精品九九99| 国产成人影院久久av| 婷婷丁香在线五月| 美女福利国产在线| 不卡一级毛片| 亚洲国产毛片av蜜桃av| 捣出白浆h1v1| 大香蕉久久成人网| 国产精品一区二区在线不卡| 一边摸一边抽搐一进一出视频| 91成人精品电影| 午夜福利免费观看在线| x7x7x7水蜜桃| 国产成人av激情在线播放| 一级黄色大片毛片| 一a级毛片在线观看| 性色av乱码一区二区三区2| 国产成人一区二区三区免费视频网站| 亚洲精品久久午夜乱码| 美女国产高潮福利片在线看| 精品亚洲成国产av| 成人免费观看视频高清| 午夜精品久久久久久毛片777| 高清在线国产一区| www日本在线高清视频| 国产精品久久久av美女十八| 久久香蕉国产精品| 后天国语完整版免费观看| 国产在线精品亚洲第一网站| 亚洲专区国产一区二区| 精品视频人人做人人爽| 免费观看人在逋| 一进一出抽搐动态| 精品国产乱子伦一区二区三区| 丰满饥渴人妻一区二区三| 国产av一区二区精品久久| 精品福利永久在线观看| 成人国语在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 国产男女超爽视频在线观看| av福利片在线| 国产精品一区二区免费欧美| 天堂俺去俺来也www色官网| 久久精品国产清高在天天线| av网站免费在线观看视频| 国产男女内射视频| 自拍欧美九色日韩亚洲蝌蚪91| av国产精品久久久久影院| 中文字幕精品免费在线观看视频| 婷婷成人精品国产| 国产午夜精品久久久久久| 国产精品成人在线| 两个人看的免费小视频| 国产一区二区三区视频了| 人人妻人人爽人人添夜夜欢视频| 久久天躁狠狠躁夜夜2o2o| 色综合欧美亚洲国产小说| 国产免费男女视频| 免费黄频网站在线观看国产| 欧美激情高清一区二区三区| 免费不卡黄色视频| 在线十欧美十亚洲十日本专区| 99re在线观看精品视频| 人成视频在线观看免费观看| 成熟少妇高潮喷水视频| 手机成人av网站| 国产精品久久久人人做人人爽| 中国美女看黄片| av网站在线播放免费| 无限看片的www在线观看| 欧美在线黄色| 天天操日日干夜夜撸| 在线av久久热| 国产免费男女视频| 精品一品国产午夜福利视频| 51午夜福利影视在线观看| 在线av久久热| 免费在线观看视频国产中文字幕亚洲| 国产日韩欧美亚洲二区| 亚洲av第一区精品v没综合| 首页视频小说图片口味搜索| 无限看片的www在线观看| 极品少妇高潮喷水抽搐| 18禁美女被吸乳视频| 校园春色视频在线观看| 岛国在线观看网站| 两个人免费观看高清视频| 一级片免费观看大全| 搡老乐熟女国产| 动漫黄色视频在线观看| 国产欧美日韩精品亚洲av| 国产精品美女特级片免费视频播放器 | 亚洲精品国产色婷婷电影| 麻豆乱淫一区二区| 日本黄色日本黄色录像| 日本一区二区免费在线视频| 欧美乱码精品一区二区三区| 精品亚洲成国产av| 最近最新中文字幕大全免费视频| 久久久国产一区二区| 777久久人妻少妇嫩草av网站| 久久国产乱子伦精品免费另类| cao死你这个sao货| 伦理电影免费视频| 妹子高潮喷水视频| 亚洲av日韩精品久久久久久密| 精品国产一区二区久久| 搡老熟女国产l中国老女人| 飞空精品影院首页| 欧美黑人欧美精品刺激| 操出白浆在线播放| 久久久水蜜桃国产精品网| 在线观看免费日韩欧美大片| 脱女人内裤的视频| 久久人人97超碰香蕉20202| 国产精品av久久久久免费| 一级片'在线观看视频| 欧美成人午夜精品| 男女床上黄色一级片免费看| 热99re8久久精品国产| av电影中文网址| 亚洲第一av免费看| 亚洲国产欧美一区二区综合| xxx96com| 丝袜在线中文字幕| 大陆偷拍与自拍| 美女扒开内裤让男人捅视频| 黄网站色视频无遮挡免费观看| 19禁男女啪啪无遮挡网站| 最新在线观看一区二区三区| 国产精品久久久久成人av| 在线观看免费视频网站a站| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美日韩高清在线视频| 成人国产一区最新在线观看| 中文字幕人妻丝袜制服| 可以免费在线观看a视频的电影网站| 欧美日韩av久久| 亚洲精品一二三| x7x7x7水蜜桃| 欧美激情久久久久久爽电影 | 国产精品免费大片| 亚洲精品一卡2卡三卡4卡5卡| 国产一卡二卡三卡精品| 久久精品亚洲av国产电影网| 黄色毛片三级朝国网站| tube8黄色片| 色精品久久人妻99蜜桃| 怎么达到女性高潮| 精品人妻1区二区| 在线看a的网站| 最近最新中文字幕大全电影3 | 97人妻天天添夜夜摸| 深夜精品福利| 很黄的视频免费| 欧美不卡视频在线免费观看 | 又黄又爽又免费观看的视频| 热99久久久久精品小说推荐| 黄片大片在线免费观看| 色播在线永久视频| 天天影视国产精品| 国产精华一区二区三区| 国产1区2区3区精品| 日韩成人在线观看一区二区三区| 亚洲一码二码三码区别大吗| 午夜免费鲁丝| 亚洲成人免费av在线播放| 十八禁人妻一区二区| 日日爽夜夜爽网站| 午夜福利视频在线观看免费| 天天操日日干夜夜撸| 国产99久久九九免费精品| 涩涩av久久男人的天堂| 在线看a的网站| 欧美精品av麻豆av| 日韩免费av在线播放| 男女下面插进去视频免费观看| 两个人看的免费小视频| 婷婷精品国产亚洲av在线 | 中文字幕精品免费在线观看视频| 亚洲精品一二三| aaaaa片日本免费| 久久国产精品影院| 天天添夜夜摸| 久久人妻av系列| 女人高潮潮喷娇喘18禁视频| 国产精品久久久久久人妻精品电影| 欧美久久黑人一区二区| tube8黄色片| 日韩有码中文字幕| 757午夜福利合集在线观看| 人人妻人人爽人人添夜夜欢视频| 国产极品粉嫩免费观看在线| 下体分泌物呈黄色| 国产一区二区三区综合在线观看| 成熟少妇高潮喷水视频| 老熟妇仑乱视频hdxx| 免费女性裸体啪啪无遮挡网站| 两个人看的免费小视频| 亚洲国产毛片av蜜桃av| 男人的好看免费观看在线视频 | 一区二区三区精品91| 成人永久免费在线观看视频| 99精品欧美一区二区三区四区| 久久草成人影院| 身体一侧抽搐| 最近最新中文字幕大全免费视频| 精品国内亚洲2022精品成人 | 男女之事视频高清在线观看| xxx96com| 天天躁狠狠躁夜夜躁狠狠躁| 国产av又大| 在线十欧美十亚洲十日本专区| 99在线人妻在线中文字幕 | 精品一品国产午夜福利视频| 亚洲人成电影免费在线| 老熟妇乱子伦视频在线观看| 亚洲精品粉嫩美女一区| 老司机福利观看| 久久久久久久国产电影| 免费人成视频x8x8入口观看| 男男h啪啪无遮挡| 另类亚洲欧美激情| 亚洲精品乱久久久久久| 久久婷婷成人综合色麻豆| 精品福利永久在线观看| 久久国产亚洲av麻豆专区| 一进一出好大好爽视频| 色婷婷av一区二区三区视频| 天堂俺去俺来也www色官网| 少妇裸体淫交视频免费看高清 | 美女高潮喷水抽搐中文字幕| 精品久久蜜臀av无| 午夜精品国产一区二区电影| 国产精品久久视频播放| av电影中文网址| 国产一区二区三区综合在线观看| 欧美日韩亚洲国产一区二区在线观看 | 建设人人有责人人尽责人人享有的| 在线观看日韩欧美| 日本一区二区免费在线视频| 欧美激情久久久久久爽电影 | 亚洲中文日韩欧美视频| 欧美午夜高清在线| 国产97色在线日韩免费| 日韩免费av在线播放| 超色免费av| 曰老女人黄片| 久久国产精品大桥未久av| 最新美女视频免费是黄的| 一级a爱片免费观看的视频| 操出白浆在线播放| 亚洲中文字幕日韩| 亚洲aⅴ乱码一区二区在线播放 | 看片在线看免费视频| 精品电影一区二区在线| 亚洲精品粉嫩美女一区| www.熟女人妻精品国产| tube8黄色片| 女警被强在线播放| 一边摸一边做爽爽视频免费| 99国产精品一区二区蜜桃av | 国产成人精品在线电影| 黄色怎么调成土黄色| 久久精品aⅴ一区二区三区四区| aaaaa片日本免费| 亚洲国产欧美日韩在线播放| 国产黄色免费在线视频| 老司机亚洲免费影院| 午夜福利在线免费观看网站| 久久热在线av| 两个人免费观看高清视频| 老司机福利观看| 日韩成人在线观看一区二区三区| 欧美人与性动交α欧美精品济南到| 成年人午夜在线观看视频| 国内毛片毛片毛片毛片毛片| 99国产精品99久久久久| 美女高潮喷水抽搐中文字幕| 免费久久久久久久精品成人欧美视频| videos熟女内射| 三级毛片av免费| 国产高清视频在线播放一区| 91国产中文字幕| 亚洲欧美精品综合一区二区三区| 精品国产乱码久久久久久男人| 日日爽夜夜爽网站| 亚洲男人天堂网一区| 搡老乐熟女国产| 久久亚洲精品不卡| 国产97色在线日韩免费| 极品少妇高潮喷水抽搐| 很黄的视频免费| 韩国精品一区二区三区| 久久性视频一级片| 国产乱人伦免费视频| 男女之事视频高清在线观看| 欧美日韩国产mv在线观看视频| 黄色视频不卡| 午夜免费鲁丝| 日本精品一区二区三区蜜桃| 搡老乐熟女国产| 亚洲视频免费观看视频| 美女扒开内裤让男人捅视频| 国产在线观看jvid| 久久久久视频综合| 色在线成人网| 夜夜夜夜夜久久久久| 女人被躁到高潮嗷嗷叫费观| 久久久久久久久久久久大奶| 在线观看免费高清a一片| 男人舔女人的私密视频| 国产蜜桃级精品一区二区三区 | 日本撒尿小便嘘嘘汇集6| 9色porny在线观看| 动漫黄色视频在线观看| 免费黄频网站在线观看国产| 午夜福利在线免费观看网站| av天堂久久9| 欧美精品高潮呻吟av久久| 欧美黄色淫秽网站| 久久国产乱子伦精品免费另类| 国产亚洲欧美98| 婷婷精品国产亚洲av在线 | 午夜两性在线视频| 91成年电影在线观看| 国产aⅴ精品一区二区三区波| 一a级毛片在线观看| 亚洲精品乱久久久久久| 波多野结衣av一区二区av| 十分钟在线观看高清视频www| 国产男靠女视频免费网站| 性少妇av在线| 乱人伦中国视频| 国产成人av教育| 涩涩av久久男人的天堂| 免费不卡黄色视频| 午夜福利视频在线观看免费| 狂野欧美激情性xxxx| 久久久国产欧美日韩av| 女人精品久久久久毛片| 精品久久久久久,| 久久久久视频综合| 亚洲美女黄片视频| 狠狠婷婷综合久久久久久88av| 久久影院123| 免费在线观看完整版高清| 国产aⅴ精品一区二区三区波| 久久热在线av| 一本一本久久a久久精品综合妖精| 国产aⅴ精品一区二区三区波| 免费在线观看完整版高清| 亚洲精品久久午夜乱码| 亚洲精品美女久久av网站| 窝窝影院91人妻| av有码第一页| 日日夜夜操网爽| 美女高潮到喷水免费观看| 脱女人内裤的视频| 成人免费观看视频高清| 99久久99久久久精品蜜桃| 久久中文字幕人妻熟女| 女人爽到高潮嗷嗷叫在线视频| 国产成人精品在线电影| 亚洲午夜理论影院| 国产有黄有色有爽视频| 色老头精品视频在线观看| 精品卡一卡二卡四卡免费| 最近最新中文字幕大全电影3 | √禁漫天堂资源中文www| av免费在线观看网站| 村上凉子中文字幕在线| av网站免费在线观看视频| 精品第一国产精品| 亚洲精品成人av观看孕妇| 亚洲熟妇熟女久久| 大码成人一级视频| 免费在线观看影片大全网站| 午夜亚洲福利在线播放| 91大片在线观看| avwww免费| 日本撒尿小便嘘嘘汇集6| 757午夜福利合集在线观看| 国产精品av久久久久免费| 亚洲精品久久午夜乱码| 日本欧美视频一区| 三级毛片av免费| 丰满饥渴人妻一区二区三| 交换朋友夫妻互换小说| netflix在线观看网站| 国产精品免费一区二区三区在线 | 大码成人一级视频| 好男人电影高清在线观看| 久久精品成人免费网站| 无限看片的www在线观看| 亚洲久久久国产精品| 午夜久久久在线观看| 午夜免费成人在线视频| 在线观看免费午夜福利视频| 成年人午夜在线观看视频| 久久人妻福利社区极品人妻图片| 日日夜夜操网爽| 人人妻人人爽人人添夜夜欢视频| 免费久久久久久久精品成人欧美视频| 日韩免费av在线播放| 18禁裸乳无遮挡免费网站照片 | 国产91精品成人一区二区三区| 熟女少妇亚洲综合色aaa.| 婷婷精品国产亚洲av在线 | 黑人巨大精品欧美一区二区mp4| 狠狠狠狠99中文字幕| 在线国产一区二区在线| 五月开心婷婷网| 国产1区2区3区精品| 91麻豆精品激情在线观看国产 | 中文字幕人妻熟女乱码| 在线观看免费视频日本深夜| 成人亚洲精品一区在线观看| 国精品久久久久久国模美| 亚洲精品美女久久av网站|