CN104459785A - Method and device for detecting multi-attribute color processing geologic body - Google Patents

Abstract

The invention relates to a method and device for detecting a multi-attribute color processing geologic body. The method comprises the steps that post-stacked earthquake data with the SEGY format are obtained; the geologic horizon explanation is carried out through the post-stacked earthquake data with the SEGY format, and time horizon data are obtained; the earthquake attribute calculation is carried out through the post-stacked earthquake data with the SEGY format, and the earthquake attribute data are obtained; the stratum attribute value is obtained from the earthquake attribute data according to the time horizon data; the normalization T transform is carried out on the stratum attribute value, and IR, IG and IB are obtained; RH, RS and RV are obtained through IR, IG and IB; edge extraction is carried out on RS, and result information SG is obtained; inverse transformation is carried out on an image formed by RH, result information SG and RV to the RGB area, and the detection of the geographic body is finished.

Description

A kind of many attributes colors countenance geologic body detection method and device

Technical field

The present invention relates to technical field of image processing, particularly one many attributes colors countenance geologic body detection method and device.

Background technology

1971, Balch geological data is carried out colored show a kind of effective means conventional seismic attribute color display technology utilizing the visual effect Direct Recognition abnormal area of coloured image to become seismic data interpretation be by certain conversion attribute value is mapped to coloured image " once show an attribute but for multiple seismic properties " this single attribute one by one color display method can not reflect that overall trend " outburst area sexual abnormality is in order to address this problem " merges skills based on many attributes of color space well, be introduced in attributive analysis.1984, RGBA color blend technology is applied to far by Onstott, in, the display of nearly geophone offset geological data, the AVO feature of geological data is described.Liu (2007) proposes a kind of according to self-defined spectral range, to the different frequency data volume after spectral factorization, adopts the RGBA color blend technology based on cosine function conversion.Principal component analysis (PCA) (PCA) technology and the technical tie-up of RGBA color blend are applied to spectral factorization data by Guo etc. (2006), improve the recognition capability of spectral factorization data to river course.2008, Henderson utilized rgb color to merge and explains the direct type geologic body of many attributes.Current RGB many attributes integration technology has been widely used in various detection.But because the selection of fusion process and attribute has direct write the two or more syllables of a word together, be easily subject to the impact of source data, cause predictably that plastid obscurity boundary is unclear.

In sum, in the prior art, usually utilize seismic properties predictably plastid distribution, its border is often difficult to determine, but contacts directly due to the selection of fusion process and attribute, is easily subject to the impact of source data, causes predictably that plastid obscurity boundary is unclear.

Summary of the invention

For solving the problem of prior art, the present invention proposes a kind of many attributes colors countenance geologic body detection method and device, solve the blurring effect that data volume produces when color bleeding, strengthen geological data image and there is obvious structural information, the importantly plastid such as meticulous depiction river course, salt dome, solution cavity.

For achieving the above object, the invention provides a kind of many attributes colors countenance geologic body detection method, the method comprises:

Obtain general format post-stack seismic data;

Utilize general format post-stack seismic data to carry out geologic horizon explanation, obtain time horizon bit data;

Utilize general format post-stack seismic data to carry out seismic properties calculating, obtain seismic attributes data;

According to described time horizon bit data, extract from described seismic attributes data along layer property value;

Be normalized conversion to described along layer property value, obtain IR, IG, IB;

IR, IG, IB is utilized to obtain RH, RS, RV;

Edge extracting is carried out to RS, obtains object information GS;

Contravariant is carried out to the image of RH, object information GS, RV formation and changes to RGB territory, obtain RGB area image, complete geologic body and detect.

Preferably, the described expression formula obtaining IR, IG, IB is:

[IR(z),IG(z),IB(z)]＝T[V1(z),V2(z),V3(z)]

T＝255×(vi–vmin)/(vmax–vmin)

