劉 焱，王 媛，辛凱強(qiáng)，孫晨紅

(1.北京中科航天人才服務(wù)有限公司陜西省分公司，陜西 西安 710000;2.河北師范大學(xué)地理科學(xué)學(xué)院，河北 石家莊 050024;3.陜西航天技術(shù)應(yīng)用研究院有限公司，陜西 西安 710000;4.自然資源陜西省衛(wèi)星應(yīng)用技術(shù)中心，陜西 西安 710050)

0 引言

石油是工業(yè)血液，經(jīng)濟(jì)發(fā)展的源動(dòng)力。但石油致害力也高，受暴雨等自然災(zāi)害影響較大，洪水、油氣滲漏等會(huì)給油氣開(kāi)采帶來(lái)潛在威脅，同時(shí)也會(huì)造成生態(tài)環(huán)境的破壞[1-2]。加強(qiáng)對(duì)石油生產(chǎn)設(shè)施的核查、監(jiān)測(cè)力度對(duì)石油勘探開(kāi)發(fā)具有重要的意義[3-4]。遙感衛(wèi)星技術(shù)作為信息獲取的一種重要技術(shù)手段，具備動(dòng)態(tài)監(jiān)測(cè)的能力，在各個(gè)領(lǐng)域發(fā)揮越來(lái)越大的作用，通過(guò)遙感圖像可以從宏觀上識(shí)別地物目標(biāo)，實(shí)現(xiàn)遙感專題信息的提取工作，例如油田油氣微滲漏異常信息、運(yùn)輸管線和交通線路等地面工程解譯，以及災(zāi)害危險(xiǎn)信息等。本文提出的油田遙感信息提取系統(tǒng)以遙感、GIS為基礎(chǔ)，以油區(qū)道路、水體以及油氣微滲漏提取為核心，能夠方便快捷地進(jìn)行數(shù)據(jù)處理、分析及專題信息提取，為油田勘探開(kāi)發(fā)提供決策支持。

1 油田遙感信息提取系統(tǒng)設(shè)計(jì)

1.1 系統(tǒng)開(kāi)發(fā)模式

在綜合比較各種遙感圖像處理系統(tǒng)開(kāi)發(fā)模式的基礎(chǔ)上，選擇開(kāi)源QGIS進(jìn)行二次開(kāi)發(fā)。QGIS作為一個(gè)專業(yè)的的桌面地理信息平臺(tái)應(yīng)用框架，實(shí)現(xiàn)了地理信息應(yīng)用領(lǐng)域通用完備底層服務(wù)，使得開(kāi)發(fā)人員可以在一個(gè)已有的通用功能的基礎(chǔ)上進(jìn)行更加專業(yè)的系統(tǒng)開(kāi)發(fā)[5-6]。QGIS應(yīng)用框架以軟件設(shè)計(jì)的重用性和系統(tǒng)的可擴(kuò)展性為原則，以縮短大型應(yīng)用軟件系統(tǒng)的開(kāi)發(fā)周期，提高開(kāi)發(fā)質(zhì)量為目標(biāo)。與傳統(tǒng)的從底層開(kāi)發(fā)的區(qū)別在于，QGIS應(yīng)用框架技術(shù)強(qiáng)調(diào)的是通用功能的可重用性，系統(tǒng)開(kāi)發(fā)人員直接可以重用它的一些底層的算法，而不必關(guān)注這些算法如何實(shí)現(xiàn)，它通過(guò)提供一些外部接口，使得從已有構(gòu)件庫(kù)中建立應(yīng)用變得簡(jiǎn)單[7-8]。同時(shí)由于框架在結(jié)構(gòu)上的可擴(kuò)展性，使得開(kāi)發(fā)人員在此基礎(chǔ)上能夠方便進(jìn)行第三方庫(kù)的引入。因此，筆者基于開(kāi)源QGIS地理信息平臺(tái)應(yīng)用框架，以Visual Studio 2015作為二次開(kāi)發(fā)平臺(tái)，利用C++、QT語(yǔ)言進(jìn)行了油田遙感信息提取系統(tǒng)的開(kāi)發(fā)。

1.2 系統(tǒng)建立流程

創(chuàng)建油田遙感信息提取系統(tǒng)首先需要收集相關(guān)地圖資料和遙感影像等數(shù)據(jù)，本系統(tǒng)中使用到衛(wèi)星影像數(shù)據(jù)包括GF-2、Landsat-7 ETM+以及Hypersion。然后，對(duì)遙感影像數(shù)據(jù)進(jìn)行預(yù)處理，包括幾何校正、正射校正、圖像增強(qiáng)、空間變換、濾波以及高光譜分析處理，并對(duì)處理后的數(shù)據(jù)歸檔及錄入數(shù)據(jù)庫(kù)。最后，從軟件開(kāi)發(fā)的角度進(jìn)行需求設(shè)計(jì)、開(kāi)發(fā)架構(gòu)設(shè)計(jì)以及功能設(shè)計(jì)，利用Microsoft Visual Studio 2015可視化集成開(kāi)發(fā)平臺(tái)完成對(duì)油田遙感專題信息提取系統(tǒng)功能的開(kāi)發(fā)，實(shí)現(xiàn)從遙感圖像中獲取與油田勘探開(kāi)發(fā)有關(guān)信息的需求。本文以中石化油田為例，主要提取道路、水體和油氣微滲漏信息，如圖1所示。

圖1 油田遙感信息提取系統(tǒng)建立流程圖

1.3 系統(tǒng)功能結(jié)構(gòu)

圖2 系統(tǒng)功能結(jié)構(gòu)圖

系統(tǒng)采用模塊化分塊設(shè)計(jì)，油田遙感信息提取系統(tǒng)按照功能不同分為多個(gè)一級(jí)模塊，在各個(gè)一級(jí)模塊下根據(jù)需要設(shè)二級(jí)分模塊，二級(jí)分模塊下即為具體的子功能模塊，通過(guò)一級(jí)模塊+二級(jí)分模塊+子功能模塊的形式，形成層次化的模塊結(jié)構(gòu)。油田遙感信息提取系統(tǒng)由工程模塊、視圖模塊、圖像校正模塊、圖像增強(qiáng)模塊、光譜分析模塊和信息提取模塊組成，如圖2所示。

1.4 數(shù)據(jù)讀取與可視化

空間數(shù)據(jù)讀取與可視化是遙感信息提取系統(tǒng)基本功能之一。空間數(shù)據(jù)從數(shù)據(jù)結(jié)構(gòu)上分為柵格文件和矢量文件，在遙感圖像中，主要包括GeoTiff、IMG、PIX等柵格格式。QGIS地理信息平臺(tái)支持一些通用遙感數(shù)據(jù)格式的讀取，對(duì)于無(wú)法讀取的數(shù)據(jù)格式，通過(guò)GDAL柵格地理數(shù)據(jù)格式庫(kù)新建數(shù)據(jù)驅(qū)動(dòng)，調(diào)用GDALAllRegister()函數(shù)來(lái)進(jìn)行擴(kuò)展。系統(tǒng)可視化界面如圖3所示。

圖3 系統(tǒng)可視化界面

1.5 油田遙感數(shù)據(jù)處理

系統(tǒng)不僅能夠?yàn)橛脩籼峁┛梢暬?wù)，而且在QGIS平臺(tái)基礎(chǔ)上進(jìn)行了數(shù)據(jù)處理應(yīng)用擴(kuò)展，油田遙感信息提取系統(tǒng)不僅開(kāi)發(fā)了基礎(chǔ)的幾何校正、投影轉(zhuǎn)換服務(wù)，而且集成了圖像分析模塊，通過(guò)圖像的拉伸、銳化、像元分析等為油田遙感專題信息提取提供輔助決策，包括圖像增強(qiáng)、光譜分析等。

圖像增強(qiáng)模塊針對(duì)遙感影像像元值，通過(guò)相應(yīng)的卷積算法進(jìn)行色調(diào)調(diào)整，突出圖像中的有用信息，使圖像中感興趣的內(nèi)容得以強(qiáng)調(diào)和清晰，從而有利于目視解譯和計(jì)算機(jī)信息提取。

光譜分析模塊從像元光譜特征出發(fā)，建立針對(duì)高光譜圖像進(jìn)行特征分析的模型，包括光譜特征、邊緣特征、紋理特征、形狀特征等數(shù)據(jù)進(jìn)行分析，以獲取相關(guān)的圖形分析結(jié)果。

1.6 油田遙感專題信息提取

油田遙感信息提取模塊在數(shù)據(jù)處理功能的基礎(chǔ)上，開(kāi)發(fā)了針對(duì)不同遙感影像專題信息提取功能，實(shí)現(xiàn)了道路、水體、油氣微滲漏等信息提取。專題信息提取模塊基于QGIS平臺(tái)框架，通過(guò)C++、QT編程語(yǔ)言實(shí)現(xiàn)了自己特有的圖形及數(shù)據(jù)管理組件，可以獨(dú)立運(yùn)行，目前只實(shí)現(xiàn)了部分關(guān)鍵信息快速提取算法。為了后續(xù)能夠擴(kuò)展信息提取模塊，平臺(tái)預(yù)留了算法擴(kuò)展接口，方便不斷完善油田遙感應(yīng)用平臺(tái)。

2 關(guān)鍵技術(shù)

2.1 道路提取

2.1.1 技術(shù)流程

道路提取算法從自動(dòng)化程度上分為全自動(dòng)和半自動(dòng)算法[9]。由于目前計(jì)算機(jī)人工智能還不能實(shí)現(xiàn)全自動(dòng)從數(shù)字影像中理解與提取地物信息，因此本系統(tǒng)在經(jīng)過(guò)反復(fù)試驗(yàn)和野外驗(yàn)證的基礎(chǔ)上，決定采用針對(duì)油田工作區(qū)遙感影像的、基于模板匹配的半自動(dòng)提取算法。

模板匹配道路提取算法是一種半自動(dòng)提取方法，它需要手動(dòng)創(chuàng)建模板，統(tǒng)計(jì)目標(biāo)區(qū)域剖面像元與模板區(qū)域剖面像元灰度值差的平方和，最小值對(duì)應(yīng)的像元點(diǎn)作為新的道路點(diǎn)，同時(shí)計(jì)算舊點(diǎn)與該點(diǎn)連接的直線與水平方向直線的夾角[10](圖4)。具體步驟如下：

圖4 道路提取算法流程圖

(1)首先通過(guò)人工交互確定道路中心線的起點(diǎn)，一般包括道路中心點(diǎn)在內(nèi)的3個(gè)以上點(diǎn)。

markMatchPoint(Point pre_pt,Point cur_pt)

(2)創(chuàng)建道路模板剖面。通過(guò)人工交互的方式獲取道路的起點(diǎn)p、道路延伸方向θ、道路寬度w, 垂直于道路延伸方向的道路剖面，剖面寬度L是道路寬度w的2倍。

cross_width=abs(pre_pt.x-cur_pt.x)

cross_height=abs(pre_pt.y-cur_pt.y)

cross_widthy=2*cross_width

Templ=getTempl(pre_pt,cur_pt)

(3)創(chuàng)建道路的目標(biāo)剖面，統(tǒng)計(jì)模板剖面像元與目標(biāo)剖面像元的灰度差平方和。從當(dāng)前道路點(diǎn)沿道路延伸方向θ前進(jìn)長(zhǎng)度s，到達(dá)點(diǎn)p，然后在道路延伸方向上實(shí)施幾何變換，形成多個(gè)目標(biāo)剖面，計(jì)算這些目標(biāo)剖面的像元灰度值與模板剖面的差的平方和，最小值對(duì)應(yīng)的像元點(diǎn)作為新的道路點(diǎn)，同時(shí)計(jì)算舊點(diǎn)與該點(diǎn)連接的直線與水平方向直線的夾角。

matchTemplate( cross_Yimg, Templ, cross_Result, CV_TM_SQDIFF )

minMaxLoc( cross_Result, &cross_minValue, &cross_maxValue, &cross_minLoc, &cross_maxLoc, mask)

(4)判斷道路提取結(jié)束條件。如果到達(dá)影像的邊界，或者延伸到之前已提取的道路，或者最小二乘的結(jié)果超過(guò)規(guī)定的閾值T，則跳到步驟(5)；否則，將該新提取的道路點(diǎn)添加到追蹤軌跡中，并跳到步驟(3)。

(5)結(jié)束道路提取，人工進(jìn)行判別、整理。

2.1.2 試驗(yàn)分析

利用中石化某油田東部區(qū)域一張GF-2全色影像對(duì)上述模版匹配道路提取算法進(jìn)行了驗(yàn)證，該圖像的大小為3 124像元×3 160像元，圖像上道路特征明顯，路面輻射比較均一，有少量樹(shù)木造成的遮蔽。試驗(yàn)統(tǒng)計(jì)發(fā)現(xiàn)在提取過(guò)程中出現(xiàn)追蹤錯(cuò)誤次數(shù)、追蹤停止次數(shù)、提示用戶手工輸入次數(shù)分別為10次、12次、69次，總耗時(shí)1 900s。

2.2 水體提取

2.2.1 技術(shù)途徑

水的光譜特征是由水的物理和化學(xué)性質(zhì)決定的。對(duì)于清水，在可見(jiàn)光波段500nm之前，水的吸收和反射率較低，透射較高；在藍(lán)—綠光波段反射率4%～5%，500nm以下的紅光部分反射率降到2%～3%，在近紅外和短波紅外波段幾乎全部吸收，反射很少，因此利用紅外波段提取水體較為可行[11]。

結(jié)合水體的光譜特征，通過(guò)特征波段構(gòu)造水體指數(shù)模型實(shí)現(xiàn)水體提取。歸一化水體指數(shù)(NDWI)模型[12]就是一種廣泛使用的算法，其表達(dá)式為:

NDWI=(band2-band4)/(band2+band4)

(1)

式中:band2、band4分別表示TM影像中的第2和第4波段，NDWI模型根據(jù)水體和植被在綠光和近紅外兩個(gè)波段的反射吸收波譜差異，使植被等干擾信息得到最大限度的抑制，而水體信息得到最大程度的增強(qiáng)。主要實(shí)現(xiàn)流程如下：

(1)計(jì)算NDWI：ndwi[i]=(green[i]-nir[i])/(nir[i]+green[i])。

(2)保存NDWI，并以文件形式輸出：SaveNDWI(rows, cols, ndwi, NdwiFile)。

2.2.2 試驗(yàn)分析

試驗(yàn)選擇中石化某油田東部區(qū)域Landsat-7 ETM+遙感影像，計(jì)算NDWI后，根據(jù)水體指數(shù)灰度直方圖，確定區(qū)分水體與非水體的閾值, 輸出水體信息分布影像。通過(guò)精度評(píng)定發(fā)現(xiàn)，Kappa系數(shù)達(dá)到0.983 3，漏提率0.15%，誤提率0.75%，總體精度達(dá)到98.65%。

2.3 油氣微滲漏信息提取

2.3.1 技術(shù)途徑

微滲漏烴類擴(kuò)散到地表后，在水和微生物的作用下，與土壤中的金屬離子反應(yīng)生成碳酸鹽，因此油氣層上方地表土壤的碳酸鹽總量顯著提高，其特征吸收波段在2 350nm左右。同時(shí)微滲漏烴類中所含的氫硫化物和碳?xì)浠锔淖兞松细矌r層的酸堿環(huán)境，pH值降低，酸性環(huán)境使得巖石中長(zhǎng)石礦物被粘土礦物所置換，改變了地表上粘土礦物的分布，致使油氣層上方地表粘土礦物豐度明顯提高[12]。試驗(yàn)發(fā)現(xiàn)粘土礦化在2.2μm處表現(xiàn)為較強(qiáng)的羥基基團(tuán)的吸收特征。

試驗(yàn)發(fā)現(xiàn)礦物含量與光譜吸收特征參數(shù)吸收深度H相關(guān)性最高，結(jié)合到具體識(shí)別上，希望構(gòu)建一個(gè)來(lái)指示油氣異常點(diǎn)[13-14]的模型，通過(guò)大量試驗(yàn)表明，有烴類成因的碳酸鹽礦物的吸收深度與粘土礦物的吸收深度基本上接近，因此根據(jù)這個(gè)原理構(gòu)建了一個(gè)指示油氣微滲漏的快速指數(shù)(FI)模型[15]。

(2)

式中:H1代表碳酸鹽蝕變特征的光譜吸收深度，H2為粘土礦物蝕變特征的光譜吸收深度。在油氣微滲漏區(qū)域，由于兩種蝕變礦物吸收深度接近，所以模型中分子較大，分母較小，F(xiàn)I值較大。具體步驟如下：

(1)打開(kāi)原始影像數(shù)據(jù)，進(jìn)行幾何校正、研究區(qū)剪裁等數(shù)據(jù)預(yù)處理工作。

(2)對(duì)預(yù)處理后的影像進(jìn)行包絡(luò)線去除。

continuumRemoval(srcFile，resultFile)；

(3)定義去除包絡(luò)線之后影像上2 200nm附近處的波段為band1，定義2 350nm附近處的波段為band2，分別計(jì)算1-band1和1-band2，即為碳酸鹽礦物和粘土礦物特征吸收位置處的吸收深度H1和H2。

(4)利用式(2)計(jì)算整幅影像的FI值。

FI[i*x_size+j]=((1-D[index_first_wavelength])+(1-D[index_second_wavelength]))/ abs((1-D[index_first_wavelength])-(1-D[index_second_wavelength]))；

(5)利用同樣的方法計(jì)算實(shí)測(cè)數(shù)據(jù)，根據(jù)實(shí)際采樣點(diǎn)的分布情況定義油氣指數(shù)的閾值T，F(xiàn)I>T為油氣微滲漏分布處。

2.3.2 試驗(yàn)分析

利用中石化某油田東部地區(qū)一張Hyperion高光譜影像進(jìn)行油氣微滲漏信息提取，該區(qū)域內(nèi)分布多個(gè)油氣地層，且伴隨有許多采油井，側(cè)面證明了該區(qū)域存在油氣微滲漏的可能。試驗(yàn)發(fā)現(xiàn)通過(guò)FI指數(shù)提取的油氣微滲漏信息與野外勘探的結(jié)果吻合性達(dá)到97%，同時(shí)由于大部分位油氣微滲漏點(diǎn)位于油氣開(kāi)采區(qū)，間接支撐了該方法的有效性。

3 結(jié)論

基于QGIS的油田遙感信息提取系統(tǒng)采用模塊化的編程方法，通過(guò)開(kāi)源QGIS平臺(tái)進(jìn)行二次開(kāi)發(fā)，集成了遙感數(shù)據(jù)基礎(chǔ)處理功能、數(shù)據(jù)分析和專題信息提取功能，為油田遙感數(shù)據(jù)處理、信息提取和分析提供一個(gè)通用平臺(tái)，輔助油田勘探開(kāi)發(fā)作業(yè)，為相關(guān)部門(mén)提供決策依據(jù)?；赒GIS應(yīng)用框架構(gòu)建的油田遙感信息提取系統(tǒng)有效縮短了軟件的開(kāi)發(fā)周期，提高了系統(tǒng)的重用性和可擴(kuò)展性。對(duì)于系統(tǒng)中涉及到的一些專業(yè)問(wèn)題，需要后續(xù)不斷研究開(kāi)發(fā)補(bǔ)充完善。隨著衛(wèi)星遙感技術(shù)的發(fā)展、油田數(shù)據(jù)量的增長(zhǎng)和系統(tǒng)的深入開(kāi)發(fā)，遙感支撐油田勘探開(kāi)發(fā)的能力將越來(lái)越強(qiáng)。