張喜林 范德江 劉 明 王 亮

(中國海洋大學(xué)海洋地球科學(xué)學(xué)院海底科學(xué)與探測技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室 山東青島 266100)

沉積物巖芯X光片圖像灰度數(shù)值及其影響因素①

張喜林 范德江 劉 明 王 亮

(中國海洋大學(xué)海洋地球科學(xué)學(xué)院海底科學(xué)與探測技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室 山東青島 266100)

沉積體的X光影像包含了沉積物密度、結(jié)構(gòu)、構(gòu)造等眾多信息,對X光影像進(jìn)行數(shù)字圖像分析可快速提取這些信息,成為進(jìn)行沉積記錄和古環(huán)境演變研究的有力手段?；陂L江水下三角洲沉積物巖芯的X光影像,利用Matlab軟件平臺把X光影像轉(zhuǎn)換成數(shù)字灰度圖像,進(jìn)而提取了巖芯灰度值;結(jié)合該巖芯的粒度等沉積學(xué)屬性,探討了灰度圖像的影響因素。結(jié)果表明:沉積物巖芯的X光圖像灰度值是沉積物組成、結(jié)構(gòu)和構(gòu)造的綜合反映,其中粒度對灰度值起到控制性的影響。平均粒徑與灰度值之間為良好的線性關(guān)系;砂粒級、粉砂粒級含量與灰度之間存在較好的正相關(guān),而黏土粒級含量與灰度存在很好的負(fù)相關(guān)。沉積物巖芯X光圖像的灰度值可以作為該巖芯的粒度的良好替代性指標(biāo)。

沉積物巖芯 X光圖像 灰度 粒度 影響因素

0 引言

沉積物巖芯的X光負(fù)片圖像是沉積物密度、結(jié)構(gòu)和構(gòu)造等屬性的綜合體現(xiàn),由于它直觀、形象和易于獲得,所以在沉積學(xué)研究中得到廣泛應(yīng)用[1]。針對沉積物巖芯而言,圖像分析具有覆蓋面廣、連續(xù)性好、信息量大、效率高和無損壞成像等特點(diǎn),這比將沉積物巖芯剖分、以一定間隔取樣、分析沉積物樣品、得到相關(guān)信息的傳統(tǒng)研究途徑具有明顯的優(yōu)勢。當(dāng)今,沉積體的高分辨記錄研究方興未艾,研究時(shí)間分辨率從千年際-百年際-十年際-年代際-季節(jié)性不斷提高[2,3],這對研究手段提出了極高的要求。很顯然,僅靠傳統(tǒng)的巖芯研究手段難以適應(yīng)當(dāng)今高分辨率沉積記錄研究的需要,而沉積物巖芯圖像的無縫連續(xù)信息正好彌補(bǔ)了傳統(tǒng)研究方法的不足。

X光圖像很早就應(yīng)用到沉積學(xué)研究,DeMaster等[4]于1985對中國東海陸架沉積作用研究中就使用了X光成像技術(shù),后來的國內(nèi)外許多學(xué)者使用X光圖像揭示沉積體中的層理構(gòu)造[5]、事件沉積層[6]、生物擾動(dòng)構(gòu)造[7,8];近期,Genty等[9]、秦小光等[10]利用數(shù)字圖像處理技術(shù)研究了石筍生長、微層特征等,孫東懷等[1]利用X圖像灰度值研究了南海珊瑚的生長,并試圖建立灰度與海水溫度聯(lián)系。迄今為止,由于對X光圖像進(jìn)行數(shù)字化處理需要專業(yè)軟件以及大型儀器的支持,僅有少數(shù)學(xué)者開展了針對結(jié)構(gòu)和成分單一的如珊瑚、石筍等特定沉積體的灰度數(shù)字圖像研究,絕大部分關(guān)于沉積物巖芯的X光圖像的研究仍然停留在對影像的直接識別上,缺乏從定量化方面進(jìn)行深入的探討。為此,本文基于長江口沉積物巖芯X光影像,利用Matlab軟件平臺把X光影像轉(zhuǎn)換成數(shù)字灰度圖像,提取了巖芯的數(shù)字灰度值,結(jié)合該巖芯的沉積學(xué)特征探討了灰度值的主要影響因素。

1 沉積物巖芯的灰度圖像分析原理

X射線又稱倫琴射線,波長介于10-2～102?之間,具有很強(qiáng)的穿透能力,并使底片曝光。由于被照射物體結(jié)構(gòu)、密度和厚度等的差別,X線被吸收的程度不同,所以到達(dá)熒屏或膠片上的X線的通量有差異,形成了黑白對比不同的影像。

通過X光片得到的圖像實(shí)質(zhì)上是顏色值或灰度值的數(shù)字矩陣,單色圖像的每個(gè)像點(diǎn)即為該點(diǎn)的灰度值。這類圖像通常顯示為從最暗的黑色到最亮的白色,共有256級灰度。圖像分析的目的是從這一矩陣中提取相應(yīng)的數(shù)字信息。

沉積物巖芯是在特定時(shí)期內(nèi)形成的,該時(shí)期內(nèi)的氣候、沉積環(huán)境的變化都可能導(dǎo)致巖芯物質(zhì)組成、結(jié)構(gòu)和構(gòu)造上的改變,并進(jìn)而導(dǎo)致其X光影像的差異。很顯然,沉積物巖芯的灰度圖像的影響因素遠(yuǎn)較成分、結(jié)構(gòu)單一的珊瑚、石筍等X圖像復(fù)雜。

2 沉積物巖芯灰度圖像處理

2.1 沉積物巖芯

本研究所使用的沉積物巖芯(18站:122°37.132' E,31°00.990'N)為2006年6月中國海洋大學(xué)“東方紅2號”科學(xué)考察船對長江口及其鄰近海域進(jìn)行調(diào)查時(shí)采用重力取樣器采集的,它位于長江水下三角洲的前緣部分,該處水深20.4 m。具體位置見圖1。

圖1 研究區(qū)及采樣站位(TWWC-臺灣暖流;YSCC-黃海沿岸流;ECSCC-浙閩沿岸流.陰影區(qū)為泥質(zhì)沉積體)Fig.1 Study area and sampling site(TWC-the TaiwanWarm Current;YSCC-the Yellow Sea Coastal Current;ECSCC-the EastChina Sea CoastalCurrent.The shaded areas denotemuddy deposition)

該站沉積物巖芯巖性以粉砂和黏土為主,兩者的含量之和在98%左右,砂含量很少,不足2%。個(gè)別層段粉砂較多,如37～38 cm處和93～95 cm處的粉砂含量明顯增多,形成透鏡狀粉砂夾層。垂向上,粉砂和黏土相對含量的波動(dòng)明顯,并因此顯示紋層構(gòu)造(圖2)?？傮w上,該站沉積物巖芯沉積物組成較穩(wěn)定,層序保存完整,反映了較為穩(wěn)定的沉積動(dòng)力環(huán)境。

2.2 沉積物巖芯X光成像

圖2 18站沉積物巖芯的粒級組成Fig.2 Grain-size compositions of Core 18

采用醫(yī)用SIEMENS-FX型X射線成像系統(tǒng)對18站柱狀沉積物進(jìn)行X射線成像,儀器工作條件為16 mA/70keV,獲得的該巖芯的X光負(fù)片圖像見圖3A。

2.3 灰度值提取技術(shù)

2.3.1 灰度圖像預(yù)處理

沉積物巖芯X光負(fù)片成像過程中受到采樣管管壁及拍片時(shí)參照物等因素的影響,會導(dǎo)致X光負(fù)片灰度值的改變,所以必須對原始X光片用進(jìn)行預(yù)處理。本研究使用Photoshop軟件進(jìn)行X光負(fù)片預(yù)處理,包括灰度轉(zhuǎn)換、對比度優(yōu)化等過程。預(yù)處理后的圖像供后續(xù)灰度值提取使用(圖3B)。

2.3.2 沉積物巖芯X光負(fù)片圖像灰度提取

出于灰度數(shù)據(jù)與實(shí)測沉積物參數(shù)的嚴(yán)格對應(yīng)的考慮,本研究選取了該柱樣0～55 cm段作為研究對象。首先利用Matlab7.0軟件中的imread函數(shù)讀入處理后的灰度圖像,然后利用圖像轉(zhuǎn)換函數(shù)rgb2gray將該圖像轉(zhuǎn)換為灰度圖像(圖4),進(jìn)而提取該灰度圖像的灰度值,它是一個(gè)550 X30的灰度值矩陣,代表的垂向分辨率0.1mm。限于篇幅,我們僅列出矩陣中的6行數(shù)據(jù),見表1。

圖3 18站沉積物巖芯X光負(fù)片圖像(A.原始X光負(fù)片圖像;B.經(jīng)預(yù)處理后的X光負(fù)片圖像)Fig.3 X-ray negative of Core 18(A.original X-ray negative;B.processed X-ray negative)

圖4 18站X光負(fù)片灰度圖像Fig.4 The gray image of X-ray radiograph negative of Core 18

