周普志,李正元,沈澤中,殷征欣,呂修亞,賀惠忠,董超,2
1. 國家海洋局南海調(diào)查技術(shù)中心,廣東廣州 510275
2. 自然資源部海洋環(huán)境探測技術(shù)與應(yīng)用重點實驗室,廣東廣州 510275
近年來,隨著信息產(chǎn)業(yè)全球一體化建設(shè)的加快,我國海域內(nèi)海底光纜建設(shè)數(shù)量均處于迅速發(fā)展階段。海底管線路由勘察是海底管線工程建設(shè)前期工作的關(guān)鍵環(huán)節(jié)之一[1],為審批路由鋪設(shè)提供了依據(jù),也為海底光纜建設(shè)的路由選擇、施工設(shè)計、鋪設(shè)施工、運營維護(hù)等提供了必需的基礎(chǔ)資料。海底管線勘察中的路由指的是海底光纜在海底從起點到終點的路徑,它的寬度一般在500~2 000 m 之間,呈條帶狀。路由勘察即是沿著預(yù)選路由專門針對海纜工程而進(jìn)行的海洋工程勘察[2]。
由于海底光纜路由勘察海域存在著已鋪設(shè)管線,特別是與之縱橫交錯的已建海底光纜。經(jīng)過了歲月的變遷, 這些海底光纜在海域中的坐標(biāo)是否有了變化,是否還存在也不明確。此外,過去敷設(shè)海底光纜時的定位儀器存在很大的誤差[3]。目前,探測海底光纜的方法有磁法探測、側(cè)掃聲吶探測、淺地層剖面探測和海纜追蹤探測等多種方法[2,4],磁法探測主要是針對掩埋較深,管徑太細(xì)的管線探測[5]。根據(jù)海底光纜的鎧裝鋼絲在地磁背景場中形成的磁信號異常對其進(jìn)行探測識別的裝備精度較高、探測效果較好[6]。隨著高精度定位技術(shù)的應(yīng)用,在深水磁力調(diào)查中引入USBL 進(jìn)行磁力拖魚的定位,能獲得可靠穩(wěn)定的高精度目標(biāo)位置[7]。
海底光纜在淺水段為避免錨害,通常采用中繼型雙層鎧裝光纜,如圖1-2所示。其由纖芯、光信號轉(zhuǎn)發(fā)器、鐵磁性包裹層、阻水保溫層、防壓防腐層等結(jié)構(gòu)組成[8]。
圖1 海底光纜截面圖Fig.1 Submarine optical cable cross section
圖2 海底光纜纖芯Fig.2 Submarine optical cable core
海底光纜產(chǎn)生磁異常的原因有兩種,一是光纜包裹層含有剛性或其它能產(chǎn)生磁異常的鐵磁性材料,光纜就相當(dāng)于一個無限延伸的圓柱體產(chǎn)生的磁異常;二是為了達(dá)到光信號遠(yuǎn)距離傳輸?shù)男枰诠饫|上都要加裝一定數(shù)量的光信號中繼器,這些中繼器必須依靠遠(yuǎn)供電源系統(tǒng)的電流才能正常工作,這個電流為直流或是低頻電流,磁力儀能采集到此類電流產(chǎn)生的磁異常。
因此,海底電纜磁場模型由鐵磁性材料自身產(chǎn)生的圓柱體磁場模型和附加電流磁場模型組成[9],即
其中△T為國際光纜磁異常;△T1為光纜包裹層磁性材料自身產(chǎn)生的磁異常;△T2為遠(yuǎn)供電源系統(tǒng)電流產(chǎn)生的磁異常。
將包裹層磁性材料理想化為無限長水平圓柱體產(chǎn)生的磁異常,如圖3所示。
圖3 海底光纜包裹層磁場模型坐標(biāo)系Fig.3 Magnetic model coordinate system of submarine optical cable integument
假設(shè)海底光纜的橫截面積S,探測高度D,光纜磁導(dǎo)率為μ0,有效磁化強(qiáng)度為Ms,有效磁化傾角為is,地磁傾角為I,測線方位角為A,磁力儀在測線上連續(xù)記錄的理想磁異常分布為
將遠(yuǎn)供系統(tǒng)電流產(chǎn)生的磁異常理想化為無限長載流導(dǎo)線在其周圍空間產(chǎn)生的磁異常,根據(jù)安培定則,線電流產(chǎn)生的磁場強(qiáng)度H為其中i為電流強(qiáng)度,r為海底光纜距離測線上任一測點的垂直距離。
由于磁力儀測得的異常值視為磁場強(qiáng)度H在地磁場T0方向上的分量,遠(yuǎn)供電源系統(tǒng)電流產(chǎn)生的磁異常為
其中α為電流磁場H與水平面的夾角,I為地磁傾角,如圖4所示。
圖4 海底光纜電磁場模型坐標(biāo)系Fig.4 Electromagnetic model coordinate system of submarine optical cable
可見,海底光纜產(chǎn)生的磁異常與探測高度、光纜外徑、電流強(qiáng)度有關(guān)。磁力儀距離海底越近、光纜外徑越大,通過電流強(qiáng)度越大的光纜產(chǎn)生的磁異常越大。
由于光纜磁異常值有一定范圍的幅寬,且因其走向、所處地理緯度以及埋深的不同,實際測量時所得到的二維圖像也是大相徑庭的。磁力儀在海底進(jìn)行光纜探測時測得的磁異常曲線是光纜本身包裹磁性材質(zhì)的異常和遠(yuǎn)供電流產(chǎn)生的異常值的疊加,在圖像上往往是非對稱圖像。
劉勝旋曾提出過對已鋪設(shè)電纜具體位置的圖像定位方法[10],如圖5 所示。其中,5(a)是關(guān)于y軸對稱的磁異常曲線,極值點在航跡線上的投影點就是海底電纜在海平面的投影點。5(b)是關(guān)于原點對稱的異常曲線,曲線的拐點在航跡線上的投影點就是光纜在海平面的投影點。5(c)是非對稱的異常曲線,則根據(jù)其與上述兩種對稱曲線的接近程度來確定光纜的定位點。如與y軸對稱的曲線較接近,則其定位點較接近極值點。5(c)如與原點對稱的曲線較接近,則其定位點較接近拐點,如圖5(d)所示。
海底光纜根據(jù)水深一般有兩種鋪設(shè)方式,在水深大于500 m 的陸坡及海盆區(qū)域,由于受到人類活動影響較少,一般采用直接敷設(shè)在海底的方式,而陸架及近岸即水深小于500 m 的區(qū)域,通常采用挖溝埋設(shè)的方式保證海底光纜安全。精確定位光纜的平面位置和走向需要根據(jù)其埋設(shè)方式和狀態(tài)選擇相應(yīng)的探測方法。
對于直接敷設(shè)的深水區(qū)光纜,可采用側(cè)掃聲吶順光纜走向布設(shè)三條測線,其中一條測線應(yīng)沿路由布置,其余測線布置在路由兩側(cè),測線間距按成圖比例尺不同可設(shè)定為25~100 m[11]。對于埋設(shè)的光纜,由于潮流或其他原因,預(yù)挖溝槽被沖蝕后難以識別,又由于光纜直徑小,側(cè)掃聲吶或淺地層剖面均難以探測到,通常采用磁法探測。首先,應(yīng)根據(jù)路由選擇依據(jù)說明材料,充分了解預(yù)選路由上可能交越的已鋪設(shè)在運營的國際光纜走向,劃定大致的交越點靶區(qū)。然后,垂直于已鋪設(shè)纜線布設(shè)測線,將多條測線測得的磁異常點位置連線,即可確定埋設(shè)光纜的平面位置和走向。
磁法探測通常用拖曳式磁力儀,而水深10 m的海底通信光纜或者動力電纜產(chǎn)生的磁異常寬度范圍通常為75~100 m[12]。測線應(yīng)垂直于已鋪設(shè)電纜走向布設(shè),在交越點探測區(qū)域,為確保磁力儀不受調(diào)查船船磁影響,線長度不小于500 m,測線間距不大于200 m[11],根據(jù)磁力儀實測信號強(qiáng)度和磁異常幅值,適當(dāng)加密測線或增加備用測線。
為保證調(diào)查區(qū)域的定位精度,通常采用星站差分GPS。市場上已經(jīng)得到廣泛應(yīng)用的星站差分系統(tǒng)有三家:美國的Starfire、荷蘭的OmniStar、挪威的VeriPos,可實現(xiàn)分米級甚至厘米級實時定位。
由于海底光纜磁法勘察過程中需要將磁力儀拖曳于船尾,為盡可能消除船磁干擾,一般選擇纜長約為三倍船長。受到流速流向的影響,如果直接使用船載GPS 計算后拖距離勢必產(chǎn)生較大的定位誤差,為獲取磁力儀拖魚的精確位置,需要使用水下導(dǎo)航定位系統(tǒng)。
水下導(dǎo)航定位系統(tǒng)主要有三種形式:超短基線(USBL)定位系統(tǒng)、短基線(SBL)定位系統(tǒng)、長基線(LBL)定位系統(tǒng)。目前在海底光纜路由勘察中,主要使用超短基線導(dǎo)航定位系統(tǒng)為水下拖體提供精確的定位信息。將超短基線聲學(xué)定位系統(tǒng)的換能器通過固定桿安裝在調(diào)查船船舷或者穿艙方式安裝(主流方式是安裝在船舷),船舷式安裝操作相對簡便。調(diào)查作業(yè)時,將水下信標(biāo)綁縛在需水下定位的目標(biāo)上即可,而且可以同時跟蹤多個水下目標(biāo),滿足多載體同時作業(yè)的定位需求。
國際上主流的超短基線主要有英國Sonardyne公司的Gyro USBL 水下聲學(xué)定位系統(tǒng),法國IXBlue公司的GAPS 型全球聲學(xué)定位系統(tǒng),挪威Kongs?berg Marinetime 公司的HiPAP 系列,英國AAE 公司的Easytrack系統(tǒng)系列等等。
光泵磁力儀建立在塞曼效應(yīng)基礎(chǔ)之上,是利用拉莫爾頻率與環(huán)境磁場間精確的比例關(guān)系來測量磁場的。T=K×f,這里f是拉莫爾頻率;K為比例系數(shù);T是地磁場,單位為nT。只要測量拉莫爾頻率f,就可以得到地磁場f的大小。光泵磁力儀靈敏度可達(dá)0.01 nT 或更高,采樣率可達(dá)10 Hz或更高,主要應(yīng)用于對靈敏度要求較高的海洋磁力異常調(diào)查[12]。
目前國內(nèi)用于工程勘察的主流光泵磁力儀為美國GEOMETRICS 公司生產(chǎn)的G-880 型或G-882型銫光泵磁力儀。
本文以2019 年底完成的灣區(qū)互聯(lián)海纜(Bto?BE)路由勘察項目中的磁法勘察和解析過程為例,對勘察過程中的測線布設(shè)、設(shè)備選用及連接、數(shù)據(jù)分析以及交越光纜的判定成圖進(jìn)行具體介紹。
BtoBE 海纜是跨太平洋的海底光纜,主干線由粵港澳大灣區(qū)直達(dá)美國舊金山大灣區(qū),是首條兩個灣區(qū)之間的直達(dá)海纜,并設(shè)有分支延伸至新加坡及菲律賓等東南亞地區(qū)。主干線中S1(HKBU1)自中國香港穿過南海北部陸架、陸坡延伸至南海北部陸坡坡底BU1 連接處,長約570 km。該段路在領(lǐng)海內(nèi)由與東亞交匯海底光纜系統(tǒng)(EAC)、亞太直達(dá)國際海底光纜(APG)、亞洲海底快線(ASE)、環(huán)球光纜(FLAG)等光纜的理論交越點為7個(如圖6所示)。
圖6 BtoBE海纜領(lǐng)海內(nèi)與已建光纜交越情況Fig.6 Intersection between BtoBE submarine cable and built optical cable
根據(jù)收集到的已鋪設(shè)光纜資料,以理論交越點為中心,半徑500 m 區(qū)域內(nèi),垂直于BtoBE 光纜路由,間隔200 m 布設(shè)磁力測線,每個交越點布設(shè)5 條測線,其中最外側(cè)的2 條測線為備用測線,測線長度為1 000 m,相鄰測線走航方向相反,如圖7所示。
圖7 磁力調(diào)查測線分布示意圖Fig.7 Magnetic survey line
為保證水下定位精度,我們采用美國C&C 公司生產(chǎn)的C-Nav3050 星載差分DGNSS 系統(tǒng)進(jìn)行時差分定位,配合使用英國Sonardyne 公司生產(chǎn)的So?nardyne Ranger 2 超短基線系統(tǒng)對磁力儀拖魚進(jìn)行水下實時定位。其中,DGPS 系統(tǒng)的平面精度優(yōu)于0.15 m。USBL其平面精度優(yōu)于斜距的0.5%,有效作用距離1 000 m,使用時將水聲換能器固定于船舷一側(cè)并浸沒海水中使之超過船底深度,水下信標(biāo)用鋼架固定在側(cè)掃拖魚上,如圖8所示。
圖8 水下信標(biāo)及側(cè)掃聲吶Fig.8 Underwater transducer and side-scan sonar
釋放同軸鎧裝纜連接側(cè)掃拖魚,纜長為2.5倍船長以消除船磁影響[13]。再使用9 m 信號纜將磁力儀拖曳于側(cè)掃拖魚尾部,如圖9所示。側(cè)掃拖魚利用自身重量下沉使得磁力儀盡可能的貼近海底光纜,穩(wěn)定拖魚采集姿態(tài)及距離海底高度,提高探測信號強(qiáng)度。
圖9 設(shè)備連接示意圖Fig.9 Magnetic survey equipments connection
磁力探測設(shè)備選用美國GEOMETRICS 公司生產(chǎn)的G-882 型銫光泵磁力儀,如圖10 所示。該型磁力儀雖受測線方向和調(diào)查地區(qū)的緯度影響,但靈敏度高,在南海500 m 深度以內(nèi)近岸海域探測鋪設(shè)光纜的效果明顯。
圖10 G-882銫光泵磁力儀Fig.10 G-882 cesium optical-pumping magnetometer
磁力儀探測時,磁力儀與水下定位系統(tǒng)同步作業(yè),水下定位系統(tǒng)實時記錄磁力拖魚位置,保持磁力儀與海床垂直距離約在10 m,盡可能使用最緩慢測量航速,以獲取精確的數(shù)據(jù)。如果有測線測得磁異常不明顯,則使用調(diào)查區(qū)域最外側(cè)的2條備用測線進(jìn)行復(fù)測。采集勘察過程中使用Ma?gLog 采集軟件采集數(shù)據(jù)(如圖11),每個交越點區(qū)域設(shè)置為一個survey工程文件,每條剖面記錄為一條測線,采樣頻率10 Hz。
圖11 交越點磁異常現(xiàn)場采集Fig.11 Collection of magnetic anomalies at the crossover point
3.4.1 光纜磁異常正演曲線特征為驗證前述圖像定位法在海底光纜磁法探測中的有效性,需要以實際鋪設(shè)的海底光纜特征來構(gòu)建磁性體模型,對其在直流通電狀態(tài)下的磁異常值進(jìn)行正演計算,以確認(rèn)光纜界面在磁異常分布圖像中所處位置。
以BtoBE 海纜勘察為例,對沿垂直于已鋪設(shè)海底光纜的測線進(jìn)行磁異常正演模擬。該海纜淺水磁異??睖y的緯度為22°N 左右,由于以近南北走向的光纜為正演模型時,A=is= 90°,正演模型可簡化為
按照近岸淺水海底海纜鋪設(shè)深度為3 m 的設(shè)計要求,如果磁力儀的探測高度為10 m,而光纜截面直徑通常僅為0.2 ~ 0.3 m。正演結(jié)果曲線如圖12所示。
圖12 走向近南北的海底光纜磁異常正演模型Fig.12 Magnetic forward model of the north-south trend cable
圖像近似對稱,異常寬度約為100 m,幅值約25 nT,光纜位置位于異常圖像幾何中心,靠近正異常幅值頂點處,即坐標(biāo)軸零點位置,此結(jié)果與前圖5(d)海底電纜的磁異常曲線相類似。包裹層磁性材料產(chǎn)生的磁場強(qiáng)度相對于海底光纜遠(yuǎn)供電流所產(chǎn)生的磁場強(qiáng)度而言,幾乎可以忽略,這也正是多數(shù)廢棄光纜不易探測到的原因。
沿BtoBE 光纜S1 段路由進(jìn)行側(cè)掃作業(yè)時,我們對與路由交越的亞太2 號(APCN2)光纜S8A 段進(jìn)行了實驗性探測,該光纜為正南北走向,設(shè)計埋深3 m,與正演模型假設(shè)條件一致,理論交越點水深約400 m,拖魚保持距離海底10 m的高度進(jìn)行探測,交越點位置如圖13所示。
圖13 BtoBE海纜S1段與APCN2海纜S8A段交越點位置Fig.13 BtoBE S1 segment and APCN2 S8A segment intersection point
測線垂直于APCN2 光纜S8A 段,沿正東西向布設(shè),經(jīng)數(shù)據(jù)處理后,磁異常曲線如圖14 所示,和正演模型的磁異常曲線相比較,曲線形態(tài)極為相似,異常寬度略窄,約為75 m,異常幅度高略小,約22 nT,這可能是因為APCN2 光纜供電電流強(qiáng)度略小,或者埋設(shè)深度略深。
圖14 東西向測線異常值曲線Fig.14 Magnetic anomaly curve of EW direction line
3.4.2 數(shù)據(jù)處理過程數(shù)據(jù)處理過程包括異常值提取、磁異常濾波、磁異常平剖圖繪制。因采用拖曳式磁力儀采集,需刪除轉(zhuǎn)彎過程中受到影響的磁力值,分析水下定位數(shù)據(jù)以及對磁力儀的航速、航向進(jìn)行計算,改正9 米纜后拖[14],消除磁力儀后拖對探測定位的影響。將實測磁力值歸算到磁力拖魚位置。計算國際地磁參考場(使用igrf-13 模型)對測線進(jìn)行地磁正常場校正。因每個交越點區(qū)域采集時間短,日變引起的磁場變化相對寬緩,對海底管線探測影響?。?5],文中未對測線磁異常值進(jìn)行日變改正。我們對校正后的測線磁異常值進(jìn)行低通濾波,以消除高頻抖動值。
以BtoBE與APG的交越點測線為例,編制濾波程序代碼,濾波效果如圖15所示。
圖15 BtoBE與APG的交越點測線磁異常值濾波Fig.15 Magnetic anomaly filtered data of the BtoBE and APG crossover point
3.4.3 交越點位置判定成果以濾波后的數(shù)據(jù)為成圖基礎(chǔ)數(shù)據(jù),采用北京金浩林公司浩拓地球物理軟件平臺(version 3.5)進(jìn)行磁力測線平剖圖繪制,坐標(biāo)系選用WGS84,投影方式選用UTM 六度帶,中央經(jīng)度115°E。
以BtoBE與APG的交越點范圍內(nèi)布設(shè)的三條測線為例,經(jīng)異常值提取及濾波后可見產(chǎn)生的磁異常變化幅度分別為28、24、26 nT(圖16)。異常值幅度的變化與磁力儀距離光纜的距離有關(guān)。繪制三條測線磁異常的平面剖面圖。用虛線勾勒連接磁異常突變的剖面幾何中點,即可得到交越光纜APG的準(zhǔn)確走向。
圖16 BtoBE與APG光纜S6段交越點1磁力測線平剖圖(WGS 1984 UTM_ZONE_50N)Fig.16 Magnetic sectional plan of the BtoBE and APG Section 6 intersection point 1(WGS 1984 UTM_ZONE_50N)
按此方法繪制其余6處交越點的已鋪設(shè)光纜異常曲線觀測平面剖面圖并推斷已鋪設(shè)光纜的走向(圖17-22)。從異常平剖圖來看,每條測線磁異常剖面完整,異常寬度在50~100 m 左右。不同的交越點,異常的正負(fù)形態(tài)有所區(qū)別,這主要是因為遠(yuǎn)供系統(tǒng)直流電的供電方向不同。將以上7處交越點判定成果進(jìn)行列表(如表1 所示),表中具體列舉了7處交越點的坐標(biāo)位置以及異常幅值等具體參數(shù)。
表1 交越點磁力探測判定成果列表Table1 Achievement of magnetic detected intersection point
圖17 BtoBE與EAC光纜D段交越點磁力測線平剖圖(WGS 1984 UTM_ZONE_50N)Fig.17 Magnetic sectional plan of the BtoBE and EAC Section D intersection point(WGS 1984 UTM_ZONE_50N)
圖18 BtoBE與APG光纜S6段交越點2磁力測線平剖圖(綠色虛線為附近疑似電纜)(WGS 1984 UTM_ZONE_50N)Fig.18 Magnetic sectional plan of the BtoBE and APG Section 6 intersection point 2(green dotted line is suspected cable)(WGS 1984 UTM_ZONE_50N)
圖19 BtoBE與ASE光纜S4段交越點磁力測線平剖圖(WGS 1984 UTM_ZONE_50N)Fig.19 Magnetic sectional plan of the BtoBE and ASE Section 4 intersection point(WGS 1984 UTM_ZONE_50N)
圖20 BtoBE與EAC光纜C段交越點磁力測線平剖圖(WGS 1984 UTM_ZONE_50N)Fig.20 Magnetic sectional plan of the BtoBE and EAC Section C intersection point(WGS 1984 UTM_ZONE_50N)
圖21 BtoBE與FLAG光纜交越點磁力測線平剖圖(WGS 1984 UTM_ZONE_50N)Fig.21 Magnetic sectional plan of the BtoBE and FLAG intersection point(WGS 1984 UTM_ZONE_50N)
圖22 BtoBE與SMW3光纜交越點磁力測線平剖圖(WGS 1984 UTM_ZONE_50N)Fig.22 Magnetic sectional plan of the BtoBE and SMW3 intersection point(WGS 1984 UTM_ZONE_50N)
本文在介紹海底光纜結(jié)構(gòu)和磁性特征的基礎(chǔ)上,建立了海底光纜的磁異常模型,并由此計算磁異常正演曲線。選擇一條南北走向海底光纜的實測磁異常曲線對正演結(jié)果加以驗證。證明了海底光纜以遠(yuǎn)供系統(tǒng)電流產(chǎn)生的磁場為主,且主要與電流強(qiáng)度和埋設(shè)深度有關(guān)的特征。由正演曲線推導(dǎo)出光纜位置,具體調(diào)查中可依據(jù)已掌握的光纜特征建立正演模型,推算光纜位置與正演曲線的對應(yīng)關(guān)系。在研究分析磁法探測過程中,有如下結(jié)論:
1)盡可能掌握預(yù)選路由上可能交越的已鋪設(shè)電纜的歷史資料是必要條件,測線垂直與電纜理論路由布設(shè)并保證測線布設(shè)密度,采集到的磁力數(shù)據(jù)至少有3 條經(jīng)過已鋪設(shè)光纜的高質(zhì)量剖面是關(guān)鍵。
2)拖曳式磁力儀用于海底光纜探測時,側(cè)掃聲吶拖魚配合磁力儀下沉,能消除船磁影響、穩(wěn)定拖魚采集姿態(tài)及距離海底高度,提高探測信號強(qiáng)度。
3)為消除后拖對磁法探測的影響,使用高精度星站差分GPS 以及水下定位技術(shù),并使用水下定位數(shù)據(jù)模擬磁力儀的航速、航向,改正9 m 纜后拖誤差,能大幅度提高光纜定位的準(zhǔn)確度。
4)判定光纜位置時,可根據(jù)已知光纜大致走向,分析光纜位置與正演曲線特征的關(guān)系,結(jié)合平剖圖勾勒出已鋪設(shè)電纜走向,進(jìn)而獲取交越點坐標(biāo)位置。