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

    西太平洋拉蒙特平頂海山富鈷結(jié)殼礦區(qū)圈定與資源量估算

    2015-10-24 07:00:08程永壽姜效典張富元孫思軍宋士吉章偉艷
    海洋學(xué)報 2015年1期
    關(guān)鍵詞:結(jié)殼礦塊海山

    程永壽,姜效典,張富元,孫思軍,宋士吉,章偉艷

    (1.中國海洋大學(xué),山東青島266003;2.國家海洋信息中心,天津300171;3.國家海洋局第二海洋研究所國家海洋局海底科學(xué)重點實驗室,浙江杭州310012;4.清華大學(xué)自功化系,北京100084)

    西太平洋拉蒙特平頂海山富鈷結(jié)殼礦區(qū)圈定與資源量估算

    程永壽1,2,姜效典1,張富元3,孫思軍2,宋士吉4,章偉艷3

    (1.中國海洋大學(xué),山東青島266003;2.國家海洋信息中心,天津300171;3.國家海洋局第二海洋研究所國家海洋局海底科學(xué)重點實驗室,浙江杭州310012;4.清華大學(xué)自功化系,北京100084)

    為能科學(xué)合理、快速量化地圈定出大洋海山鈷結(jié)殼優(yōu)質(zhì)礦區(qū),本文基于國際海底管理局提出的礦區(qū)選取模型,利用我國西太平洋海山鈷結(jié)殼資源調(diào)查的拖網(wǎng)采樣資料,結(jié)合西太平洋海山鈷結(jié)殼分布特征,提出人機交互式的礦區(qū)圈定方法,對西太平洋馬爾庫斯-威克海山區(qū)拉蒙特平頂海山圈定出鈷結(jié)殼資源遠景區(qū)和符合國際海底管理局規(guī)章要求的51個鈷結(jié)殼礦塊,估算出拉蒙特平頂海山和51個礦塊的鈷結(jié)殼資源量及金屬量,同時通過對拉蒙特平頂海山日本申請礦區(qū)的結(jié)殼資源分布特征分析,表明采用人機交互式的礦區(qū)圈定方法得出的結(jié)果,不僅能快速定量地圈定出鈷結(jié)殼資源前景較好的礦區(qū),也為大洋海山鈷結(jié)殼礦區(qū)圈定和資源量估算提供了新方法。

    鈷結(jié)殼;資源量;礦區(qū)圈定方法;海洋地質(zhì)

    程永壽,姜效典,張富元,等.西太平洋拉蒙特平頂海山富鈷結(jié)殼礦區(qū)圈定與資源量估算[J].海洋學(xué)報,2015,37(1):106—114,doi:10.3969/j.issn.0253-4193.2015.01.011

    Cheng Yongshou,Jiang Xiaodian,Zhang Fuyuan,et al.Del ineation of cobalt crust blocks and estimation of Co-rich crust resource of Lamont Guyotin the Western Pacific[J].Haiyang Xuebao,2015,37(1):106—114,doi:10.3969/j.issn.0253-4193.2015.01.011

    1 引言

    拉蒙特平頂海山(Lamont Guyot)位于西太平洋馬爾庫斯-威克海山區(qū)(Marcus-Wake Mountains),南面為馬紹爾群島(Marshall Islands),西南為麥哲倫海山區(qū)(Maggellan Seamounts),東部為中太平洋海山區(qū)(Mid-Pacific Mountains)。馬爾庫斯-威克海山區(qū)屬于多熱點成因的板內(nèi)巖漿作用的產(chǎn)物,是富鈷結(jié)殼的一個重要成礦帶,具有“平頂星形海山”的典型形貌特點,平坦山頂以及山脊、山谷相間的山坡形貌[1—3]。富鈷結(jié)殼主要分布于海山、島嶼、海臺和洋脊等水下高地的頂部和斜坡地帶,水深范圍為400~4 000 m,鈷富集地帶一般位于800~2 500 m[1]。富鈷結(jié)殼富含鈷、鎳、鉑、稀土等金屬,其中鈷含量最高達1.2%,是多金屬結(jié)核鈷含量的4倍,比陸地原生鈷礦高幾十倍。1998年8月,在國際海底管理局第4屆會議期間,俄羅斯代表團向會議提出國際海底區(qū)域除多金屬結(jié)核外,還存在著海底熱液硫化物和富鈷鐵錳結(jié)殼資源,根據(jù)《聯(lián)合國海洋法公約》有關(guān)條款(第2章第162條),要求就另兩種資源的勘探制訂規(guī)章規(guī)定。2001年以來,在國際海底管理局一系列研討會和會議中,均把《“區(qū)域”內(nèi)多金屬硫化物和富鈷結(jié)殼探礦和勘探規(guī)章草案》的準備作為重點內(nèi)容。2012年7月,國際海底管理局大會通過《“區(qū)域”內(nèi)富鈷鐵錳結(jié)殼探礦和勘探規(guī)章(ISBA/18/A/11)》(以下簡稱規(guī)章),中國大洋協(xié)會及時提交了西太平洋面積約3 000 k m2鈷結(jié)殼礦區(qū)申請[4]。1997年開始由大洋1號和海洋4號分別對中太平洋、馬紹爾海山區(qū)及麥哲倫海山區(qū)進行了鈷結(jié)殼系統(tǒng)調(diào)查,已調(diào)查了大約20余座海山,分布在10°~25°N的太平洋海山區(qū)[5],并進行了海山鈷結(jié)殼資源量計算和前景評估,優(yōu)選出部分重點海山進行詳細調(diào)查,對鈷結(jié)殼資源潛力不理想的海山暫時放棄調(diào)查,其中包括拉蒙特平頂海山。富鈷結(jié)殼礦區(qū)圈定和資源評價旨在求取評價區(qū)域內(nèi)的資源總量和有用金屬含量,所涉及的參數(shù)主要包括:鈷結(jié)殼的厚度、豐度、覆蓋率、見礦率、含水率,Mn、Cu、Co、Ni等元素含量,礦區(qū)的面積、坡度和水深。這些重要參數(shù)和結(jié)殼資源評價方法,直接影響到富鈷結(jié)殼礦區(qū)的范圍、資源量和礦石質(zhì)量[6—14]。本文在圈定鈷結(jié)殼礦區(qū)的傳統(tǒng)方法基礎(chǔ)上,結(jié)合規(guī)章對礦塊的有關(guān)條款要求,提出了人機交互式的礦區(qū)圈定方法。2013年7月國際海底管理局審議通過了日本提交的3 000 k m2富鈷結(jié)殼勘探申請區(qū),其中包括馬爾庫斯-威克海山區(qū)拉蒙特平頂海山的51個礦塊約1 020 k m2礦區(qū)(ISBA/ 19/LTC/4)。基于收集的拉蒙特平頂海山鈷結(jié)殼測站資料,圈定出面積為1 020 k m2礦區(qū),并估算了拉蒙特平頂海山和礦區(qū)的結(jié)殼資源量。

    2 資料和方法

    本文鈷結(jié)殼樣品來自于西太平洋馬爾庫斯-威克海山區(qū)(21°~22°N,159°~161°E,水深1 600~6 000 m),共23個地質(zhì)采樣站位。拉蒙特平頂海山地形資料取自2011年航次調(diào)查,多波束水深測量設(shè)備型號為SeaBeam2112.360。全覆蓋的水深數(shù)據(jù)用MB-system軟件經(jīng)過各種改正、人工編輯去除異常等后處理,生成DTM,再用Fledermaus三維軟件生成海底地形三維圖像。多波束調(diào)查數(shù)據(jù)可靠,并且精度較高,保證了海底地形模型的真實性。海山坡度為水平面與局部地表之間的正切值。筆者編制軟件模塊利用DTM網(wǎng)格數(shù)據(jù)生成海山的坡度網(wǎng)格數(shù)據(jù),用于后續(xù)計算海底表面積和鈷結(jié)殼資源量。

    多波束水深網(wǎng)格數(shù)據(jù)分析:拉蒙特平頂海山沿NW-SE向延伸,長約123 k m,寬約40 k m,海山頂部水深為1 200 m。水深范圍為1 200~5 000 m,面積約7 000 k m2。拉蒙特平頂海山的山頂可以分成東、西兩部分。西部山頂長59 k m,寬11.1 k m,面積為864 k m2(1 500 m水深等值線圈閉),山頂內(nèi)分布著許多凹凸不平的小海丘。東部山頂長37 k m,寬22 k m,面積為658 k m2(按2 000 m水深等值線圈閉),山頂處存在兩座面積約10 k m2的孤立小海山。東、西兩個山頂形成一個明顯的階梯平臺,高差約500 m,山頂平均坡度小于2°。平坦山頂以下至2 500 m等值線內(nèi)是一明顯的陡坡帶,一般為20°左右。2 500 m以深坡度相對較小,平均值為15°。海山南部坡度較陡,北部坡度漸緩。

    收集的23個測站全部為拖網(wǎng)采樣,站位主要分布在水深1 700~3 300 m斜坡上。拖網(wǎng)設(shè)備不具備水下定位功能,常用船體的位置來代替拖網(wǎng)所在的海底位置,導(dǎo)致測站定位不準確,測站位置與相應(yīng)的水深記錄均存在較大差異。拖網(wǎng)采樣作業(yè)時,均記錄了拖網(wǎng)開始和拖網(wǎng)離開海底時的經(jīng)度、緯度、水深以及鋼纜與大滑輪垂向之間的夾角θ,利用坐標校正公式計算得到的測站站位與船體的經(jīng)度、緯度之差,即可計算出拖網(wǎng)實際所處的地理位置[15]。富鈷結(jié)殼分為板狀結(jié)殼、礫狀結(jié)殼和鈷結(jié)核3種類型[12,16]。拖網(wǎng)采樣使幾種類型的結(jié)殼彼此混雜,如何計算各站結(jié)殼的平均厚度將影響對結(jié)殼分布規(guī)律的認識和結(jié)殼資源量。調(diào)查資料分析認為不同類型的結(jié)殼在空間上出現(xiàn)重疊分布的概率較小,若僅把不同類型結(jié)殼厚度累加作為測站結(jié)殼的平均厚度顯然與自然事實不符。拖網(wǎng)采樣中多種類型結(jié)殼共存時,選擇樣品中具有代表性類型的結(jié)殼厚度作為測站的平均厚度,或?qū)Σ煌愋偷慕Y(jié)殼厚度進行加權(quán)平均。本文采用的測站鈷結(jié)殼厚度計算方法為:測站有板狀結(jié)殼時采用板狀結(jié)殼厚度;沒有板狀結(jié)殼,有礫狀結(jié)殼和鈷結(jié)核時,采用礫狀結(jié)殼厚度;只有鈷結(jié)核測站采用鈷結(jié)核厚度,即取一種結(jié)殼類型的厚度作為測站厚度,3種結(jié)殼厚度不累加計算[12,17]。23個測站拖網(wǎng)采樣有3個測站為結(jié)膜,其厚度小于0.1 cm,結(jié)殼豐度視為0。20個測站結(jié)殼樣品的主要化學(xué)元素含量數(shù)據(jù)均來自于X射線熒光分析法。經(jīng)過處理后,得到了用于結(jié)殼資源量計算的數(shù)據(jù)集,拉蒙特平頂海山鈷結(jié)殼地質(zhì)采樣數(shù)據(jù)包括:站號、經(jīng)度、緯度、水深、基巖、結(jié)殼含水率、結(jié)殼密度、厚度、豐度、和Fe、Mn、Cu、Co、Ni等主要化學(xué)元素百分含量。

    站位數(shù)據(jù)統(tǒng)計結(jié)果表明,拉蒙特平頂海山在淺水區(qū)(1 700~1 900 m)比較多見碳酸鹽巖、磷塊巖、泥巖、燧石、礁灰?guī)r等沉積巖,出現(xiàn)率為57%,深水區(qū)(2 500~3 300 m)則富集火山巖,包括玄武巖、火山角礫巖、火山碎屑巖,出現(xiàn)率為100%。1 900~2 500 m之間缺乏取樣資料,基巖狀況不明。據(jù)現(xiàn)有資料統(tǒng)計,拉蒙特平頂海山火山巖和沉積巖上生長的結(jié)殼厚度基本一致,平均厚度均為2.65 cm,屬薄層型結(jié)殼。拉蒙特平頂海山結(jié)殼常見中厚層單至雙層板狀,結(jié)殼厚度的分布大體以159°50′E為界分東、西兩區(qū),兩區(qū)的結(jié)殼厚度有著明顯的差異(圖1)。159°50′E以西14個測站中板狀結(jié)殼平均厚度為2 cm;以東的9個測站板狀結(jié)殼平均厚度為4 cm,位于坡度相對較陡的上斜坡的2個測站的結(jié)殼厚度大于6 cm,屬厚層板狀結(jié)殼。

    拉蒙特平頂海山鈷結(jié)殼的Co含量為0.50%~ 0.86%,平均值為0.63%(表1),與太平洋其他海山Co的平均百分含量基本一致。Co含量的頻率分布曲線呈單峰式,峰值位于0.50%~0.70%。Co含量大于0.50%和大于0.60%的樣品比例分別為91%和61%,其他主要金屬元素平均值為Fe 16.00%,Mn 21.63%,Ni 0.42%和Cu 0.13%。

    圖1 拉蒙特平頂海山地質(zhì)采樣鈷結(jié)殼厚度直方圖Fig.1 Histogram of cobalt crust sample thickness of Lamont Guyot

    表1 拉蒙特平頂海山富鈷結(jié)殼主要金屬元素含量統(tǒng)計表Tab.1 Main metal contents of Co-rich crusts of Lamont Guyot

    3 遠景區(qū)圈定

    富鈷結(jié)殼礦體在海底呈片狀、環(huán)狀和斑狀3種分布,說明鈷結(jié)殼厚度分布具有顯著的“斑點”特征(不是連續(xù)分布)[18]。利用有效數(shù)據(jù)圈定合理的資源遠景區(qū),是礦區(qū)圈定和資源評價的基礎(chǔ),其準確度也決定了礦區(qū)圈定的合理性和資源評價的精度。遠景區(qū)圈定的任務(wù)是在作為研究對象的結(jié)殼礦體的范圍內(nèi),按照一定的標準參數(shù)劃分出若干單獨的礦體,利用外推法將低至采樣站含礦參數(shù)數(shù)據(jù)推廣到其范圍之外。

    地形對富鈷結(jié)殼覆蓋率、產(chǎn)狀、形態(tài)具有明顯的控制作用,坡度太緩和太陡不利于結(jié)殼的生長發(fā)育[17]。拉蒙特平頂海山的頂面一般較大且平坦,一般認為,山頂由穩(wěn)定的沉積層(有孔蟲砂)構(gòu)成,早期形成的結(jié)殼往往被后期的沉積物覆蓋,因此,結(jié)殼分布的上界通常與沉積物分布的下界(尖滅處)有關(guān)[19]。礦體的上部邊界為平頂海山頂面松散沉積物發(fā)育的下部邊界,礦體的下部邊界為3 000 m等深線。潘家華1))潘家華,等.大洋“十一五”研究課題內(nèi)部報告——大洋富鈷結(jié)殼生長與富集特征研究[R].認為海山斜坡只是結(jié)殼的分布場所,近山頂?shù)摹碍h(huán)狀”分布模式僅為結(jié)殼富集的一般規(guī)律,高品質(zhì)結(jié)殼在海山上主要呈現(xiàn)為“隆起帶效應(yīng)”、“窄脊效應(yīng)”和“坡折區(qū)效應(yīng)”等的富集模式和富集特征。在拉蒙特平頂海山的西北方向緊鄰的所謂尖頂海山上發(fā)現(xiàn)的結(jié)殼類型主要為連續(xù)性好的中厚層板狀結(jié)殼,結(jié)殼的豐度、覆蓋率等也均優(yōu)于此平頂海山。23個測站基本分布于拉蒙特平頂海山的斜坡上,面積多達1 520 k m2山頂上數(shù)據(jù)幾乎為空白,特別是沒有淺鉆站位,無法證實山頂結(jié)殼的真實分布情況。

    本文將拉蒙特平頂海山4 000 m水深作為遠景區(qū)的外邊界,把劃分出若干子區(qū)域分為3種類型(見圖2),即Ⅰ級區(qū)(符合標準查明的資源潛力大的區(qū)域)、Ⅱ級區(qū)(有資源潛力需進一步調(diào)查區(qū))和結(jié)殼空白區(qū)(查明的無結(jié)殼資源潛力的區(qū)域)。站位外推法有2種方案:符合標準的采樣站與不符合標準的采樣站毗鄰,要考慮鈷結(jié)殼參數(shù)的穩(wěn)定性及其有用組分的均勻分布,分析每個采樣站位含礦參數(shù)的影響,確定遠景區(qū)邊界;符合標準的采樣站與無礦段毗鄰,遠景區(qū)邊界通過采樣站之間的中間線。在Arcgis軟件中,疊加地形DTM、水深等值線、站位等圖層,對站位圖層取距站點5 k m為半徑的緩沖區(qū)分析,在4 000 m水深圈閉區(qū)內(nèi)圈出Ⅰ級區(qū)4個,面積為1 531 k m2,2個空白區(qū)面積為275 k m2,Ⅱ級區(qū)面積為2 532 k m2。

    4 基于勘探規(guī)章的礦塊圈定方法

    規(guī)章對結(jié)殼申請涵蓋區(qū)塊的基本要求包括:(1)“鈷結(jié)殼區(qū)塊”是指國際海底管理局規(guī)定的一個或多個網(wǎng)格單元,可以是正方形或長方形,面積不超過20 k m2;(2)每一份請求核準勘探鈷結(jié)殼工作計劃的申請書所涉區(qū)域由不超過150個鈷結(jié)殼區(qū)塊組成,這些區(qū)塊應(yīng)排列為組群;(3)5個毗連鈷結(jié)殼區(qū)塊構(gòu)成一個鈷結(jié)殼區(qū)組群。在任何一點相接觸的2個鈷結(jié)殼區(qū)塊視為毗連區(qū)塊,鈷結(jié)殼區(qū)塊組群不一定毗連但須鄰近,且完全局限在一個不超過550 k m ×550 k m的地理區(qū)域內(nèi)。如何在拉蒙特平頂海山參照上述遠景區(qū)圈定符合規(guī)章要求的若干結(jié)殼區(qū)塊的方法,本文稱為網(wǎng)格微分法和人機交互式礦區(qū)圈定法。

    所謂網(wǎng)格微分法是將相關(guān)圖層設(shè)為WGS84坐標系下的等面積投影,生成覆蓋所有遠景區(qū)的由若干劃分為20 k m2大小的網(wǎng)格組成的目標區(qū)圖層。調(diào)用“createfishnet”和“Feature to Polygon”功能,生成邊長為4 472.13 m的正方形網(wǎng)格礦塊區(qū)域圖層,再和分類遠景區(qū)疊加分析,刪除與遠景區(qū)以外的和完全在結(jié)殼空白區(qū)的區(qū)塊,保留的區(qū)塊逐個分析對比排序,優(yōu)先保留資源潛力大的區(qū)塊,直至最后保留的區(qū)塊數(shù)滿足海山的面積要求為止。網(wǎng)格微分方法操作方便,能夠最大限度地反映眾多地質(zhì)變量包含的礦產(chǎn)資源信息,利于判斷含礦與非礦的判斷。結(jié)殼分布的邊界往往為長寬不定的條帶狀曲線而且很不規(guī)則,特別是隨地形和空間變化劇烈的區(qū)域如資源分布較好區(qū)域與結(jié)殼空白區(qū)交替出現(xiàn)的區(qū)域,對于礦塊面積為不超過20 k m2固定大小礦塊,邊界附近存在的“鋸齒狀”礦塊要滿足規(guī)章條款的選區(qū)模型要求,不可避免地造成落在遠景區(qū)等目標區(qū)以外的不愿意保留的空白區(qū)域面積ΔS過大,從而帶來較大的誤差。連片的區(qū)塊局部整體移動能減小ΔS,但又可能“動一處而引發(fā)全身”。

    如何使ΔS盡量小,使圈定的礦區(qū)的資源好的面積盡量大,本文提出了人機交互式礦區(qū)圈定法(見圖2)。此法借鑒了計算機圖形學(xué)中多邊形的種子填充算法,即把遠景區(qū)視為待填充的多邊形,假設(shè)在多邊形區(qū)域內(nèi)已知一個像素(區(qū)塊),可以通過左、右、上、下4個方向移動的四聯(lián)通區(qū)域,從4個方向?qū)ふ蚁乱粋€區(qū)塊,又稱為四向算法。如圖3,輸入的矩形區(qū)塊A,可選擇區(qū)塊A的左、右、上、下邊的任意一邊為邊界生成下一個矩形區(qū)塊Ai(i=1,3為左右;i=2,4為上下),以區(qū)塊A的左邊界Bl為邊生成區(qū)塊A1,再對區(qū)塊A1在Bl兩端的范圍內(nèi)依據(jù)遠景區(qū)上下移動;區(qū)塊A2、A3和A4類似。本方法的實現(xiàn)還有一個關(guān)鍵就是面積,規(guī)章要求單個區(qū)塊的面積不超過20 k m2。本方法在圖層設(shè)為地理坐標系時,首個區(qū)塊的輸入是選擇矩形的一條邊的邊長為前提,其鄰邊逐步增加若干適當大小的Δb,計算區(qū)塊的地球橢球體表面積[19],直到區(qū)塊面積滿足要求為止,此時得到首個區(qū)塊的坐標。采用此法能保證區(qū)塊面積最大又不超過20 k m2,在地理坐標系下若區(qū)塊采用固定長寬,海山范圍較大時,必帶來區(qū)塊的面積變大或變小。與其他區(qū)塊毗連的區(qū)塊輸入方法說明:設(shè)區(qū)塊A為已知,要以A的左邊Bl為共邊輸入?yún)^(qū)塊A1,選擇Bl為區(qū)塊A1的已知右邊(已知x2,y1,y2,求x1),x1=x2-n×Δb,計算判斷區(qū)塊A1的橢球體面積,直到區(qū)塊A1面積滿足要求(S>20-eps)為止。區(qū)塊A1在y1~y2范圍內(nèi)或向上或向下移動時,仍以A1的右邊為已知邊,計算其他邊的過程與前述類似。

    圖2 拉蒙特平頂海山鈷結(jié)殼遠景區(qū)示意圖Fig.2 Sketch of resource perspective areas distribution of Lamont Guyot

    圖3 人機交互礦區(qū)圈定法示意圖Fig.3 Sketch of model of man-machine interactive mining del ineation resource area

    本文采用人機交互礦區(qū)圈定法在拉蒙特平頂海山根據(jù)遠景區(qū),圈出4個礦塊群和51個礦塊(見圖4b),單個礦塊的地球橢球體面積為20 k m2(長寬均約為4 472 m),總面積為1 020 k m2,51個區(qū)塊平均水深為2 292 m,最大水深為3 500 m;平均坡度為12.3°,最大坡度為21°。拉蒙特平頂海山日本申請區(qū)的單個區(qū)塊東西長為5 k m,南北寬為4 k m,共51個礦塊(見圖4a),總面積為1 020 k m2,平均水深為1 544 m,最大水深為1 990 m,礦塊平均坡度為5.7°,最大坡度為16°,此海山上日本申請區(qū)主要分布在水深1 700 m以淺的平頂海山山頂。左圖中的礦塊編號順序是先從下至上再從左至右排列,其圈礦方法同樣是采用上述網(wǎng)格微分法得到規(guī)則網(wǎng)格后,刪除部分礦塊。

    5 鈷結(jié)殼資源量估算

    鈷結(jié)殼資源量估算和資源評價的方法有多種[6],主要有算術(shù)平均法[7]、地質(zhì)塊段法、加權(quán)平均法、克立格法[7]、最近區(qū)域法[9]、神經(jīng)網(wǎng)絡(luò)法[20]、分形理論法[21—22]和基于網(wǎng)格剖分積分計算富鈷結(jié)殼資源[14]。鈷結(jié)殼資源量估算和資源評價的主要指標有濕結(jié)殼資源量、干結(jié)殼資源量、金屬資源量以及鈷等量金屬量等。結(jié)殼資源評價方法的選取與調(diào)查程度有關(guān),在預(yù)查和普查階段測站數(shù)據(jù)少、網(wǎng)度大,且分布不均勻,應(yīng)用克立格法、鄰近區(qū)域法等地質(zhì)統(tǒng)計法誤差太大,而算術(shù)平均法可滿足評價精度,通常采用算術(shù)平均法和地質(zhì)塊段法;詳查和勘探階段,隨著數(shù)據(jù)量的增加和調(diào)查網(wǎng)度縮小,應(yīng)用克立格法等地質(zhì)統(tǒng)計法可達到較佳效果,獲得最優(yōu)估計結(jié)果并給出估計誤差。從調(diào)查資料分析,拉蒙特平頂海山的結(jié)殼調(diào)查尚處于普查階段,資源評價宜采用算術(shù)平均法,能較好地滿足階段評估精度要求。

    圖4 拉蒙特平頂海山礦塊分布圖Fig.4 Distribution of cobalt blocks of Lamont Guyot

    地質(zhì)塊段法是在算術(shù)平均法的基礎(chǔ)上加以改進的資源量計算方法,方法的原理是將一個礦體投影到一個平面上,根據(jù)礦石種類、品位和儲量級別等地質(zhì)特征的不同,將一個礦體劃分為若干個不同厚度的理想塊體,即塊段,然后使用算術(shù)平均法求出每個塊段的資源量,各塊段資源量的總和即為整個礦體的資源量。按照水深范圍計算結(jié)殼資源量是采用地質(zhì)塊段法對海山資源的初步評價方法之一。資源量計算的算法基本流程:選擇待計算的區(qū)域圖層,添加其面積和坡度屬性,調(diào)用查詢站位圖層或數(shù)據(jù)庫,遍歷區(qū)域圖層的每個區(qū)塊或水深段區(qū)域計算其結(jié)殼厚度、豐度和金屬元素百分含量等要素的平均值,最終獲得礦區(qū)的資源量[23]。對區(qū)域內(nèi)分布有采樣站位數(shù)據(jù)的取其平均值,部分沒有站位數(shù)據(jù)的區(qū)域,暫取海山范圍內(nèi)結(jié)殼站位數(shù)據(jù)的平均值。

    拉蒙特平頂海山1 200~1 600 m和3 000~4 000 m水深段面積較大(表2),沒有結(jié)殼采樣資料的水深段取海山結(jié)殼的均值,該水深段濕結(jié)殼量為4 810×104t和7 661×104t,分別占海山資源量的25%和40%。拉蒙特平頂海山濕結(jié)殼資源量為19 340.4×104t(見表3),干結(jié)殼資源量為13 209.4 ×104t,錳金屬量為2 691.1×104t,銅金屬量為15.7 ×104t,鈷金屬量為80.4×104t,鎳金屬量為49.5× 104t,鈷等量金屬量為158.8×104t。拉蒙特平頂海山圈定的51個礦塊平均坡度為12.3°,平均厚度為2.90 cm,濕結(jié)殼資源量為5 490.97×104t,干結(jié)殼資源量為3 766.41×104t,錳金屬量為775.5×104t,銅金屬量為4.51×104t,鈷金屬量為23.07×104t,鎳金屬量為14.15×104t,鈷等量金屬量為45.66×104t。拉蒙特平頂海山日本申請的51個區(qū)塊平均坡度為5.7°,表明其主要分布在海山山頂,區(qū)塊內(nèi)僅有2個拖網(wǎng)采樣站,結(jié)殼厚度不足1 cm,取海山結(jié)殼站位數(shù)據(jù)的平均值計算其結(jié)殼資源量為5 071.57×104t。

    表2 拉蒙特平頂海山結(jié)殼資源量(按水深)Tab.2 Estimation of cobalt crust resource amount of Lamont Guyot(based on water depth)

    續(xù)表2

    表3 拉蒙特平頂海山圈定礦區(qū)和日本申請區(qū)結(jié)殼資源量Tab.3 Estimation of cobalt crust resource amount of delineation blocks and Japan application blocks in Lamont Guyot

    6 結(jié)論與討論

    (1)基于初查階段有限的調(diào)查資料,以水深4 000 m作為遠景區(qū)的外邊界,將拉蒙特平頂海山劃分為3種類型區(qū)域,即Ⅰ級區(qū)(符合標準查明的資源潛力大的區(qū)域)、Ⅱ級區(qū)(有資源潛力需進一步調(diào)查分析區(qū))和結(jié)殼空白區(qū)(查明的無結(jié)殼資源潛力的區(qū)域)。

    (2)拉蒙特平頂海山平均坡度為11°,結(jié)殼濕資源量為19 340.4×104t,主要分布在1 700~3 000 m斜坡上,西部結(jié)殼資源分布不如東部;基于拖網(wǎng)站位和遠景區(qū)圈定的51個長寬均約為4 472 m、面積為20 k m2區(qū)塊,符合規(guī)章條款要求,其面積為1 020 k m2,占海山面積的14.5%,結(jié)殼平均厚度為2.30~2.90 cm,濕結(jié)殼資源量為19 340.4×104t,資源量基本正確。山頂上資料少,無法判斷結(jié)殼分布上邊界,日本在此海山上的申請面積約為1 020 k m2,分布在水深1 700 m以淺山頂,所計算的資料量的結(jié)果可信度難以判定;日本申請區(qū)位于拉蒙特海山山頂,無法確定為何其沒有全部包含此海山的結(jié)殼資源潛力大的東部斜坡。結(jié)殼分布的上界與沉積物分布的下界(尖滅處)有關(guān)和近山頂?shù)摹碍h(huán)狀”分布模式僅為結(jié)殼富集的一般規(guī)律,需要結(jié)合海山實際具體分析。

    (3)與日本采用的網(wǎng)格微分法得到的結(jié)殼申請區(qū)塊礦區(qū)相比,本文的圈定結(jié)果表明所提出的人機交互式礦區(qū)圈定方法具有明顯優(yōu)勢。基于GIS開發(fā)的國際海域金屬礦產(chǎn)資源礦區(qū)申請輔助軟件實現(xiàn)了人機交互礦區(qū)圈定方法,已經(jīng)應(yīng)用于我國大洋熱液硫化物和富鈷結(jié)殼礦區(qū)申請工作中。

    致謝:何高文、馬維林、宋成兵、楊克紅、任向文等專家都為本文的完成給予了幫助,在此表示由衷感謝!

    [1]Hein JR.Cobalt-rich ferromanganese crusts:global distribution,composition,origin and research activities[C]//International Seabed Authority. Polymetic Massive Sulphides and Cobalt-rich Ferromanganese Crusts.Kingston:Status and Prospects,2002:36-89.

    [2]Hein JR,Conrad TA,Dunham RE.Seamount Characteristics and Mine-Site Model Appl ied to Exploration-and Mining-Lease-Block Selection for Cobalt-Rich Ferromanganese Crusts[J].Marine Georesources&Geotechnology,2009,27(2):160-176.

    [3]任向文,劉季花,石學(xué)法,等.西太平洋Lamont海山中新世以來富鈷結(jié)殼成礦環(huán)境的演化[J].海洋科學(xué)進展,2006,24(1):17-29.

    Ren Xiangwen,Liu Jihua,Shi Xuefa,et al.Evolution of ore-forming condition of co-rich crustsfrom Lamont Guyotin the Western Pacific since the miocene[J].Journal of marine sciences,2006,24(1):17-29.

    [4]中國大洋協(xié)會向國際海底管理局提交首份富鈷結(jié)殼礦區(qū)勘探申請[J].國土資源-科技成果,2012,50(2):32.

    COMRAsubmitted its first cobalt-rich crust mining exploration appl ication to ISBA.Land and resources-scientific and technological achievements,2012,50(2):32.

    [5]崔迎春.中太平洋海區(qū)富鈷結(jié)殼地球化學(xué)特征及成因機制[D].青島:中國科學(xué)院海洋研究所,2008:1-4.

    Cui Yingchun.Geochemical characteristics and origin of cobalt-rich crusts from Central Pacific ocean[D].Qingdao:Institute of oceanology,Chinese Academy of Sciences,2008:1-4.

    [6]周林立,胡光道.富鈷結(jié)殼資源量評估方法[J].海洋地質(zhì)動態(tài),2005,21(6):29-31.

    Zhou Linl i,Hu Guangdao.Methods to evaluate cobalt-rich crust resource[J].Marine geology letters,2005,21(6):29-31.

    [7]張富元.大洋多金屬結(jié)核資源量計算和資源評價方法探討[J].海洋學(xué)報,1995,17(4):85-93.

    Zhang Fuyuan.Research of polymetal l ic nodule resource quantity estimation and resource evaluation method[J].Haiyang Xuebao,1995,17(4):85 -93.

    [8]張富元.克立格法在大洋多金屬結(jié)核資源量計算和資源評價中應(yīng)用[J].海洋學(xué)報,1995,17(4):94-100.

    Zhang Fuyuan.Appl ication of krigging method of polymetal l ic nodule resource quantity estimation and resource evaluation[J].Haiyang Xuebao,1995,17(4):94-100.

    [9]武光海,周懷陽,楊樹鋒.最近區(qū)域法與地質(zhì)塊段法在富鈷結(jié)殼資源量評估中的綜合應(yīng)用[J].海洋地質(zhì)與第四紀地質(zhì),2000,20(4):87-92.

    Wu Guanghai,Zhou Huaiyang,Yang Shufeng.Combined appl ication ofthe methods ofthe nearest domain and geological block to resource evaluation of cobalt-rich crust on a seamountint the Pacific[J].Marine geology&quaternary geology,2000,20(4):87-92.

    [10]馬維林,初鳳友,金翔龍.富鈷結(jié)殼資源評價和圈礦方法探討[J].海洋學(xué)報,2007,29(2):67-73.

    Ma Wei l in,Chu Fengyou,Jin Xianglong.Method approach of resource assessment and ore del ineation for cobalt-rich crust[J].Haiyang Xuebao,2007,29(2):67-73.

    [11]武光海,馬維林,劉捷紅,等.海山富鈷結(jié)殼的圈礦方法[J].海洋學(xué)研究,2005,2(4):15-19.

    Wu Guanghai,Ma Wei l in,Liu Jiehong,et al.Amethod forfinding the ore block of Co-rich crusts on seamounts[J].Journal of marine sciences,2005,2(4):15-19.

    [12]張富元,章偉艷,朱克超,等.鈷結(jié)殼礦區(qū)圈定和資源評價的參數(shù)指標[J].地球科學(xué):中國地質(zhì)大學(xué)學(xué)報,2008,33(2):251-258.

    Zhang Fuyuan,Zhang Weiyan,Zhu Kechao,et al.Parameter and index for del ineation and evaluation of Co-rich crust resources[J].Earth Science-Journal of China University of Geosciences,2008,33(2):251-258.

    [13]楊克紅,初鳳友,馬維林,等.富鈷結(jié)殼多維數(shù)據(jù)模型及礦區(qū)圈定研究[J].海洋學(xué)報,2007,29(6):58-63.

    Yang Kehong,Chu Fengyou,Ma Wei l in,et al.Research on multi-dimension data models and ore area del ineation of cobalt-rich crust[J].Haiyang Xuebao,2007,29(6):58-63.

    [14]吳自銀,初鳳友,馬維林,等.網(wǎng)格剖分積分法計算富鈷結(jié)殼資源量[J].海洋學(xué)報,2007,29(5):81-87.

    Wu Ziyin,Chu Fengyou,Ma Wei l in,et al.Calculation on the resource quantity of cobalt-rich crust using the grid subdivision andintegral method[J].Haiyang Xuebao,2007,29(5):81-87.

    [15]武光海,程永壽,劉捷紅,等.拖網(wǎng)及深海攝像站位位置與水深的初步校正[J].海洋學(xué)研究,2008,26(3):90-97.

    Wu Guanghai,Cheng Yongshou,Liu Jiehong,et al.Position and water depth revisesforthe station of dredge and video equipment[J].Joural of Marine Science,2008,26(3):90-97.

    [16]張海生,趙鵬大,陳守余,等.中太平洋海山多金屬結(jié)殼的成礦特征[J].地球科學(xué):中國地質(zhì)大學(xué)學(xué)報,2001,26(2):205-209.

    Zhang Haisheng,Zhao Pengda,Chen Shouyu,et al.Mineral ization characters of cobalt-rich ferromanganese nodule and crustin central pacific ocean seamount[J].Earth Science-Journal of China University of Geosciences,2001,26(2):205-209.

    [17]Hein JR,Morgenson LA,Clague DA,et al.Cobalt-rich ferromanganese crusts from the Exclusive Economic Zone of the United States and nodules from the oceanic Paci fic,in Geology and Resource Potential of the Continental Margin of Western North America and Adjacent Ocean Basins-Beaufort Sea to Baja Cal i fornia[M]//Circum-Pacif ic Counci lfor Energy and Mineral Resources.Earth Science Series V.6.Houston,1987:753-772.

    [18]楊克紅,初鳳友,馬維林,等.富鈷結(jié)殼資源遠景區(qū)圈定方法[J].海洋地質(zhì)與第四紀地質(zhì),2007,27(6):125-130.

    Yang Kehong,Chu Fengyou,Ma Wei l in,et al.Methods for del ineating the resource perspective area of cobalt-rich crusts[J].Marine Geology and Quaternary Geology,2007,27(6):125-130.

    [19]何高文,梁東紅,宋成兵,等.淺地層剖面測量和海底攝像聯(lián)合應(yīng)用確定平頂海山富鈷結(jié)殼分布界線[J].地球科學(xué):中國地質(zhì)大學(xué)學(xué)報,2005,30(4):509-512.

    He Gaowen,Liang Donghong,Song Chengbing,et al.Determining the distribution boundary of cobalt-rich crusts of guyot by synchronous appl ication of sub-bottom profi l ing and deep-sea video recording[J].Earth Science-Journal of China University of Geosciences,2005,30(4):509-512.

    [20]李裕偉,李林松.鄰域法人工神經(jīng)元網(wǎng)絡(luò)礦產(chǎn)儲量計算[M].北京:地質(zhì)出版社,1997:29-34.

    Li Yuwei,Zhao Jingman,Li Chenyang,et al.Potential mineralresources evaluation methods based on GMS、DSSand GIS[M].Bei j ing:Seismological Press,2007:291-358.

    [21]楊克紅,章偉艷,胡光道,等.分形理論在富鈷結(jié)殼資源量評價方面的應(yīng)用研究[J].地質(zhì)與勘探,2004,40(6):61-64.

    Yang Kehong,Zhang Weiyan,Hu Guangdao,et al.Appl ication offractaltheory in evaluating cobalt-rich incrusttation resource quantity[J].Geology and Prospecting,2004,40(6):61-64.

    [22]何高文.富鈷結(jié)殼分布的分形特征[J].海洋地質(zhì)與第四紀地質(zhì),2001,21(1):89-92.

    He Gaowen.Prel iminary fractal analysis of cobalt-rich crust distribution[J].Marine Geology&Quaternary Geology,2001,21(1):89-92.

    [23]程永壽,馮志勇,殷汝廣,等.基于GIS的富鈷結(jié)殼資源量計算系統(tǒng)開發(fā)[J].海洋通報,2007,26(6):85-90.

    Cheng Yongshou,F(xiàn)eng Zhiyong,Yin Ruguang,et al.Development of the evaluation system for the resource quantity of the cobalt-rich crust based on GIS[J].Marine Science Bulletin,2007,26(6):85-90.

    [24]張富元,章偉艷,朱克超,等.太平洋海山鈷結(jié)殼資源評價[M].北京:海洋出版社,2011:120-121.

    Zhang Fuyuan,Zhang Weiyan,Zhu Kechao,et al.Evaluation offerromanganese crust resources on seamountsin the Pacific Ocean[M].Bei j ing:China Ocean Press,2011:120-121.

    Delineation of cobalt crust blocks and estimation of Co-rich crust resource of Lamont Guyotin the western Pacific

    Cheng Yongshou1,2,Jiang Xiaodian1,Zhang Fuyuan3,Sun Si jun2,Song Shi j i4,Zhang Weiyan3

    (1.Ocean University of China,Qingdao266003,China;2.National Marine Data and Information Service,Tianjin300171,China;3.Key Laboratory of Submarine Geosciences,State Oceanic Ad ministration,The Second Institute of Oceanography,State Oceanic Administration,Hangzhou310012,China;4.Automation Department of Tsinghua University,Beijing100084,China)

    In order to define reasonable and preferable mine area for oceanic Co-rich crust mine appl ication,a new method of man-machine interactive mining del ineation is putforward based on the model of mining selection which was proposed through the International Seabed Authority(ISBA),as well as the crusts distribution of Pacific seamount and the parameter index of mine del ineation.Besides,the prospective areas of Lamont Guyot crusts resources and 51 crust blocks which meet regulations about blocks required by ISBAare del ineated based on the surveying data of Co-rich crust resources on Lamont Guyot by means of dredge haul ing,the amount of Co-rich crust and metal are estimated.Finally,by analyzing the location of areas 51 blocks in this guyot appl ied by Japan in 2013,it is indicated that the results of the man-machine interactive mining del ineation method is more unanimous with the goal of mining del ineation of crusts resource.It provides a new method for oceanic crust cobalt-rich mine appl ication.

    Co-rich crusts;resources amount;del ineating method of cobalt crust resources;marine geology

    P744

    A

    0253-4193(2015)01-0106-09

    2013-07-22;

    2014-09-05。

    中國大洋協(xié)會資助項目(DYXM-125-25-02);海洋公益項目(201005029)。

    程永壽(1971—),男,河南省嵩縣人,博士研究生,主要從事海洋地球物理研究。E-mai l:yshcheng@163.com

    猜你喜歡
    結(jié)殼礦塊海山
    西太平洋鐵錳結(jié)殼中兩類不同成因磷酸鹽的元素特征、形成機制及指示意義
    On the heating mechanism of electron cyclotron resonance thruster immerged in a non-uniform magnetic field
    我就是我,富鈷結(jié)殼
    夢里鮮花開放
    37°女人(2019年4期)2019-04-20 02:17:42
    稀土儲量估算研究
    佐藤海山的詩
    作品(2016年10期)2016-12-06 10:41:54
    迎春花
    大江南北(2016年3期)2016-11-22 07:37:20
    全球三大洋海山鈷結(jié)殼資源量估算
    富鈷結(jié)殼中礦物組成對稀土元素分布的制約
    充填采礦法的礦塊布置與構(gòu)造形式
    在线天堂最新版资源| 亚洲av二区三区四区| 窝窝影院91人妻| 精品99又大又爽又粗少妇毛片 | 欧美一级a爱片免费观看看| 久久久国产成人免费| 国产伦人伦偷精品视频| 丰满人妻一区二区三区视频av| 此物有八面人人有两片| 免费看光身美女| 欧美最黄视频在线播放免费| 国产 一区 欧美 日韩| 成人鲁丝片一二三区免费| 很黄的视频免费| 少妇人妻精品综合一区二区 | 性插视频无遮挡在线免费观看| 国产精品三级大全| 琪琪午夜伦伦电影理论片6080| 亚洲自偷自拍三级| 国产精品av视频在线免费观看| 亚洲av熟女| 成年女人永久免费观看视频| 亚洲综合色惰| 我的女老师完整版在线观看| 亚州av有码| 91麻豆精品激情在线观看国产| 黄片wwwwww| 亚洲成人中文字幕在线播放| 国产av在哪里看| 久久欧美精品欧美久久欧美| 欧美一区二区精品小视频在线| av在线蜜桃| 老熟妇仑乱视频hdxx| .国产精品久久| 婷婷亚洲欧美| 成人国产综合亚洲| 国产免费男女视频| 婷婷亚洲欧美| 国产女主播在线喷水免费视频网站 | 乱系列少妇在线播放| 赤兔流量卡办理| 两个人视频免费观看高清| 在线免费观看的www视频| 国产真实伦视频高清在线观看 | 亚洲内射少妇av| 久久久久久久久大av| 精品国产三级普通话版| 亚洲成人精品中文字幕电影| 欧美又色又爽又黄视频| 三级国产精品欧美在线观看| 在线观看66精品国产| 内射极品少妇av片p| 黄色欧美视频在线观看| 午夜久久久久精精品| 欧美另类亚洲清纯唯美| 国产单亲对白刺激| 免费看日本二区| 美女 人体艺术 gogo| 可以在线观看的亚洲视频| 国产成人影院久久av| 99国产精品一区二区蜜桃av| 观看美女的网站| 女生性感内裤真人,穿戴方法视频| 99久久中文字幕三级久久日本| 久久亚洲精品不卡| 深夜a级毛片| 2021天堂中文幕一二区在线观| 两个人的视频大全免费| 在线观看舔阴道视频| 免费一级毛片在线播放高清视频| 亚洲成人久久爱视频| 成人无遮挡网站| 欧美色欧美亚洲另类二区| xxxwww97欧美| 国产日本99.免费观看| 在线免费十八禁| 麻豆成人午夜福利视频| 日韩亚洲欧美综合| 久久午夜亚洲精品久久| 久久国内精品自在自线图片| 精品欧美国产一区二区三| 高清毛片免费观看视频网站| 国产精品久久视频播放| 免费大片18禁| 国产一区二区三区在线臀色熟女| 亚洲av一区综合| 日本爱情动作片www.在线观看 | 精品久久久久久,| 日韩一本色道免费dvd| 免费在线观看影片大全网站| 国产女主播在线喷水免费视频网站 | 国产三级中文精品| 麻豆久久精品国产亚洲av| 日本黄大片高清| 亚洲第一区二区三区不卡| 欧美日韩黄片免| 久久精品人妻少妇| ponron亚洲| 男人舔奶头视频| 日韩欧美精品v在线| 黄色女人牲交| 哪里可以看免费的av片| 成人特级av手机在线观看| 丰满乱子伦码专区| 亚洲av成人av| 日韩欧美一区二区三区在线观看| 天堂动漫精品| 我要搜黄色片| 欧美成人a在线观看| 真人一进一出gif抽搐免费| 亚洲四区av| 午夜免费激情av| 亚洲四区av| 精品一区二区三区视频在线| 亚洲人成伊人成综合网2020| 夜夜夜夜夜久久久久| 别揉我奶头 嗯啊视频| 色av中文字幕| 日韩欧美国产在线观看| 人妻丰满熟妇av一区二区三区| 搡老熟女国产l中国老女人| 特级一级黄色大片| 亚洲精华国产精华液的使用体验 | 国产av不卡久久| 有码 亚洲区| 国产男人的电影天堂91| 亚洲av日韩精品久久久久久密| 午夜亚洲福利在线播放| 欧美成人性av电影在线观看| 国产高清有码在线观看视频| 国产黄色小视频在线观看| 校园人妻丝袜中文字幕| 欧美在线一区亚洲| 嫩草影院新地址| 一个人观看的视频www高清免费观看| 成人午夜高清在线视频| 国产精品久久久久久av不卡| 一边摸一边抽搐一进一小说| 久久久久九九精品影院| 在现免费观看毛片| 国产人妻一区二区三区在| 啦啦啦啦在线视频资源| 成年女人毛片免费观看观看9| 精品无人区乱码1区二区| x7x7x7水蜜桃| 国内精品久久久久久久电影| 一级av片app| 久久久久久久久久久丰满 | 久久久久国产精品人妻aⅴ院| .国产精品久久| 三级国产精品欧美在线观看| 亚洲欧美日韩无卡精品| 日韩中字成人| 国产精品一及| 69av精品久久久久久| 欧美丝袜亚洲另类 | 日韩欧美在线乱码| 中文字幕人妻熟人妻熟丝袜美| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av天美| 日本一本二区三区精品| 精品久久久久久久人妻蜜臀av| 亚洲性久久影院| 欧美人与善性xxx| 日韩中文字幕欧美一区二区| 国产精品久久久久久亚洲av鲁大| 日韩欧美精品v在线| 亚洲av成人av| 成人午夜高清在线视频| 久久人人爽人人爽人人片va| 久久精品综合一区二区三区| 一个人观看的视频www高清免费观看| 国产精品无大码| 日本色播在线视频| 亚洲人成网站在线播放欧美日韩| 国产女主播在线喷水免费视频网站 | 国产精品一区二区三区四区久久| 欧美精品国产亚洲| 欧美性猛交╳xxx乱大交人| 日日干狠狠操夜夜爽| 午夜福利在线观看吧| 午夜激情福利司机影院| 在线看三级毛片| 99久久成人亚洲精品观看| 精品一区二区三区av网在线观看| 赤兔流量卡办理| 亚洲专区国产一区二区| 在线观看av片永久免费下载| 亚洲自拍偷在线| 国产高清视频在线观看网站| 亚洲av一区综合| 嫁个100分男人电影在线观看| 亚洲欧美激情综合另类| 又黄又爽又刺激的免费视频.| 国产精品不卡视频一区二区| 亚洲av二区三区四区| 成人国产一区最新在线观看| 丰满的人妻完整版| 国产精品三级大全| 97热精品久久久久久| 18禁黄网站禁片午夜丰满| 人妻久久中文字幕网| 精品久久久噜噜| 联通29元200g的流量卡| 亚洲精华国产精华液的使用体验 | 日本 欧美在线| 欧美xxxx性猛交bbbb| 亚洲久久久久久中文字幕| 在线观看舔阴道视频| 亚洲av二区三区四区| 午夜福利在线在线| 日本精品一区二区三区蜜桃| а√天堂www在线а√下载| 熟女人妻精品中文字幕| 直男gayav资源| 国产成人一区二区在线| 在线观看av片永久免费下载| 亚洲综合色惰| 欧美一区二区亚洲| 给我免费播放毛片高清在线观看| 男人狂女人下面高潮的视频| 黄色欧美视频在线观看| 国产免费男女视频| 亚洲美女黄片视频| 久久久久久久午夜电影| 三级毛片av免费| 亚洲第一电影网av| 伦理电影大哥的女人| 国产乱人视频| 国产黄a三级三级三级人| 一进一出抽搐动态| 欧美激情国产日韩精品一区| 一a级毛片在线观看| 亚洲精品久久国产高清桃花| 在线播放国产精品三级| 变态另类成人亚洲欧美熟女| 少妇丰满av| 国产黄色小视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产老妇女一区| 在线免费观看的www视频| 日本欧美国产在线视频| 麻豆一二三区av精品| 最好的美女福利视频网| 欧美日韩中文字幕国产精品一区二区三区| 国产人妻一区二区三区在| 赤兔流量卡办理| av中文乱码字幕在线| 偷拍熟女少妇极品色| 又爽又黄无遮挡网站| 少妇熟女aⅴ在线视频| 五月伊人婷婷丁香| 国国产精品蜜臀av免费| 精品一区二区三区视频在线| 日本一本二区三区精品| 亚洲国产高清在线一区二区三| 国产精品,欧美在线| 久久天躁狠狠躁夜夜2o2o| 91麻豆av在线| 欧美+日韩+精品| 精品人妻视频免费看| 99在线人妻在线中文字幕| 精品欧美国产一区二区三| 精品一区二区三区av网在线观看| 国产熟女欧美一区二区| 亚洲精品日韩av片在线观看| 91在线精品国自产拍蜜月| 别揉我奶头 嗯啊视频| 99久久精品一区二区三区| 亚洲av二区三区四区| 亚洲中文字幕日韩| 亚洲综合色惰| 日韩高清综合在线| 久久久久久久久久久丰满 | 动漫黄色视频在线观看| 九九久久精品国产亚洲av麻豆| 欧美成人性av电影在线观看| 午夜福利成人在线免费观看| 国产精品国产高清国产av| 欧美又色又爽又黄视频| or卡值多少钱| 最近在线观看免费完整版| 亚洲人成网站高清观看| 国内精品美女久久久久久| 亚洲18禁久久av| 99久久精品国产国产毛片| 人妻丰满熟妇av一区二区三区| 一区二区三区激情视频| 国产免费av片在线观看野外av| 欧美黑人巨大hd| 蜜桃久久精品国产亚洲av| 国产精品99久久久久久久久| 深夜a级毛片| 一夜夜www| 久久99热6这里只有精品| 中文字幕高清在线视频| 日日啪夜夜撸| 久久久精品欧美日韩精品| 无人区码免费观看不卡| 91狼人影院| ponron亚洲| 免费高清视频大片| 真人一进一出gif抽搐免费| 一本精品99久久精品77| 麻豆成人av在线观看| 久久这里只有精品中国| 精品久久久久久久末码| 国产精品久久电影中文字幕| 熟妇人妻久久中文字幕3abv| 村上凉子中文字幕在线| 美女xxoo啪啪120秒动态图| 桃色一区二区三区在线观看| 99久久精品一区二区三区| 免费观看在线日韩| 天堂动漫精品| 日韩 亚洲 欧美在线| 小蜜桃在线观看免费完整版高清| 午夜福利在线观看免费完整高清在 | 欧美性猛交黑人性爽| 国产麻豆成人av免费视频| 小蜜桃在线观看免费完整版高清| 此物有八面人人有两片| 国产一区二区三区在线臀色熟女| 亚洲欧美激情综合另类| 亚洲,欧美,日韩| 亚洲专区国产一区二区| 韩国av一区二区三区四区| 国产av在哪里看| 国产伦人伦偷精品视频| 成人二区视频| 国产成人一区二区在线| 波多野结衣高清无吗| av在线天堂中文字幕| 免费看av在线观看网站| 特级一级黄色大片| 国产精品无大码| 99九九线精品视频在线观看视频| 午夜精品久久久久久毛片777| 国产色婷婷99| 又粗又爽又猛毛片免费看| 色综合站精品国产| 在线看三级毛片| 国产精品电影一区二区三区| 悠悠久久av| a级毛片a级免费在线| 夜夜夜夜夜久久久久| 搞女人的毛片| 国内揄拍国产精品人妻在线| 久久亚洲真实| 天堂影院成人在线观看| 国产精品1区2区在线观看.| 黄色女人牲交| 欧美色视频一区免费| 国产激情偷乱视频一区二区| 国产极品精品免费视频能看的| АⅤ资源中文在线天堂| 精品一区二区三区视频在线| 亚洲欧美日韩高清在线视频| 国产午夜福利久久久久久| 亚洲av成人av| 国产真实乱freesex| 免费人成在线观看视频色| 我的老师免费观看完整版| 精品午夜福利在线看| 我要看日韩黄色一级片| 91av网一区二区| 亚洲va日本ⅴa欧美va伊人久久| 一本精品99久久精品77| 欧美在线一区亚洲| 一区福利在线观看| 中文亚洲av片在线观看爽| 2021天堂中文幕一二区在线观| 亚洲无线在线观看| 日本撒尿小便嘘嘘汇集6| 久久久久国产精品人妻aⅴ院| 国内揄拍国产精品人妻在线| 午夜福利在线观看吧| 久久国产精品人妻蜜桃| 国产精品一区二区性色av| 嫩草影院精品99| 97超级碰碰碰精品色视频在线观看| 91狼人影院| 男人舔女人下体高潮全视频| 午夜福利在线观看吧| 啦啦啦观看免费观看视频高清| 免费在线观看影片大全网站| 久久99热这里只有精品18| 亚洲aⅴ乱码一区二区在线播放| 国产av麻豆久久久久久久| 女同久久另类99精品国产91| 亚洲综合色惰| 久久久久久久精品吃奶| 极品教师在线视频| 亚洲av五月六月丁香网| 亚州av有码| 亚洲人成伊人成综合网2020| 亚洲va日本ⅴa欧美va伊人久久| 久久精品夜夜夜夜夜久久蜜豆| 美女cb高潮喷水在线观看| 国产伦精品一区二区三区四那| 亚洲国产精品sss在线观看| 精品久久久久久久末码| 国产av麻豆久久久久久久| 国产精品无大码| www.www免费av| 国产精品人妻久久久久久| 99国产极品粉嫩在线观看| 亚洲av熟女| 精品一区二区三区av网在线观看| 国产男靠女视频免费网站| 一边摸一边抽搐一进一小说| 国产精品人妻久久久久久| 久久精品久久久久久噜噜老黄 | 人妻少妇偷人精品九色| 国产精品99久久久久久久久| 美女黄网站色视频| 97超级碰碰碰精品色视频在线观看| 一级黄色大片毛片| 女生性感内裤真人,穿戴方法视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲成人免费电影在线观看| 久久久久久大精品| 亚洲国产精品久久男人天堂| 色哟哟哟哟哟哟| 欧美在线一区亚洲| 成人精品一区二区免费| 国内揄拍国产精品人妻在线| 一进一出抽搐动态| 国产精品1区2区在线观看.| 亚洲av日韩精品久久久久久密| 人妻久久中文字幕网| 欧美三级亚洲精品| 免费看美女性在线毛片视频| 国产 一区 欧美 日韩| 日韩欧美国产一区二区入口| 最近最新免费中文字幕在线| 精品国内亚洲2022精品成人| 国产一区二区在线观看日韩| 一本一本综合久久| 内射极品少妇av片p| 日日撸夜夜添| 乱系列少妇在线播放| 国产探花极品一区二区| 91麻豆精品激情在线观看国产| 美女高潮的动态| 国产免费一级a男人的天堂| 丝袜美腿在线中文| 国产激情偷乱视频一区二区| 丰满的人妻完整版| 桃色一区二区三区在线观看| 我的老师免费观看完整版| 免费搜索国产男女视频| 国产激情偷乱视频一区二区| 香蕉av资源在线| 麻豆av噜噜一区二区三区| 三级国产精品欧美在线观看| 欧美黑人巨大hd| 精品一区二区三区视频在线观看免费| 国产不卡一卡二| 一本一本综合久久| 伦理电影大哥的女人| 美女黄网站色视频| 欧美另类亚洲清纯唯美| 中文字幕av成人在线电影| 嫩草影院入口| 久久久久免费精品人妻一区二区| 久久中文看片网| 亚洲av不卡在线观看| 久久久午夜欧美精品| 给我免费播放毛片高清在线观看| 亚洲va日本ⅴa欧美va伊人久久| 有码 亚洲区| 观看免费一级毛片| 少妇被粗大猛烈的视频| 美女 人体艺术 gogo| 免费观看在线日韩| 一区二区三区四区激情视频 | 精华霜和精华液先用哪个| 一边摸一边抽搐一进一小说| 欧美激情国产日韩精品一区| 亚洲av美国av| 欧美高清性xxxxhd video| 男女啪啪激烈高潮av片| 亚洲欧美清纯卡通| 婷婷亚洲欧美| 嫩草影院精品99| 精品久久久噜噜| 搡女人真爽免费视频火全软件 | av在线天堂中文字幕| 成人综合一区亚洲| 日日撸夜夜添| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av在线| 真人一进一出gif抽搐免费| 久久香蕉精品热| 欧美成人性av电影在线观看| 韩国av一区二区三区四区| 美女高潮喷水抽搐中文字幕| 最新在线观看一区二区三区| 精品人妻偷拍中文字幕| 97热精品久久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲无线在线观看| 成人美女网站在线观看视频| 亚洲综合色惰| 九九爱精品视频在线观看| 色在线成人网| 韩国av在线不卡| 日韩人妻高清精品专区| 午夜免费激情av| 精品福利观看| 校园人妻丝袜中文字幕| 18+在线观看网站| 日本精品一区二区三区蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 波多野结衣高清无吗| 99精品久久久久人妻精品| 国产一级毛片七仙女欲春2| 男女视频在线观看网站免费| 国产精品亚洲美女久久久| 噜噜噜噜噜久久久久久91| 成人一区二区视频在线观看| 亚州av有码| 国内久久婷婷六月综合欲色啪| 精品久久久久久久久亚洲 | 亚洲成人久久性| 1024手机看黄色片| 成人特级黄色片久久久久久久| 亚洲国产精品成人综合色| 尾随美女入室| 国产一区二区亚洲精品在线观看| 国产高清三级在线| 免费av毛片视频| 特大巨黑吊av在线直播| 精品人妻1区二区| 日韩 亚洲 欧美在线| 日韩欧美 国产精品| 性插视频无遮挡在线免费观看| 日日摸夜夜添夜夜添av毛片 | 天堂av国产一区二区熟女人妻| 国产精品免费一区二区三区在线| 亚洲美女搞黄在线观看 | 黄片wwwwww| 亚洲欧美精品综合久久99| 97人妻精品一区二区三区麻豆| bbb黄色大片| 一进一出抽搐gif免费好疼| 麻豆国产av国片精品| 三级国产精品欧美在线观看| 午夜福利视频1000在线观看| 亚洲色图av天堂| 毛片一级片免费看久久久久 | 日本 av在线| 婷婷丁香在线五月| av在线观看视频网站免费| 97超级碰碰碰精品色视频在线观看| 久久热精品热| 在线免费观看不下载黄p国产 | 在线免费十八禁| 身体一侧抽搐| 日韩欧美精品免费久久| 久久香蕉精品热| 亚洲五月天丁香| 黄色日韩在线| 色精品久久人妻99蜜桃| 小说图片视频综合网站| www.www免费av| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品久久国产高清桃花| 三级男女做爰猛烈吃奶摸视频| 成人无遮挡网站| 中文亚洲av片在线观看爽| 久久九九热精品免费| 99久久久亚洲精品蜜臀av| 成人性生交大片免费视频hd| 一进一出抽搐gif免费好疼| 在线免费观看不下载黄p国产 | bbb黄色大片| 此物有八面人人有两片| 99久久成人亚洲精品观看| 国产精品免费一区二区三区在线| 淫妇啪啪啪对白视频| 精品久久久久久久久av| 99久久精品一区二区三区| or卡值多少钱| 熟妇人妻久久中文字幕3abv| 91在线精品国自产拍蜜月| 白带黄色成豆腐渣| 久久久久九九精品影院| 1000部很黄的大片| 无人区码免费观看不卡| 国产av在哪里看| 最新在线观看一区二区三区| 亚洲成人中文字幕在线播放| 日本在线视频免费播放| 亚洲av美国av| 看片在线看免费视频| 综合色av麻豆| 老司机深夜福利视频在线观看| 嫩草影院新地址| 黄色日韩在线| 日本一二三区视频观看| 国产爱豆传媒在线观看| 两人在一起打扑克的视频| 91av网一区二区| 又爽又黄a免费视频| 哪里可以看免费的av片| 午夜免费男女啪啪视频观看 | 观看美女的网站| 国产黄片美女视频| 俺也久久电影网| 中文字幕av成人在线电影| 久久久国产成人精品二区| 日本 av在线| 亚洲欧美日韩高清在线视频|