表1 18站X光負(fù)片灰度值矩陣Table1 Thematrix of the gray image of X-ray radiograph negative of Core 18(該矩陣的行按照沉積物巖芯的長度進(jìn)行排列,列按照巖芯的橫截面進(jìn)行排列)

選取灰度矩陣中的第5、10、15、20、25、30等六列做圖,得到6條隨深度變化的灰度曲線(圖5)。它們表明:該巖芯X光負(fù)片灰度值大部分集中在70～100之間,且這6條曲線總體變化趨勢非常一致,它們在11 cm、26 cm、38 cm等層位都出現(xiàn)明顯的灰度峰值。

圖5 18站沉積物巖芯灰度曲線(A-F分別代表矩陣中的第5、10、15、20、25、30六列灰度變化曲線)Fig.5 The gray values curves of the Core 18

相關(guān)性分析表明這6列灰度曲線之間的相關(guān)性很強(qiáng),僅在巖芯邊緣(30列)相關(guān)性變差(表2)。巖芯邊緣相關(guān)性變差可能與管壁兩側(cè)的巖芯變形或者與X光成像時(shí)的邊緣效應(yīng)有關(guān)。因此,我們選用位于巖芯中央的C列的灰度值作為該巖芯灰度垂向分布的代表,并據(jù)此探討影響巖芯灰度值制約因素。

表2 巖芯中6列灰度值的相關(guān)系數(shù)表(n=550)Table2 Correlations of the 6 rows of the gray values(n=550)

3 影響沉積物巖芯X光負(fù)片影像灰度的因素

通過X射線成像得到的負(fù)片,其灰度值的變化實(shí)際上就是X射線穿透強(qiáng)度的表現(xiàn)。影響該圖像灰度的因素十分復(fù)雜,除了采樣、樣品保存以及灰度成像過程的人為影響外,主要受到沉積物屬性的影響。我們主要從以下幾個(gè)方面進(jìn)行探討。

3.1 沉積物粒度組成的影響

該巖芯在室內(nèi)以0.25 cm間隔取樣并使用激光粒度分析儀進(jìn)行了粒度分析,其平均粒徑見圖6。通過該巖芯的粒徑和灰度的對比分析發(fā)現(xiàn)兩者之間的變化趨勢基本一致,比如在24 cm和38 cm處都出現(xiàn)了峰值,并且在38 cm處為最大的峰值。

圖6 18站沉積物巖芯平均粒徑和灰度曲線(A.灰度;B.平均粒徑)Fig.6 The curves ofmean grain-size and gray value in Core 18(A.gray value;B.mean grain-size)

為了進(jìn)一步說明平均粒徑和灰度之間的關(guān)系,對不同層位沉積物的平均粒徑及其對應(yīng)的灰度值作相關(guān)性分析,結(jié)果表明兩者之間存在良好的線性相關(guān)(圖7),有如下關(guān)系:

圖7 灰度值與平均粒徑相關(guān)圖Fig.7 The correlation between gray value and mean grain-size

根據(jù)X光負(fù)片成像原理,當(dāng)物質(zhì)密度大時(shí)顏色較淺、偏白;反之則較深、較暗。平均粒徑指示了沉積物粒徑頻率分布的中心趨向,平均粒徑越大,沉積物總體組成越粗,相對密度大,吸收的X光就越多,能透過的X光就越少,反映在X光負(fù)片上顏色淺,灰度值相應(yīng)的增加;反之亦然。由于X光灰度值與沉積物平均粒徑具有非常好的線性相關(guān)性,完全有可能使用該關(guān)系式計(jì)算粒度數(shù)據(jù)。

接下來進(jìn)一步分析沉積物灰度與沉積物中的砂、粉砂、黏土各個(gè)組分之間的相關(guān)關(guān)系,如圖8所示。

從圖8可以看到與灰度值與粉砂、砂粒級都為正相關(guān)關(guān)系,與粉砂的相關(guān)性優(yōu)于砂,相關(guān)系數(shù)平方分別為0.69、0.39?；叶戎蹬c粘土粒級則為良好的負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)平方為0.68。這進(jìn)一步說明沉積物中砂、粉砂吸收的X光多,在X光負(fù)片上體現(xiàn)的顏色淺;相反地,黏土吸收的X光少,在X光負(fù)片上所體現(xiàn)出的顏色深。

考慮到平均粒度是沉積物粒級組成的綜合反映,據(jù)此計(jì)算該孔粒度導(dǎo)致的灰度波動(dòng)方差貢獻(xiàn)為47.67,占總方差貢獻(xiàn)的94.13%。表明粒度是控制X光負(fù)片灰度的決定性因素。

3.2 其它影響因素分析

雖然平均粒徑對灰度值的變化有著非常重要的影響,但是并未完全控制灰度值的變化。通過粒度歸一化處理,可以消除粒度對灰度的影響,即:以18站巖芯的總平均粒徑為標(biāo)準(zhǔn),把不同層位的平均粒徑換算成總平均粒徑,該層位的歸一化灰度則依據(jù)上述得到的灰度與平均粒徑的關(guān)系式進(jìn)行換算得到。該方法也稱為“當(dāng)量粒徑法”[11]。剔除粒度這一主要的影響后,得到灰度值隨深度變化的曲線圖(如圖9)。

圖8 灰度值與砂、粉砂、黏土粒級的相關(guān)圖(Ⅰ.灰度值和砂的相關(guān)圖,Ⅱ.灰度值和粉砂的相關(guān)圖,Ⅲ.灰度值和黏土的相關(guān)圖)Fig.8 The correlations between gray value and sand,silt and clay(Ⅰ.gray value and sand content,Ⅱ.gray value and silt content,Ⅲ.gray value and clay content)

從圖9可以清晰的看出,灰度值的變化存在明顯的高頻和低頻的周期性波動(dòng),高頻波動(dòng)的間隔在1～ 3 cm左右,低頻波動(dòng)間隔在10～12 cm左右,根據(jù)該站的沉積速率2.1 cm/a[3],該巖芯高頻、低頻波動(dòng)對應(yīng)的時(shí)間約為1年及4年。1年的波動(dòng)與長江口沉積作用的季節(jié)性過程相符合,而4年的波動(dòng)則可能與ENSO事件有關(guān)。

圖9 18站粒度歸一化后的灰度曲線Fig.9 The grain-size normalized gray curve of core 18

無論是何種原因造成的周期性波動(dòng),都是通過沉積物屬性的差異體現(xiàn)出來的。以季節(jié)性變化為例,夏秋季節(jié),長江輸入的水量、砂量大,沉積物顆粒細(xì),同時(shí)該季節(jié)生物生產(chǎn)力高,沉積物中硅藻等生物組分含量高;冬春季節(jié),長江輸入的水量、砂量少,沉積物顆粒較粗,以硅藻為代表的生物生產(chǎn)力小,沉積物中硅藻等生物組分含量低[3]。正是這種季節(jié)性過程導(dǎo)致了沉積物巖芯組成的改變,并影響到X光影像。4年的周期變化應(yīng)該也與周期性的沉積物組成的改變相聯(lián)系。受到目前所掌握的數(shù)據(jù)的限制,現(xiàn)在尚無法確定是哪種具體的組分影響到X光圖像灰度的1年及4年的周期性波動(dòng)。這有待下一步深入研究。

4 結(jié)論

(1)沉積物巖芯的X射線灰度圖像信息量大、信息綜合性強(qiáng),可以通過該圖像直接獲取高分辨率沉積記錄信息。

(2)沉積物粒度與X射線負(fù)片圖像的灰度值有著密切的關(guān)系,其中平均粒徑和灰度具有良好的線性關(guān)系,砂、粉砂粒級與灰度值之間為正相關(guān),黏土粒級與灰度值存在較好的負(fù)相關(guān)。

(3)除了粒度的影響之外,沉積物X射線負(fù)片圖像灰度值可能還受到沉積物中礦物、生物組分等因素的影響。

References)

1 孫東懷,劉禹,譚明。古環(huán)境記錄的數(shù)字圖像分析及應(yīng)用[J]?？茖W(xué)通報(bào),2002,47(21):1613-1619[Sun Donghuai,Liu Yu,Tan Ming.The application and analysis of the digital image of paleoenvironmental records[J].hinese Science Bulletin,2002,47(21):1613-1619]

2 楊作升,陳曉輝。百年來長江口泥質(zhì)區(qū)高分辨率沉積粒度變化及影響因素探討[J]。第四紀(jì)研究,2007,27(5):690-699[Yang Zuosheng,Cheng Xiao.Centurial high resolution records of sediment grain-size variation in themud area off the Changjiang(Yangtze River) estuary and its influencial factors[J].uaternary Science,2007,27 (5):690-699]

3 范德江,齊紅艷,孫曉霞,等。長江水下三角洲現(xiàn)代沉積速率厘定的新途徑:季節(jié)性Bio-Si記錄[J]。中國海洋大學(xué)報(bào),2009,39 (5):1025-1028[Fan Dejiang,Qi Hongyan,Sun Xiaoxia,et al.A new method to evaluating sedimentary rate on the Yangtze River subaqueous delta:Seasonal sedimentary record of Bio-Si[J].eriodical of Ocean University of China,2009,39(5):1025-1028]

4 DeMaster D J,McKee B A,Nittrouer C A,et al.Rates of sediment accumulation and particle reworking based on radiochemical measurements from continental shelf deposits in the East China Sea[J].ontinental Shelf Research,1985,4(1-2):143-158

5 Kniskern T A,Kuehl S A.Spatial and temporal variability of seabed disturbance in the York River subestuary[J].stuarine,Coastal and Shelf Science,2003,58:37-55

6 Bentley S J,Sheremet A,JaegerJM.Event sedimentation,bioturbation,and preserved sedimentary fabric:Field andmodel comparisons in three contrastingmarine settings[J].arine Geoglogy,2006,26(17-18):2108-2124

7 范德江,徐琳,齊紅艷。長江水下三角洲淺表沉積層中的生物擾動(dòng)構(gòu)造[J]。海洋與湖沼,2008.39(6):578-584[Fan Dejiang,Xu Lin, Qi Hongyan.Bioturbation in surficial sediment of subaqueous Changjiang river delta[J].ceanologia et Limnologia Sinica,2008, 39(6):578-584]

8 楊群慧,周懷陽,季福武,等。海底生物擾動(dòng)作用及其對沉積過程和記錄的影響[J]。地球科學(xué)進(jìn)展,2008,23(9):932-941[Yang Qunhui,Zhou Huaiyang,Ji Fuwu,et al.Bioturbation in seabed sediments and its effects on marine sedimentary processes and records[J].dvances in Earth Sciences,2008,23(9):932-941]

9 Genty D,Baker A,Maire R.Comparison of annual luminescent and visible laminae in stalagmites,Comptes Rendus DeL'Academie des Sciences,SerⅡ.ascicule A[J].ciences Dela Terreet des Planetes,1997,325:193-200

10 秦小光,劉東生,譚明,等。北京石花洞石筍微層灰度變化特征及其氣候意義[J]。中國科學(xué):D輯,2000,30(3):240-248[Qin Xiaoguang,Liu Dongsheng,Tan Ming,et al.The micro level gray scale characteristics of Beijing Shihuadong stalagmite and its climatic Significance[J].cience in China:Series D,2000,30(3):240-248]

11 焦春文。工程類比法探討巖基沖刷模擬試驗(yàn)[J]。長江科學(xué)院院報(bào),1998,15(2):42-45[Jiao Chunwen.Study on anti scour test in rock river bed by project comparison method[J].ournal of Yangtze River Scientific Research Institute,1998,15(2):42-45]

The Gray Value of X-ray Radiograph in Sediment Core and Its Influencing Factors

ZHANG Xi-lin FAN De-jiang LIU Ming WANG Liang
(College of M arine Geosciences,Key Laboratory of Submarine Geosciences and Technology of M inistry of Education, Ocean University of China,Qingdao Shandong 266100)

X-ray radiograph of sedimentary core containsmuch information,including sedimentary texture,density, structure,etc.We can quickly extract the information of climate and environment changes from sedimentary core by using digital image processingmethod,witchmay be a powerful research tools to exploring high-resolution sedimentary records and palaeo-environmentevolution.Based on thematlab software platform,we obtained the X-ray images of the sediment core taken from the subaqueous Yangtze River delta and changed them into digital gray images then for extracting their gray values.Combined with the sediment properties in the core such as grain size parameters,influencing factors of the digital gray valueswere discussed.The results indicate that the gray values of sediment core are influenced comprehensively by sediment composition,texture and structure,with being dominantly influenced by sediment grain-size.There is a good linear relationship between mean grain-size and gray value.Furthermore,there is a positive correlation between sand,silt and gray value.In contrary,it is of a good negative relationship between clay and gray value.Generally,the X-ray image gray value of sediment coremay be an ideal proxy of comprehensive sedimentary records.

sediment core;X-ray radiograph;gray value;grain-size;influencing factor

張喜林 男 1984出生 碩士研究生 海洋沉積學(xué)

范德江 E-mail:difan@ouc.edu.cn

P512.2

A

1000-0550(2012)02-0318-07

①國家自然科學(xué)基金項(xiàng)目(批準(zhǔn)號:40976020,41030856)和國家重大基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目(編號:2010CB951202)聯(lián)合資助。

2010-08-01;收修改稿日期:2011-05-05