Wherein, z represents the planimetric position of a certain property value, and T conversion adopts linear transformation, and vi comprises V1, and V2, V3, V1, V2, V3 are all expressed as data component; Vmin and vmax is respectively V1, the minimum value of data component and maximal value in V2, V3.

Preferably, the described expression formula obtaining RH, RS, RV is:

RH＝60×(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)。

Preferably, the expression formula of described object information GS is:

GS＝RS*Sobel

Wherein, $\mathrm{Sobel}=\left[\begin{array}{ccc}-1& 0& 1\\ -2& 0& 2\\ -1& 0& 1\end{array}\right].$

Preferably, the expression formula of described RGB area image is:

IR＝RV×(1-GS)

IG＝RV×(1-f×GS)

IB＝RV×(1-(1-f)×GS)。

For achieving the above object, present invention also offers a kind of many attributes colors countenance geologic body pick-up unit, this device comprises:

Post-stack seismic data acquiring unit, for obtaining general format post-stack seismic data;

Time horizon bit data acquiring unit, for utilizing general format post-stack seismic data to carry out geologic horizon explanation, obtains time horizon bit data;

Seismic attributes data acquiring unit, for utilizing general format post-stack seismic data to carry out seismic properties calculating, obtains seismic attributes data;

Along layer property value acquiring unit, for according to described time horizon bit data, extract from described seismic attributes data along layer property value;

Converter unit, for being normalized conversion to described along layer property value, obtains IR, IG, IB;

HSV information acquisition unit, obtains RH, RS, RV for utilizing IR, IG, IB;

Edge extracting unit, for carrying out edge extracting to RS, obtains object information GS;

Inverse transformation unit, the image for being formed RH, object information GS, RV carries out contravariant and changes to RGB territory, obtains RGB area image, completes geologic body and detects.

Preferably, the expression formula of IR, IG, IB of described converter unit acquisition is:

[IR(z),IG(z),IB(z)]＝T[V1(z),V2(z),V3(z0]

T＝255×(vi–vmin)/(vmax–vmin)

Wherein, z represents the planimetric position of a certain property value, and T conversion adopts linear transformation, and vi comprises V1, and V2, V3, V1, V2, V3 are all expressed as data component; Vmin and vmax is respectively V1, the minimum value of data component and maximal value in V2, V3.

Preferably, the expression formula of RH, RS, RV of described HSV information acquisition unit acquisition is:

RH＝60×(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)。

Preferably, the expression formula of the object information GS of described edge extracting unit acquisition is:

GS＝RS*Sobel

Wherein, $\mathrm{Sobel}=\left[\begin{array}{ccc}-1& 0& 1\\ -2& 0& 2\\ -1& 0& 1\end{array}\right].$

Preferably, the expression formula of RGB area image that described inverse transformation unit obtains is:

IR＝RV×(1-RS)

IG＝RV×(1-f×RS)

IB＝RV×(1-(1-f)×RS)。

Technique scheme has following beneficial effect: the technical program to identification river course, salt dome, solution cavity etc. importantly plastid have good recognition capability, particularly have outstanding effect to identification branch channel.On the basis of precise arrangements imaging, boundary information can be obtained from seismic data volume and seismic image, improve local geometric features structure in seismic properties, more careful portraying is carried out to target geological body, thus obtain more real geology, geophysics understanding, the reservoir prediction research work in auxiliary reservoir exploration stage.

Accompanying drawing explanation

In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.

Fig. 1 is one many attributes colors countenance geologic body detection method process flow diagram that the present invention proposes;

Fig. 2 is one many attributes colors countenance geologic body pick-up unit block diagram that the present invention proposes;

Fig. 3 is the time horizon bit data schematic diagram that the present embodiment obtains;

Fig. 4 is the 10hz frequency geological data section attribute data schematic diagram of the present embodiment;

Fig. 5 is the 20hz frequency geological data section attribute data schematic diagram of the present embodiment;

Fig. 6 is the 30hz frequency geological data section attribute data schematic diagram of the present embodiment;

Fig. 7 is the schematic diagram after the gray processing process of the chromaticity diagram that the present embodiment normalization conversion obtains;

Fig. 8 is the schematic diagram of HSV chromaticity diagram after gray proces that the present embodiment obtains;

Fig. 9 is the RS component section schematic diagram of the present embodiment based on Fig. 8;

Figure 10 is the section schematic diagram of the marginal date GS of the present embodiment;

Figure 11 is the schematic diagram after the gray processing process of the RGB area image of the present embodiment.

Embodiment

Below in conjunction with the accompanying drawing in the embodiment of the present invention, be clearly and completely described the technical scheme in the embodiment of the present invention, obviously, described embodiment is only the present invention's part embodiment, instead of whole embodiments.Based on the embodiment in the present invention, those of ordinary skill in the art, not making the every other embodiment obtained under creative work prerequisite, belong to the scope of protection of the invention.

The principle of work of the technical program is: based on many attributes convergence analysis technology of color space, and the uncontinuity feature of outstanding geological data, portrays geological objects boundary.Successively experience: RGB attribute merges, RGB to HSV changes, S component extracts, rim detection becomes figure with structural form RGB.Boundary information can be obtained from seismic data volume and seismic image, improve local geometric features structure in seismic properties, because geological data image exists obvious structural information, this technology has stronger applicability to geological data, in estimation partial discontinuous, present stronger architectural feature.Therefore, be suitable for portraying the target geological bodies such as tomography, crack, solution cavity and river course, be specially adapted to geologic body identification in seismic prospecting, thus improve the precision of reservoir prediction.Seismic survey work discloses the Main Means of underground structure, for obtaining geologic body image more clearly.The present invention effectively improves seismic prospecting prediction and discloses the precision of geologic body, oil-gas exploration, coal exploration, engineering exploration, etc. field have broad application prospects.

Based on above-mentioned principle of work, the present invention proposes a kind of many attributes colors countenance geologic body detection method, as shown in Figure 1.The method comprises:

Step 101): obtain general format post-stack seismic data;

Step 102): utilize general format post-stack seismic data to carry out geologic horizon explanation, obtain time horizon bit data;

Step 103): utilize general format post-stack seismic data to carry out seismic properties calculating, obtain seismic attributes data;

Step 104): according to described time horizon bit data, extract from described seismic attributes data along layer property value;

Step 105): be normalized conversion to described along layer property value, obtain IR, IG, IB;

Step 106): utilize IR, IG, IB to obtain RH, RS, RV;

Step 107): edge extracting is carried out to RS, obtains object information GS;

Step 108): contravariant is carried out to the image of RH, object information GS, RV formation and changes to RGB territory, obtain RGB area image, complete geologic body and detect.

Preferably, the described expression formula obtaining IR, IG, IB is:

[IR(z),IG(z),IB(z)]＝T[V1(z),V2(z),V3(z0](1)

T＝255×(vi–vmin)/(vmax–vmin)(2)

Wherein, z represents the planimetric position of a certain property value, and T conversion adopts linear transformation, and vi comprises V1, and V2, V3, V1, V2, V3 are all expressed as data component; Vmin and vmax is respectively V1, the minimum value of data component and maximal value in V2, V3.

Preferably, the described expression formula obtaining RH, RS, RV is:

RH＝60×(IG-IB)/(max(IR,IG,IB)-min(IR,IG,IB))(3)

RS＝(max(IR,IG,IB)-min(IR,IG,IB))/max(IR,IG,IB)(4)

RV＝max(IR,IG,IB)(5)

Preferably, the expression formula of described object information GS is:

GS＝RS*Sobel

Wherein, $\mathrm{Sobel}=\left[\begin{array}{ccc}-1& 0& 1\\ -2& 0& 2\\ -1& 0& 1\end{array}\right].$

Preferably, the expression formula of RGB area image that described inverse transformation obtains is:

IR＝RV×(1-GS)

IG＝RV×(1-f×GS)

IB＝RV×(1-(1-f)×GS)

Based on above-mentioned principle of work, present invention also offers a kind of many attributes colors countenance geologic body pick-up unit.As shown in Figure 2, this device comprises:

Post-stack seismic data acquiring unit 201, for obtaining general format post-stack seismic data;

Time horizon bit data acquiring unit 202, for utilizing general format post-stack seismic data to carry out geologic horizon explanation, obtains time horizon bit data;

Seismic attributes data acquiring unit 203, for utilizing general format post-stack seismic data to carry out seismic properties calculating, obtains seismic attributes data;

Along layer property value acquiring unit 204, for according to described time horizon bit data, extract from described seismic attributes data along layer property value;

Converter unit 205, for being normalized conversion to described along layer property value, obtains IR, IG, IB;

HSV information acquisition unit 206, obtains RH, RS, RV for utilizing IR, IG, IB;

Edge extracting unit 207, for carrying out edge extracting to RS, obtains object information GS;

Inverse transformation unit 208, the image for being formed RH, object information GS, RV carries out contravariant and changes to RGB territory, obtains RGB area image, completes geologic body and detects.

Below in conjunction with embodiment and the accompanying drawing of certain oil field A block.The embodiment of the present invention is described in further detail.

1), utilize conventional poststack seismologic record to carry out geologic horizon explanation, obtain time horizon bit data.As shown in Figure 3, be the time horizon bit data schematic diagram of the present embodiment acquisition.

2), utilize conventional poststack seismologic record to carry out seismic properties calculating, obtain seismic attributes data.As shown in Figure 4, be the 10hz frequency geological data section attribute data schematic diagram of the present embodiment.As shown in Figure 5, be the 20hz frequency geological data section attribute data schematic diagram of the present embodiment.As shown in Figure 6, be the 30hz frequency geological data section attribute data schematic diagram of the present embodiment.

3), to described be normalized conversion along layer property value, obtain IR, IG, IB and synthesize chromaticity diagram.Obtain chromaticity diagram after gray processing process as shown in Figure 7.Wherein,

[IR(z),IG(z),IB(z)]＝T[V1(z),V2(z),V3(z0]

T＝255×(vi–vmin)/(vmax–vmin)

4), utilize IR, IG, IB to obtain RH, RS, RV, be converted to HSV chromaticity diagram.Based on Fig. 7 chromaticity diagram conversion obtain HSV chromaticity diagram, HSV chromaticity diagram after gray proces as shown in Figure 8.Wherein

RH＝60×(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)

5), adopt Sobel operator to carry out convolution to RS and obtain marginal date GS.Wherein, Fig. 9 is the RS component section schematic diagram based on Fig. 8, and Figure 10 is the section schematic diagram of the marginal date GS of the present embodiment.Wherein, the expression formula of object information GS is:

GS＝RS*Sobel

Wherein, $\mathrm{Sobel}=\left[\begin{array}{ccc}-1& 0& 1\\ -2& 0& 2\\ -1& 0& 1\end{array}\right].$

6) the RGB area image, utilizing inverse transformation to obtain.As shown in figure 11, be the schematic diagram after the gray processing process of the RGB area image of the present embodiment.

Fig. 3 and Figure 11 compares known, after adopting the technical program process, good imaging is obtained to the trunk of paleo-streams, branch channel, to identification river course, salt dome, solution cavity etc. importantly plastid have good recognition capability, particularly have outstanding effect to identification branch channel.

Above-described embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only the specific embodiment of the present invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (10)

1. the colors countenance of attribute more than a geologic body detection method, it is characterized in that, the method comprises:

Obtain general format post-stack seismic data;

Utilize general format post-stack seismic data to carry out geologic horizon explanation, obtain time horizon bit data;

Utilize general format post-stack seismic data to carry out seismic properties calculating, obtain seismic attributes data;

According to described time horizon bit data, extract from described seismic attributes data along layer property value;

Be normalized conversion to described along layer property value, obtain IR, IG, IB;

IR, IG, IB is utilized to obtain RH, RS, RV;

Edge extracting is carried out to RS, obtains object information GS;

Contravariant is carried out to the image of RH, object information GS, RV formation and changes to RGB territory, obtain RGB area image, complete geologic body and detect.

2. the method for claim 1, is characterized in that, the described expression formula obtaining IR, IG, IB is:

[IR(z),IG(z),IB(z)]＝T[V1(z),V2(z),V3(z)]

T＝255×(vi–vmin)/(vmax–vmin)

Wherein, z represents the planimetric position of a certain property value, and T conversion adopts linear transformation, and vi comprises V1, and V2, V3, V1, V2, V3 are all expressed as data component; Vmin and vmax is respectively V1, the minimum value of data component and maximal value in V2, V3.

3. the method for claim 1, is characterized in that, the described expression formula obtaining RH, RS, RV is:

RH＝60×(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. the method for claim 1, is characterized in that, the expression formula of described object information GS is:

GS＝RS*Sobel

Wherein, $\mathrm{Sobel}\left[\begin{array}{ccc}-1& 0& 1\\ -2& 0& 2\\ -1& 0& 1\end{array}\right].$

5. the method for claim 1, is characterized in that, the expression formula of described RGB area image is:

IR＝RV×(1-GS)

IG＝RV×(1-f×GS)

IB＝RV×(1-(1-f)×GS) 。

6. the colors countenance of attribute more than a geologic body pick-up unit, it is characterized in that, this device comprises:

Post-stack seismic data acquiring unit, for obtaining general format post-stack seismic data;

Time horizon bit data acquiring unit, for utilizing general format post-stack seismic data to carry out geologic horizon explanation, obtains time horizon bit data;

Seismic attributes data acquiring unit, for utilizing general format post-stack seismic data to carry out seismic properties calculating, obtains seismic attributes data;

Along layer property value acquiring unit, for according to described time horizon bit data, extract from described seismic attributes data along layer property value;

Converter unit, for being normalized conversion to described along layer property value, obtains IR, IG, IB;

HSV information acquisition unit, obtains RH, RS, RV for utilizing IR, IG, IB;

Edge extracting unit, for carrying out edge extracting to RS, obtains object information GS;

Inverse transformation unit, the image for being formed RH, object information GS, RV carries out contravariant and changes to RGB territory, obtains RGB area image, completes geologic body and detects.

7. device as claimed in claim 6, is characterized in that, the expression formula of IR, IG, IB that described converter unit obtains is:

[IR(z),IG(z),IB(z)]＝T[V1(z),V2(z),V3(z0]

T＝255×(vi–vmin)/(vmax–vmin)

Wherein, z represents the planimetric position of a certain property value, and T conversion adopts linear transformation, and vi comprises V1, and V2, V3, V1, V2, V3 are all expressed as data component; Vmin and vmax is respectively V1, the minimum value of data component and maximal value in V2, V3.

8. device as claimed in claim 6, is characterized in that, the expression formula of RH, RS, RV that described HSV information acquisition unit obtains is:

RH＝60×(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. device as claimed in claim 6, is characterized in that, the expression formula of the object information GS that described edge extracting unit obtains is:

GS＝RS*Sobel

Wherein, $\mathrm{Sobel}\left[\begin{array}{ccc}-1& 0& 1\\ -2& 0& 2\\ -1& 0& 1\end{array}\right].$

10. device as claimed in claim 6, it is characterized in that, the expression formula of the RGB area image that described inverse transformation unit obtains is:

IR＝RV×(1-RS)

IG＝RV×(1-f×RS)

IB＝RV×(1-(1-f)×RS) 。