CN119024438B - Logging remote detection imaging domain coherent imaging method, device and equipment - Google Patents
Logging remote detection imaging domain coherent imaging method, device and equipment Download PDFInfo
- Publication number
- CN119024438B CN119024438B CN202411513967.XA CN202411513967A CN119024438B CN 119024438 B CN119024438 B CN 119024438B CN 202411513967 A CN202411513967 A CN 202411513967A CN 119024438 B CN119024438 B CN 119024438B
- Authority
- CN
- China
- Prior art keywords
- imaging
- data
- offset
- bootstrap
- coherence
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 198
- 238000001514 detection method Methods 0.000 title claims abstract description 82
- 230000001427 coherent effect Effects 0.000 title claims abstract description 46
- 238000000034 method Methods 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims description 9
- 238000010586 diagram Methods 0.000 description 7
- 238000005070 sampling Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 230000006870 function Effects 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
- G01V1/305—Travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供一种测井远探测成像域相干成像方法、装置及设备，涉及油气勘探与开发技术领域，方法包括：获取测井远探测分偏移距数据；对测井远探测分偏移距数据进行分偏移距成像，得到分偏移距成像结果；对分偏移距成像结果沿着偏移距方向水平叠加，得到初始叠加剖面；对初始叠加剖面进行深度方向自举相干增强，得到初始叠加剖面的每个位置点对应的初步增强成像数据；针对每个初步增强成像数据对应的位置点，计算该位置点处的第二高分辨率自举相干；第二高分辨率自举相干为分偏移距成像结果在偏移距方向的第二高分辨率自举相干；根据第二高分辨率自举相干和初步增强成像数据，生成相干成像结果。本发明能够提高测井远探测数据井周裂缝成像精度。
The present invention provides a method, device and equipment for coherent imaging in a logging remote detection imaging domain, and relates to the technical field of oil and gas exploration and development. The method comprises: obtaining remote detection offset data of logging; performing offset imaging on the remote detection offset data of logging to obtain offset imaging results; horizontally stacking the offset imaging results along the offset direction to obtain an initial stacking section; performing bootstrap coherent enhancement in the depth direction on the initial stacking section to obtain preliminary enhanced imaging data corresponding to each position point of the initial stacking section; for each position point corresponding to the preliminary enhanced imaging data, calculating the second high-resolution bootstrap coherence at the position point; the second high-resolution bootstrap coherence is the second high-resolution bootstrap coherence of the offset imaging result in the offset direction; generating a coherent imaging result according to the second high-resolution bootstrap coherence and the preliminary enhanced imaging data. The present invention can improve the imaging accuracy of wellbore fractures in remote detection data of logging.
Description
Technical Field
The invention relates to the technical field of oil and gas exploration and development, in particular to a method, a device and equipment for coherent imaging of a logging remote detection imaging domain.
Background
The oil and gas resources are stored in underground rock cracks, and the detection of the cracks of stratum around the well plays an important role in oil and gas exploration and development. The logging acoustic remote detection technology realizes the detection of stratum conditions within the range of tens of meters around a well by exciting acoustic waves of thousands of hertz and receiving acoustic reflection signals in the stratum through a limited number of receivers. Because of its high detection frequency (e.g., 3000Hz dominant frequency), and its relatively close proximity to the reservoir location in the wellbore, it has a natural advantage in imaging small-scale fractures of the reservoir.
The logging remote detection essence is to accurately position the spatial position of the underground geologic body by utilizing the acoustic wave reflection signals of the underground geologic body, and the received signals need to be imaged. Because the well logging remote detection data and the seismic data have similar physical bases, the imaging of the well logging remote detection on the small-scale fracture mainly refers to the imaging methods of the seismic imaging, such as F-K offset imaging, kirchhoff integral imaging and the like. However, the seismic data is usually observed with extremely large aperture (such as an offset distance of >8 km), so the seismic imaging method usually assumes large aperture data observation and uses intra-aperture imaging domain data superposition to achieve the purpose of suppressing imaging artifacts. However, the well logging remote detection generally adopts a few receivers to receive the data, and the observation aperture is very limited, so the well logging remote detection data does not meet the related requirements, and the imaging accuracy is reduced, imaging artifacts exist, and the interpretability of fracture imaging results is affected.
Based on this, how to improve imaging accuracy of the logging remote detection data well Zhou Liefeng becomes a technical problem to be solved urgently.
Disclosure of Invention
In view of the above, the present invention provides a method, apparatus and device for coherent imaging in a logging remote detection imaging domain.
The invention adopts the following technical scheme:
in a first aspect, the present invention provides a method for coherent imaging in a well logging remote detection imaging domain, comprising:
acquiring well logging far detection offset data, wherein the well logging far detection offset data is three-dimensional data with three dimensions of depth, time and offset;
Performing sub-offset imaging on the well logging remote detection sub-offset data to obtain a sub-offset imaging result;
Horizontally superposing the sub-offset imaging result along the offset direction to obtain an initial superposition profile;
Performing depth-direction bootstrap coherent enhancement on the initial superposition profile to obtain preliminary enhancement imaging data corresponding to each position point of the initial superposition profile;
Calculating second high-resolution bootstrap coherence at a position point corresponding to each preliminary enhanced imaging data, wherein the second high-resolution bootstrap coherence is second high-resolution bootstrap coherence of the sub-offset imaging result in the offset direction;
Generating a final enhanced coherent imaging result from the second high resolution bootstrapping coherence and the preliminary enhanced imaging data.
Optionally, the acquiring the data of the well logging far detection sub offset distance specifically includes:
Acquiring logging far detection pre-stack data;
And recombining and arranging the logging remote detection pre-stack data to obtain the logging remote detection sub-offset data.
Optionally, performing depth-direction bootstrapping coherence enhancement on the initial superposition profile to obtain preliminary enhancement imaging data corresponding to each position point of the initial superposition profile, which specifically includes:
Determining each scanning slope based on a preset maximum scanning slope, a preset minimum scanning slope and a preset scanning slope number in the depth direction;
For each position point of the initial superposition profile, respectively acquiring first partial window data of the position point relative to a first partial window along the direction of each scanning slope, wherein the first partial window is a window with depth direction length and radial distance direction length;
Calculating a first high resolution bootstrap coherence at each of the scan slopes using the first local window data at the respective scan slope;
and extracting the preliminary enhancement effective signals along the direction of the scanning slope corresponding to the maximum first high-resolution bootstrap coherence to obtain the preliminary enhancement imaging data corresponding to the position points.
Optionally, the calculating the second high-resolution bootstrap coherence at the location point specifically includes:
Based on a second local window, obtaining second local window data of the partial offset imaging result in the offset direction at the position point, wherein the second local window is a window with a radial distance direction length and an offset direction length;
using the second local window data, a second high resolution bootstrap coherence at the location point is calculated.
Optionally, generating a final enhanced coherent imaging result according to the second high-resolution bootstrap coherence and the preliminary enhanced imaging data, specifically including:
normalizing all the second high-resolution bootstrap coherence;
converting the normalized second high-resolution bootstrap coherence into a coherence coefficient;
generating the final enhanced coherent imaging result using the coherence coefficient and the preliminary enhanced imaging data.
In a second aspect, the present invention provides a well logging remote detection imaging domain coherent imaging apparatus comprising:
the acquisition module is used for acquiring the well logging far detection partial offset data, wherein the well logging far detection partial offset data is three-dimensional data with three dimensions of depth, time and offset;
The sub offset imaging module is used for performing sub offset imaging on the well logging far detection sub offset data to obtain a sub offset imaging result;
the superposition module is used for horizontally superposing the sub-offset imaging result along the offset direction to obtain an initial superposition section;
The initial enhancement module is used for carrying out depth-direction bootstrap coherent enhancement on the initial superposition profile to obtain initial enhancement imaging data corresponding to each position point of the initial superposition profile;
The calculation module is used for calculating second high-resolution bootstrap coherence at the position point corresponding to each initial enhancement imaging data, wherein the second high-resolution bootstrap coherence is second high-resolution bootstrap coherence of the sub-offset imaging result in the offset direction;
And the generation module is used for generating a final enhanced coherent imaging result according to the second high-resolution bootstrap coherence and the preliminary enhanced imaging data.
In a third aspect, the present invention provides a computer device comprising:
at least one processor, and
A memory communicatively coupled to the at least one processor, wherein,
The memory stores instructions executable by the at least one processor to enable the at least one processor to implement a well logging far detection imaging domain coherence imaging method as described above.
The method comprises the steps of obtaining well logging far detection offset data, wherein the well logging far detection offset data are three-dimensional data with depth, time and offset, performing offset imaging on the well logging far detection offset data to obtain offset imaging results, horizontally superposing the offset imaging results along the offset direction to obtain an initial superposition section, performing depth-direction bootstrap coherence enhancement on the initial superposition section to obtain preliminary enhancement imaging data corresponding to each position point of the initial superposition section, calculating second high-resolution bootstrap coherence at the position point corresponding to each position point of the preliminary enhancement imaging data, wherein the second high-resolution bootstrap coherence is second high-resolution bootstrap coherence of the offset imaging results in the offset direction, and generating final enhancement coherent imaging results according to the second high-resolution bootstrap coherence and the preliminary enhancement imaging data. Based on the method, the problem of imaging artifacts generated by imaging the well-logging far detection data with limited aperture by a conventional imaging method is solved by introducing a high-precision coherence measure method, the imaging precision of the well-logging far detection data well Zhou Liefeng is effectively improved, and technical support is provided for high-precision crack imaging in oil and gas exploration and development.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic flow chart of a method for coherent imaging in a well logging remote detection imaging domain according to an embodiment of the present invention;
FIG. 2 is a diagram of well log far survey pre-stack data acquired by an embodiment of the present invention;
FIG. 3 is a conventional well logging far detection imaging result;
FIG. 4 is a plot of the result of far-detection imaging field coherent imaging using the far-detection imaging field coherent imaging method of the present invention corresponding to FIG. 3;
FIG. 5 is a schematic structural diagram of a coherent imaging device of a logging remote detection imaging domain according to an embodiment of the present invention;
fig. 6 is a schematic structural diagram of a computer device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail below. It will be apparent that the described embodiments are only some, but not all, embodiments of the invention. All other embodiments, based on the examples herein, which are within the scope of the invention as defined by the claims, will be within the scope of the invention as defined by the claims.
Examples
Fig. 1 is a schematic flow chart of a method for coherent imaging in a logging remote detection imaging domain according to an embodiment of the present invention. As shown in fig. 1, the present process includes:
step 101, acquiring well logging far detection sub offset data. The well logging far detection sub offset data is three-dimensional data with three dimensions of depth, time and offset.
Specifically, step 101 of acquiring logging distance detection offset data may include:
Firstly, acquiring logging remote detection pre-stack data, and then, recombining and arranging the logging remote detection pre-stack data to obtain logging remote detection sub-offset data . Wherein, In order to be of depth,In order to be able to take time,Is the offset distance.
In a specific example, FIG. 2 is a diagram of well log far survey pre-stack data acquired by an embodiment of the present invention. As shown in fig. 2, the direction of the coordinates Nr represents the number of receivers, which is the number of offset distances, and is the number of sampling points in the direction of offset distances, depth is Depth, and Time is Time.
In addition, in step 101, the well logging remote detection offset data can be obtained, or the well logging remote detection offset data stored in advance or input by a user can be directly obtained.
Step 102, carrying out partial offset imaging on the well logging remote detection partial offset data to obtain a partial offset imaging result.
Specifically, using formula (1) to perform offset imaging on the well logging far detection offset data to obtain an offset imaging resultEquation (1) is as follows:
......(1)
wherein, In order to be a radial distance from each other,For the imaging weighting coefficients to be used,Is offset fromIs a double-pass travel time of the computer program.
Partial offset imaging resultsCan be expressed in discrete mathematical form as. Wherein, Expressed in depthThe number of sampling points in the direction isAt radial distance ofThe number of sampling points in the direction isAt an offset distanceThe number of sampling points in the direction is. And, the starting sample point numbers for all directions are from 1,Reference numerals representing sampling pointsA result of the partial offset imaging at which, among other things,,,。
In a specific example, the values of the relevant parameters may be:,, Depth sampling interval Radial distance sampling interval。
And 103, horizontally superposing the bisection offset imaging result along the offset direction to obtain an initial superposition section.
Specifically, the imaging results of the bisection offset are horizontally overlapped along the offset direction by using the formula (2) to obtain an initial overlapped sectionEquation (2) is as follows:
......(2)
Step 104, performing depth-direction Bootstrap (Bootstrap) coherent enhancement on the initial superposition profile to obtain preliminary enhancement imaging data corresponding to each position point of the initial superposition profile.
In the embodiment of the present invention, step 104, performing depth-direction bootstrapping coherence enhancement on the initial superposition profile to obtain preliminary enhancement imaging data corresponding to each position point of the initial superposition profile, which may specifically include:
(1041) Each scanning slope is determined based on a preset maximum scanning slope, a preset minimum scanning slope and a preset number of scanning slopes in the depth direction.
Specifically, the invention sets the maximum scanning slope in the depth directionMinimum scan slopeNumber of scan slopesSo that the in-phase axial slope of the effective signal in the initial superposition section is included in the intervalIn (3), a series of scan slopes at equal intervals are calculated by using the formula (3)Equation (3) is as follows:
......(3)
wherein, Is marked with the reference numberIs set up to the scanning slope of (c),。
In a specific example, the related parameters can be the maximum scanning slopeMinimum scan slopeNumber of scan slopes。
(1042) And acquiring first local window data of the position points relative to a first local window along the direction of each scanning slope for each position point of the initial superposition profile, wherein the first local window is a window with a depth direction length and a radial distance direction length.
In particular, the first partial window may be designed to be a depth-wise long windowLong in the radial distance directionThe window size of the partial window, i.e. the first partial window, isWherein, the method comprises the steps of, wherein,AndAre all non-negative odd numbers.
It should be noted that, the acquisition process of the preliminary enhanced imaging data corresponding to each position point is the same or similar, so that the position points are usedFor example, a procedure of acquiring preliminary enhanced imaging data corresponding to one position point is described. And in the computer array index, the distance between adjacent position points in the depth direction is an index unit 1 under the same radial distance, and the distance between adjacent position points in the radial distance direction is also an index unit 1 under the same depth.
At the point of acquisitionIn the corresponding process of preliminarily enhancing imaging data, after a first local window is obtained, enabling a window center point and a position point of the first local windowOverlap along each scan slopeIs to acquire a position pointFirst partial window data about a first partial windowCan obtainFirst partial window dataFirst partial Window dataThe calculation formula of (2) is as follows:
......(4)
wherein, ,AndThe windows of the first partial window in the depth direction and the radial distance direction are half-long,,,,,,,,In order to take the whole downwards,Is rounded upward.
(1043) First high resolution bootstrap coherence at each scan slope is calculated using the first local window data at the corresponding scan slope, respectively.
Specifically, the first high-resolution bootstrap coherence under any scan slope is calculated by first, applying a first local window data under the scan slopeSubstituting the first similarity coefficient into the formula (5) to obtain the first similarity coefficientEquation (5) is as follows:
......(5)
then, for the first partial window data Performing non-repeated bootstrap for a first preset number of times to obtain first local window bootstrap data of the first preset number. Here, the first partial window dataIs common in depth directionThe data of each track is taken as a whole, and the data of the track and the track are rearranged randomly to obtain first local window bootstrap data without repeated bootstrapFirst partial Window dataIs shared in depth directionLanes, altogether, can be generatedThe local window bootstraps the data, wherein,,Representation ofIs the nth partial window bootstrapping data.
Finally, based on equation (6), bootstrap data according to the first partial windowAnd a first similarity coefficientCalculating to obtain a first high-resolution bootstrap coherence under the current scanning slopeEquation (6) is as follows:
......(6)
wherein, ,Indicating nth first partial window bootstrap dataIs used for the correlation coefficient of (a),Expressed in the scanning slopeDirection bootstrapping data with n first partial windowsThe first high resolution bootstrap coherence calculated.
In the present invention,The subscript n of (c) may be set according to the actual situation, and the present invention is not limited thereto, for example, n may be set to 3, i.e., the first preset number is equal to 3.
(1044) And extracting the preliminary enhancement effective signals along the direction of the scanning slope corresponding to the maximum first high-resolution bootstrap coherence to obtain preliminary enhancement imaging data corresponding to the position points.
Specifically, determining the scanning slope corresponding to the first high-resolution bootstrap coherence with the maximum valueThen, along the scan slopeExtracting the primary enhancement effective signal in the direction of (2) to obtain the position pointCorresponding preliminary enhanced imaging dataPreliminary enhancement of imaging dataThe calculation formula of (2) is as follows:
......(7)
wherein, To follow the scanning slopeIs a first partial window data acquired in the direction of (a), here the number of the elements is the number,。Is the gaussian weight coefficient of the model,In this embodiment, take。
Step 105, calculating second high-resolution bootstrap coherence at a position point corresponding to each preliminary enhanced imaging data, wherein the second high-resolution bootstrap coherence is second high-resolution bootstrap coherence of a sub-offset imaging result in an offset direction.
Here, calculating the second high resolution bootstrap coherence at the location point may specifically include:
(1) And obtaining second partial window data of the partial offset imaging result in the offset direction at the position point based on a second partial window, wherein the second partial window is a window with a radial distance direction length and an offset direction length.
In particular, the second partial window may be designed to be a window that is long in the radial distance directionIn the direction of offset distance is longThe window size of the partial window, i.e. the second partial window, is. At the position point corresponding to the preliminary enhanced imaging dataObtaining second partial window data of the partial offset imaging result in the offset directionThe formula used is as follows:
......(8)
wherein, ,Is an integer in the range of,。
(2) Second high resolution bootstrap coherence at the location point is calculated using the second local window data.
Specifically, first, the second partial window dataSubstituting formula (9) to obtain a second similarity coefficientEquation (9) is as follows:
......(9)
Then, for the second partial window data Performing non-repeated bootstrap for a second preset number of times to obtain bootstrap data of a second local window of the second preset number. Here, the second partial window dataIs common in the offset directionThe data of each track is taken as a whole, and the data of the track and the track are rearranged randomly to obtain second local window bootstrap data without repeated bootstrap. Because of the second partial window dataIs common in the offset directionLanes, so that a total ofSecond partial Window bootstrapping data that is not repeatedWherein, the method comprises the steps of, wherein,。Representation ofIs the r second partial window bootstrapping data.
Finally, based on equation (10), bootstrap data according to the second partial windowAnd a second coefficient of similarityCalculating to obtain position points corresponding to the preliminary enhanced imaging dataSecond high resolution bootstrap coherence atEquation (10) is as follows:
......(10)
wherein, ,Indicating the r second partial window bootstrap dataIs used for the correlation coefficient of (a),Is expressed at a position pointBootstrap data with r second partial windowsThe calculated second high resolution bootstrap coherence.
In an embodiment of the present invention,Subscript of (2)The value of (2) can be set according to the actual situation, the invention is not particularly limited, and for example, can be setThe value of (2) is 4, i.e. the second preset number is equal to 4.
Step 106, generating a final enhanced coherent imaging result according to the second high-resolution bootstrap coherent and preliminary enhanced imaging data.
In the embodiment of the present invention, step 106, generating a final enhanced coherent imaging result according to the second high-resolution bootstrap coherent and preliminary enhanced imaging data, may specifically include:
(1061) Normalization processing is performed on all second high-resolution bootstrap coherence.
Specifically, through the above process, a second high-resolution bootstrap coherence at the corresponding position of each preliminary enhanced imaging data point is obtained. Dividing each second high-resolution bootstrap coherence by the maximum value of all second high-resolution bootstrap coherence (namely the maximum second high-resolution bootstrap coherence) to obtain normalized high-resolution bootstrap coherence with the element value range of the second high-resolution bootstrap coherence between 0 and 1。
(1062) And converting the normalized second high-resolution bootstrap coherence into a coherence coefficient.
In particular, to suppress offset-generated artifacts and interference signals in the data, to enhance the effective signal, it is necessary to bootstrap coherent for a normalized high resolutionAnd performing coherence coefficient conversion. Specifically, a cutoff threshold is set,The coherence factor is calculated as follows:
......(11)
Wherein, Is constant to control the coherence coefficientThe weight values from 1 to 0 decay the transition.
In the present invention, the term "a" is used to refer to,、、The value of (2) may be otherwise as long as it satisfies,,And (3) obtaining the product.
(1063) The final enhanced coherent imaging result is generated using the coherence coefficient and the preliminary enhanced imaging data.
Specifically, the correlation coefficientAnd preliminary enhancement of imaging dataSubstituting formula (12) to suppress offset-induced artifacts and extract effective imaging data to obtain final enhanced coherent imaging resultsEquation (12) is as follows:
......(12)
The method comprises the steps of obtaining well logging far detection offset data, carrying out offset imaging on the well logging far detection offset data to obtain an offset imaging result, horizontally superposing the offset imaging result along the offset direction to obtain an initial superposition section, carrying out depth-direction bootstrap coherence enhancement on the initial superposition section to obtain preliminary enhancement imaging data corresponding to each position point of the initial superposition section, calculating second high-resolution bootstrap coherence at the position point corresponding to each position point of the preliminary enhancement imaging data, wherein the second high-resolution bootstrap coherence is second high-resolution bootstrap coherence of the offset imaging result in the offset direction, and generating a final enhancement coherence imaging result according to the second high-resolution bootstrap coherence and the preliminary enhancement imaging data. Based on the method, the problem of imaging artifacts generated by imaging the well-logging far detection data with limited aperture by a conventional imaging method is solved by introducing a high-precision coherence measure method, the imaging precision of the well-logging far detection data well Zhou Liefeng is effectively improved, and technical support is provided for high-precision crack imaging in oil and gas exploration and development.
In a specific example, fig. 3 is a conventional well-logging far-detection imaging result, and fig. 4 is a well-logging far-detection imaging field coherent imaging result obtained using the well-logging far-detection imaging field coherent imaging method of the present invention corresponding to fig. 3. Comparing fig. 3 and fig. 4, it can be found that the conventional well logging remote detection imaging result has low precision, imaging false images exist, the crack is not clearly delineated, the imaging false images are remarkably suppressed after the method is applied, the crack is clearly delineated, and the imaging precision is improved.
Based on one general inventive concept, the invention also provides a well logging remote detection imaging domain coherent imaging device. Fig. 5 is a schematic structural diagram of a coherent imaging device of a logging remote detection imaging domain according to an embodiment of the present invention. As shown in fig. 5, the present apparatus includes:
The acquisition module 51 is configured to acquire the far detection offset data. The well logging far detection sub offset data is three-dimensional data with three dimensions of depth, time and offset.
The offset imaging module 52 is configured to perform offset imaging on the far log detection offset data to obtain an offset imaging result.
And the superposition module 53 is configured to horizontally superimpose the split offset imaging results along the offset direction, so as to obtain an initial superimposed profile.
The preliminary enhancement module 54 is configured to perform depth-direction bootstrapping coherence enhancement on the initial superimposed profile, so as to obtain preliminary enhanced imaging data corresponding to each position point of the initial superimposed profile.
The calculation module 55 is configured to calculate, for each location point corresponding to the preliminary enhanced imaging data, a second high-resolution bootstrap coherence at the location point, where the second high-resolution bootstrap coherence is a second high-resolution bootstrap coherence of the sub-offset imaging result in the offset direction.
A generation module 56 for generating a final enhanced coherent imaging result from the second high resolution bootstrapping coherence and the preliminary enhanced imaging data.
Optionally, the acquiring module 51 may specifically be configured to:
(1) And acquiring logging remote detection pre-stack data.
(2) And recombining and arranging the logging remote detection pre-stack data to obtain the logging remote detection sub-offset data.
Optionally, the preliminary enhancement module 54 may specifically be configured to:
(1) Each scanning slope is determined based on a preset maximum scanning slope, a preset minimum scanning slope and a preset number of scanning slopes in the depth direction.
(2) And aiming at each position point of the initial superposition profile, acquiring first partial window data of the position point relative to a first partial window along the direction of each scanning slope, wherein the first partial window is a window with depth direction length and radial distance direction length.
(3) First high-resolution bootstrap coherence at the corresponding scan slope is calculated using the first local window data at each of the scan slopes, respectively.
(4) And extracting the preliminary enhancement effective signals along the direction of the scanning slope corresponding to the maximum first high-resolution bootstrap coherence to obtain the preliminary enhancement imaging data corresponding to the position points.
Optionally, the calculating module 55 may specifically be configured to:
(1) And obtaining second partial window data of the partial offset imaging result in the offset direction at the position point based on a second partial window, wherein the second partial window is a window with a radial distance direction length and an offset direction length.
(2) Using the second local window data, a second high resolution bootstrap coherence at the location point is calculated.
Optionally, the generating module 56 may specifically be configured to:
(1) And normalizing all the second high-resolution bootstrap coherence.
(2) And converting the normalized second high-resolution bootstrap coherence into a coherence coefficient.
(3) Generating the final enhanced coherent imaging result using the coherence coefficient and the preliminary enhanced imaging data.
Based on a general inventive concept, the present invention also provides a computer apparatus. Fig. 6 is a schematic structural diagram of a computer device according to an embodiment of the present invention. As shown in fig. 6, the apparatus 600 includes:
at least one processor 610, and
A memory 630 communicatively coupled to the at least one processor 610, wherein,
The memory 630 stores instructions 620 executable by the at least one processor 610, the instructions 620 being executable by the at least one processor 610 to enable the at least one processor 610 to implement a well logging far detection imaging domain coherence imaging method as described above.
It is to be understood that the same or similar parts in the above embodiments may be referred to each other, and that in some embodiments, the same or similar parts in other embodiments may be referred to.
It should be noted that in the description of the present invention, the terms "first," "second," and the like are used for descriptive purposes only and are not to be construed as indicating or implying relative importance. Furthermore, in the description of the present invention, unless otherwise indicated, the meaning of "plurality" means at least two.
Any process or method descriptions in flow diagrams or otherwise described herein may be understood as representing modules, segments, or portions of code which include one or more executable instructions for implementing specific logical functions or steps of the process, and additional implementations are included within the scope of the preferred embodiment of the present invention in which functions may be executed out of order from that shown or discussed, including substantially concurrently or in reverse order from that shown or discussed, depending on the functionality involved, as would be understood by those reasonably skilled in the art of the embodiments of the present invention.
It is to be understood that portions of the present invention may be implemented in hardware, software, firmware, or a combination thereof. In the above-described embodiments, the various steps or methods may be implemented in software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented in hardware, as in another embodiment, may be implemented using any one or combination of techniques known in the art, discrete logic circuits with logic gates for implementing logic functions on data signals, application specific integrated circuits with appropriate combinational logic gates, programmable Gate Arrays (PGAs), field Programmable Gate Arrays (FPGAs), and the like.
Those of ordinary skill in the art will appreciate that all or a portion of the steps carried out in the method of the above-described embodiments may be implemented by a program to instruct related hardware, where the program may be stored in a computer readable storage medium, and where the program, when executed, includes one or a combination of the steps of the method embodiments.
In addition, each functional unit in the embodiments of the present invention may be integrated in one processing module, or each unit may exist alone physically, or two or more units may be integrated in one module. The integrated modules may be implemented in hardware or in software functional modules. The integrated modules may also be stored in a computer readable storage medium if implemented in the form of software functional modules and sold or used as a stand-alone product.
The above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, or the like.
In the description of the present specification, a description referring to terms "one embodiment," "some embodiments," "examples," "specific examples," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiments or examples. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.
While embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and not to be construed as limiting the invention, and that variations, modifications, alternatives and variations may be made to the above embodiments by one of ordinary skill in the art within the scope of the invention.
Claims (6)
1. A method of well logging far detection imaging domain coherent imaging, comprising:
acquiring well logging far detection offset data, wherein the well logging far detection offset data is three-dimensional data with three dimensions of depth, time and offset;
Performing sub-offset imaging on the well logging remote detection sub-offset data to obtain a sub-offset imaging result;
Horizontally superposing the sub-offset imaging result along the offset direction to obtain an initial superposition profile;
Performing depth-direction bootstrap coherent enhancement on the initial superposition profile to obtain preliminary enhancement imaging data corresponding to each position point of the initial superposition profile;
Calculating second high-resolution bootstrap coherence at a position point corresponding to each preliminary enhanced imaging data, wherein the second high-resolution bootstrap coherence is second high-resolution bootstrap coherence of the sub-offset imaging result in the offset direction;
Generating a final enhanced coherent imaging result according to the second high resolution bootstrapping coherence and the preliminary enhanced imaging data;
The calculating the second high-resolution bootstrap coherence at the location point specifically includes:
Based on a second local window, obtaining second local window data of the partial offset imaging result in the offset direction at the position point, wherein the second local window is a window with a radial distance direction length and an offset direction length;
calculating a second high-resolution bootstrap coherence at the location point using the second local window data based on a second high-resolution bootstrap coherence calculation formula, the second high-resolution bootstrap coherence calculation formula being as follows:
wherein, Is expressed at a position pointBootstrap data with r second partial windowsThe second high-resolution bootstrap coherence obtained by calculation; By applying to the second partial window data Obtained by non-repeated bootstrap, meansIs the r second partial window bootstrapping data;, indicating the r second partial window bootstrap data Is a coefficient of coherence of (2); representing the length of the corresponding local window in the offset direction; representing the length of the corresponding partial window in the radial distance direction; Representing data according to a second partial window And calculating a second similarity coefficient.
2. The method for coherent imaging in the far-ranging imaging domain according to claim 1, wherein the acquiring the far-ranging offset data specifically comprises:
Acquiring logging far detection pre-stack data;
And recombining and arranging the logging remote detection pre-stack data to obtain the logging remote detection sub-offset data.
3. The method for coherent imaging of a well logging far detection imaging domain according to claim 1, wherein the step of performing depth-wise bootstrap coherent enhancement on the initial superimposed profile to obtain preliminary enhanced imaging data corresponding to each position point of the initial superimposed profile specifically comprises:
Determining each scanning slope based on a preset maximum scanning slope, a preset minimum scanning slope and a preset scanning slope number in the depth direction;
For each position point of the initial superposition profile, respectively acquiring first partial window data of the position point relative to a first partial window along the direction of each scanning slope, wherein the first partial window is a window with depth direction length and radial distance direction length;
calculating a first high-resolution bootstrap coherence under the corresponding scanning slope by using the first local window data under each scanning slope, wherein the calculation formula of the first high-resolution bootstrap coherence is as follows:
wherein, Expressed in the scanning slopeDirection bootstrapping data with n first partial windowsA first high resolution bootstrap coherence calculated; By applying to the first partial window data Obtained by non-repeated bootstrap, meansN first partial window bootstrapping data;, indicating nth first partial window bootstrap data Is a coefficient of coherence of (2); representing the length of the corresponding partial window in the radial distance direction; Representing the length of the corresponding partial window in the depth direction; representing data according to a first partial window Calculating a first similarity coefficient;
and extracting the preliminary enhancement effective signals along the direction of the scanning slope corresponding to the maximum first high-resolution bootstrap coherence to obtain the preliminary enhancement imaging data corresponding to the position points.
4. The method of imaging domain coherence imaging for well logging of claim 1, wherein generating final enhanced coherence imaging results from the second high resolution bootstrap coherence and the preliminary enhanced imaging data, comprises:
normalizing all the second high-resolution bootstrap coherence;
converting the normalized second high-resolution bootstrap coherence into a coherence coefficient;
generating the final enhanced coherent imaging result using the coherence coefficient and the preliminary enhanced imaging data.
5. A far-detection imaging field coherent imaging apparatus, applied to the far-detection imaging field coherent imaging method according to any one of claims 1 to 4, comprising:
the acquisition module is used for acquiring the well logging far detection partial offset data, wherein the well logging far detection partial offset data is three-dimensional data with three dimensions of depth, time and offset;
The sub offset imaging module is used for performing sub offset imaging on the well logging far detection sub offset data to obtain a sub offset imaging result;
the superposition module is used for horizontally superposing the sub-offset imaging result along the offset direction to obtain an initial superposition section;
The initial enhancement module is used for carrying out depth-direction bootstrap coherent enhancement on the initial superposition profile to obtain initial enhancement imaging data corresponding to each position point of the initial superposition profile;
The calculation module is used for calculating second high-resolution bootstrap coherence at the position point corresponding to each initial enhancement imaging data, wherein the second high-resolution bootstrap coherence is second high-resolution bootstrap coherence of the sub-offset imaging result in the offset direction;
And the generation module is used for generating a final enhanced coherent imaging result according to the second high-resolution bootstrap coherence and the preliminary enhanced imaging data.
6. A computer device, comprising:
at least one processor, and
A memory communicatively coupled to the at least one processor, wherein,
The memory stores instructions executable by the at least one processor to enable the at least one processor to implement the well logging far ranging imaging field coherence imaging method of any one of claims 1 to 4.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411513967.XA CN119024438B (en) | 2024-10-29 | 2024-10-29 | Logging remote detection imaging domain coherent imaging method, device and equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202411513967.XA CN119024438B (en) | 2024-10-29 | 2024-10-29 | Logging remote detection imaging domain coherent imaging method, device and equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN119024438A CN119024438A (en) | 2024-11-26 |
| CN119024438B true CN119024438B (en) | 2025-01-24 |
Family
ID=93535682
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202411513967.XA Active CN119024438B (en) | 2024-10-29 | 2024-10-29 | Logging remote detection imaging domain coherent imaging method, device and equipment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN119024438B (en) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113866825A (en) * | 2021-08-23 | 2021-12-31 | 中国石油大学(华东) | Angle-domain least-squares reflectivity inversion method based on coherent superposition |
| CN114114417A (en) * | 2020-08-25 | 2022-03-01 | 中国石油天然气股份有限公司 | Seismic imaging method and device |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100118654A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling migration method |
| US9841518B2 (en) * | 2014-02-26 | 2017-12-12 | Schlumberger Technology Corporation | Noise attenuation |
| CN116307133A (en) * | 2023-02-28 | 2023-06-23 | 成都理工大学 | Tight oil and gas reservoir natural fracture prediction method, system, medium, equipment and terminal |
| CN117890904B (en) * | 2024-03-14 | 2024-06-21 | 中南大学 | Method for extracting diffraction attribute field of ground penetrating radar, storage medium and equipment |
| CN118566986A (en) * | 2024-05-20 | 2024-08-30 | 广东海洋大学 | Submarine cold spring plume seismic imaging processing method |
-
2024
- 2024-10-29 CN CN202411513967.XA patent/CN119024438B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114114417A (en) * | 2020-08-25 | 2022-03-01 | 中国石油天然气股份有限公司 | Seismic imaging method and device |
| CN113866825A (en) * | 2021-08-23 | 2021-12-31 | 中国石油大学(华东) | Angle-domain least-squares reflectivity inversion method based on coherent superposition |
Also Published As
| Publication number | Publication date |
|---|---|
| CN119024438A (en) | 2024-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10520619B2 (en) | FWI model domain angle stacks with amplitude preservation | |
| CN104237940B (en) | Diffracted wave imaging method and device based on dynamic characteristics | |
| Spica et al. | Shallow VS imaging of the Groningen area from joint inversion of multimode surface waves and H/V spectral ratios | |
| CN105116448B (en) | A kind of converted wave azimuthal anisotropy bearing calibration and device | |
| Qiu et al. | Eikonal tomography of the Southern California plate boundary region | |
| CN105182408A (en) | Manufacturing method and device for synthesizing earthquake record | |
| CN116165707B (en) | A tunnel engineering detection method and device based on viscoelastic wave field simulation | |
| Delorey et al. | Basin shear-wave velocities beneath Seattle, Washington, from noise-correlation Rayleigh waves | |
| CN112180444A (en) | Detection method and device for stratum speed structure and storage medium | |
| CN111880223B (en) | Attribute analysis method for well-side sound wave remote detection imaging fracture-cave carving | |
| CN104570090B (en) | The extraction of full waveform inversion noise filter operator and the method filtered using its noise | |
| CN119024438B (en) | Logging remote detection imaging domain coherent imaging method, device and equipment | |
| CN111158050A (en) | Data acquisition system and method and tunnel seismic wave advanced prediction method | |
| CN109307887B (en) | Earthquake weak reflector identification method and system | |
| Eken et al. | The use of direct shear waves in quantifying seismic anisotropy: exploiting regional arrays | |
| Kong et al. | Applicability of the multiple-event stacking technique for shear-wave splitting analysis | |
| Yang et al. | Corner frequency ratios of P and S waves and strain drops of earthquakes recorded by a tight network around the Karadere segment of the North Anatolian Fault Zone: evidence for non-classical source processes | |
| Yang et al. | A frequency-Hankel transform method to extract multimodal Rayleigh wave dispersion spectra from active and passive source surface wave data | |
| CN114488346B (en) | Underground space abnormal body detection method, device, equipment and medium | |
| US20240069230A1 (en) | Method and system for processing seismic images to progressively enhance an rgt image of a geological formation | |
| US20170336523A1 (en) | Seismic signal processing method, apparatus and system | |
| CN115421195A (en) | Method, device and equipment for generating velocity field in seismic exploration and storage medium | |
| KR102806861B1 (en) | Automated method to calculate surface-wave velocities from beamforming surface-wave tests | |
| CN109490952B (en) | Seismic coherence analysis method and system | |
| Kinoshita et al. | Tilt motions recorded at two WISE sites for the 2003 Tokachi-Oki earthquake (M 8.3) |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |