王欣凱, 夏小明, 程 林
?
海島界定數(shù)量的變維分形特征——以嵊泗縣為例
王欣凱1, 2, 夏小明2, 程 林3
(1. 浙江大學(xué) 地球科學(xué)學(xué)院, 浙江 杭州 310012; 2. 國(guó)家海洋局第二海洋研究所 杭州 310012; 3. 河北省科學(xué)院地理科學(xué)研究所, 河北 石家莊 050051)
我國(guó)多次海島數(shù)量普查因界定標(biāo)準(zhǔn)不一而結(jié)果相差巨大, 為了滿足海洋主管部門對(duì)海島的管理要求, 亟需研究合理的方法界定海島數(shù)量。利用高分辨率遙感影像解譯嵊泗縣海島2 875個(gè), 以分形理論討論海島數(shù)量與海島面積、海島間距的關(guān)聯(lián)。結(jié)果表明, 嵊泗縣50%的海島面積小于50 m2或島間距不足6 m, 海島數(shù)量與界定標(biāo)準(zhǔn)間呈分維數(shù)先增大后穩(wěn)定的變維分形特征, 分維數(shù)在220 m2和28.5 m之后穩(wěn)定在1.80和2.25。以200~500 m2為界區(qū)分大小島, 小島數(shù)量和海島間距也存在變維分形特征, 且分維數(shù)在間距大于50 m后趨于穩(wěn)定。故嵊泗海島最佳界定方法為以220 m2篩選出計(jì)數(shù)海島后再補(bǔ)充間距大于50 m的較小海島, 該方法下界定嵊泗海島681個(gè)。
海島數(shù)量; 海島面積; 海島間距; 變維分形
我國(guó)是海洋大國(guó), 沿海密布海島, 然而海島數(shù)量卻在數(shù)次全國(guó)性海島普查項(xiàng)目中相去甚遠(yuǎn)。例如: 1976年, 確認(rèn)我國(guó)海島6 536個(gè); 20世紀(jì)80年代末, 確認(rèn)我國(guó)不含港、澳、臺(tái)地區(qū)的海島6 971個(gè)[1]; 最新一輪的“908專項(xiàng)”調(diào)查, 確認(rèn)我國(guó)海島10 312個(gè)[2]。除客觀的海島滅失及新島被發(fā)現(xiàn)外, 海島界定方法的不同是主要原因。當(dāng)前我國(guó)海島數(shù)量的界定方法主要依據(jù)海島面積和海島間距兩個(gè)界定參數(shù), 界定參數(shù)取值越精細(xì)統(tǒng)計(jì)所得海島數(shù)量越多。該特征類似于Mandelbrot提出的分形, 分形理論起源于對(duì)地學(xué)形態(tài)的研究, 其指出分形集中對(duì)象的測(cè)量結(jié)果隨所用尺度變化而不同, 并利用分維數(shù)來描述此類變化特征[3]。以往對(duì)分形特征的研究, 多集中于有具體形態(tài)的對(duì)象, 如岸線、地形等。那么海島數(shù)量這一抽象對(duì)象是否也與界定參數(shù)構(gòu)成分形特征?本文由此入手, 以浙江嵊泗縣海島為例, 利用浙江省航空攝影測(cè)量成果獲取嵊泗縣可識(shí)別海島的地理信息, 依據(jù)海島數(shù)量和界定參數(shù)之間的數(shù)量關(guān)系討論其分形特征, 并提出了依靠分維數(shù)變化規(guī)律界定海島數(shù)量的方法。本文將分形理論推廣到了抽象對(duì)象, 研究結(jié)果為我國(guó)海島數(shù)量的管理提出了科學(xué)的參考依據(jù)。
本文的研究區(qū)域嵊泗海域位于東海的杭州灣外側(cè),所轄海島眾多且距大陸較遠(yuǎn)(圖1), 海岸線人類活動(dòng)較少干預(yù), 多處于自然狀態(tài)。本文依托“浙江908專項(xiàng)”海島岸線專題的成果, 參考2007—2008年的航空遙感影像, 通過人工解譯的方法提取嵊泗縣所有可判識(shí)海島的岸線信息。所用影像主要源于2008年拍攝的地面分辨率0.8 m的浙江省沿海低潮時(shí)航空攝影成果, 并以Google Earth Pro軟件截取的拍攝時(shí)間為2007年3—7月的地面分辨率0.5~0.7 m的影像作為補(bǔ)充。
海岸線為大潮高潮位時(shí)海面所淹沒的水陸界線[4], 由于所用影像中的潮時(shí)各異, 故本文參考國(guó)家海洋局所編規(guī)程和前人的研究[4-6], 依靠高空間分辨率影像的優(yōu)勢(shì), 從海岸線的野外地貌特征進(jìn)行判識(shí)解譯。即人工型海岸線位置位于人工構(gòu)筑物的外緣, 基巖型和沙礫質(zhì)岸線的位置位于巖石和沙灘上顏色由深到淺的分界線。利用ArcGIS 軟件將所用影像和提取岸線統(tǒng)一到CGCS2000坐標(biāo)系并以121°E中央經(jīng)線進(jìn)行高斯-克呂格投影, 再對(duì)海島面積、間距等各類地理信息進(jìn)行計(jì)算。
圖1 研究區(qū)域簡(jiǎn)圖
分形理論指出, 分形對(duì)象的測(cè)量結(jié)果隨所用尺度呈指數(shù)變化[3], 并滿足式(1)中所示關(guān)系。
式中:()為被測(cè)對(duì)象結(jié)果,為測(cè)量尺度,為待定常數(shù),為被測(cè)對(duì)象分維數(shù)。對(duì)式(1)兩邊同時(shí)取對(duì)數(shù), 可得:
即理想狀態(tài)下, 具有分析特征對(duì)象的測(cè)量結(jié)果和尺度將在雙對(duì)數(shù)坐標(biāo)系中投射出一條直線, 且斜率與分維數(shù)相關(guān)。本文將界定參數(shù)和取值對(duì)應(yīng)為測(cè)量尺度, 相應(yīng)海島數(shù)量N和N對(duì)應(yīng)為測(cè)量結(jié)果(), 同樣將結(jié)果繪制于雙對(duì)數(shù)坐標(biāo)系中分析海島數(shù)量的變化趨勢(shì)。同時(shí)依界定參數(shù)升序, 將前30個(gè)數(shù)據(jù)點(diǎn)作為基礎(chǔ)樣本, 利用最小二乘法計(jì)算樣本的分維數(shù), 隨后逐個(gè)將下一數(shù)據(jù)點(diǎn)加入樣本空間計(jì)算分維數(shù), 至樣本包含全部數(shù)據(jù)點(diǎn)。分析上述所得與界定參數(shù)相關(guān)的分維數(shù)序列, 研究海島數(shù)量的分形特征。
本文共識(shí)別嵊泗縣海島2 875個(gè), 相比于“浙江908專項(xiàng)”的508個(gè)[2], 海島數(shù)量增至5.6倍。其中面積大于4 m2的海島2 113個(gè), 最大的為泗礁山島的22.6 km2; 海島間距大于2 m的海島2 586個(gè), 最遠(yuǎn)的為半洋礁的3.9 km。> 2的海島中52.6%位于<7的區(qū)間里,<100區(qū)間中的海島占94.4%。以2個(gè)單位為步長(zhǎng)放大<100的區(qū)間(圖2), 在∈(2, 4)時(shí)海島數(shù)量的分布密度最大, 為443個(gè)。隨著增大, 海島數(shù)量的分布密度迅速減小, 至=20附近時(shí)減小趨勢(shì)停止, 并大致在∈(20, 34)的區(qū)間內(nèi)穩(wěn)定在25附近。當(dāng)>34時(shí), 海島數(shù)量的分布密度進(jìn)一步減小至零星分布。
圖2 基于面積的海島數(shù)量密度分布
海島間距方面的分布情況與之類似,>2的海島中55.8%位于<6的區(qū)間里,<35區(qū)間中的海島占98%。以1個(gè)單位為步長(zhǎng)放大<35的區(qū)間(圖3), 在∈(2, 3)時(shí)海島數(shù)量的分布密度最大, 為447個(gè)。隨著增大, 海島數(shù)量的分布密度迅速減小, 至=20附近時(shí), 減小的速率明顯降低, 隨后進(jìn)一步減小至零星分布。
圖3 基于間距的海島數(shù)量密度分布
依據(jù)海島數(shù)量的統(tǒng)計(jì)結(jié)果, 從海島面積和海島間距兩個(gè)方面, 將界定值和與相應(yīng)海島數(shù)量的自然對(duì)數(shù)繪制成圖4。圖4顯示隨著界定參數(shù)值增大, 界定的海島數(shù)量迅速減小。兩條曲線形狀和變化趨勢(shì)大致相同, 但海島面積所示圖像直線相關(guān)性更好, 斜率變化不大, 而海島間距所指示圖像直線相關(guān)性稍差, 斜率有著明顯的先變陡再變緩過程。
圖4 海島數(shù)量與界定參數(shù)的雙對(duì)數(shù)關(guān)系
依據(jù)本文構(gòu)建的海島面積和海島間距兩個(gè)樣本組, 以樣本組中每個(gè)樣本空間內(nèi)和最大值的自然對(duì)數(shù)為橫坐標(biāo), 樣本分維數(shù)為縱坐標(biāo)繪制曲線(圖5)。圖5顯示代表海島面積和海島間距的ln和ln兩條曲線都表現(xiàn)出分維數(shù)先增大后迅速趨于平穩(wěn)的變維分形特征。在海島面積方面, 拐點(diǎn)為ln=2.69, 對(duì)應(yīng)海島面積為220 m2。2<220 m2時(shí), 樣本分維數(shù)從1.25逐漸增大至1.70;2>220 m2時(shí), 樣本分維數(shù)穩(wěn)定于1.70至1.80之間, 全樣本分維數(shù)為1.80。在海島間距方面, 拐點(diǎn)為ln=3.35, 對(duì)應(yīng)海島間距為28.5 m。即<28.5 m時(shí), 分維數(shù)從1.35逐漸增大至2.25; 當(dāng)>28.5 m時(shí), 樣本分維數(shù)穩(wěn)定于2.25附近, 全樣本分維數(shù)為2.25。
圖5 海島數(shù)量的分維數(shù)與界定參數(shù)的關(guān)系
對(duì)比該兩條曲線, 初始階段海島間距樣本的分維數(shù)比海島面積的約大0.25, 兩者分維數(shù)增大的速率相近, 接著海島面積方面分維數(shù)穩(wěn)定, 海島間距方面分維數(shù)繼續(xù)增大, 兩組分維數(shù)都穩(wěn)定后, 海島間距的分維數(shù)比海島面積約大0.5。
依據(jù)本文計(jì)算結(jié)果, 海島數(shù)量在海島面積和海島間距兩方面都存在著明顯的分形特征, 但分維數(shù)在所研究的尺度范圍內(nèi)并非一直穩(wěn)定, 而是存在著先增大后穩(wěn)定的過程。汪富泉等[8]在研究大自然事物的分形特征時(shí)也發(fā)現(xiàn)了類似現(xiàn)象, 分維數(shù)在一段尺度區(qū)間內(nèi)時(shí)是穩(wěn)定的, 但超過該區(qū)間后分維數(shù)將明顯變化。本文將這種分維數(shù)隨尺度范圍不同而發(fā)生改變的特征稱為變維分形特征。由于同類特征的分維數(shù)相似, 異類特征的分維數(shù)不同[9-10]。本文推斷分維數(shù)隨尺度變化的原因在于不同尺度下事物的特征發(fā)生了改變。這種現(xiàn)象在自然界中是存在的, 如圖6所示, 沙灘表面在米級(jí)的尺度下可視作平面, 放大到毫米級(jí)尺度則是延綿起伏曲面, 放大到微米尺度則脫離二維屬性表現(xiàn)離散特征。在描述沙灘面積時(shí), 以上三個(gè)尺度下的對(duì)象特征和分維數(shù)顯然都會(huì)差異巨大, 表現(xiàn)出變維分形特征。
變維分形特征在分維數(shù)變化的同時(shí)又包含了分維數(shù)的穩(wěn)定部分。正如本文海島數(shù)量的分維數(shù), 在中后段也是穩(wěn)定的, 擁有一致的分形特征。本文將分維數(shù)穩(wěn)定的區(qū)間稱為該分形特征的特征區(qū)間。在特征區(qū)間內(nèi)所有尺度描述的對(duì)象特征是一致的, 結(jié)果隨尺度的變化可由分維數(shù)推算, 屬量變范疇。尺度越過特征區(qū)間, 分維數(shù)發(fā)生改變, 該尺度轉(zhuǎn)而描述對(duì)象新的特征, 結(jié)果隨尺度的變化不符合原分維數(shù)的推算, 發(fā)生質(zhì)變。因此若描述具有變維分形特征的對(duì)象, 并非尺度越精細(xì)越準(zhǔn)確, 因?yàn)槌叨瘸^對(duì)應(yīng)特征區(qū)間的左側(cè)邊界時(shí)對(duì)該特征的描述必然失真。綜上, 特征區(qū)間左邊界是描述對(duì)應(yīng)特征最精確的尺度, 而右邊界則是能有效描述該特征的最粗略尺度。
圖6 不同尺度下的沙灘形態(tài)特征
海島數(shù)量界定目的在于篩選出利用價(jià)值高、海洋權(quán)益大且數(shù)量穩(wěn)定的海島, 以提高海島資源管理的能效, 而界定方法是否合理的關(guān)鍵在于如何確定界定參數(shù)的取值。界定值過大, 納入管理體系的海島數(shù)量偏少, 將造成海島資源的浪費(fèi); 界定值過小, 納入管理體系的低價(jià)值海島偏多, 又會(huì)影響海島的有效管理。
一般來說, 海島的面積越大, 越易于開發(fā)利用; 間距越遠(yuǎn), 給予的海洋權(quán)益越大。依據(jù)基巖海島的地質(zhì)成因, 將本文的海島分為構(gòu)造類和海蝕類兩類。構(gòu)造運(yùn)動(dòng)時(shí)間尺度和空間尺度都較大, 故所形成的海島在短時(shí)期內(nèi)是穩(wěn)定的, 且面積大、島間距遠(yuǎn), 利用價(jià)值也高。反觀海蝕作用形成的海島, 依附于海岸淺灘存在, 海島新增和滅失現(xiàn)象頻發(fā), 故所形成海島數(shù)量不穩(wěn)定, 且面積小、島間距短, 利用價(jià)值較低。因此合理的海島數(shù)量界定方法應(yīng)將構(gòu)造類海島納入海島管理而排出海蝕類。依本文計(jì)算結(jié)果, 海島面積和海島間距兩方面的分維數(shù)都在較大尺度端穩(wěn)定, 其特征區(qū)間的左邊界分別為ln=2.69和ln=3.35, 即=220 m2、=28.5 m。由于構(gòu)造類海島剛好對(duì)應(yīng)著大面積遠(yuǎn)間距的特征, 故可假設(shè)分維數(shù)穩(wěn)定的區(qū)間為構(gòu)造類的特征區(qū)間。也就是說,>220 m2的海島屬于構(gòu)造類;>28.5 m的海島屬于構(gòu)造類。使用小于上述界定值的尺度描述海島, 將漸脫離構(gòu)造類特征進(jìn)入海蝕類特征, 致使海島數(shù)量冗多, 對(duì)海島管理的資源造成浪費(fèi)。因此, 在界定嵊泗縣海島的數(shù)量時(shí), 單獨(dú)使用海島面積合理界定值應(yīng)為220 m2, 單獨(dú)使用海島間距合理界定值應(yīng)為28.5 m。
在實(shí)際的海島數(shù)量界定中, 往往利用多種界定參數(shù)組合判斷, 而非單一的參數(shù)。海島面積因其直觀且與利用價(jià)值關(guān)聯(lián)緊密常作為優(yōu)先指標(biāo)。例如, “浙江908專項(xiàng)”先將滿足面積大于500 m2的較大海島納入數(shù)量統(tǒng)計(jì), 而后再利用50 m間距對(duì)剩余海島進(jìn)行篩選補(bǔ)充。那么優(yōu)先界定參數(shù)的取值大小對(duì)其他輔助界定參數(shù)取值的影響如何呢?本文同樣以海島面積為主海島間距為輔, 以200、300、400和500 m2為面積界值, 分4種情況篩選出面積大于面積界值的海島作為大島, 并以余下小島與大島的間距為界定參數(shù), 計(jì)算小海島的分維數(shù)序列, 結(jié)果如圖7所示。結(jié)果表明余下的小海島數(shù)量與海島間距同樣存在著變維分形特征, 且小海島的分維數(shù)顯著小于全體海島。隨面積界值的增大, 同間距尺度下的分維數(shù)降低, 其中200~300 m2的降幅較大, 300、400和500 m2之間的降幅較小。分維數(shù)穩(wěn)定的區(qū)間方面, 余下海島數(shù)量的特征區(qū)間較之全體海島向右偏移, 且上述4種情況下, 均偏移至ln=3.9(=50 m)附近。結(jié)合上文的討論, 經(jīng)過200~500 m2的面積篩選后, 海島間距大于50 m的海島歸屬構(gòu)造類分形特征。在嵊泗縣海島范圍內(nèi), 使用200~500 m2的面積作為優(yōu)先界定參數(shù)時(shí), 余下海島間距的合理界定值應(yīng)為50 m, 該值恰巧與“浙江908”中的界定值相同。
綜上, 在分形理論下, 所得最佳界定方法為優(yōu)先對(duì)面積大于220 m2的海島計(jì)數(shù), 再將海島間距大于50 m的余下海島補(bǔ)充計(jì)數(shù), 所得嵊泗縣海島界定數(shù)量為681個(gè)。
圖7 面積篩選后海島數(shù)量與海島間距的分維數(shù)
本文通過高分辨率遙感影像提取嵊泗縣海島地理信息, 使用分形理論的研究方法, 以海島面積和海島間距為界定參數(shù), 探討海島界定數(shù)量與界定參數(shù)取值之間的聯(lián)系。隨后, 利用變維分形的思想, 討論合理的界定值大小, 得到了以下結(jié)論。(1)2008年嵊泗縣可影像識(shí)別的海島數(shù)量為2 875個(gè), 其中面積大于4 m2的海島2 113個(gè), 海島間距大于2 m的海島2 586個(gè)。(2)絕大多數(shù)的海島都分布在小面積、短間距的區(qū)間之中, 超過50%的海島面積不超過50 m2, 間距不足6 m, 且隨著面積或間距的增大, 海島分布密度持續(xù)降低。(3)海島界定數(shù)量在海島面積和海島間距兩方面都存在分維數(shù)先增加后穩(wěn)定的變維分形特征, 在海島面積大于220 m2, 海島間距大于28.5m時(shí), 海島界定數(shù)量的分維數(shù)趨于穩(wěn)定。(4)變維分形特征顯示嵊泗縣海島合理的界定方法為先對(duì)面積大于220 m2的海島計(jì)數(shù), 再將余下海島中間距大于50 m的海島補(bǔ)充計(jì)數(shù), 所得數(shù)量為681個(gè)。
本文的完成得益于賈建軍、時(shí)連強(qiáng)、劉毅飛、蔡廷祿、王寧、曹一峰等參與“浙江908專項(xiàng)”海島調(diào)查同志們的辛苦工作, 華東師范大學(xué)賈建軍教授對(duì)文章的修改提供了幫助。在此一并致謝。
[1] 宋小棣. 浙江海島資源綜合調(diào)查與研究[M]. 杭州: 浙江科學(xué)技術(shù)出版社, 1995. Song Xiaodi. Survey and Research of The Island Resource in Zhejiang[M]. Hangzhou: Zhejiang Science and Technology Presss, 1995.
[2] 夏小明, 賈建軍, 陳勇, 等. 中國(guó)海島(礁)名錄[M]. 北京: 海洋出版社, 2012. Xia Xiaoming, Jia Jianjun, Chen Yong, et al. China Islands Directories[M]. Beijing: China Ocean Press, 2012.
[3] Mandelbrot B B. How long is the coast of Britain?[J]. Science, 1967, 156: 636-638.
[4] 國(guó)家海洋局908專項(xiàng)辦公室. 我國(guó)近海海洋綜合調(diào)查與評(píng)價(jià)專項(xiàng)海岸線修測(cè)技術(shù)規(guī)程[M]. 試行本. 北京: 海洋出版社, 2007. State Oceanic Administration, People’s Republic of China, 908 Project Office .Chinese Offshore Investigation and Assessment Project Coastline Measurement Technical Regulation[M]. Beijing: China Ocean Press, 2007.
[5] 孫偉富, 馬毅, 張杰, 等. 不同類型海岸線遙感解譯標(biāo)志建立和提取方法研究[J]. 測(cè)繪通報(bào), 2011, 3: 41- 44. Sun Weifu, Ma Yi, Zhang Jie, et al. Study of remote sensing interpretation keys and extraction technique of different types of shoreline[J]. Bulletin of Surveying and Mapping, 2011, 3: 41-44.
[6] 李靜, 張鷹. 基于遙感測(cè)量的海岸線變化與分析[J]. 河海大學(xué)學(xué)報(bào)(自然科學(xué)版), 2012, 40(2): 224-228. Li Jing, Zhang Ying. Analysis of coastline change based on remote sensing measurement[J]. Journal of Hohai University (Natural Sciences), 2012, 40(2): 224-228.
[7] 國(guó)家海洋局908專項(xiàng)辦公室. 海島界定技術(shù)規(guī)程[M]. 北京: 海洋出版社, 2011. State Oceanic Administration, People’s Republic of China, 908 Project Office. Specifications for Island Definition Techniques[M]. Beijing: China Ocean Press, 2011.
[8] 汪富泉, 李后強(qiáng). 分形——大自然的藝術(shù)構(gòu)造[M]. 濟(jì)南: 山東教育出版社, 1996: 7-36. Wang Fuquan, Li Houqiang. Fractal: the Art Structure of Nature[M]. Ji’nan: Shangdong Education Press, 1996: 7-36.
[9] 蘇媛媛, 李欽幫, 侯敏. 黃河三角洲海岸線分形研究[J].資源與產(chǎn)業(yè), 2008, 10(6): 103-107. Su Yuanyuan, Li Qinbang, Hou Ming. Stydy on coastal fractal of Yellow River delta[J]. Reaources & Industries, 2008, 10(6): 103-107.
[10] 王欣凱, 賈建軍. 浙江海島岸線分維數(shù)與其人工岸線比例的關(guān)系[J]. 海洋學(xué)研究, 2011 , 29(4): 25-37. Wang Xinkai, Jia Jianjun. Relationship between the fractal dimension and percentage of artificial coastlines of islands in Zhejiang Province[J]. Journal of Marine Sciences, 2011, 29(4): 25-37.
(本文編輯: 劉珊珊)
The variable fractal dimensions characterizing of the statistical number of islands — a case study of Shengsi County
WANG Xin-kai1, 2, XIA Xiao-ming2, CHENG Lin3
(1. School of Earth Sciences, Zhejiang University, Hangzhou 310012, China; 2. The Second Institution of Oceanography, State Oceanic Administration, Hangzhou 310012, China; 3. Institute Geographical Sciences, Hebei Academy of Sciences, Shijiazhuang 050051, China)
It is imperative that the oceanic authority who manages the islands count them. However, because of the inconsistencies in statistical standards, different numbers of the islands in China can be attained. This paper examines the islands of Shengsi County as a case study and uses the fractal theory to determine the most suitable standard by which to count islands. From high resolution remote-sensing images, 2 875 islands were detected visually in Shengsi County; more than half of these are smaller than 50 m2in area or with a separation distance shorter than 6 m. In both aspects of island area and separation distance, the relevant statistics for island numbers display the characteristics of variable fractal dimensions. At first the fractal dimensions increase then stabilize at 1.80 and 2.25 for areas larger than 220 m2and distances greater than 28.5 m, respectively. For islands separated by 200, 300, 400 and 500 m2, the number of smaller islands and their separation distances from the larger islands also display the characteristics of variable fractal dimensions. And for separation distances greater than 50 m, the fractal dimensions are stable for each of the four cases previously discussed. Therefore, the rational standard for counting islands in Shengsi is by screening out those islands larger than 220 m2and then by supplementing them with those islands whose separation distance is larger than 50 m. Using this method, the number of islands in Shengsi is 681.
number of islands; area of islands; distance of islands; variable fractal dimension
Dec. 23, 2016
王欣凱(1986-), 男, 湖北武漢人, 博士在讀, 主要從事海岸帶、海島調(diào)查與管理研究, 電話: 0571-81963500, E-mail: carlinwxk@163.com
P714
A
1000-3096(2017)09-0130-06
10.11759/hykx20161223001
2016-12-23;
2017-03-19
國(guó)家自然科學(xué)基金(41576095); 我國(guó)近海海洋綜合調(diào)查與評(píng)價(jià)專項(xiàng)(908-ZC- I-09, ZJ908-01-02)
[National Science Foundation of China, No. 41576095; Chinese Offshore Investigation and Assessment Program, No. 908-ZC-I-09, No. ZJ908-01-02]