CN114255174A - Three-dimensional fault information identification method and device - Google Patents
Three-dimensional fault information identification method and device Download PDFInfo
- Publication number
- CN114255174A CN114255174A CN202011024395.0A CN202011024395A CN114255174A CN 114255174 A CN114255174 A CN 114255174A CN 202011024395 A CN202011024395 A CN 202011024395A CN 114255174 A CN114255174 A CN 114255174A
- Authority
- CN
- China
- Prior art keywords
- data
- initial
- difference
- dimensional volume
- dimensional
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 52
- 238000012545 processing Methods 0.000 claims abstract description 81
- 238000001914 filtration Methods 0.000 claims abstract description 40
- 238000010606 normalization Methods 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 14
- 238000004590 computer program Methods 0.000 claims description 13
- 230000009466 transformation Effects 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000004422 calculation algorithm Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 14
- 230000001427 coherent effect Effects 0.000 description 13
- 230000006870 function Effects 0.000 description 13
- 230000008569 process Effects 0.000 description 11
- 238000004891 communication Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 5
- 238000011161 development Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000002708 enhancing effect Effects 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/40—Image enhancement or restoration using histogram techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
Abstract
The invention discloses a method and a device for identifying three-dimensional fault information, wherein the method comprises the following steps: acquiring initial three-dimensional volume data of a third generation coherence attribute of a three-dimensional seismic data volume of a predetermined area, wherein the initial three-dimensional volume data has fault information; performing median filtering processing on the initial three-dimensional volume data, and determining difference data between the three-dimensional volume data subjected to the median filtering processing and the initial three-dimensional volume data; performing histogram equalization processing and threshold adjustment processing on the difference data to enhance fault information; and identifying fault information in the initial three-dimensional volume data according to the data after fault information enhancement. By the method and the device, the seismic data quality can be improved, the difference between the fault area and the background value is enhanced, noise pollution is eliminated, and therefore the fault identification precision can be improved.
Description
Technical Field
The invention relates to the field of seismic data processing, in particular to a method and a device for identifying three-dimensional fault information.
Background
Faults, fractures, karsts and the like are important components of fracture type oil and gas reservoirs such as carbonate heterogeneous reservoirs, and play an important role in the oil and gas reservoir forming process. The fault development condition is one of important factors for controlling the generation, storage, covering, looping, transportation and maintenance of reservoir oil and gas, so that many large oil and gas reservoirs in the world are closely related to fault development. Meanwhile, the fault is an important oil and gas migration channel, and the fault can communicate with dispersed storage spaces such as a karst cave system to form a large-scale storage layer and can provide necessary permeability for the interior of each storage space. Therefore, the high-precision fault detection work is carried out to support subsequent exploration and development researches such as modeling and oil reservoir simulation, and the method has important significance.
Therefore, fault identification is an important ring of seismic interpretation work in oil and gas exploration. However, conventional post-stack fault detection approaches suffer from the effects of seismic data quality, which has been plagued for a long time in both of these respects. On one hand, in a stratum broken region, cloud-fog-shaped fuzzy information shields a fault due to low offset imaging quality, and the real position of fault development is not accurately described; on the other hand, the presence of formation dip causes cross talk in the lateral direction of the fault signature.
That is, fault information in seismic data cannot be effectively identified due to poor quality of the seismic data.
Disclosure of Invention
In view of the above, the present invention provides a method and an apparatus for identifying three-dimensional fault information, so as to solve at least one of the above-mentioned problems.
According to a first aspect of the present invention, there is provided a method of identifying three-dimensional tomographic information, the method comprising:
acquiring initial three-dimensional volume data of a third generation coherence attribute of a three-dimensional seismic data volume of a predetermined area, wherein the initial three-dimensional volume data has fault information;
performing median filtering processing on the initial three-dimensional volume data, and determining difference data between the three-dimensional volume data subjected to the median filtering processing and the initial three-dimensional volume data;
performing histogram equalization processing and threshold adjustment processing on the difference data to enhance the fault information;
and identifying fault information in the initial three-dimensional data according to the data after the fault information enhancement.
According to a second aspect of the present invention, there is provided an apparatus for identifying three-dimensional tomographic information, the apparatus comprising:
a data acquisition unit for acquiring initial three-dimensional volume data of a third generation coherence property of a three-dimensional seismic data volume of a predetermined region, the initial three-dimensional volume data having fault information;
a median filtering unit, configured to perform median filtering processing on the initial three-dimensional volume data;
a difference data determination unit for determining difference data between the three-dimensional volume data after median filtering processing and the initial three-dimensional volume data;
an equalization processing unit configured to perform histogram equalization processing and threshold adjustment processing on the difference data to enhance the tomographic information;
and the identification unit is used for identifying the fault information in the initial three-dimensional volume data according to the data enhanced by the fault information.
According to a third aspect of the present invention, there is provided an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the method when executing the program.
According to a fourth aspect of the invention, a computer-readable storage medium is provided, on which a computer program is stored which, when being executed by a processor, carries out the steps of the above-mentioned method.
According to the technical scheme, the acquired initial three-dimensional volume data is subjected to median filtering, difference data between the three-dimensional volume data subjected to the median filtering and the initial three-dimensional volume data is determined, then histogram equalization processing and threshold value adjustment processing are performed on the difference data, so that fault information is enhanced, fault information in the initial three-dimensional volume data can be identified according to the data enhanced by the fault information, and the fault identification precision can be improved by improving the seismic data quality, enhancing the difference between a fault area and a background value, and eliminating noise pollution.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
Fig. 1 is a flowchart of a three-dimensional tomographic information identification method according to an embodiment of the present invention;
FIG. 2 is a cross-sectional view of the original coherence properties according to an embodiment of the present invention;
FIG. 3 is a graph of enhancement results for an original coherence property profile according to an embodiment of the present invention;
FIG. 4 is a raw coherence properties slice diagram according to an embodiment of the present invention;
FIG. 5 is a graph of enhancement results for slicing of original coherence properties according to an embodiment of the present invention;
fig. 6 is a block diagram of the structure of a three-dimensional tomographic information recognition apparatus according to an embodiment of the present invention;
fig. 7 is a schematic block diagram of a system configuration of an electronic apparatus 600 according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In carrying out the present invention, the applicant has found the following related art:
a typical representative of the tomographic technique is a coherent series technique. Coherent technology has originated in the last 90 th century, and three generations of algorithms have been developed. The first generation algorithm based on cross-correlation (abbreviated C1 algorithm) was proposed by Bahorich and Frmer in 1995, the second generation algorithm using multi-channel similarity (abbreviated C2 algorithm) was proposed by marcfort equal to 1998, and the third generation coherent algorithm based on feature structure (abbreviated C3 algorithm) was proposed by Gersztenkorn and marcfort. The C1 algorithm has harsh application precondition, the C2 algorithm has the defects of being sensitive to waveform and insensitive to the change of transverse amplitude, and in comparison, the C3 algorithm makes up the defects of the two algorithms, and obtains the coherence between each sample point and surrounding data in the shifted three-dimensional data body to form a three-dimensional data body representing the coherence, namely the data coherence in a calculation time window. Therefore, continuity can be suppressed, discontinuity can be highlighted, lateral change of seismic features can be quantitatively reflected, and the obtained result is more visual than geological interpretation of a seismic horizontal slice. The method is mainly applied to more objective and more detailed fault interpretation and river channel, sand body and fracture prediction. However, overall, the noise immunity of the coherent attribute technology is not good, and the coherent attribute technology is easily polluted by noise and is also easily subjected to crosstalk caused by an inclined stratum, so that the contrast between a fault and the surrounding background is often too low, which causes the result image to describe the fault less clearly, and even influences the judgment of an interpreter on fault interpretation. For example, in a fault development area formed by multiple stages of construction motion, improper setting of offset parameters causes energy dispersion of fault boundaries, causes fault information to present cloud distribution, and the like.
In view of the above problems, embodiments of the present invention provide a three-dimensional fault information identification scheme, which is based on third-generation coherent attribute three-dimensional volume data, and can improve fault identification accuracy of a research region by enhancing a difference between a fault region and a background value, eliminating noise pollution, and so on. Embodiments of the present invention are described in detail below with reference to the accompanying drawings.
Fig. 1 is a flowchart of a three-dimensional tomographic information recognition method according to an embodiment of the present invention, as shown in fig. 1, the flowchart including:
And 102, performing three-dimensional median filtering processing on the initial three-dimensional volume data, and determining difference data between the three-dimensional volume data subjected to median filtering processing and the initial three-dimensional volume data.
For example, for a point on the three-dimensional volume whose raw data is 0.56, the value of the point after median filtering is 0.64, and the difference value is 0.08 (0.64-0.56). The difference calculation method is to directly calculate the difference between two three-dimensional volume data point to point.
In an embodiment, initial difference data may be determined according to the three-dimensional volume data after the median filtering processing and the initial three-dimensional volume data, and then the initial difference data may be normalized to obtain difference data in a gray domain.
Specifically, the normalizing the initial difference data includes: determining a difference maximum value and a difference minimum value according to the initial difference data; and carrying out normalization processing on the initial difference data according to the difference maximum value and the difference minimum value.
In one example, the data volume is u, the range of values is [0,1], and the median filtering result of u is calculated to obtain um; calculating a difference du between the data volume u and the median filtering result um, wherein du is um-u; then, the transformation formula is provided for du normalized between [0,255 ]. For example, if the original data at a certain point is 0.56, the point value after median filtering is 0.64, and the difference value is 0.08(0.64-0.56), the point value becomes 23 after normalization.
Specifically, firstly, carrying out matrix transformation operation on the difference data of the gray domain, and carrying out histogram equalization processing after the matrix transformation operation; then, the processed difference data (U) is equalized to the histogram according to a predetermined ruleenhanced) The threshold value adjustment processing is performed for each pixel in (1).
In one embodiment, the matrix transformation algorithm is a global algorithm, and the histogram equalization method is not dimensional, and may be one-dimensional, two-dimensional, or three-dimensional. To achieve a fast application histogram algorithm, the difference data volume (three-dimensional, say a 5 x 6 x 7 "block") is here directly converted into a one-dimensional vector (5 x 6 x 7 in length).
In the case of wells, if there are wells in the work area, such as the inner A-well, the fractures of the target zone of the A-well develop, provided that the U is at the A-wellenhancedThe pixel value is presented as 159, then 159 is set here as the threshold, above which value is set to 255, and below which value the original value remains.
Under the condition of no well, according to geological experience, a large fault has a certain influence range, U, in a transverse rangeenhancedWill become progressively larger in the direction laterally away from the fault, at which time the U at the appropriate lateral distance is selectedenhancedThe value is used as a threshold value, and the following method refers to the case of a well.
The predetermined rule here may be that the threshold value of each pixel is adjusted according to the difference between the tomographic section and the background information. The "background" here is a value of the coherence property when the formation is flat and no fracture occurs, and is generally close to 1.
And 104, identifying fault information in the initial three-dimensional volume data according to the data after the fault information enhancement.
Specifically, fault information in the initial three-dimensional volume data is identified according to the time slice data and the section data after the fault information is enhanced.
The method comprises the steps of carrying out median filtering processing on the obtained initial three-dimensional volume data, determining difference data between the three-dimensional volume data subjected to the median filtering processing and the initial three-dimensional volume data, and then carrying out histogram equalization processing and threshold value adjustment processing on the difference data to enhance fault information, so that fault information in the initial three-dimensional volume data can be identified according to the data enhanced by the fault information.
For better understanding of the present invention, a specific flow example of three-dimensional fault information identification of the embodiment of the present invention is given below:
1) inputting data: the input data is three-dimensional volume data U of a third generation coherence attribute extracted from a three-dimensional seismic data volume, and the value range floating range of each pixel in the U is [0,1 ].
2) Preprocessing three-dimensional volume data of coherent attributes:
a) performing median filtering on the attribute body U to obtain a result U after the median filteringmAssuming that the dimension of Um is nmime, Nxline, and Ninline, where nmime (z axis, time axis), Nxline (x axis, xline direction), and Ninline (y axis, inline direction) are the number of time sampling points, xline lines, and inline lines, respectively, and the size of the filter window is set according to the size of the fault scale, for example, 3, 5, and the like;
b) calculating attribute U and value filtering result UmThe three-dimensional data volume difference value Δ U of (1), wherein Δ U is U-Um；
c) Normalizing the value range of delta U to [0,255] of the gray scale range]Within the range, the value range can be adjusted according to the value range distribution condition of U, and the calculation formula is shown as (1), wherein round is the rounding operator, and delta UminIs the minimum value of Δ U, Δ UmaxIs the maximum value of Δ U
ΔUnormal＝round((ΔU-ΔUmin)/(ΔUmax-ΔUmin)*255) (1)
d) For delta UnormalA matrix transformation is performed to generate a vector of length ntitime Nxline Ninline.
The preprocessing method plays a key role in fault boundary feature retention and strengthening.
3) The generated vector is equalized by using a conventional global histogram, and the obtained result is Uenhanced。
The conventional global histogram equalization processing has the characteristic of global optimization, and the risk of local distortion is low.
4) Adjusting U according to the difference between fault and background informationenhancedObtaining a fault enhancement result.
By the enhancement processing, a clouded tomographic blurring region in the image can be eliminated, and a linear tomographic structure in the region can be highlighted.
5) Extracting the enhancement results U separatelyenhancedThe time slice data, the profile data, and the raw data are compared to identify fault information.
According to the embodiment of the invention, the three-dimensional volume data based on the third-generation coherence attribute result is subjected to enhancement processing by a targeted preprocessing method and means such as global histogram equalization and threshold setting, so that the linear structure of the fault is strengthened, and the fault identification precision is improved.
To further understand the embodiments of the present invention, an example is given below. In this example, the main fault distribution trend in the X-ray work area is known to be clear, but the boundary of the secondary fault is blurred due to formation interference.
According to the embodiment of the invention, the process for identifying the break layer information in the X work area specifically comprises the following steps:
1) extracting a third generation coherent attribute body of the X work area three-dimensional post-stack seismic data body;
2) the coherent attribute body is enhanced by using the technical solution provided by the embodiment of the present invention to obtain results, and as shown in the drawings, fig. 2 to fig. 5 can be seen, where fig. 2 is a cross-sectional view of the original coherent attribute, fig. 3 is a view of an enhancement result of the cross-sectional view of the original coherent attribute, fig. 4 is a view of a slice of the original coherent attribute, and fig. 5 is a view of an enhancement result of the slice of the original coherent attribute.
As can be seen from the comparison and analysis of the imaging results, after the enhancement method provided by the embodiment of the invention is applied, the fault identification precision of the research region is integrally improved, the enhancement result can clearly show the linear structure of the fault, and the embodiment of the invention improves the fault identification precision.
Based on similar inventive concepts, the embodiment of the present invention further provides an apparatus for identifying three-dimensional fault information, and preferably, the apparatus may be used to implement the process in the above method embodiment.
Fig. 6 is a block diagram showing the structure of the three-dimensional tomographic information recognition apparatus, which includes, as shown in fig. 6: a data acquisition unit 61, a median filtering unit 62, a difference data determination unit 63, an equalization processing unit 64, and a recognition unit 65, wherein:
a data acquisition unit 61 for acquiring initial three-dimensional volume data of a third generation coherence property of a three-dimensional seismic data volume of a predetermined area, the initial three-dimensional volume data having tomographic information;
a median filtering unit 62, configured to perform median filtering processing on the initial three-dimensional volume data;
a difference data determining unit 63 configured to determine difference data between the three-dimensional volume data after the median filtering processing and the initial three-dimensional volume data;
an equalization processing unit 64 for performing histogram equalization processing and threshold adjustment processing on the difference data to enhance the tomographic information;
and the identification unit 65 is used for identifying the fault information in the initial three-dimensional volume data according to the data after the fault information enhancement.
Specifically, the identifying unit 65 identifies the tomographic information in the initial three-dimensional volume data from the slice data and the section data after enhancement of the tomographic information.
The initial three-dimensional volume data acquired by the data acquisition unit 61 is subjected to median filtering processing by the median filtering unit 62, the difference data determination unit 63 determines difference data between the three-dimensional volume data subjected to the median filtering processing and the initial three-dimensional volume data, and then the equalization processing unit 64 performs histogram equalization processing and threshold adjustment processing on the difference data to enhance fault information, so that the identification unit 65 can identify fault information in the initial three-dimensional volume data according to the data subjected to fault information enhancement.
Specifically, the above-described difference data determination unit 63 includes: an initial difference determination module and a normalization module, wherein:
the initial difference determining module is used for determining initial difference data according to the three-dimensional volume data after median filtering processing and the initial three-dimensional volume data;
and the normalization module is used for performing normalization processing on the initial difference data to obtain difference data of a gray scale domain.
In one embodiment, the normalization module specifically includes: a difference value determining submodule for determining a difference maximum value and a difference minimum value according to the initial difference data; and the normalization submodule is used for performing normalization processing on the initial difference data according to the difference maximum value and the difference minimum value.
In a specific implementation process, the equalization processing unit 64 specifically includes: a histogram equalization module and a threshold adjustment module, wherein:
the histogram equalization module is used for carrying out matrix transformation operation on the difference data of the gray scale domain and carrying out histogram equalization processing after the matrix transformation operation;
and the threshold adjusting module is used for carrying out threshold adjusting processing on each pixel in the difference data after the histogram equalization processing according to a preset rule.
For specific execution processes of the units, the modules, and the sub-modules, reference may be made to the description in the foregoing method embodiments, and details are not described here again.
In practical operation, the units, the modules and the sub-modules may be combined or may be arranged singly, and the present invention is not limited thereto.
The present embodiment also provides an electronic device, which may be a desktop computer, a tablet computer, a mobile terminal, and the like, but is not limited thereto. In this embodiment, the electronic device may be implemented by referring to the above method embodiment and the embodiment of the three-dimensional fault information identification device, and the contents thereof are incorporated herein, and repeated descriptions are omitted.
Fig. 7 is a schematic block diagram of a system configuration of an electronic apparatus 600 according to an embodiment of the present invention. As shown in fig. 7, the electronic device 600 may include a central processor 100 and a memory 140; the memory 140 is coupled to the central processor 100. Notably, this diagram is exemplary; other types of structures may also be used in addition to or in place of the structure to implement telecommunications or other functions.
In one embodiment, the three-dimensional fault information recognition function may be integrated into the central processor 100. The central processor 100 may be configured to control as follows:
acquiring initial three-dimensional volume data of a third generation coherence attribute of a three-dimensional seismic data volume of a predetermined area, wherein the initial three-dimensional volume data has fault information;
performing median filtering processing on the initial three-dimensional volume data, and determining difference data between the three-dimensional volume data subjected to the median filtering processing and the initial three-dimensional volume data;
performing histogram equalization processing and threshold adjustment processing on the difference data to enhance the fault information;
and identifying fault information in the initial three-dimensional data according to the data after the fault information enhancement.
As can be seen from the above description, in the electronic device provided in the embodiment of the present application, the median filtering process is performed on the acquired initial three-dimensional volume data, the difference data between the three-dimensional volume data after the median filtering process and the initial three-dimensional volume data is determined, and then, the histogram equalization process and the threshold adjustment process are performed on the difference data to enhance the fault information, so that the fault information in the initial three-dimensional volume data can be identified according to the data enhanced by the fault information.
In another embodiment, the three-dimensional tomographic information recognition apparatus may be configured separately from the central processor 100, for example, the three-dimensional tomographic information recognition apparatus may be configured as a chip connected to the central processor 100, and the three-dimensional tomographic information recognition function may be realized by the control of the central processor.
As shown in fig. 7, the electronic device 600 may further include: communication module 110, input unit 120, audio processing unit 130, display 160, power supply 170. It is noted that the electronic device 600 does not necessarily include all of the components shown in fig. 7; furthermore, the electronic device 600 may also comprise components not shown in fig. 7, which may be referred to in the prior art.
As shown in fig. 7, the central processor 100, sometimes referred to as a controller or operational control, may include a microprocessor or other processor device and/or logic device, the central processor 100 receiving input and controlling the operation of the various components of the electronic device 600.
The memory 140 may be, for example, one or more of a buffer, a flash memory, a hard drive, a removable media, a volatile memory, a non-volatile memory, or other suitable device. The information relating to the failure may be stored, and a program for executing the information may be stored. And the central processing unit 100 may execute the program stored in the memory 140 to realize information storage or processing, etc.
The input unit 120 provides input to the cpu 100. The input unit 120 is, for example, a key or a touch input device. The power supply 170 is used to provide power to the electronic device 600. The display 160 is used to display an object to be displayed, such as an image or a character. The display may be, for example, an LCD display, but is not limited thereto.
The memory 140 may be a solid state memory such as Read Only Memory (ROM), Random Access Memory (RAM), a SIM card, or the like. There may also be a memory that holds information even when power is off, can be selectively erased, and is provided with more data, an example of which is sometimes called an EPROM or the like. The memory 140 may also be some other type of device. Memory 140 includes buffer memory 141 (sometimes referred to as a buffer). The memory 140 may include an application/function storage section 142, and the application/function storage section 142 is used to store application programs and function programs or a flow for executing the operation of the electronic device 600 by the central processing unit 100.
The memory 140 may also include a data store 143, the data store 143 for storing data, such as contacts, digital data, pictures, sounds, and/or any other data used by the electronic device. The driver storage portion 144 of the memory 140 may include various drivers of the electronic device for communication functions and/or for performing other functions of the electronic device (e.g., messaging application, address book application, etc.).
The communication module 110 is a transmitter/receiver 110 that transmits and receives signals via an antenna 111. The communication module (transmitter/receiver) 110 is coupled to the central processor 100 to provide an input signal and receive an output signal, which may be the same as in the case of a conventional mobile communication terminal.
Based on different communication technologies, a plurality of communication modules 110, such as a cellular network module, a bluetooth module, and/or a wireless local area network module, may be provided in the same electronic device. The communication module (transmitter/receiver) 110 is also coupled to a speaker 131 and a microphone 132 via an audio processor 130 to provide audio output via the speaker 131 and receive audio input from the microphone 132 to implement general telecommunications functions. Audio processor 130 may include any suitable buffers, decoders, amplifiers and so forth. In addition, an audio processor 130 is also coupled to the central processor 100, so that recording on the local can be enabled through a microphone 132, and so that sound stored on the local can be played through a speaker 131.
Embodiments of the present invention also provide a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to implement the steps of the three-dimensional fault information identification method.
In summary, the embodiment of the present invention performs enhancement processing by using three-dimensional volume data of the third-generation coherence attribute result and combining with conventional global histogram equalization, threshold setting, and other means through a targeted preprocessing method, so as to strengthen the linear structure of the fault, and thus the method is an effective fault identification technology. The advantages of the embodiment of the invention include the following aspects:
(1) the enhancement method can eliminate a cloudy fault fuzzy region in the image and highlight a linear fault structure in the region; (2) the calculation process is simple, easy to realize and high in calculation efficiency; (3) the capability of mining partial hidden fault information is provided; (4) the conventional global histogram equalization processing method has the characteristic of global optimization, and the risk of local distortion is low.
The preferred embodiments of the present invention have been described above with reference to the accompanying drawings. The many features and advantages of the embodiments are apparent from the detailed specification, and thus, it is intended by the appended claims to cover all such features and advantages of the embodiments which fall within the true spirit and scope thereof. Further, since numerous modifications and changes will readily occur to those skilled in the art, it is not desired to limit the embodiments of the invention to the exact construction and operation illustrated and described, and accordingly, all suitable modifications and equivalents may be resorted to, falling within the scope thereof.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The principle and the implementation mode of the invention are explained by applying specific embodiments in the invention, and the description of the embodiments is only used for helping to understand the method and the core idea of the invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, there may be variations in the specific embodiments and the application scope, and in summary, the content of the present specification should not be construed as a limitation to the present invention.
Claims (10)
1. A method for identifying three-dimensional fault information, the method comprising:
acquiring initial three-dimensional volume data of a third generation coherence attribute of a three-dimensional seismic data volume of a predetermined area, wherein the initial three-dimensional volume data has fault information;
performing median filtering processing on the initial three-dimensional volume data, and determining difference data between the three-dimensional volume data subjected to the median filtering processing and the initial three-dimensional volume data;
performing histogram equalization processing and threshold adjustment processing on the difference data to enhance the fault information;
and identifying fault information in the initial three-dimensional data according to the data after the fault information enhancement.
2. The method of claim 1, wherein determining difference data between the median filtered three-dimensional volume data and the initial three-dimensional volume data comprises:
determining initial difference data according to the three-dimensional volume data after median filtering and the initial three-dimensional volume data;
and carrying out normalization processing on the initial difference data to obtain difference data of a gray scale domain.
3. The method of claim 2, wherein normalizing the initial difference data comprises:
determining a difference maximum value and a difference minimum value according to the initial difference data;
and carrying out normalization processing on the initial difference data according to the difference maximum value and the difference minimum value.
4. The method of claim 2, wherein performing histogram equalization processing and threshold adjustment processing on the difference data comprises:
performing matrix transformation operation on the difference data of the gray scale domain, and performing histogram equalization processing after the matrix transformation operation;
and performing threshold value adjustment processing on each pixel in the difference data after the histogram equalization processing according to a preset rule.
5. The method of claim 1, wherein identifying fault information in the initial three-dimensional volume data from the fault information enhanced data comprises:
and identifying fault information in the initial three-dimensional volume data according to the time slice data and the profile data after the fault information enhancement.
6. An apparatus for recognizing three-dimensional tomographic information, the apparatus comprising:
a data acquisition unit for acquiring initial three-dimensional volume data of a third generation coherence property of a three-dimensional seismic data volume of a predetermined region, the initial three-dimensional volume data having fault information;
a median filtering unit, configured to perform median filtering processing on the initial three-dimensional volume data;
a difference data determination unit for determining difference data between the three-dimensional volume data after median filtering processing and the initial three-dimensional volume data;
an equalization processing unit configured to perform histogram equalization processing and threshold adjustment processing on the difference data to enhance the tomographic information;
and the identification unit is used for identifying the fault information in the initial three-dimensional volume data according to the data enhanced by the fault information.
7. The apparatus of claim 6, wherein the difference data determining unit comprises:
the initial difference determining module is used for determining initial difference data according to the three-dimensional volume data after median filtering processing and the initial three-dimensional volume data;
and the normalization module is used for performing normalization processing on the initial difference data to obtain difference data of a gray scale domain.
8. The apparatus of claim 7, wherein the normalization module comprises:
a difference value determining submodule for determining a difference maximum value and a difference minimum value according to the initial difference data;
and the normalization submodule is used for performing normalization processing on the initial difference data according to the difference maximum value and the difference minimum value.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the steps of the method according to any of claims 1 to 5 are implemented when the processor executes the program.
10. A computer-readable storage medium, on which a computer program is stored which, when being executed by a processor, carries out the steps of the method of any one of claims 1 to 5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011024395.0A CN114255174A (en) | 2020-09-25 | 2020-09-25 | Three-dimensional fault information identification method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011024395.0A CN114255174A (en) | 2020-09-25 | 2020-09-25 | Three-dimensional fault information identification method and device |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114255174A true CN114255174A (en) | 2022-03-29 |
Family
ID=80789204
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011024395.0A Pending CN114255174A (en) | 2020-09-25 | 2020-09-25 | Three-dimensional fault information identification method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114255174A (en) |
-
2020
- 2020-09-25 CN CN202011024395.0A patent/CN114255174A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107688201B (en) | RBM-based seismic prestack signal clustering method | |
| CN103792576A (en) | Reservoir heterogeneous detection method and device based on gradient structure tensor | |
| CN112381845B (en) | Rock core image generation method, model training method and device | |
| CN108847252B (en) | Acoustic feature extraction method based on acoustic signal spectrogram texture distribution | |
| CN108830140B (en) | Volcanic lithology identification method based on electric imaging logging fractal dimension | |
| CN112415591A (en) | Diffracted wave imaging method and device, electronic equipment and storage medium | |
| CN111582114B (en) | Earthquake fault identification method, device, equipment and storage medium | |
| CN112418345A (en) | Method and device for quickly identifying fine-grained small target | |
| CN115184998A (en) | Rayleigh wave frequency dispersion curve automatic extraction method based on improved U-net neural network | |
| CN114255174A (en) | Three-dimensional fault information identification method and device | |
| CN109948291A (en) | A kind of discontinuous boundary line direction-adaptive recognition methods of sand body | |
| CN114252913A (en) | Method and device for identifying plane fault information | |
| CN114114421B (en) | Deep learning-based guided self-learning seismic data denoising method and device | |
| CN113743193B (en) | Pre-stack seismic data linear interference suppression method and system | |
| WO2022250801A1 (en) | Seismic first breaks onset times determination by interface tracking using level-sets method | |
| CN105301655B (en) | A kind of common imaging gather line noise minimizing technology and device | |
| CN112147687B (en) | Reservoir gas content prediction method and prediction system | |
| CN114755740A (en) | Rock distribution determination method, device, equipment and medium | |
| CN114139229A (en) | Fault enhancement method, fault development interpretation method, storage medium, and electronic apparatus | |
| CN113075747B (en) | Reservoir fracture development area prediction method and device | |
| CN113900141A (en) | Oil gas distribution prediction method and device | |
| CN111290020A (en) | Fracture detection method and device based on structural filtering processing and frequency division attribute fusion | |
| CN113325474B (en) | Method for discriminating biological reef | |
| CN113960657B (en) | Seismic data characteristic acquisition method and system | |
| CN110873901B (en) | Pseudo well curve frequency increasing method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |