楊達(dá)昌, 陳 潔,2, 高子弘, 韓亞超
(1.中國國土資源航空物探遙感中心,北京 100083; 2.中國科學(xué)院遙感與數(shù)字地球研究所,北京 100101)
地質(zhì)是高光譜遙感最早應(yīng)用的一個(gè)領(lǐng)域,也是高光譜遙感應(yīng)用中最成功的領(lǐng)域之一。特別是由于高光譜數(shù)據(jù)具有的連續(xù)窄波段特性,能夠精細(xì)地反映出礦物和巖石的診斷性光譜特征(特征波段寬度和位置等),光譜分辨率越高,礦物的吸收、反射診斷性光譜特征就越明顯,這從根本上提升了從遙感圖像上提取地質(zhì)信息的質(zhì)量和數(shù)量。
作為一種直接的油氣調(diào)查手段,烴類微滲漏遙感探測(cè)技術(shù)已在國內(nèi)外受到廣泛關(guān)注,許多學(xué)者對(duì)此進(jìn)行了研究。起初,受限于光譜分辨率,提取方法主要是基于油氣微滲漏會(huì)引起上覆蓋層的色調(diào)異常理論,通過圖像增強(qiáng)的方法對(duì)多光譜遙感數(shù)據(jù)進(jìn)行處理,進(jìn)而實(shí)現(xiàn)油氣遠(yuǎn)景區(qū)預(yù)測(cè)。Segal等[1]通過Landsat TM數(shù)據(jù)的波段比值和彩色空間變換方法,證實(shí)了由于烴類微滲漏造成的地表巖石蝕變褪色現(xiàn)象。但是,因?yàn)槎喙庾V遙感數(shù)據(jù)的光譜分辨率不高,并且引起地表巖石和礦物的色調(diào)變化因素較多、機(jī)理復(fù)雜,僅依據(jù)該現(xiàn)象解釋油氣微滲漏會(huì)造成一定的誤判,使得多光譜遙感油氣微滲漏信息探測(cè)結(jié)果在未知油氣區(qū)存在多解性。隨著近年來高光譜遙感技術(shù)的發(fā)展,油氣微滲漏遙感探測(cè)手段進(jìn)一步豐富。Ellis等[2]建立了油氣滲漏區(qū)土壤光譜庫,為其他地質(zhì)研究提供了依據(jù); 周密[3]利用星載高光譜Hyperion數(shù)據(jù)提取了東勝地區(qū)油氣微滲漏的空間分布信息; 李娜等[4]基于ASTER數(shù)據(jù),開展了基于波段運(yùn)算和主成分分析的烴類微滲漏蝕變信息提??; 陳圣波等[5]通過對(duì)甘肅慶陽油氣區(qū)地面黃土樣品實(shí)測(cè)波譜曲線的分析,明確了油氣微滲漏引起的蝕變礦物含量與光譜吸收特征之間的關(guān)系。
本文利用天宮一號(hào)(TG-1)高光譜數(shù)據(jù)進(jìn)行油氣微滲漏信息提取,研究油氣微滲漏地表蝕變的光譜響應(yīng)機(jī)理,實(shí)現(xiàn)油氣微滲漏共生異常的分類識(shí)別,完成油氣微滲漏地表共生異常區(qū)的初步圈定,生成油氣勘探綜合異??臻g分布圖,為油氣勘探靶區(qū)的圈定和新區(qū)油氣勘探工作部署提供依據(jù)。
埋藏在地下深部、具有動(dòng)態(tài)平衡的油氣藏,其內(nèi)部具有相較于上覆蓋層更高的壓力,這使得油氣藏中的烴類物質(zhì)通過滲透、水動(dòng)力和擴(kuò)散的運(yùn)移方式,穿過上覆蓋層的孔隙、斷裂和其他不整合面到達(dá)地表,甚至擴(kuò)散至近地表的空氣中。這種微滲漏的烴類物質(zhì)與土壤和巖石發(fā)生蝕變作用,使油氣藏上方地表在遙感影像中具有某種特殊的光譜異常表現(xiàn),通過對(duì)烴類微滲漏信息的提取即可實(shí)現(xiàn)油氣資源的遙感探測(cè)。遙感探測(cè)微滲漏的指示標(biāo)志有: 熱異常暈、紅層褪色暈、低價(jià)鐵富集暈、粘土礦化暈和烴類異常暈。其中熱異常暈是由于油氣藏上方的地表溫度升高引起的異?,F(xiàn)象,可利用熱紅外波段進(jìn)行溫度異常提?。?而紅層褪色暈是一種非蝕變巖層的褪色現(xiàn)象,可利用可見光波段色彩增強(qiáng)進(jìn)行信息提取。鑒于本文主要討論烴類微滲漏引起的礦物和烴類物質(zhì)的光譜異常信息,故針對(duì)后3種異常指示標(biāo)志進(jìn)行詳細(xì)闡述。
當(dāng)三價(jià)鐵離子轉(zhuǎn)化為二價(jià)鐵離子時(shí),在形成紅層褪色暈的同時(shí),也形成了二價(jià)鐵離子的富集。由于二價(jià)鐵離子與三價(jià)鐵離子有著不同的波譜特性,因而出現(xiàn)波譜異常。這種由三價(jià)鐵離子轉(zhuǎn)化為二價(jià)鐵離子所形成的遙感影像波譜異常稱為低價(jià)鐵富集暈。由波譜分析知,二價(jià)鐵離子礦物(如黃鐵礦)在1 000 nm處有強(qiáng)吸收帶,在1 000~1 500 nm處有寬吸收帶; 而含有三價(jià)鐵離子的礦物(如褐鐵礦、黃鉀鐵礬、赤鐵礦和針鐵礦等)則在800~900 nm處表現(xiàn)出強(qiáng)吸收帶,并在400~500 nm處表現(xiàn)為吸收特征。利用該波譜特征差異,可以進(jìn)行高價(jià)鐵離子和低價(jià)鐵離子的區(qū)分,提取低價(jià)鐵富集暈信息,從而達(dá)到識(shí)別烴類微滲漏的目的。
在發(fā)生烴類微滲漏時(shí),由于微滲漏的液體或氣體中所含的氫硫化物和碳?xì)浠锔淖兞松细矌r石的氧化-還原環(huán)境,使得巖石被酸化,而酸性溶液可使某些在堿性環(huán)境中穩(wěn)定的礦物(如長(zhǎng)石和粘土)被蝕變或被其他粘土礦物(如高嶺石、綠泥石和蒙脫石等)代替。而這一系列蝕變作用的結(jié)果是,在烴類微滲漏上方,粘土礦物的分布與構(gòu)成被改變,導(dǎo)致在該區(qū)域存在粘土礦物的高豐度異常。這種由于烴類微滲漏導(dǎo)致的粘土富集所出現(xiàn)的遙感影像波譜異常稱為粘土礦化暈,它們往往具有較強(qiáng)的羥基基團(tuán)的吸收特征,其波譜特征為: 在短波紅外波段2 200 nm處有很強(qiáng)的羥基吸收帶; 在2 000~2 500 nm處較硅酸鹽巖、碳酸鹽巖呈低反射率。利用這一波譜特征,可以實(shí)現(xiàn)粘土礦化信息的提取,識(shí)別烴類微滲漏。
在烴類微滲漏過程中,由于烴類物質(zhì)揮發(fā)或滲漏到地表或近地表,通常會(huì)吸附或溶解在土壤孔隙間,使得油氣富集區(qū)地表形成烴類含量異常。也稱為烴類異常暈。因此在高光譜遙感影像上往往出現(xiàn)烴類異常,可以用于標(biāo)識(shí)烴類物質(zhì)的滲漏范圍。
首先,對(duì)烴類物質(zhì)微滲漏現(xiàn)象以及由此引起的地表蝕變進(jìn)行梳理和歸納; 然后,以烴類物質(zhì)微滲漏引起的地表土壤及巖石地球化學(xué)異常、地表土壤吸附烴異常等特征建立遙感指示標(biāo)志; 最后,利用所建立的指示標(biāo)志實(shí)現(xiàn)油氣微滲漏的高光譜信息提取技術(shù)和方法研究,并使用TG-1高光譜數(shù)據(jù)進(jìn)行實(shí)驗(yàn)驗(yàn)證。技術(shù)路線如圖1所示。
圖1 油氣微滲漏高光譜信息提取技術(shù)路線
航天高光譜以衛(wèi)星作為遙感平臺(tái),具有全球覆蓋、可獲取任意目標(biāo)光譜信息的優(yōu)點(diǎn),美國EO-1衛(wèi)星搭載的Hyperion和歐空局Proba衛(wèi)星搭載的CHRIS,已經(jīng)展示了航天高光譜成像技術(shù)的信息獲取能力和應(yīng)用潛力。于2008年發(fā)射的環(huán)境一號(hào)衛(wèi)星(HJ-1A),是我國首個(gè)搭載高光譜儀的衛(wèi)星,但其主要面向環(huán)境與災(zāi)害監(jiān)測(cè)預(yù)報(bào),其工作譜段(459~956 nm)和空間分辨率(100 m)都嚴(yán)重制約了其在地質(zhì)和巖礦識(shí)別領(lǐng)域的應(yīng)用。而TG-1目標(biāo)飛行器,搭載了我國自主研制的、空間和光譜綜合指標(biāo)最高的高光譜成像設(shè)備,具有空間分辨率高、波段多、光譜分辨率高的特點(diǎn),真正意義上實(shí)現(xiàn)了“圖譜合一”。其主要技術(shù)指標(biāo)及與國際上在軌運(yùn)行的航天高光譜傳感器的參數(shù)對(duì)比如表1所示。
表1 TG-1高光譜成像儀及同類型航天傳感器技術(shù)指標(biāo)對(duì)比Tab.1 Comparison of TG-1 hyperspectral imaging system and the other space sensors
通過對(duì)比分析,TG-1搭載的高光譜成像儀可見光—近紅外、短波紅外高光譜數(shù)據(jù)的各項(xiàng)指標(biāo)(如空間分辨率、光譜分辨率)在國內(nèi)外均處于領(lǐng)先水平,某些性能甚至優(yōu)于國外同類型數(shù)據(jù)。另外,由于烴類微滲漏產(chǎn)生的蝕變礦物光譜特征基本都位于1 000~ 2 500 nm波段范圍,故本次研究選用甘肅環(huán)縣—慶陽地區(qū)TG-1短波紅外數(shù)據(jù)作為研究對(duì)象。數(shù)據(jù)獲取時(shí)間為2012年08月07日13: 57,為L(zhǎng)2級(jí)圖像數(shù)據(jù)產(chǎn)品,格式為標(biāo)準(zhǔn)GeoTIFF(6.0)。
2.3.1 冗余波段去除
TG-1高光譜短波紅外數(shù)據(jù)共有75個(gè)波段,其中1~6波段為與可見光—近紅外波段重合,予以剔除; 然后再剔除受水汽影響的波段,波譜范圍為1 356.35~1 416.94 nm,1 820.48~ 1 931.47 nm和2 395.5~2 565 nm,分別對(duì)應(yīng)波段為20~21、36~42、69~75。剩余53個(gè)波段將用于后續(xù)數(shù)據(jù)分析。
2.3.2 壞線修復(fù)及條帶去除
采用基于多光譜影像的光譜重構(gòu)方法進(jìn)行高光譜數(shù)據(jù)的壞線修復(fù)處理,并使用全局去條紋法進(jìn)行條帶去除,即通過像元的列平均值、標(biāo)準(zhǔn)差和波段平均值、標(biāo)準(zhǔn)差之間的差異對(duì)像元進(jìn)行分波段線性化修正,消除垂直條紋的影響。修復(fù)結(jié)果有較好的目視效果和良好的光譜保真度,基本滿足專題制圖和信息反演的需要。
2.3.3 大氣校正
為了消除因大氣散射和吸收引起的地物光譜曲線失真,采用基于大氣輻射傳輸模型的FLAASH法對(duì)TG-1數(shù)據(jù)進(jìn)行大氣校正。該方法根據(jù)高光譜數(shù)據(jù)獲取的時(shí)間和地理位置,選擇相應(yīng)的模型完成,模型嚴(yán)密、解算精度高,能較真實(shí)地還原地物光譜特征。為了對(duì)校正后的光譜數(shù)據(jù)進(jìn)行檢驗(yàn)和評(píng)價(jià),利用研究區(qū)內(nèi)的植被光譜與標(biāo)準(zhǔn)光譜進(jìn)行吸收特征比對(duì)。結(jié)果顯示,雖然由于兩者的光譜分辨率差異較大,重建光譜與標(biāo)準(zhǔn)波譜存在一定的波長(zhǎng)位置偏移,但偏移量均在合理的光譜區(qū)間內(nèi),曲線特征吻合較好,基本滿足后續(xù)研究要求。
2.4.1 粘土異常信息提取
烴類微滲漏的液體或氣體中所含的氫硫化物和碳?xì)浠锔淖兞松细矌r石的氧化-還原環(huán)境,導(dǎo)致粘土礦物的集中分布。將美國地質(zhì)調(diào)查局發(fā)布的典型粘土礦物的標(biāo)準(zhǔn)波譜重采樣至TG-1高光譜短波紅外譜段范圍,同時(shí)對(duì)其波譜曲線進(jìn)行包絡(luò)線去除,對(duì)比分析其光譜特征(圖2和表2)。
圖2 典型粘土礦物去除包絡(luò)線后的波譜曲線
礦物名稱英文吸收位置反射位置1反射位置2反射位置3反射位置4高嶺石Kaolinite2 205753.5~1 273.51 503.5~2 0252 225~2 265—蒙脫石Montmorillonite2 215和1 915753.5~1 323.51 546.0~1 8352 055~2 1752 255~2 295白云母Muscovite2 205750.5~1 368.51 503.5~2 1352 255~2 295—
從圖2和表2中可以得出如下結(jié)論: ①3種典型粘土礦物的共同吸收波段位于2 200 nm附近; ②共同高反射位置為1 700~1 800 nm和2 100 nm附近; ③由于1 700~1 800 nm波段處與水汽影響波段接近,可能會(huì)造成錯(cuò)誤的信息提取,故而使用2 100 nm附近的波段作為高反射響應(yīng)波段。
采用能突出所需指示標(biāo)志并壓制其他地物特征的波段進(jìn)行信息提取,選定大氣校正后短波紅外波段的第47波段(中心波長(zhǎng)為2 044.827 nm)與第55波段(中心波長(zhǎng)為2 210.329 nm)進(jìn)行比值處理。粘土富集異常信息的提取公式為
B47/55=R47/R55,
(1)
式中:R47和R55分別為第47和55波段的反射率值;B47/55為這2個(gè)波段的比值。
受限于空間分辨率和混合礦物等因素的影響,實(shí)測(cè)的波段反射率比值與標(biāo)準(zhǔn)波譜所得的數(shù)值之間肯定存在一定的誤差,因此在進(jìn)行異常信息提取時(shí),對(duì)通過標(biāo)準(zhǔn)光譜曲線計(jì)算出的范圍進(jìn)行適當(dāng)?shù)臄U(kuò)展,使其更符合實(shí)際情況。高嶺石、蒙脫石和白云母的標(biāo)準(zhǔn)波譜曲線計(jì)算得出的B47/55值分別為1.834,1.268和1.814,將此標(biāo)準(zhǔn)范圍擴(kuò)展20%,得到實(shí)際提取的變化范圍為1.01~2.20。
2.4.2 低價(jià)鐵離子異常信息提取
二價(jià)鐵的富集是由于烴類微滲漏使三價(jià)鐵離子轉(zhuǎn)化為二價(jià)鐵離子而形成的,所以在進(jìn)行低價(jià)鐵離子的異常信息提取時(shí),除了要按照上述方法尋找二價(jià)鐵離子代表性礦物的高反射和吸收波段外,還需要考慮二價(jià)、三價(jià)鐵離子及混合礦物的共存和區(qū)分問題。這3種物質(zhì)的代表性礦物光譜特征如圖3和表3所示。同樣,對(duì)這些典型礦物的波譜進(jìn)行了包絡(luò)線去除,能進(jìn)一步突出吸收和反射特征。
從圖3和表3中可以看出: ①含二價(jià)鐵礦物吸收位置波長(zhǎng)要大于含三價(jià)鐵礦物及其混合礦物,且屬于TG-1短波紅外數(shù)據(jù)的譜段內(nèi); ②含二價(jià)鐵礦物反射率在1 200—1 800 nm處,急速上升,而含三價(jià)鐵礦物反射率幾乎無變化,混合礦物的反射率雖然經(jīng)歷了上升再下降的過程,但變化幅度不大。
為了突出二價(jià)鐵礦物并壓制其他,選定短波紅外波段的第34波段(中心波長(zhǎng)為1 742.99 nm)和第14波段(中心波長(zhǎng)為1 196.036 nm)進(jìn)行比值處理,低價(jià)鐵離子富集異常信息提取公式為
B34/14=R34/R14,
(2)
式中:R34和R14分別為第34和14波段的反射率值;B34/14為這2個(gè)波段的比值。
按照式(2)計(jì)算表3中的5種代表性礦物的B34/14值分別為1.028,1.151,1.360,4.839和1.520。低價(jià)鐵的比值最大,混合礦物的比值介于兩者之間。故將實(shí)際提取低價(jià)鐵異常信息的取值范圍定為1.60~5.80。
2.4.3 烴類異常信息提取
物質(zhì)烴類運(yùn)移到上覆蓋面時(shí),通常會(huì)吸附或溶解在土壤孔隙間,導(dǎo)致地表烴類含量異常。烴類物質(zhì)在1 725 nm和2 270~2 460 nm處有明顯的吸收特征,如圖4所示。
圖4 典型烴類物質(zhì)反射光譜曲線
但由于2 400 nm附近位于水汽吸收波段,不予采用,且Kühn等[6]也指出1 730 nm是烴類探測(cè)的關(guān)鍵性特征,故本文采用該處作為烴類的吸收位置對(duì)其進(jìn)行提取,同時(shí)選取位于2個(gè)吸收谷之間的高值位置(2 070 nm)作為烴類的高反射波段。
選定大氣校正后短波紅外波段的第48波段(中心波長(zhǎng)為2 066.381 nm)和第33波段(中心波長(zhǎng)為1 717.600 nm)進(jìn)行比值處理。烴類物質(zhì)富集異常信息提取公式為
B48/33=R48/R33,
(3)
式中:R48和R33分別為第48和33波段的反射率值;B48/33為這2個(gè)波段的比值。
按照式(3)計(jì)算的B48/33值約為1.50,故將實(shí)際提取烴類異常信息的取值范圍定為1.40~ 1.70。
2.4.4 油氣微滲漏綜合信息提取
根據(jù)粘土礦物蝕變異常、二價(jià)鐵離子蝕變異常和烴類物質(zhì)蝕變異常對(duì)油氣微滲漏的影響程度,建立綜合評(píng)價(jià)模型,得出研究區(qū)油氣微滲漏情況的綜合評(píng)價(jià)結(jié)果。由已有研究成果[4]確定各蝕變礦物權(quán)重值分別為: 粘土礦物為0.97,二價(jià)鐵離子為0.58,烴類為0.70。油氣微滲漏綜合信息Y提取公式為
Y=0.97B47/55+0.58B34/14+0.70B48/33。
(4)
利用經(jīng)數(shù)據(jù)預(yù)處理后的TG-1短波紅外高光譜數(shù)據(jù),選取具有典型反射率特征的波段,采用波段比值法進(jìn)行粘土礦物蝕變、低價(jià)鐵蝕變和烴類異常信息提取,其信息專題圖如圖5所示?;谏鲜鲂畔?,利用油氣微滲漏綜合信息提取法,獲得研究區(qū)局部提取結(jié)果如圖6所示。底圖均為全色波段圖像。
(a) 粘土礦物蝕變信息 (b) 低價(jià)鐵蝕變信息 (c) 烴類物質(zhì)異常信息
圖5油氣微滲漏信息提取專題圖
Fig.5Extractionthematicmapsofoilandgasmicro-seepage
圖6 慶陽地區(qū)油氣微滲漏信息提取專題圖
從圖5中可以發(fā)現(xiàn),3種不同蝕變信息的分布情況較為一致,驗(yàn)證了烴類微滲漏會(huì)引起地表物質(zhì)發(fā)生蝕變而造成粘土、低價(jià)鐵礦物和烴類物質(zhì)富集理論的可靠性。在圖6中,將3種指示標(biāo)志均有出現(xiàn)的區(qū)域視為油氣富集靶區(qū),可作進(jìn)一步研究分析的依據(jù)。
根據(jù)研究區(qū)地質(zhì)資料,該區(qū)域位于鄂爾多斯盆地的天環(huán)凹陷中部,緊鄰伊陜斜坡,構(gòu)造特征為傾角不足1° 的西傾單斜,鼻狀隆起發(fā)育。三疊紀(jì)時(shí)受西南物源控制,沉積物以碎屑巖為主,縱向上巖性變化大,形成局部封閉圈,成為油氣聚集有利場(chǎng)所。王建民等[7]通過對(duì)鉆井及試采資料的分析,利用密井網(wǎng)繪圖方法,結(jié)合地質(zhì)構(gòu)造和演化特征指出了慶陽地區(qū)具有油氣富集特征; 李振華等[8]以油氣成藏地質(zhì)學(xué)理論為基礎(chǔ),通過成藏條件對(duì)比,認(rèn)為慶陽地區(qū)具有良好的油氣儲(chǔ)藏潛力; 楊華等[9]通過模擬實(shí)驗(yàn)和資料分析,得出慶陽地區(qū)具有致密氣藏形成的有利條件。由此可見,研究區(qū)內(nèi)烴源巖發(fā)生過大規(guī)模生烴階段,油氣資源富集,存在烴類微滲漏遙感異常指示標(biāo)志的可能性很大,可利用上述遙感方法進(jìn)行信息提取。
2015年,張?zhí)煊頪10]曾應(yīng)用快速指數(shù)法和光譜角匹配法對(duì)慶城地區(qū)的Hyperion高光譜數(shù)據(jù)進(jìn)行油氣微滲漏信息提取,初步圈定了油氣空間分布,與實(shí)地的采油井分布及已知油田區(qū)一致性較好。并且其提取方法同樣是基于油氣微滲漏理論,遙感指示標(biāo)志與本文相同,印證了本文方法的可行性。從圖6提取結(jié)果的異常信息分布情況看,本文提取的異常信息分布與前人的地質(zhì)分析結(jié)果和提取結(jié)果也都吻合較好。
提取結(jié)果表明,數(shù)據(jù)覆蓋區(qū)具有明顯的烴類微滲漏遙感異常指示標(biāo)志存在。本文在3種異常信息同時(shí)存在的地區(qū),劃定了一個(gè)15像素×15像素的區(qū)域,進(jìn)行平均波譜曲線計(jì)算。該曲線在2 200 nm和2 300 nm處存在明顯的雙吸收峰,并且后者還具有一定的吸收寬度。2 200 nm處的吸收特征,主要受粘土礦物和烴類礦物影響,而2 300 nm處,粘土、碳酸鹽巖和烴類物質(zhì)均表現(xiàn)有吸收特征,其中碳酸鹽巖礦物的吸收強(qiáng)度最大,故而造成該波段處存在一定的吸收寬度。
基于油氣微滲漏理論,選取了3種具有顯著波譜特性的烴類微滲漏遙感指示標(biāo)志作為遙感識(shí)別目標(biāo),在分析了烴類對(duì)其上覆蓋層蝕變物質(zhì)光譜吸收特征的基礎(chǔ)上,提出運(yùn)用波段比值法突出異常信息,成功提取了甘肅環(huán)縣—慶陽地區(qū)的3種烴類微滲漏蝕變信息,與該地區(qū)已有地質(zhì)調(diào)查資料和油氣儲(chǔ)藏分析結(jié)果吻合; 對(duì)異常信息區(qū)的波譜解譯也表明,烴類微滲漏的蝕變特征均有所體現(xiàn),驗(yàn)證了利用比值法進(jìn)行天宮一號(hào)高光譜數(shù)據(jù)油氣微滲漏信息提取方法的準(zhǔn)確性和可行性。
本研究在提取烴類微滲漏高光譜異常信息過程中,由于僅采用了普遍的壞線處理、大氣校正等方法,地物真實(shí)光譜的恢復(fù)程度一般; 遙感指示標(biāo)志均只應(yīng)用比值法進(jìn)行突出,也給信息提取帶來了一定誤差。在下一階段的研究中,應(yīng)針對(duì)天宮一號(hào)數(shù)據(jù)特點(diǎn)和識(shí)別目標(biāo)的波譜特性,綜合運(yùn)用沙漏法、監(jiān)督/非監(jiān)督法等提升遙感指示標(biāo)志的提取效果,探索更實(shí)用、更有效的數(shù)據(jù)預(yù)處理技術(shù)和方法。
參考文獻(xiàn)(References):
[1] Segal D B,Merin I S.Successful use of Landsat thematic mapper data for mapping hydrocarbon microseepage-induced mineralogic alteration, Lisbon Valley,Utah[J].Photogrammetric Engineering and Remote Sensing,1989,55(7):1137-1145.
[2] Ellis J M,Davis H H,Zamudio A.Exploring for onshore oil seeps with hyperspectral imaging[J].Oil and Gas Journal,2001,99(37):49-58.
[3] 周 密.油氣微滲漏的高光譜信息提取研究[D].北京:中國地質(zhì)大學(xué)(北京),2007.
Zhou M.Study on the Hyperspectral Information Extraction of Oil and Gas Microseepage[D].Beijing:China University of Geosciences (Beijing),2007.
[4] 李 娜,周 萍.基于ASTER數(shù)據(jù)和反射光譜分析的烴類微滲漏信息提取[J].地質(zhì)力學(xué)學(xué)報(bào),2015,21(2):218-227.
Li N,Zhou P.Hydrocarbon microleakage information extraction based on ASTER data and reflectance spectroscopy[J].Journal of Geomechanics,2015,21(2):218-227.
[5] 陳圣波,趙 靚,王晉年.黃土覆蓋區(qū)油氣微滲漏地表蝕變高光譜特征響應(yīng)機(jī)理研究[J].光譜學(xué)與光譜分析,2012,32(6):1616-1619.
Chen S B,Zhao L,Wang J N.Study on hyperspectral response to surface alteration by oil and gas microseepage under loess area[J].Spectroscopy and Spectral Analysis,2012,32(6):1616-1619.
[6] Kühn F,Oppermann K,H?rig B.Hydrocarbon index:An algorithm for hyperspectral detection of hydrocarbons[J].International Journal of Remote Sensing,2004,25(12):2467-2473.
[7] 王建民,王佳媛.鄂爾多斯盆地伊陜斜坡上的低幅度構(gòu)造與油氣富集[J].石油勘探與開發(fā),2013,40(1):49-57.
Wang J M,Wang J Y.Low-amplitude structures and oil-gas enrichment on the Yishaan Slope, Ordos Basin[J].Petroleum Exploration and Development,2013,40(1):49-57.
[8] 李振華,陳占軍,余 嵐,等.慶陽地區(qū)上古生界天然氣資源潛力評(píng)估[J].地下水,2015,37(6):245-247.
Li Z H,Chen Z J,Yu L,et al.Study on the potential of forming reservoir in Upper Paleozoic layer,Qing-yang area[J].Ground Water,2015,37(6):245-247.
[9] 楊 華,付金華,劉新社,等.鄂爾多斯盆地上古生界致密氣成藏條件與勘探開發(fā)[J].石油勘探與開發(fā),2012,39(3):295-303.
Yang H,Fu J H,Liu X S,et al.Accumulation conditions and exploration and development of tight gas in the Upper Paleozoic of the Ordos Basin[J].Petroleum Exploration and Development,2012,39(3):295-303.
[10] 張?zhí)煊?基于高光譜遙感的黃土覆蓋區(qū)油氣微滲漏信息提取[D].長(zhǎng)春:吉林大學(xué),2015.
Zhang T Y.Extraction of Micro-Seepage of Oil and Gas Information in Loess Covered Area Based on Hyperspectral Remote Sensing[D].Changchun:Jilin University,2015.
[11] 易 歡,李健強(qiáng),韓海輝,等.遙感技術(shù)在阿爾金貝殼灘地區(qū)礦產(chǎn)資源綜合調(diào)查中的應(yīng)用[J].中國地質(zhì)調(diào)查,2016,3(4):1-5.
Yi H,Li J Q,Han H H,et al.Application of remote sensing in integrated survey on mineral exploration in Beketan,Altyn[J].Geological Survey of China,2016,3(4):1-5.
[12] 閆柏琨,董新豐,王 喆,等.航空高光譜遙感礦物信息提取技術(shù)及其應(yīng)用進(jìn)展——以中國西部成礦帶調(diào)查為例[J].中國地質(zhì)調(diào)查,2016,3(4):55-62.
Yan B K,Dong X F,Wang Z,et al.Mineral information extraction technology by airborne hyperspectral remote sensing and its application progress:An example of mineralization belts of western China[J].Geological Survey of China,2016,3(4):55-62.