一種多屬性色彩處理地質體檢測方法及裝置制造方法
【專利摘要】本發明涉及一種多屬性色彩處理地質體檢測方法及裝置，其中，該方法包括：獲取格式為SEGY的疊后地震數據；利用格式為SEGY的疊后地震數據進行地質層位解釋，獲得時間層位數據；利用格式為SEGY的疊后地震數據進行地震屬性計算，獲得地震屬性數據；根據所述時間層位數據，從所述地震屬性數據中提取沿層屬性值；對所述沿層屬性值進行歸一化T變換，獲得IR、IG、IB；利用IR、IG、IB獲得RH、RS、RV；對RS進行邊緣提取，獲得結果信息SG；對RH、結果信息SG、RV形成的圖像進行反變換到RGB域，完成地質體檢測。
【專利說明】一種多屬性色彩處理地質體檢測方法及裝置
【技術領域】
[0001] 本發明涉及圖像處理【技術領域】，特別涉及一種多屬性色彩處理地質體檢測方法及 裝置。
【背景技術】
[0002] 1971年，Balch將地震數據進行彩色顯示利用彩色圖像的視覺效果直接識別異常 區域成為地震資料解釋的一種有效手段常規地震屬性彩色顯示技術是通過某種變換將屬 性數值映射成彩色圖像〃一次顯示一種屬性但對于多個地震屬性〃這種單個屬性逐一彩色 顯示方法不能很好地反映整體趨勢"突出區域性異常為了解決這一問題"基于顏色空間的 多屬性融合技，被引入屬性分析中。1984年，Onstott將RGBA顏色融合技術應用于遠、中、 近炮檢距地震數據的顯示，來描述地震數據的AVO特征。Liu(2007)提出一種根據自定義頻 譜范圍，對譜分解后的不同頻率數據體，采用基于余弦函數變換的RGBA顏色融合技術。Guo 等（2006)將主成分分析（PCA)技術和RGBA顏色融合技術聯合應用于譜分解數據，提高了 譜分解數據對河道的識別能力。2008年，Henderson利用RGB色彩融合對多屬性直接機型 地質體解釋。目前RGB多屬性融合技術已經廣泛應用于各種檢測。但由于融合過程與屬性 的選擇有直接連寫，容易受到源數據的影響，導致預測地質體邊界模糊不清。
[0003] 綜上所述，在現有技術中，通常利用地震屬性預測地質體分布，其邊界往往難以確 定，但由于融合過程與屬性的選擇有直接聯系，容易受到源數據的影響，導致預測地質體邊 界模糊不清。
【發明內容】
[0004] 為解決現有技術的問題，本發明提出一種多屬性色彩處理地質體檢測方法及裝 置，解決了數據體在色彩融合時產生的模糊效應，增強地震數據圖像存在明顯的結構信息， 精細刻畫河道、鹽丘、溶洞等重要地質體。
[0005] 為實現上述目的，本發明提供了一種多屬性色彩處理地質體檢測方法，該方法包 括：
[0006] 獲取常規格式疊后地震數據；
[0007] 利用常規格式疊后地震數據進行地質層位解釋，獲得時間層位數據；
[0008] 利用常規格式疊后地震數據進行地震屬性計算，獲得地震屬性數據；
[0009] 根據所述時間層位數據，從所述地震屬性數據中提取沿層屬性值；
[0010] 對所述沿層屬性值進行歸一化變換，獲得IR、IG、IB ;
[0011] 利用 IR、IG、IB 獲得 RH、RS、RV;
[0012] 對RS進行邊緣提取，獲得結果信息GS ;
[0013] 對RH、結果信息GS、RV形成的圖像進行反變換到RGB域，獲得RGB域圖像，完成地 質體檢測。
[0014] 優選地，所述獲得IR、IG、IB的表達式為：
[0015] [IR(z), IG(z), IB(z)] = T [VI (z), V2 (z), V3 (z)]
[0016] T = 255 X (vi - vmin)/ (vmax - vmin)
[0017] 其中，z表示某一屬性值的平面位置，T變換采用線性變換，vi包括VI，V2，V3, VI， V2, V3均表示為數據分量；vmin和vmax分別為VI，V2, V3中數據分量的最小值和最大值。
[0018] 優選地，所述獲得RH、RS、RV的表達式為：
[0019] RH = 60X (IG-IB)/(max(IR, IG, IB)-min(IR, IG, IB))
[0020] RS = (max (IR, IG, IB)-min (IR, IG, IB))/max (IR, IG, IB)
[0021] RV = max (IR, IG, IB) 〇
[0022] 優選地，所述結果信息GS的表達式為：
[0023] GS = RS*Sobel
[0024] 其中，
【權利要求】
1. 一種多屬性色彩處理地質體檢測方法，其特征在于，該方法包括： 獲取常規格式疊后地震數據； 利用常規格式疊后地震數據進行地質層位解釋，獲得時間層位數據； 利用常規格式疊后地震數據進行地震屬性計算，獲得地震屬性數據； 根據所述時間層位數據，從所述地震屬性數據中提取沿層屬性值； 對所述沿層屬性值進行歸一化變換，獲得IR、IG、IB; 利用IR、IG、IB獲得RH、RS、RV; 對RS進行邊緣提取，獲得結果信息GS; 對RH、結果信息GS、RV形成的圖像進行反變換到RGB域，獲得RGB域圖像，完成地質體 檢測。
2. 如權利要求1所述的方法，其特征在于，所述獲得IR、IG、IB的表達式為： [IR(z),IG(z),IB(z)] =T[Vl(z),V2(z),V3(z)] T= 255X(vi-vmin)/(vmax-vmin) 其中，z表示某一屬性值的平面位置，T變換采用線性變換，vi包括VI，V2,V3,VI，V2，V3均表示為數據分量；vmin和vmax分別為VI，V2,V3中數據分量的最小值和最大值。
3. 如權利要求1所述的方法，其特征在于，所述獲得RH、RS、RV的表達式為： RH= 60X(IG-IB)/(max(IR,IG,IB)-min(IR,IG,IB)) RS=(max(IR,IG,IB)-min(IR,IG,IB))/max(IR,IG,IB) RV=max(IR,IG,IB)。
4. 如權利要求I所述的方法，其特征在于，所述結果信息GS的表達式為： GS=RS*Sobel '-1οΓ 其中，&?/= -2 0 2 -1 0 1 *- -1〇
5. 如權利要求1所述的方法，其特征在于，所述RGB域圖像的表達式為： IR=RVX(I-GS) IG=RVX(1-fXGS) IB=RVX(l-(l-f)XGS)。
6. -種多屬性色彩處理地質體檢測裝置，其特征在于，該裝置包括： 疊后地震數據獲取單元，用于獲取常規格式疊后地震數據； 時間層位數據獲取單元，用于利用常規格式疊后地震數據進行地質層位解釋，獲得時 間層位數據； 地震屬性數據獲取單元，用于利用常規格式疊后地震數據進行地震屬性計算，獲得地 震屬性數據； 沿層屬性值獲取單元，用于根據所述時間層位數據，從所述地震屬性數據中提取沿層 屬性值； 變換單元，用于對所述沿層屬性值進行歸一化變換，獲得IR、IG、IB; HSV信息獲取單元，用于利用IR、IG、IB獲得RH、RS、RV; 邊緣提取單元，用于對RS進行邊緣提取，獲得結果信息GS; 反變換單元，用于對RH、結果信息GS、RV形成的圖像進行反變換到RGB域，獲得RGB域 圖像，完成地質體檢測。
7. 如權利要求6所述的裝置，其特征在于，所述變換單元獲得的IR、IG、IB的表達式 為： [IR(z)，IG(z)，IB(z)] =T[V1 (z)，V2 (z)，V3(zO] T= 255X(vi-vmin)/(vmax-vmin) 其中，z表示某一屬性值的平面位置，T變換采用線性變換，vi包括VI，V2,V3,VI，V2，V3均表示為數據分量；vmin和vmax分別為VI，V2,V3中數據分量的最小值和最大值。
8. 如權利要求6所述的裝置，其特征在于，所述HSV信息獲取單元獲得的RH、RS、RV的 表達式為： RH= 60X(IG-IB)/(max(IR,IG,IB)-min(IR,IG,IB)) RS=(max(IR,IG,IB)-min(IR,IG,IB))/max(IR,IG,IB) RV=max(IR,IG,IB) 。
9. 如權利要求6所述的裝置，其特征在于，所述邊緣提取單元獲得的結果信息GS的表 達式為： GS=RS*Sobel "-1οΓ 其中，= -2 0 2 一1 0 1 - O
10. 如權利要求6所述的裝置，其特征在于，所述反變換單元得到的RGB域圖像的表達 式為： IR=RVX(I-RS) IG=RVX(1-fXRS) IB=RVX(l-(l-f)XRS)。
【文檔編號】G06T7/00GK104459785SQ201410852343
【公開日】2015年3月25日 申請日期:2014年12月31日 優先權日:2014年12月31日
【發明者】楊慶, 楊午陽, 高建虎, 周春雷, 魏新建, 何欣, 王恩利, 李琳, 鄢高韓, 王萬里, 祿娟, 李三家, 陳德武 申請人:中國石油天然氣股份有限公司