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

    聯(lián)合重力異常和重力垂直梯度異常數(shù)據(jù)反演皇帝山海域海底地形

    2022-02-04 09:20:12陽(yáng)凡林沈瑞杰梅賽屠澤杰辛明真
    海洋學(xué)報(bào) 2022年12期
    關(guān)鍵詞:水深賦權(quán)重力

    陽(yáng)凡林,沈瑞杰,梅賽,屠澤杰,辛明真

    (1.山東科技大學(xué) 測(cè)繪與空間信息學(xué)院,山東 青島 266590;2.自然資源部海洋測(cè)繪重點(diǎn)實(shí)驗(yàn)室,山東 青島 266590;3.中國(guó)地質(zhì)調(diào)查局 青島海洋地質(zhì)研究所,山東 青島 266237)

    1 引言

    海洋約占地球表面積的71%,蘊(yùn)藏著極其豐富的礦產(chǎn)和生物資源。海底地形數(shù)據(jù)是研究海洋生物、海洋化學(xué)及海洋地球物理學(xué)等相關(guān)學(xué)科的基礎(chǔ)。自20 世紀(jì)50 年代回波測(cè)深技術(shù)廣泛應(yīng)用以來(lái),搭載單波束和多波束的研究船可以開(kāi)展高精度、高分辨率的海底地形勘測(cè)。然而船載水深測(cè)量存在效率低、測(cè)量結(jié)果分布不均勻的問(wèn)題,使用回聲測(cè)深技術(shù)對(duì)全球海底地形進(jìn)行完全覆蓋式勘測(cè)仍需要很長(zhǎng)時(shí)間[1]。

    衛(wèi)星測(cè)高技術(shù)的出現(xiàn),可以全天候、連續(xù)準(zhǔn)確地在全球范圍內(nèi)直接獲取海面高度的觀測(cè)數(shù)據(jù),其測(cè)量效率和覆蓋范圍是船載水深測(cè)量無(wú)法比擬的。隨著衛(wèi)星測(cè)高資料的積累,利用測(cè)高數(shù)據(jù)獲取的海洋重力信息已成為反演海底地形的主要數(shù)據(jù)來(lái)源。目前,利用重力數(shù)據(jù)進(jìn)行海底地形反演的研究方法主要有:重力地質(zhì)法(Gravity Geologic Method,GGM)[2–3]、導(dǎo)納函數(shù)法(Admittance Method)[4]、最小二乘配置法(Least Squares Collocation,LSC)[5–6]和Smith and Sandwell(SAS)法[7]。其中,GGM 和SAS 法使用最為廣泛,兩者大多以重力異常數(shù)據(jù)作為輸入。Kim 等[8]將GGM 與SAS法構(gòu)建的海底地形模型進(jìn)行比較,發(fā)現(xiàn)GGM 在預(yù)測(cè)短波長(zhǎng)(≤12 km)分量方面具有優(yōu)勢(shì),而SAS 方法在預(yù)測(cè)長(zhǎng)波長(zhǎng)(≥25 km)分量方面具有優(yōu)勢(shì);李倩倩和鮑李峰[9],彭聰?shù)萚10]均對(duì)SAS 方法和GGM 進(jìn)行了比較分析,后者給出了兩種方法的水深和海底地形復(fù)雜程度的適用條件;范雕等[11]應(yīng)用導(dǎo)納函數(shù)法建立了不同均衡補(bǔ)償模式下的海底地形模型,并對(duì)各項(xiàng)地球物理參數(shù)的影響進(jìn)行了詳細(xì)分析;郭金運(yùn)等[12]利用迭代延拓法,提升了重力異常向下延拓的精度。

    除重力異常外,重力垂直梯度異常(Vertical Gravity Gradient Anomalies,VGG)也是反演海底地形的重要數(shù)據(jù),其反映了重力場(chǎng)在垂直方向上的變化情況,且重力垂直梯度異常相比于重力異常含有更多的高頻信號(hào)成分,能夠更好地描述異常特征[13–14];Wang[15]指出大洋地殼的質(zhì)量補(bǔ)償對(duì)重力梯度的影響很小,提出可以利用重力垂直梯度異常反演水深的思路;Wessel 和Lyons[16]通過(guò)對(duì)高斯型海山模型的研究,指出重力垂直梯度異常具有放大短波信號(hào)、抑制長(zhǎng)波信號(hào)的作用;吳云孫等[17]推導(dǎo)了重力垂直梯度異常與海底地形之間的關(guān)系并反演了南海海域的海底地形模型;胡敏章等[18]通過(guò)疊加重力異常(波長(zhǎng)<100 km)和重力垂直梯度異常(波長(zhǎng)100~200 km)反演的海底地形構(gòu)建了中國(guó)海域的海底地形模型,其精度與V15.1 相近;范雕等[19]利用多元回歸分析技術(shù),通過(guò)解算重力異常和重力垂直梯度異常的偏回歸系數(shù)和常數(shù)項(xiàng)來(lái)反演海底地形。

    目前,聯(lián)合重力異常和重力垂直梯度異常數(shù)據(jù)進(jìn)行海底地形反演的研究較少。本文針對(duì)重力異常、重力垂直梯度異常和船測(cè)水深數(shù)據(jù)各自的特性,通過(guò)穩(wěn)健回歸分析降低了船測(cè)數(shù)據(jù)和重力數(shù)據(jù)中異常點(diǎn)的偏差影響,分別使用重力異常和重力垂直梯度異常數(shù)據(jù)反演了太平洋皇帝山海域(35°~45°N,165°~175°E)的海底地形模型,采用滑動(dòng)窗口賦權(quán)的方法,選擇滑動(dòng)窗口內(nèi)的最優(yōu)權(quán)值,聯(lián)合兩種重力數(shù)據(jù)構(gòu)建的模型。最后,以船測(cè)水深數(shù)據(jù)為外部檢核,與DTU18 和SIO V23.1 全球海底地形模型以及地形剖面處的實(shí)測(cè)水深數(shù)據(jù)進(jìn)行比較分析。

    2 研究區(qū)域和數(shù)據(jù)來(lái)源

    夏威夷皇帝海山鏈位于太平洋西北側(cè),該區(qū)域存在海底山脈且兩側(cè)海隆聚集,海底深度跨度大(200~8 000 m),是開(kāi)展海底地形反演的熱點(diǎn)區(qū)域。本文選取的研究區(qū)域位于其南部(35°~45°N,165°~175°E),考慮到邊緣效應(yīng),將研究區(qū)域的經(jīng)緯度均向外擴(kuò)充1°,研究區(qū)域如圖1 所示。

    圖1 研究區(qū)域與船載測(cè)深軌跡分布Fig.1 Study area and shipborne bathymetry trace distribution map

    2.1 船載測(cè)深數(shù)據(jù)

    本文采用的船載測(cè)深數(shù)據(jù)來(lái)自于美國(guó)國(guó)家地球物理數(shù)據(jù)中心(National Geophysical Data Center,NGDC)發(fā)布的實(shí)測(cè)數(shù)據(jù),其中單波束數(shù)據(jù)時(shí)間跨度為42 a(1961–2002 年),多波束數(shù)據(jù)時(shí)間跨度為16 a(2004–2019 年)。由于早期數(shù)據(jù)可能存在測(cè)量誤差,以全球地形模型ETOPO1 作為約束條件進(jìn)行互差比較,依照“3σ”準(zhǔn)則,采用10′×10′的窗口大小,5′的移動(dòng)步長(zhǎng)滑動(dòng)剔除單波束數(shù)據(jù)中的粗差點(diǎn),得到61 169 個(gè)單波束船測(cè)點(diǎn)。均勻選取其中4/5 的船測(cè)點(diǎn)作為控制點(diǎn),剩余船測(cè)單波束數(shù)據(jù)和濾波抽稀后的多波束數(shù)據(jù)用于后續(xù)模型的精度評(píng)定及地形剖面對(duì)比。

    2.2 重力數(shù)據(jù)

    重力異常數(shù)據(jù)(圖2a)和重力垂直梯度異常數(shù)據(jù)(圖2b)均來(lái)源于加州大學(xué)圣地亞哥分校(University of California,San Diego,UCSD)的斯克里普斯海洋研究所(Scripps Institution of Oceanography,SIO)于2021年8 月發(fā)布的V31.1 版本,格網(wǎng)分辨率為1′×1′。對(duì)比先前版本,V31.1 模型補(bǔ)充了新的Altika、Cryosat LRM、Cryosat SAR 和Sentinel-3A/B 衛(wèi)星測(cè)高數(shù)據(jù)。為了方便說(shuō)明,下文將重力異常和重力垂直梯度異常數(shù)據(jù)合稱為重力數(shù)據(jù)。

    圖2 重力模型示意圖Fig.2 Gravimetric model diagram

    2.3 海底地形模型

    本文采用目前廣泛使用的ETOPO1、DTU18 和SIO V23.1 海底地形模型(1′×1′)作為參考模型。

    ETOPO1 由美國(guó)國(guó)家地球物理中心(National Geophysical Data Center,NGDC)和美國(guó)國(guó)家海洋和大氣管理局(National Oceanic and Atmospheric Administration,NOAA)于2008 年8 月發(fā)布,該模型集成了來(lái)自區(qū)域和全球數(shù)據(jù)集的地形、測(cè)深和海岸線數(shù)據(jù)。

    DTU18 由丹麥科技大學(xué)(Technical University of Denmark,DTU)于2019 年9 月發(fā)布,該模型補(bǔ)充使用了3 a 的Sentinel-3A 衛(wèi)星和7 a 的Cryosat-2 LRM 衛(wèi)星數(shù)據(jù),并使用FES2014 海潮模型進(jìn)行潮汐改正。

    SIO V23.1 由美國(guó)SIO 的Sandwell 團(tuán)隊(duì)于2021 年11 月發(fā)布,該模型結(jié)合了最新的重力模型數(shù)據(jù)、衛(wèi)星測(cè)高數(shù)據(jù)以及優(yōu)化的船測(cè)水深數(shù)據(jù),是公認(rèn)精度最高的全球海底地形模型之一。

    3 海底地形反演方法

    Smith 和Sandwell[7]于1994 年提出了依據(jù)線性回歸分析技術(shù)反演海底地形的方法。其基本原理是利用重力數(shù)據(jù)和海底地形在強(qiáng)相關(guān)性波段內(nèi)的線性關(guān)系,反演一定波段的海底地形,并疊加其他波段的船測(cè)水深數(shù)據(jù)得到全波段海底地形模型。由于線性回歸分析技術(shù)會(huì)受到重力數(shù)據(jù)或船測(cè)數(shù)據(jù)中一定異常點(diǎn)的影響,因此采用穩(wěn)健回歸分析法,提高線性回歸分析中比例因子的魯棒性;并提出滑動(dòng)窗口賦權(quán),聯(lián)合由重力異常和重力垂直梯度異常構(gòu)建的模型,通過(guò)對(duì)窗口內(nèi)的兩種模型分配最優(yōu)權(quán)值,構(gòu)建研究海域最終的海底地形模型。下面對(duì)使用的方法進(jìn)行具體介紹。

    3.1 反演波段的選擇

    Parker[20]于1973 年將快速傅里葉變換[21]應(yīng)用到位理論,根據(jù)Parker 公式,在不顧及地殼撓曲均衡影響下,海底地形與重力異常G0(k)之間的關(guān)系可以描述為

    同理,海底地形與重力垂直梯度異常G1(k)之間的關(guān)系可以描述為[18]

    式中,G為地球引力常數(shù);k為徑向頻率(kx,ky)=(1/λx,1/λy) ),其中 (kx,ky)和(λx,λy)分別為x和y方向的頻率和波長(zhǎng);d為平均海深值;ρc和ρw分別為地殼和海水密度;Z(k)為導(dǎo)納函數(shù)。

    由Parker 公式可知,在頻率域內(nèi)重力數(shù)據(jù)和海底地形在一定波段范圍內(nèi)具有相干性,若將重力數(shù)據(jù)和水深看作兩種不同的信號(hào),兩者的相干性可以表示為

    式中,r2(k)為互譜相干函數(shù),當(dāng)r2(k)趨近于1 時(shí),說(shuō)明兩種信號(hào)呈線性相關(guān),趨近于0 時(shí),說(shuō)明兩種信號(hào)完全不相干;G(k)、H(k)分別代表重力信號(hào)和地形信號(hào)的傅里葉變換;G*(k)、H*(k)表示為G(k)、H(k)的復(fù)數(shù)共軛。研究海域內(nèi)重力數(shù)據(jù)與海底地形相干性關(guān)系如圖3 所示。

    圖3 重力數(shù)據(jù)與海底地形的相干性Fig.3 Coherence between gravimetric data and seafloor topography

    重力垂直梯度異常在50~200 km 波長(zhǎng)范圍內(nèi)與海底地形相干性大于0.5,而重力異常數(shù)據(jù)在50~400 km 波長(zhǎng)范圍內(nèi)相干性大于0.5,考慮到在波長(zhǎng)大于200 km 時(shí),會(huì)受到地殼均衡補(bǔ)償影響[22],海底地形產(chǎn)生很少或沒(méi)有重力場(chǎng)變化,因此本文選擇的反演波段為50~200 km。

    3.2 穩(wěn)健回歸分析

    穩(wěn)健回歸分析的原理是基于最小二乘估計(jì)的迭代加權(quán)法,通過(guò)對(duì)權(quán)的優(yōu)選以提高擬合精度與可靠性,使參數(shù)的估值盡可能避免粗差的影響,得到正常模式下的最佳估值[23],對(duì)于線性回歸模型

    Huber 函數(shù)將殘差分為兩組:一組殘差超過(guò)初始標(biāo)準(zhǔn)差的某個(gè)因子c,另一組殘差不超過(guò)這個(gè)限制。u代表標(biāo)準(zhǔn)化的殘差指標(biāo)(ui=vi/σ0),σ0為單位權(quán)中誤差,c表示調(diào)和系數(shù),本文取c=2。參數(shù)a的穩(wěn)健估計(jì)可用迭代加權(quán)最小二乘求解,陳艷國(guó)[25]指出了具體求解步驟。反演波段的重力數(shù)據(jù)和船測(cè)水深數(shù)據(jù)線性擬合結(jié)果如圖4 所示。

    圖4 反演波段重力數(shù)據(jù)與殘余海深線性擬合結(jié)果Fig.4 Results of linear fitting between gravity data of inversion band and residual sea depth

    可以看出,在研究海域,重力數(shù)據(jù)與殘余海深在反演波段表現(xiàn)出了良好的線性關(guān)系。文獻(xiàn)[26]指出,傳統(tǒng)的最小二乘分析相比于穩(wěn)健回歸分析更容易受到離群點(diǎn)的影響,導(dǎo)致其斜率(即比例因子)傾向于離群點(diǎn),而利用Huber 損失函數(shù)的穩(wěn)健回歸分析通過(guò)給予離群點(diǎn)較小的權(quán)值,能夠減小對(duì)離群點(diǎn)的敏感度問(wèn)題,以此提高比例因子的魯棒性。

    3.3 滑動(dòng)窗口賦權(quán)

    滑動(dòng)窗口賦權(quán)是建立在最小二乘賦權(quán)[27]的基礎(chǔ)上,將重力異常反演的模型G0(i)和重力垂直梯度異常反演的模型G1(i)與船測(cè)水深值H(i)進(jìn)行計(jì)算。假設(shè)滑動(dòng)窗口內(nèi)G0(i)所占的權(quán)重為xi,則G1(i)所占的權(quán)重為1?xi,融合得到的G′(i)可表示為

    將融合后的G′(i)和H(i)進(jìn)行比較,以其差值的標(biāo)準(zhǔn)差最小為準(zhǔn)則分配最優(yōu)權(quán)值,計(jì)算公式為

    以此對(duì)滑動(dòng)窗口內(nèi)的G0(i)和G1(i)分配最優(yōu)權(quán)值,進(jìn)而融合得到區(qū)域不同地形起伏條件下的權(quán)值分配。研究海域內(nèi),重力垂直梯度異常模型的權(quán)值分配如圖5 所示。

    圖5 重力垂直梯度異常模型權(quán)值xi 分配示意圖,相應(yīng)的,重力異常模型的權(quán)值分配為1?xiFig.5 The schematic diagram of weight distribution xi of vertical gravity gradient anomalies model.Accordingly,the weight distribution of gravity anomalies model is 1?xi

    基于上述分析,本文反演海底地形模型的基本流程如圖6 所示。

    圖6 海底地形反演流程Fig.6 Flowchart of seafloor topography inversion

    ①通過(guò)相干性分析,獲取重力數(shù)據(jù)和海底地形之間相干性高的波段作為反演波段,將重力數(shù)據(jù)進(jìn)行帶通濾波和向下延拓處理,獲得反演波段的重力數(shù)據(jù)G(x)并插值到船測(cè)控制點(diǎn)處;

    ② 將格網(wǎng)化的船載測(cè)深數(shù)據(jù)進(jìn)行低通濾波得到長(zhǎng)波海深模型,內(nèi)插至船測(cè)控制點(diǎn)處與海深值做差得到控制點(diǎn)處殘余海深;

    ③將步驟①的中波重力數(shù)據(jù)與步驟②的殘余海深以滑動(dòng)窗口方式進(jìn)行穩(wěn)健回歸分析,解算格網(wǎng)點(diǎn)處比例因子S(x)和常數(shù)項(xiàng);

    ④ 將比例因子S(x)和重力數(shù)據(jù)G(x)相乘并疊加常數(shù)項(xiàng)得到反演波段海底地形,與長(zhǎng)波海深模型相加得到中長(zhǎng)波海深模型;

    ⑤ 將中長(zhǎng)波海深模型內(nèi)插至船測(cè)控制點(diǎn)得到控制點(diǎn)處中長(zhǎng)波水深值,與控制點(diǎn)處水深測(cè)量值做差得到控制點(diǎn)處短波水深并格網(wǎng)化,最后與中長(zhǎng)波海底地形模型相加得到全波段海底地形模型;

    ⑥ 采用滑動(dòng)窗口賦權(quán)融合重力異常和重力垂直梯度異常數(shù)據(jù)構(gòu)建的海底地形模型,得到最終的海底地形模型。

    4 模型構(gòu)建結(jié)果及精度評(píng)價(jià)

    4.1 多源海底地形模型的構(gòu)建

    以20′×20′的窗口大小[28],通過(guò)穩(wěn)健回歸分析擬合格網(wǎng)處比例因子和常數(shù)項(xiàng)。分別將重力異常和重力垂直梯度異常構(gòu)建的模型稱為SGA(Single Gravity Anomalies)和SVG(Single Vertical Gradient)。以SGA和SVG 模型與船測(cè)檢核點(diǎn)間標(biāo)準(zhǔn)差最小為準(zhǔn)則,采用滑動(dòng)窗口賦權(quán)法以30′×30′的窗口大小,將SGA 和SVG 模型分配窗口內(nèi)的最優(yōu)權(quán)值,融合構(gòu)建的最終模型稱為MGM(Multiple Gravity Model),反演結(jié)果如圖7 所示。

    圖7 MGM 海底地形模型Fig.7 The seafloor topography model of multiple gravity model

    4.2 與單一重力數(shù)據(jù)源模型比較分析

    為驗(yàn)證本文方法的優(yōu)越性,對(duì)SGA、SVG、傳統(tǒng)最小二乘賦權(quán)法(MGM_LS)和MGM 4 種模型進(jìn)行比較分析,考慮到構(gòu)建以上4 種模型時(shí)均未使用多波束數(shù)據(jù)進(jìn)行約束,且多波束數(shù)據(jù)集中分布在地形起伏變化劇烈海域,因此使用該數(shù)據(jù)作為外部檢核條件,得到的統(tǒng)計(jì)結(jié)果如表1 所示。

    表1 海底地形模型與多波束數(shù)據(jù)差值統(tǒng)計(jì)結(jié)果Table 1 Difference statistics of seafloor topographic model with multibeam data

    以標(biāo)準(zhǔn)差作為質(zhì)量評(píng)估依據(jù),SVG 模型相較于SGA 模型,其精度提高了13.11%,表明在地形起伏變化劇烈的區(qū)域,重力垂直梯度異常能更好的反映地形起伏變化;使用傳統(tǒng)最小二乘賦權(quán)法構(gòu)建的模型MGM_LS,其精度低于SVG 模型,而通過(guò)滑動(dòng)窗口賦權(quán)構(gòu)建的模型MGM 相較于僅使用單一重力數(shù)據(jù)構(gòu)建的SGA 和SVG 模型,精度分別提高了14.92%和2.08%左右,表明本文采用的滑動(dòng)窗口賦權(quán)融合了兩種重力數(shù)據(jù)在不同地形起伏變化下的反演優(yōu)勢(shì),在復(fù)雜海域的水深反演表現(xiàn)更為突出。

    4.3 MGM 精度評(píng)價(jià)

    船測(cè)水深數(shù)據(jù)在預(yù)處理后,可作為研究海域內(nèi)的真值對(duì)反演模型進(jìn)行外部檢核。將MGM 內(nèi)插至船測(cè)水深處與船測(cè)水深值做差,并定義模型與檢核點(diǎn)海深差值與檢核點(diǎn)海深值之比作為相對(duì)誤差,得到的差值分布如圖8 所示。

    結(jié)合圖8a 和圖8b,MGM 與單波束測(cè)深數(shù)據(jù)在絕大部分海域差值較小,與多波束測(cè)深數(shù)據(jù)差值較大的區(qū)域分布在地形起伏變化劇烈的海山鏈、海隆與海溝交界處;根據(jù)直方圖統(tǒng)計(jì),差值大于1 000 m 的僅占4.48%,可能是由于多波束實(shí)測(cè)數(shù)據(jù)集中分布在地形變化幅度大的海域,存在粗差未剔除;由圖8c 可知,在兩側(cè)相對(duì)平坦海域,MGM 相對(duì)誤差保持在3%以下,精度較高;在中心海山鏈區(qū)域,相對(duì)誤差在10%~50%范圍內(nèi)占比約為7%,是模型的主要誤差源。

    圖8 MGM 差值分布圖Fig.8 Difference distribution diagram of multiple gravity model

    為進(jìn)一步驗(yàn)證MGM 的可靠性,引入DTU18 和V23.1 模型進(jìn)行比較分析,得到的統(tǒng)計(jì)結(jié)果如表2所示。

    表2 海底地形模型在檢核點(diǎn)處差值統(tǒng)計(jì)結(jié)果Table 2 Difference statistics of seafloor topographic model at check point

    由表2 知,在研究海域,本文聯(lián)合重力異常和重力垂直梯度異常構(gòu)建的模型MGM,其誤差極值更小,相較于DTU18 和V23.1 模型,其精度提高了約36%。其次,MGM 與船測(cè)數(shù)據(jù)間相關(guān)系數(shù)達(dá)0.998 0,表明MGM 與實(shí)測(cè)水深數(shù)據(jù)更為貼近。由于3 種模型中均存在與船測(cè)水深數(shù)據(jù)差值較大的點(diǎn),但僅占檢核點(diǎn)總數(shù)的極小部分,為直觀看出其分布情況,取相對(duì)誤差大于6%的點(diǎn)位作為精度較低的異常點(diǎn)。分別統(tǒng)計(jì)模型與檢核點(diǎn)差值在400 m 內(nèi)的頻率分布直方圖(圖9)以及精度較低的異常點(diǎn)位分布圖(圖10)。

    圖9 檢核點(diǎn)處差值結(jié)果統(tǒng)計(jì)直方圖Fig.9 Statistical histogram of the difference results at the check point

    圖10 異常點(diǎn)空間位置分布圖(黑色圓點(diǎn)代表異常點(diǎn)位)Fig.10 Spatial distribution map of outliers (black dots represent outliers)

    根據(jù)直方圖可看出,MGM 與船測(cè)水深數(shù)據(jù)一致性良好,兩者差值在300 m 以內(nèi)的占比達(dá)99.42%,優(yōu)于DTU18 和V23.1 模型的98.12%和98.66%,相較于以上兩種模型,MGM 的異常點(diǎn)數(shù)量分別減少了21.18%和17.32%,表明在地形起伏變化劇烈海域,MGM 仍與船測(cè)水深數(shù)據(jù)間有較好的擬合效果?;谇拔姆治觯? 種模型在異常點(diǎn)位的空間位置分布上具有一致性,主要集中在中心海山鏈附近,該部分海域地形起伏變化劇烈,對(duì)反演地形容易造成較大的影響。

    地形剖面圖可以直觀反映出剖面線上的地形起伏情況。為進(jìn)一步探究皇帝海山鏈區(qū)域?qū)Ψ囱菥鹊挠绊?,提取其地形剖面進(jìn)行比較分析[29],結(jié)果如圖11所示。其中藍(lán)色曲線代表MGM 繪制的地形剖面圖,紅色曲線代表船測(cè)水深數(shù)據(jù),可作為該地形剖面的實(shí)際地形。

    圖11 皇帝海山鏈地形剖面圖Fig.11 The topographic profile of the Emperor Seamount Chain

    由圖11 知,皇帝海山鏈區(qū)域的整體地形走勢(shì)清晰明顯,該剖面的海深范圍為984~6 475 m,平均水深約為4 428 m。MGM 與實(shí)測(cè)數(shù)據(jù)的差異主要分布在42°~44°N 海域附近,該區(qū)域地形起伏變化劇烈,海山海溝連續(xù)交錯(cuò),對(duì)反演結(jié)果造成了一定程度的影響。整體上,MGM 與實(shí)測(cè)水深數(shù)據(jù)在描述海山地形走勢(shì)上表現(xiàn)出了良好的一致性,兩者相關(guān)性達(dá)0.996 3。

    綜上所述,文中利用滑動(dòng)窗口賦權(quán)聯(lián)合重力異常與重力垂直梯度異常反演海底地形的方法具有較強(qiáng)的可行性,尤其在地形起伏變化劇烈的海域,可有效提高海底地形反演精度。

    5 結(jié)論

    本文利用重力異常和重力垂直梯度異常數(shù)據(jù),采用滑動(dòng)窗口賦權(quán)和穩(wěn)健回歸分析反演了太平洋皇帝山海域的海底地形模型,將反演結(jié)果與船測(cè)水深數(shù)據(jù)進(jìn)行比較分析,得出以下結(jié)論:

    (1)在海底地形起伏變化劇烈的區(qū)域,依據(jù)重力垂直梯度異常反演的海底地形模型更能準(zhǔn)確地反映地形起伏變化。

    (2)采用滑動(dòng)窗口賦權(quán)聯(lián)合重力異常和重力垂直梯度異常數(shù)據(jù),相比于單一數(shù)據(jù)源能有效提高精度,但在地形起伏變化劇烈區(qū)域的細(xì)節(jié)表現(xiàn)仍有待提高。

    (3)在與船測(cè)檢核點(diǎn)比較中發(fā)現(xiàn),MGM 中精度較低的異常點(diǎn)數(shù)量相較于DTU18 和V23.1 模型有所減少,反映出在地形起伏變化劇烈的區(qū)域,其檢核精度略優(yōu)于DTU18 和V23.1 模型。

    采用回歸分析法反演海底地形,對(duì)船測(cè)數(shù)據(jù)的分布和數(shù)量要求較高,后續(xù)可以針對(duì)不同船測(cè)數(shù)據(jù)的分布情況以及重力異常和重力垂直梯度異常數(shù)據(jù)在不同地形起伏條件下的適用性開(kāi)展進(jìn)一步研究。

    猜你喜歡
    水深賦權(quán)重力
    瘋狂過(guò)山車——重力是什么
    書法靜水深流
    河北水利(2022年10期)2022-12-29 11:48:12
    基于水深分段選擇因子的多光譜影像反演水深
    論鄉(xiāng)村治理的有效賦權(quán)——以A縣扶貧項(xiàng)目為例
    企業(yè)數(shù)據(jù)賦權(quán)保護(hù)的反思與求解
    試論新媒體賦權(quán)
    活力(2019年15期)2019-09-25 07:22:12
    基于改進(jìn)AHP熵博弈賦權(quán)的輸變電工程評(píng)價(jià)
    仰斜式重力擋土墻穩(wěn)定計(jì)算復(fù)核
    一張紙的承重力有多大?
    GPS RTK技術(shù)在水深測(cè)量中的應(yīng)用
    嫁个100分男人电影在线观看| 精品亚洲成a人片在线观看| 99久久99久久久精品蜜桃| 桃花免费在线播放| 亚洲精华国产精华精| 大香蕉久久网| 男人添女人高潮全过程视频| 亚洲精品在线美女| 狂野欧美激情性bbbbbb| 国产亚洲精品一区二区www | 免费一级毛片在线播放高清视频 | 波多野结衣av一区二区av| 久久久久网色| 亚洲va日本ⅴa欧美va伊人久久 | 亚洲精品成人av观看孕妇| 在线十欧美十亚洲十日本专区| 久久人人爽av亚洲精品天堂| 在线av久久热| 亚洲五月婷婷丁香| 捣出白浆h1v1| 美女午夜性视频免费| 亚洲熟女毛片儿| 一区二区三区乱码不卡18| 免费久久久久久久精品成人欧美视频| 亚洲国产av影院在线观看| 后天国语完整版免费观看| 交换朋友夫妻互换小说| 国产精品.久久久| 男女午夜视频在线观看| 国产亚洲一区二区精品| 丰满饥渴人妻一区二区三| 我要看黄色一级片免费的| 国产精品欧美亚洲77777| 欧美成狂野欧美在线观看| 亚洲色图 男人天堂 中文字幕| 大型av网站在线播放| 欧美日韩亚洲综合一区二区三区_| 亚洲精品国产精品久久久不卡| a级片在线免费高清观看视频| 日韩欧美一区视频在线观看| 青春草视频在线免费观看| 女人高潮潮喷娇喘18禁视频| 自线自在国产av| 国产精品99久久99久久久不卡| 青春草亚洲视频在线观看| 婷婷丁香在线五月| 国产精品秋霞免费鲁丝片| 欧美人与性动交α欧美精品济南到| 视频在线观看一区二区三区| 国产成人精品无人区| 捣出白浆h1v1| 蜜桃在线观看..| 好男人电影高清在线观看| 亚洲五月色婷婷综合| 热99国产精品久久久久久7| 久久 成人 亚洲| 国产亚洲午夜精品一区二区久久| 亚洲伊人久久精品综合| 一区二区三区精品91| 人人妻人人澡人人看| 亚洲黑人精品在线| 大陆偷拍与自拍| 欧美激情极品国产一区二区三区| 久热这里只有精品99| 精品一品国产午夜福利视频| 亚洲成人免费av在线播放| 午夜精品国产一区二区电影| 亚洲精品国产区一区二| 水蜜桃什么品种好| 欧美 亚洲 国产 日韩一| 91九色精品人成在线观看| 午夜91福利影院| av国产精品久久久久影院| 高清在线国产一区| 欧美另类亚洲清纯唯美| 99久久精品国产亚洲精品| 91国产中文字幕| a级毛片在线看网站| 99精国产麻豆久久婷婷| 久久久水蜜桃国产精品网| 在线精品无人区一区二区三| a级片在线免费高清观看视频| 国产91精品成人一区二区三区 | 悠悠久久av| 欧美亚洲日本最大视频资源| 手机成人av网站| 黄色视频在线播放观看不卡| 不卡av一区二区三区| 国产免费视频播放在线视频| 脱女人内裤的视频| 欧美日韩国产mv在线观看视频| 亚洲av成人一区二区三| 国产av国产精品国产| xxxhd国产人妻xxx| av在线播放精品| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区在线不卡| av天堂久久9| 欧美久久黑人一区二区| 精品一区在线观看国产| av天堂久久9| 中文字幕av电影在线播放| 国产免费一区二区三区四区乱码| netflix在线观看网站| 女人爽到高潮嗷嗷叫在线视频| 王馨瑶露胸无遮挡在线观看| 无限看片的www在线观看| 国产精品久久久久久精品电影小说| 免费在线观看视频国产中文字幕亚洲 | 欧美大码av| 欧美另类一区| 国产精品欧美亚洲77777| 搡老乐熟女国产| www.熟女人妻精品国产| 免费久久久久久久精品成人欧美视频| 国产亚洲午夜精品一区二区久久| 日韩中文字幕视频在线看片| 99久久精品国产亚洲精品| 亚洲免费av在线视频| 一区二区av电影网| 黄色片一级片一级黄色片| 18禁黄网站禁片午夜丰满| 三级毛片av免费| 视频区图区小说| 亚洲免费av在线视频| 久久久精品94久久精品| 亚洲国产中文字幕在线视频| 菩萨蛮人人尽说江南好唐韦庄| 国产熟女午夜一区二区三区| 99香蕉大伊视频| 精品欧美一区二区三区在线| 国产精品亚洲av一区麻豆| 日韩欧美免费精品| 90打野战视频偷拍视频| 久久久久久久大尺度免费视频| 热re99久久精品国产66热6| 亚洲精品成人av观看孕妇| 欧美成人午夜精品| 999精品在线视频| 纯流量卡能插随身wifi吗| 久久人人97超碰香蕉20202| 欧美+亚洲+日韩+国产| 久久精品国产a三级三级三级| av在线老鸭窝| 午夜影院在线不卡| 精品一品国产午夜福利视频| 日韩欧美一区二区三区在线观看 | 99国产精品一区二区蜜桃av | 狂野欧美激情性bbbbbb| 一区二区三区激情视频| 妹子高潮喷水视频| 国产亚洲精品第一综合不卡| 久久中文看片网| 成人18禁高潮啪啪吃奶动态图| 日本猛色少妇xxxxx猛交久久| 黑人猛操日本美女一级片| 久久午夜综合久久蜜桃| 黄色a级毛片大全视频| 国精品久久久久久国模美| 王馨瑶露胸无遮挡在线观看| 99久久99久久久精品蜜桃| 在线十欧美十亚洲十日本专区| 色老头精品视频在线观看| 老司机靠b影院| 黄片小视频在线播放| 日韩欧美免费精品| 亚洲国产精品一区二区三区在线| 制服诱惑二区| 成人18禁高潮啪啪吃奶动态图| 色播在线永久视频| 狂野欧美激情性bbbbbb| 自线自在国产av| 久久中文看片网| 亚洲精品久久久久久婷婷小说| 国产高清videossex| 丰满少妇做爰视频| 成人av一区二区三区在线看 | 黄片大片在线免费观看| 精品一区二区三区四区五区乱码| 亚洲av电影在线观看一区二区三区| 国产免费视频播放在线视频| 欧美日韩av久久| 久久精品aⅴ一区二区三区四区| a级片在线免费高清观看视频| 日本黄色日本黄色录像| 国产免费福利视频在线观看| 99热网站在线观看| 美女福利国产在线| 动漫黄色视频在线观看| 一区二区日韩欧美中文字幕| 久久女婷五月综合色啪小说| 国产精品.久久久| 欧美日韩亚洲综合一区二区三区_| 天天操日日干夜夜撸| 国产精品二区激情视频| 捣出白浆h1v1| 精品国产一区二区三区四区第35| 蜜桃在线观看..| 亚洲中文av在线| 久久国产精品大桥未久av| 老鸭窝网址在线观看| a在线观看视频网站| 在线观看免费午夜福利视频| 美女主播在线视频| 精品久久久久久电影网| av片东京热男人的天堂| 成人国产av品久久久| 美女大奶头黄色视频| 美女主播在线视频| svipshipincom国产片| 美女福利国产在线| 麻豆乱淫一区二区| 国产日韩欧美亚洲二区| av视频免费观看在线观看| 在线天堂中文资源库| 日韩熟女老妇一区二区性免费视频| 黑丝袜美女国产一区| 免费看十八禁软件| 久久香蕉激情| 亚洲情色 制服丝袜| 一边摸一边做爽爽视频免费| 成年动漫av网址| 国产精品一区二区在线观看99| 久久免费观看电影| 99热国产这里只有精品6| av网站在线播放免费| 午夜成年电影在线免费观看| 久久99热这里只频精品6学生| 99国产精品99久久久久| 亚洲va日本ⅴa欧美va伊人久久 | 巨乳人妻的诱惑在线观看| 国产免费av片在线观看野外av| 久久av网站| 国产精品 国内视频| 亚洲一卡2卡3卡4卡5卡精品中文| 脱女人内裤的视频| 日韩制服丝袜自拍偷拍| 一本一本久久a久久精品综合妖精| 黑丝袜美女国产一区| 亚洲欧美精品自产自拍| 男女下面插进去视频免费观看| 人人妻人人澡人人看| 精品一区在线观看国产| 久久ye,这里只有精品| 久久国产亚洲av麻豆专区| 免费在线观看视频国产中文字幕亚洲 | a 毛片基地| 动漫黄色视频在线观看| 亚洲第一av免费看| 午夜免费观看性视频| 19禁男女啪啪无遮挡网站| 中文字幕人妻丝袜制服| 18在线观看网站| 三级毛片av免费| 亚洲色图综合在线观看| 王馨瑶露胸无遮挡在线观看| a 毛片基地| 在线亚洲精品国产二区图片欧美| 免费黄频网站在线观看国产| a级片在线免费高清观看视频| 午夜福利在线观看吧| 精品欧美一区二区三区在线| 老司机靠b影院| www.自偷自拍.com| 日韩精品免费视频一区二区三区| 中文字幕人妻丝袜一区二区| 亚洲九九香蕉| 亚洲精品久久久久久婷婷小说| 国产一区二区在线观看av| 欧美中文综合在线视频| 免费一级毛片在线播放高清视频 | 99国产精品免费福利视频| 国产免费福利视频在线观看| 免费人妻精品一区二区三区视频| 色婷婷av一区二区三区视频| 男女国产视频网站| 精品久久久精品久久久| 日韩一区二区三区影片| 在线观看人妻少妇| 久久久久精品人妻al黑| 交换朋友夫妻互换小说| 国产精品秋霞免费鲁丝片| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩高清在线视频 | 亚洲精品久久午夜乱码| 欧美日韩av久久| 亚洲精品一区蜜桃| 十八禁网站网址无遮挡| 久久天躁狠狠躁夜夜2o2o| 国产区一区二久久| 啦啦啦在线免费观看视频4| 国产精品久久久久久精品古装| 精品视频人人做人人爽| 一个人免费在线观看的高清视频 | 日本猛色少妇xxxxx猛交久久| 亚洲av片天天在线观看| 亚洲精品在线美女| 国产欧美日韩一区二区三 | 日韩制服骚丝袜av| 人人妻人人添人人爽欧美一区卜| 国产老妇伦熟女老妇高清| 十八禁高潮呻吟视频| 深夜精品福利| 国产黄色免费在线视频| 18禁裸乳无遮挡动漫免费视频| 黄色毛片三级朝国网站| 亚洲中文字幕日韩| 精品一品国产午夜福利视频| 丝袜美腿诱惑在线| 国产一区二区在线观看av| 国精品久久久久久国模美| 亚洲三区欧美一区| 久久狼人影院| 中文字幕另类日韩欧美亚洲嫩草| 制服诱惑二区| 欧美精品av麻豆av| 一级黄色大片毛片| 久久精品熟女亚洲av麻豆精品| 亚洲中文av在线| 久久精品亚洲av国产电影网| 黑人猛操日本美女一级片| 老司机亚洲免费影院| 丝瓜视频免费看黄片| 90打野战视频偷拍视频| 久久性视频一级片| 国产男人的电影天堂91| 久久综合国产亚洲精品| 亚洲国产精品一区三区| 亚洲欧美一区二区三区久久| 老司机深夜福利视频在线观看 | 欧美精品人与动牲交sv欧美| 一进一出抽搐动态| 久久人人爽av亚洲精品天堂| 咕卡用的链子| 超碰成人久久| 女人高潮潮喷娇喘18禁视频| 日本黄色日本黄色录像| 欧美亚洲 丝袜 人妻 在线| 欧美日韩亚洲国产一区二区在线观看 | 黄片大片在线免费观看| 欧美另类一区| 欧美国产精品va在线观看不卡| 自线自在国产av| 免费在线观看完整版高清| 日韩人妻精品一区2区三区| 制服诱惑二区| 午夜免费鲁丝| 一本大道久久a久久精品| 咕卡用的链子| 精品福利观看| 日本五十路高清| 精品人妻熟女毛片av久久网站| 欧美xxⅹ黑人| 精品欧美一区二区三区在线| 丝袜美足系列| 精品少妇一区二区三区视频日本电影| 国产在线一区二区三区精| 80岁老熟妇乱子伦牲交| 免费观看av网站的网址| 女警被强在线播放| av在线播放精品| 国产精品一二三区在线看| 考比视频在线观看| 亚洲av欧美aⅴ国产| 免费久久久久久久精品成人欧美视频| 国产成人a∨麻豆精品| 午夜精品国产一区二区电影| 精品少妇黑人巨大在线播放| 久久精品亚洲av国产电影网| 国产高清videossex| 777久久人妻少妇嫩草av网站| 99国产精品免费福利视频| 视频在线观看一区二区三区| av一本久久久久| 黑人巨大精品欧美一区二区mp4| 国产欧美日韩一区二区三 | 视频区图区小说| 黑人操中国人逼视频| 欧美 日韩 精品 国产| 亚洲国产日韩一区二区| 亚洲黑人精品在线| 欧美日韩亚洲国产一区二区在线观看 | 叶爱在线成人免费视频播放| 亚洲国产av新网站| 一区二区三区四区激情视频| 一区二区三区乱码不卡18| 国产高清视频在线播放一区 | 亚洲综合色网址| 亚洲av成人一区二区三| 亚洲少妇的诱惑av| 美女高潮喷水抽搐中文字幕| 亚洲精品中文字幕一二三四区 | 免费观看人在逋| 一区二区三区精品91| 青青草视频在线视频观看| 国产亚洲精品一区二区www | 亚洲av欧美aⅴ国产| 十分钟在线观看高清视频www| av不卡在线播放| 国产淫语在线视频| 国产精品成人在线| 日韩中文字幕视频在线看片| 欧美精品啪啪一区二区三区 | 久久热在线av| 曰老女人黄片| 午夜福利免费观看在线| 男女国产视频网站| 国产免费一区二区三区四区乱码| 少妇精品久久久久久久| 午夜福利免费观看在线| 中文字幕色久视频| 亚洲av成人一区二区三| 巨乳人妻的诱惑在线观看| 母亲3免费完整高清在线观看| 窝窝影院91人妻| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av国产av综合av卡| 制服人妻中文乱码| 午夜影院在线不卡| av超薄肉色丝袜交足视频| 各种免费的搞黄视频| 久久99热这里只频精品6学生| 亚洲国产毛片av蜜桃av| 97在线人人人人妻| 成人免费观看视频高清| 免费人妻精品一区二区三区视频| 黄片播放在线免费| 激情视频va一区二区三区| 成在线人永久免费视频| 国产精品亚洲av一区麻豆| 每晚都被弄得嗷嗷叫到高潮| 色94色欧美一区二区| 熟女少妇亚洲综合色aaa.| 99久久综合免费| 亚洲欧美激情在线| 最近中文字幕2019免费版| 夜夜夜夜夜久久久久| 97人妻天天添夜夜摸| 国产一卡二卡三卡精品| 夜夜夜夜夜久久久久| 国产亚洲精品一区二区www | 极品少妇高潮喷水抽搐| cao死你这个sao货| 制服人妻中文乱码| 丝袜美足系列| 涩涩av久久男人的天堂| 视频在线观看一区二区三区| 欧美精品av麻豆av| 国产深夜福利视频在线观看| 亚洲av日韩在线播放| 国产1区2区3区精品| 国产xxxxx性猛交| 别揉我奶头~嗯~啊~动态视频 | 国产精品国产av在线观看| 一级,二级,三级黄色视频| 菩萨蛮人人尽说江南好唐韦庄| 水蜜桃什么品种好| 国产亚洲欧美在线一区二区| 欧美黄色淫秽网站| 免费在线观看影片大全网站| 午夜福利免费观看在线| 国产精品欧美亚洲77777| 每晚都被弄得嗷嗷叫到高潮| 天堂中文最新版在线下载| 国产高清国产精品国产三级| 亚洲av成人不卡在线观看播放网 | 亚洲欧洲日产国产| www.自偷自拍.com| 99久久国产精品久久久| 精品高清国产在线一区| 久久精品aⅴ一区二区三区四区| 精品熟女少妇八av免费久了| a 毛片基地| 国产av一区二区精品久久| 色精品久久人妻99蜜桃| 一区二区三区四区激情视频| 欧美+亚洲+日韩+国产| 狠狠婷婷综合久久久久久88av| 欧美黄色片欧美黄色片| 桃花免费在线播放| 精品国产超薄肉色丝袜足j| 国产亚洲av片在线观看秒播厂| 亚洲av国产av综合av卡| 50天的宝宝边吃奶边哭怎么回事| 天天躁狠狠躁夜夜躁狠狠躁| 美女大奶头黄色视频| 国产成人精品久久二区二区免费| 日韩欧美一区二区三区在线观看 | 超色免费av| 91大片在线观看| 新久久久久国产一级毛片| 欧美亚洲 丝袜 人妻 在线| 天堂8中文在线网| 搡老岳熟女国产| a在线观看视频网站| 老司机深夜福利视频在线观看 | tocl精华| 国产成人免费无遮挡视频| www.自偷自拍.com| 老汉色av国产亚洲站长工具| 成人国产av品久久久| 久久这里只有精品19| 午夜福利影视在线免费观看| 人妻一区二区av| 欧美精品一区二区免费开放| 超碰成人久久| 美女中出高潮动态图| 69精品国产乱码久久久| 精品卡一卡二卡四卡免费| 亚洲av成人不卡在线观看播放网 | 亚洲成国产人片在线观看| 一本综合久久免费| 深夜精品福利| 少妇的丰满在线观看| 激情视频va一区二区三区| av有码第一页| 多毛熟女@视频| 亚洲中文字幕日韩| 男女免费视频国产| 久久综合国产亚洲精品| 欧美精品亚洲一区二区| 国产av一区二区精品久久| 99精品久久久久人妻精品| 黄色毛片三级朝国网站| 乱人伦中国视频| 日韩视频一区二区在线观看| www.精华液| 欧美乱码精品一区二区三区| 精品一区二区三卡| av片东京热男人的天堂| 日韩免费高清中文字幕av| 久久久久精品人妻al黑| 色婷婷av一区二区三区视频| 国产精品一二三区在线看| www.av在线官网国产| 国产极品粉嫩免费观看在线| 老熟妇乱子伦视频在线观看 | 国产淫语在线视频| 母亲3免费完整高清在线观看| 狂野欧美激情性xxxx| 黄色视频在线播放观看不卡| 精品人妻熟女毛片av久久网站| 国产成人欧美在线观看 | 人人妻,人人澡人人爽秒播| 满18在线观看网站| 91九色精品人成在线观看| 午夜福利一区二区在线看| 中文字幕色久视频| 精品少妇一区二区三区视频日本电影| 日韩熟女老妇一区二区性免费视频| 一本—道久久a久久精品蜜桃钙片| 精品福利观看| 成年美女黄网站色视频大全免费| 欧美大码av| 精品国产一区二区三区四区第35| a在线观看视频网站| 丝袜脚勾引网站| 亚洲中文av在线| 色视频在线一区二区三区| 国产精品麻豆人妻色哟哟久久| 中文字幕人妻熟女乱码| 亚洲欧美精品综合一区二区三区| 亚洲精品中文字幕一二三四区 | 99精品久久久久人妻精品| 久久中文看片网| 69精品国产乱码久久久| 老司机在亚洲福利影院| 90打野战视频偷拍视频| 午夜免费鲁丝| 欧美日韩亚洲高清精品| 成人手机av| 老司机深夜福利视频在线观看 | 超色免费av| 一二三四在线观看免费中文在| 国内毛片毛片毛片毛片毛片| 国产成人av激情在线播放| 久久精品aⅴ一区二区三区四区| 精品国产乱码久久久久久小说| 欧美国产精品va在线观看不卡| 人人妻人人澡人人爽人人夜夜| 亚洲 欧美一区二区三区| 巨乳人妻的诱惑在线观看| 国产成人影院久久av| 夜夜骑夜夜射夜夜干| 99精品欧美一区二区三区四区| 国产亚洲欧美精品永久| 天天躁夜夜躁狠狠躁躁| 电影成人av| 一区二区三区乱码不卡18| 黑丝袜美女国产一区| 男男h啪啪无遮挡| 宅男免费午夜| av天堂久久9| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频| 狠狠精品人妻久久久久久综合| 别揉我奶头~嗯~啊~动态视频 | 日本精品一区二区三区蜜桃| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av成人不卡在线观看播放网 | 欧美乱码精品一区二区三区| 午夜精品国产一区二区电影| 一级,二级,三级黄色视频| 脱女人内裤的视频| 捣出白浆h1v1| 少妇的丰满在线观看| 精品国产一区二区三区四区第35| 亚洲免费av在线视频| 在线亚洲精品国产二区图片欧美| 一进一出抽搐动态| 一本大道久久a久久精品| 51午夜福利影视在线观看| 中文字幕精品免费在线观看视频| 久久国产亚洲av麻豆专区| 捣出白浆h1v1|