CN110501757B - Point-line-surface three-dimensional point-by-point well depth design method - Google Patents
Point-line-surface three-dimensional point-by-point well depth design method Download PDFInfo
- Publication number
- CN110501757B CN110501757B CN201910572832.3A CN201910572832A CN110501757B CN 110501757 B CN110501757 B CN 110501757B CN 201910572832 A CN201910572832 A CN 201910572832A CN 110501757 B CN110501757 B CN 110501757B
- Authority
- CN
- China
- Prior art keywords
- point
- well depth
- line
- data
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a point-line-surface three-dimensional point-by-point well depth design method, which comprises the following steps: carrying out analytic inversion on different offset distance sections by using old data of a work area to obtain a near-surface forward model; designing forward modeling simulation shallow seismic line acquisition parameters; seeking the speed corresponding relation of an inversion model of the shallow seismic line for correction; searching a speed interface of an excited lithologic layer on the basis of the corrected shallow seismic line chromatography inversion model; and carrying out 1 km-1 km micro-logging survey on the work area, layering the lithology of the work area one by one, and carrying out three-dimensional interpolation to form a well depth design drawing. Performing same-parameter chromatographic inversion on old data passing through the three-dimensional work area, and constraining and correcting the result of the step 5; and carrying out three-dimensional point-by-point excitation well depth design and mapping of old data of a work area and shallow seismic line chromatography inversion constraints. The method can effectively guarantee the rationality and accuracy of the mountain land data excitation well depth design, thereby improving the quality of the collected seismic data and laying a good foundation for subsequently improving the imaging precision of prestack depth migration.
Description
Technical Field
The invention relates to the technical field of oilfield development, in particular to a point-line-surface three-dimensional point-by-point well depth design method.
Background
The surface relief and surface structure change of the mountain land are large, due to the influence of complex conditions of the surface, the near surface and the underground, interference waves develop in areas with low signal to noise ratio, the quality of seismic data is poor, and the static correction problem influencing the seismic data processing effect is very prominent, so the importance of surface structure investigation work is undoubted, and meanwhile, accurate investigation results can provide accurate basis for well depth design. At present, the mainstream mountain surface layer investigation method is to adopt a encountering observation method or a chasing blasting observation method in a small refraction method according to the conventional density, and only to use a micro-logging method to verify at a line-measuring intersection. However, the small refraction method is suitable for being carried out under the condition that the whole arrangement height difference is less than 2m, objective conditions are obviously absent in mountainous areas, and meanwhile, the method can cause the density of micro-logging survey points with higher accuracy to be too low, so that the surface layer survey result is difficult to ensure.
The earthquake excitation technology is a key technology for improving the signal-to-noise ratio of earthquake data in a low signal-to-noise ratio area, and the patent is an in-depth research and analysis of the earthquake excitation technology by combining theory and practice. The basic idea is to utilize comprehensive surface structure survey of various methods combined with forward simulation, actual test, point-to-surface type to clarify the surface structure of the complex earth surface and to utilize a dynamic design technology based on the terrain and the surface structure to design the well depth, thereby effectively improving the seismic data quality in the area with low signal-to-noise ratio.
Similar implementation is mentioned in the chinese patent application with application number 201310356187.4 in a method for designing cannon well depth based on micro-logging and non-seismic data: (1) carrying out exploration area micro-logging well drilling and measurement, obtaining a micro-logging oscillogram and an interpretation result through processing and interpretation, and carrying out preliminary well depth design to obtain excitation depth values on each excitation point of the earthquake; (2) adopting a non-seismic method to carry out surface structure investigation, carrying out geophysical inversion on the obtained data, and controlling by using micro-logging data to obtain a surface structure explanation sectional diagram; (3) displaying each excitation depth value obtained in the step (1) in the corresponding surface structure explanation sectional view obtained in the step (2), and if the excitation point depth is not located in the favorable excitation layer, adjusting the excitation point depth to the favorable excitation layer; if the depth of the excitation point is in the beneficial layer, the depth of the excitation point is adjusted to a shallower depth; (4) and adding the adjusted excitation depth to the length of the explosive column to obtain the comprehensively designed cannon well depth. Wherein the preliminary well depth design of (1) comprises: determining a layer with good excitation effect, and determining the optimal excitation depth on each micro-logging oscillogram; and obtaining the excitation depth value of each excitation point of the earthquake according to a linear interpolation method by using the optimal excitation depth of all the micro-logging in the research area. (2) The non-seismic method comprises exploration methods such as gravity, magnetic force, resistivity method and electromagnetic method, and the density of the measuring points is higher than that of the micro-logging method. The application starts from model establishment and forward modeling, starts with data analysis of experimental points and experimental lines, takes chromatographic inversion and micro-logging constraint as means, can reflect complex transverse change of the earth surface more finely, better clears the surface structure, and is more suitable for seismic data acquisition well depth design of complex earth surface mountain front zones.
Similar implementation processes are mentioned in a point-by-point excitation well depth design method based on a consistent near-surface model in a Chinese patent application with the application number of 201010011424. X: (1) performing joint interpretation on test data obtained by three methods of lithology coring, static sounding and micro-logging lithology calibration to obtain an optimal lithology layering interface; (2) jointly interpreting test data obtained by micro-logging lithology calibration and small refraction test data to obtain optimal speed hierarchical data; (3) classifying and editing the obtained speed hierarchical data, lithologic hierarchical interface data, speed, absorption attenuation Q value, diving surface position data and well data; converting the speed hierarchical data, the lithologic hierarchical interface data and the diving surface position data into a hierarchical data format for loading; loading well data in a self-defined data format; compiling an attribute load for the well data for the velocity and absorption attenuation Q value data; (4) gridding lithologic layered interface data, speed layered data and diving surface position data obtained after joint interpretation as scatter data to generate gridded data; (5) generating a framework model based on block modeling by using a binary tree method; (6) editing the bedding surface according to the lithologic hierarchical interface data, the speed hierarchical data and the relative position relationship and the mutual constraint relationship of the diving surface position data, and eliminating singular points and smoothing the bedding surface; (7) loading well data, speed data and an absorption attenuation Q value, and finally generating a consistent near-surface model based on multivariate parameters; (8) calculating the boundary range of the optimal excitation lithology, diving surface and speed layered interface as the optimal excitation position of the excitation well depth; (9) processing the calculated optimal excitation position as layer data, and dispersing the layer data according to the field construction design requirement, wherein the dispersed data is the excitation well depth of each excitation point in the field; (10) writing a data interface program, and outputting the excitation well depth value of each excitation point in a table form; (11) and reloading the generated excitation well depth scatter data into the system to generate a contour map in a CGM format, and drawing by using a plotter. The optimal stimulated lithologic layer is obtained by analyzing the quality of data on an actual test point, then is gradually extrapolated to a whole work area from point to line, the recursion of each step is supported by the actually collected data, the shallow seismic line and the chromatographic inversion data of old data are utilized for constraint, the depth of the optimal stimulated lithologic layer is corrected, and the precision and the reliability of well depth design are improved.
The western region has severe surface fluctuation and interference wave development, and the seismic data quality is poor due to the influence of complex conditions of the surface, the near surface and the underground. The thickness and the speed of the low and deceleration zone are complex and changeable, the precision of surface layer structure investigation carried out by the traditional small refraction and micro-logging method is low, an accurate near-surface model of the mountain front zone is difficult to establish, and the excitation well designed on the basis of the accurate near-surface model is inaccurate, so that the seismic data quality and subsequent near-surface correction are influenced. Therefore, a novel point-line-surface three-dimensional point-by-point well depth design method is invented, and the technical problems are solved.
Disclosure of Invention
The invention aims to provide a point-line-surface three-dimensional point-by-point well depth design method capable of establishing an accurate near-surface velocity model.
The object of the invention can be achieved by the following technical measures: the point-line-surface three-dimensional point-by-point well depth design method comprises the following steps: step 1, carrying out analytic inversion on different offset sections by using old data of a work area to obtain a near-surface forward model; step 2, designing forward modeling simulation shallow seismic line acquisition parameters; step 3, comparing the corresponding relation between the optimal stimulated lithologic layer speed of the test point and the data quality, and seeking the speed of the shallow seismic line inversion model for correction; step 4, searching a speed interface for exciting a lithologic layer on the basis of the corrected shallow seismic line chromatographic inversion model; step 5, carrying out 1 km-1 km micro-logging survey on the work area, layering lithology port by port, and carrying out three-dimensional interpolation to form a well depth design drawing; step 6, performing same-parameter chromatographic inversion on old data passing through the three-dimensional work area, and constraining and correcting the result of the step 5; and 7, carrying out three-dimensional point-by-point excitation well depth design and mapping of old data of the work area and shallow seismic line chromatography inversion constraints.
The object of the invention can also be achieved by the following technical measures:
the point-line-surface three-dimensional point-by-point well depth design method further comprises the step of carrying out lithologic coarsening and partitioning on the earth surface through field earth surface reconnaissance and micro-logging investigation before the step 1.
In step 1, performing chromatographic inversion on different migration distance sections by using old seismic data in a work area range, analyzing a series of chromatographic inversion models, and selecting migration distance sections capable of stably depicting changes of a near-surface structure to establish an entity forward model.
In step 2, high-density and small-track-pitch forward modeling is carried out, parameter degradation chromatography inversion analysis is carried out, an observation system of shallow seismic lines is determined, and seismic data are collected.
The point-line-surface three-dimensional point-by-point well depth design method further comprises the step 2 of selecting test points and researching excitation receiving factors to determine the relation between the speed of the optimal excitation lithology layer and the data quality.
In the step of selecting test points and researching the excitation receiving factors, the relation between the speed of an excitation lithologic layer and the quality of seismic data is researched and determined by researching the seismic data of different well depths, dosages, well combinations, detector combinations and detector embedding factors and comparing single-hole micro-logging tomography inversion curves.
In step 3, lithology logging and micro-logging data are analyzed, and the thickness errors of the low-speed layer and the deceleration layer of the shallow seismic line chromatographic inversion model are corrected by combining a micro-logging lithology layering speed curve.
In step 4, the velocity of the stimulated lithologic layer is utilized to search an equal velocity interface from the shallow seismic line tomography inversion model, thickness constraint correction is carried out in step 3, and the well depth corresponding to the velocity interface is selected as the stimulated well depth.
The point-line-surface three-dimensional point-by-point well depth design method further comprises the following steps of after the step 4, collecting seismic data on a test line according to the depth determined in the step 4 as an excitation well depth; and (3) performing parameter degradation analysis on the shallow seismic lines acquired in the step (2) and laying the crossed shallow seismic lines in the work area.
And in the step of acquiring seismic data on the test line according to the depth determined in the step 4 as the excitation well depth, performing data trial acquisition and processing on the test line, and analyzing the single shot signal-to-noise ratio, the coverage times and the superposition effect of the data, thereby improving the feasibility and the accuracy of the well depth design method.
In the step of carrying out parameter degradation analysis and carrying out the layout of the criss-cross shallow seismic lines in the work area, carrying out chromatography inversion degradation on the shot distance and the track distance parameters of the shallow seismic line data designed in the step 2, selecting stable shot distance and track distance parameters capable of reflecting the change of the near-surface structure, and laying typical criss-cross shallow seismic lines in the three-dimensional work area as constraint conditions for carrying out three-dimensional well depth design.
In step 5, performing chromatographic inversion on all micro-logging data in the work area range, performing lithology layering one by one, designing the well depth according to the speed of the optimal stimulated lithology layer, and performing Kriging interpolation imaging in a three-dimensional space.
In step 6, all old seismic data in the work area are subjected to same-parameter chromatographic inversion, effective closure of three shallow seismic lines newly acquired in the three-dimensional work area is guaranteed, the three shallow seismic lines are used as constraint conditions, the designed well depth is corrected, and three-dimensional point-by-point excitation well depth design is completed.
In step 7, the calculated three-dimensional point-by-point well depth design data and plane distribution are projected to a near-surface model of shallow seismic line chromatography inversion in the work area, the depth error condition is compared, and the constraint condition in step 5 is modified until the precision requirement is met.
According to the point-line-surface three-dimensional point-by-point well depth design method, diversified information is comprehensively utilized, an accurate near-surface speed model can be established, the lithology change relation of the near-surface can be clarified, the point-by-point excitation well depth design of the targeted space variation is developed by gradually extrapolating in the test points, the test lines and the three-dimensional work area plane, the rationality and the accuracy of the mountain land data excitation well depth design can be effectively guaranteed, the quality of collected seismic data is improved, and a good foundation is laid for the follow-up improvement of the imaging precision of pre-stack depth migration.
The invention provides a point-line-surface three-dimensional point-by-point well depth design method, which aims at the situation of double complex work areas, and the traditional method establishes the situation of low precision of a near-surface model and inaccurate excitation well depth. Compared with the prior well depth design method, the method has three main advantages:
(1) the conventional investigation method of the near-surface model in the mountain front zone mainly utilizes small refraction and micro-logging, and the precision is not enough, shallow seismic lines with small track spacing and small shot spacing are designed according to the actual geological condition, and the modeling precision of the near-surface model is effectively improved by utilizing the chromatographic inversion of the shallow seismic lines;
(2) according to the method, a near-surface structure is checked on the basis of (1) according to a gradually-propelled technical idea, the relation between the excitation well depth and the excitation speed is implemented by combining a point test and section data, the well depth is designed to be expanded to the whole three-dimensional work area from a point and a surface, and a plurality of series of constraint conditions are utilized in the middle of the well depth, so that the method has strong systematicness and practicability;
(3) the implementation process can be monitored completely, the field collecting and well drilling efficiency is effectively improved, the field collecting is practically guided, and a good effect is achieved according to the actual implementation of the design scheme.
Drawings
FIG. 1 is a flow chart of an embodiment of a point-line-surface three-dimensional point-by-point well depth design method of the present invention;
FIG. 2 is a diagram illustrating the effect of combining a shallow seismic line tomography inversion model on surface coarsening zoning in accordance with an embodiment of the present invention;
FIG. 3 is a graph of a three-dimensional micro-log of discrete points in accordance with an embodiment of the present invention;
FIG. 4 is a lithology distribution and tomographic inversion velocity profile of a three-dimensional micro-log in an embodiment of the present invention;
FIG. 5 is a schematic diagram of the inverse model closure results for all old two-dimensional lines and shallow seismic lines in an embodiment of the invention;
FIG. 6 is a three-dimensional well depth design depth plan without old data constraints in an embodiment of the present invention;
FIG. 7 is a three-dimensional well depth design depth plan utilizing old data constraints in an embodiment of the present invention;
FIG. 8 is a projection of a planned excitation well depth onto a shallow seismic line in an embodiment of the invention.
Detailed Description
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
As shown in fig. 1, fig. 1 is a flow chart of the point-line-surface three-dimensional point-by-point well depth design method of the present invention.
102, establishing a near-surface forward model for the offset chromatographic inversion of old data in a work area;
103, forward modeling the design of acquisition parameters of the simulated shallow seismic lines; performing high-density and small-track-pitch forward modeling, performing parametric regression analysis, determining an observation system of a shallow seismic line and collecting seismic data;
104, selecting test points and exciting receiving factor research to determine the relation between the speed of the optimal excitation lithology layer and the data quality;
106, searching a speed interface for exciting a lithologic layer on the basis of the corrected shallow seismic line tomography inversion model;
and 108, performing parameter degradation analysis on the shallow seismic lines acquired in the step 103, and performing layout of the crossed shallow seismic lines in the work area by combining the step 101.
And step 109, carrying out 1km x 1km micro-logging survey on the work area, layering lithology port by port, exciting a well depth control point by combining with the dispersion of 101 and 104, and carrying out Krigin interpolation mapping in a three-dimensional space.
And step 110, performing isoparametric tomographic inversion on the old data passing through the three-dimensional work area, ensuring that the old data and a newly acquired shallow seismic line tomographic inversion model are effectively closed, and correcting the result of the step 109 by taking the closed shallow seismic line tomographic inversion model as a constraint.
And step 111, designing and mapping the three-dimensional point-by-point excitation well depth of the old data of the work area and shallow seismic line chromatography inversion constraints.
In one embodiment of the present invention, the method comprises the following steps:
(1) performing coarsening and zoning on the lithology of the earth surface through field earth surface reconnaissance and micro-logging investigation;
(2) carrying out chromatographic inversion on different migration distance sections by using old seismic data in a work area range, analyzing a series of chromatographic inversion models, and selecting migration distance sections capable of stably depicting the change of a near-surface structure to establish an entity forward model;
(3) forward modeling of the shot spacing 2m and the track spacing 2m is carried out on the established near-surface model, chromatographic inversion of degradation of different shot spacings and track spacings is carried out, and shot spacing and track spacing parameters of near-surface changes can be stably reflected through contrastive analysis, so that the design of an observation system of a shallow seismic line is completed.
(4) By researching seismic data of different factors such as well depth, dosage, well combination, detector combination and detector embedding, and comparing single-hole micro-logging chromatographic inversion curves, the relation between the speed of the stimulated lithologic layer and the quality of the seismic data is researched and determined.
(5) And analyzing lithology logging and micro-logging data, and correcting the thickness errors of the low-speed layer (loess) and the low-speed layer (gravel) of the shallow seismic line chromatographic inversion model by combining a micro-logging lithology layering speed curve.
(6) And (5) searching an equal-velocity interface from the shallow seismic line tomography inversion model by using the stimulated lithologic layer velocity determined in the step (4), performing thickness constraint correction, and selecting the well depth corresponding to the velocity interface as the stimulated well depth.
(7) Data trial collection and processing are carried out on a test line, the single shot signal-to-noise ratio, the covering times, the overlapping effect and the like of the data are analyzed, and the feasibility and the accuracy of the well depth design method are improved.
(8) And (4) carrying out chromatography inversion degradation on the shot distance and track distance parameters of the shallow seismic line data designed in the step (3), selecting stable shot distance and track distance parameters capable of inverting the change of the near-surface structure, and laying typical crisscross shallow seismic lines in the three-dimensional work area as constraint conditions for developing three-dimensional well depth design.
(9) And (3) carrying out chromatographic inversion on all micro-logging data in the range of the work area, carrying out lithology layering one by one, designing the well depth according to the speed of the optimal excitation lithology layer in the step (4), and then carrying out Kriging interpolation imaging in a three-dimensional space.
(10) And (3) performing chromatographic inversion with the same parameters on all old seismic data in the work area, ensuring that three shallow seismic lines newly acquired in the three-dimensional work area are effectively closed, and correcting the well depth designed in the step (9) by taking the three shallow seismic lines as constraint conditions to complete the three-dimensional point-by-point excitation well depth design.
(11) And (3) projecting the calculated three-dimensional point-by-point well depth design data and plane distribution to a shallow seismic line chromatography inversion near-surface model in the work area, comparing the depth error condition, and modifying the constraint condition in the step (10) until the precision requirement is met.
Fig. 2 is an effect diagram of a shallow seismic line chromatography inversion model on an earth surface coarsening subarea in an embodiment of the invention, and the earth surface of a work area is sequentially divided into a loess gravel coverage area and a loess gravel coverage Jurassic exposure area from north to south by the aid of the inversion model, so that a basis and a basis are provided for selection of a subsequent well depth correction and interpolation method;
FIG. 3 is a diagram of a three-dimensional micro-log distribution of discrete points, wherein the micro-log lines are arranged according to the area of the work area, the grid is 1km x 1km, and the local position is changed according to the actual surface condition;
FIG. 4 is a lithology distribution and tomographic inversion velocity profile of a three-dimensional micro-log in an embodiment of the invention. The lower half part of the graph is a curve diagram of time velocity pairs of a certain micro-logging, the upper half part is a chromatography inversion curve, the mutation point of the velocity, namely the depth of a well penetrating through a loess layer, is about 15m, and the corresponding lithologic layer is a gravel layer, so that the designed well depth is 15 m.
Fig. 5 is a schematic diagram of the closure results of the inversion models of all old two-dimensional lines and shallow seismic lines in an embodiment of the present invention, and it can be seen from the intersection points of different transverse and longitudinal lines that the thicknesses of the low-velocity layer and the low-velocity layer are consistent, and the closure of the inverted near-surface model is good.
Fig. 6 is a three-dimensional well depth design depth plan without using old data constraint in an embodiment of the present invention, where wells with larger depths are mainly concentrated in loess and gravel areas in the northeast of a work area, and elevation changes in the southwest of the work area are larger, a sampling grid of a micro-logging cannot well control a horizontal line change trend of a near-surface velocity model, and a well depth error is larger;
FIG. 7 is a three-dimensional well depth design depth plan constrained by old data according to an embodiment of the present invention, wherein the well depth plan after constraint is more reasonable, and particularly, the well depth in the southwest of the work area is effectively corrected;
FIG. 8 is a projection of a planned excitation well depth onto a shallow seismic line in an embodiment of the invention. As can be seen from the projection diagram, 90% of the designed well depth is excited below the optimal excitation lithology layer, and the acquisition quality of seismic data can be ensured.
The point-line-surface three-dimensional point-by-point well depth design method disclosed by the invention is a method for gradually promoting and designing complex mountain front zone point-by-point excitation well depth design by applying multi-information such as shallow seismic line chromatographic inversion, micro-logging, old seismic data and the like, and the precision of well depth design and the acquisition quality of seismic data are improved. The method utilizes the designed shallow seismic line and the constraint condition to improve the precision of the near-surface model, and other methods can also be realized by adopting an encryption micro-logging control point method, but the cost is higher. In addition, a plurality of interpolation methods are available for interpolating the well depth on the three-dimensional control points into a three-dimensional body, and the kriging interpolation method adopted by the invention can also be realized by other mathematical algorithms.
Claims (10)
1. The point-line-surface three-dimensional point-by-point well depth design method is characterized by comprising the following steps of:
step 1, carrying out analytic inversion on different offset sections by using old data of a work area to obtain a near-surface forward model;
step 2, designing forward modeling simulation shallow seismic line acquisition parameters;
step 3, seeking a speed corresponding relation of the shallow seismic line chromatography inversion model for correction;
step 4, searching a speed interface for exciting a lithologic layer on the basis of the corrected shallow seismic line chromatographic inversion model;
step 5, carrying out 1 km-1 km micro-logging survey on the work area, layering lithology port by port, and carrying out three-dimensional interpolation to form a well depth design drawing;
step 6, performing same-parameter chromatographic inversion on old data passing through the three-dimensional work area, and constraining and correcting the result of the step 5;
step 7, carrying out three-dimensional point-by-point excitation well depth design and mapping of old data of a work area and shallow seismic line chromatographic inversion constraints;
in step 4, searching an equal-velocity interface from the shallow seismic line tomography inversion model by using the stimulated lithologic layer velocity, performing thickness constraint correction in step 3, and selecting the well depth corresponding to the velocity interface as the stimulated well depth;
after the step 4, acquiring seismic data on a test line according to the depth determined in the step 4 as an excitation well depth; performing parameter degradation analysis on the shallow seismic lines acquired in the step 2, and laying the crossed shallow seismic lines in the work area; in the step of acquiring seismic data on the test line according to the depth determined in the step 4 as the excitation well depth, performing data trial acquisition and processing on the test line, analyzing the single shot signal-to-noise ratio, the coverage times and the superposition effect of the data, and improving the feasibility and the accuracy of the well depth design method; in the step of carrying out parameter degradation analysis and carrying out the layout of the criss-cross shallow seismic lines in the work area, carrying out chromatography inversion degradation on the shot distance and the track distance parameters of the shallow seismic line data designed in the step 2, selecting stable shot distance and track distance parameters capable of reflecting the change of the near-surface structure, and laying typical criss-cross shallow seismic lines in the three-dimensional work area as constraint conditions for carrying out three-dimensional well depth design.
2. The point-line-surface three-dimensional point-by-point well depth design method according to claim 1, further comprising, before step 1, performing coarsening and zoning on the surface lithology through field surface exploration and micro-logging investigation.
3. The point-line-surface three-dimensional point-by-point well depth design method according to claim 1, characterized in that in step 1, the old seismic data in the work area range are used for carrying out chromatography inversion on different migration distance sections, a series of chromatography inversion models are analyzed, and the migration distance sections capable of stably depicting the change of the near-surface structure are selected to establish the entity forward model.
4. The point-line-surface three-dimensional point-by-point well depth design method according to claim 1, characterized in that in step 2, high-density and small-track-pitch forward modeling is performed, parametric degeneration tomography inversion analysis is performed, an observation system of a shallow seismic line is determined, and seismic data are collected.
5. The point-line-surface three-dimensional point-by-point well depth design method according to claim 1, further comprising, after the step 2, performing test point selection and excited reception factor research to determine a relationship between a speed and data quality of an optimal excited lithologic layer.
6. The point-line-surface three-dimensional point-by-point well depth design method according to claim 5, characterized in that in the step of performing test point selection and excitation receiving factor research, the relation between the speed of exciting the lithologic layer and the quality of seismic data is researched and determined by researching seismic data of different well depths, dosages, well combinations, detector combinations and detector embedding factors and comparing single-hole micro-logging chromatographic inversion curves.
7. The method for designing the three-dimensional point-by-point well depth of the point-line surface according to claim 1, wherein in the step 3, lithology logging and micro logging data are analyzed, and the thickness errors of a low-speed layer and a low-speed layer of the shallow seismic line chromatography inversion model are corrected by combining a micro logging lithology layering speed curve.
8. The point-line-surface three-dimensional point-by-point well depth design method according to claim 1, characterized in that in step 5, all micro-logging data in a work area range are subjected to chromatographic inversion, lithology layering is performed one by one, the well depth is designed according to the speed of an optimal excitation lithology layer, and then kriging interpolation imaging is performed in a three-dimensional space.
9. The point-line-surface three-dimensional point-by-point well depth design method according to claim 1, characterized in that in step 6, all old seismic data in the work area are subjected to same-parameter tomographic inversion to ensure effective closure with a newly acquired shallow seismic line in the three-dimensional work area, and the designed well depth is corrected by taking the closed shallow seismic line as a constraint condition to complete three-dimensional point-by-point excitation well depth design.
10. The point-line-surface three-dimensional point-by-point well depth design method according to claims 1 and 9, characterized in that in step 7, the calculated three-dimensional point-by-point well depth design data and plane distribution are projected onto a shallow layer seismic tomography inversion near-surface model in a work area, and the constraint conditions in steps 4 and 6 are modified to compare the depth error condition until the precision requirement is met.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910572832.3A CN110501757B (en) | 2019-06-27 | 2019-06-27 | Point-line-surface three-dimensional point-by-point well depth design method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910572832.3A CN110501757B (en) | 2019-06-27 | 2019-06-27 | Point-line-surface three-dimensional point-by-point well depth design method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110501757A CN110501757A (en) | 2019-11-26 |
| CN110501757B true CN110501757B (en) | 2021-02-26 |
Family
ID=68585782
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910572832.3A Active CN110501757B (en) | 2019-06-27 | 2019-06-27 | Point-line-surface three-dimensional point-by-point well depth design method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110501757B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113267810B (en) * | 2020-02-17 | 2022-08-30 | 中国石油天然气集团有限公司 | Seismic exploration full-depth velocity modeling method and device |
| CN113835120B (en) * | 2020-06-24 | 2024-02-09 | 中国石油化工集团有限公司 | Forward modeling analysis method for different excitation energy magnitudes |
| CN114966826B (en) * | 2021-02-18 | 2023-06-30 | 中国石油化工股份有限公司 | High-low frequency fusion near-surface modeling method based on field near-surface survey data |
| CN113031069A (en) * | 2021-03-02 | 2021-06-25 | 吉林大学 | Multi-information constraint intelligent chromatography static correction method for karst area |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8670288B2 (en) * | 2009-02-04 | 2014-03-11 | Schlumberger Technology Corporation | Velocity model for well time-depth conversion |
| CN102176054B (en) * | 2011-02-18 | 2012-11-28 | 中国石油化工股份有限公司 | Near-surface comprehensive information processing explanation method |
| CN103869368B (en) * | 2012-12-10 | 2016-05-11 | 中国石油天然气集团公司 | A kind of artillery first arrival comprehensive modeling static correcting method without the constraint of near-surface investigation data |
| US10823868B2 (en) * | 2015-10-21 | 2020-11-03 | Baker Hughes Holdings Llc | Estimating depth-dependent lateral tectonic strain profiles |
| CN106569282B (en) * | 2016-10-17 | 2019-03-12 | 中国石油化工股份有限公司 | A kind of earthquake-capturing excitation well depth designing method |
| CN106990438B (en) * | 2017-03-28 | 2018-12-28 | 中国石油天然气集团公司 | Complicated earth surface region micro logging point and well depth determine method and device |
| CN109884710B (en) * | 2019-03-20 | 2021-02-26 | 中国石油化工股份有限公司 | Micro-logging tomography method aiming at excitation well depth design |
| CN109884707B (en) * | 2019-03-20 | 2021-09-07 | 中国石油化工股份有限公司 | Static correction method for near-surface layered time-depth curve |
-
2019
- 2019-06-27 CN CN201910572832.3A patent/CN110501757B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN110501757A (en) | 2019-11-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110501757B (en) | Point-line-surface three-dimensional point-by-point well depth design method | |
| US5835882A (en) | Method for determining barriers to reservoir flow | |
| Procter et al. | Spatial and layer-controlled variability in fracture networks | |
| US8793113B2 (en) | Method and apparatus for near well structural modeling based on borehole dips | |
| CN109884710B (en) | Micro-logging tomography method aiming at excitation well depth design | |
| Hosseini et al. | Geostatistical modeling and spatial distribution analysis of porosity and permeability in the Shurijeh-B reservoir of Khangiran gas field in Iran | |
| CN109799540B (en) | Volcanic rock type uranium deposit magnetic susceptibility inversion method based on geological information constraint | |
| CN111352172B (en) | Method for acquiring spatial distribution position of uranium anomaly in sand body by well-seismic combination method | |
| Zhang et al. | Application of multi-point geostatistics in deep-water turbidity channel simulation: A case study of Plutonio oilfield in Angola | |
| Allo et al. | Characterization of a carbonate geothermal reservoir using rock-physics-guided deep neural networks | |
| CN110095811B (en) | Method and device for constructing and processing velocity model of paste rock stratum | |
| Trippetta et al. | Carbonate-ramp reservoirs modelling best solutions: Insights from a dense shallow well database in Central Italy | |
| CN107526108B (en) | Carbonate rock fractured cave volume correction method | |
| Noroozi et al. | 3D Geometrical-Stochastical modeling of rock mass joint networks: case study of the right bank of Rudbar Lorestan Dam plant | |
| Liu et al. | Geostatistical modeling for fine reservoir description of Wei2 block of Weicheng oilfield, Dongpu depression, China | |
| CN107945271B (en) | Three-dimensional pressure field modeling method based on geological block tracking | |
| CN109188522A (en) | Velocity field construction method and device | |
| Dehni et al. | Implicit modeling of salinity reconstruction by using 3D combined models | |
| CN112305594A (en) | Oil-gas distribution determination method and system for heterogeneous reservoir | |
| Masoud et al. | Reservoir Characterization and Geostatistical Model of the Cretaceous and Cambrian-Ordovician Reservoir Intervals, Meghil Field, Sirte Basin, Libya | |
| CN112180443A (en) | Shale gas two-dimensional seismic sweet spot optimization method and device | |
| Li et al. | Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China | |
| Dou et al. | Reservoir Modeling of Braided River Reservoirs Based on Geological Knowledge Database: A Case Study of P1x Formation of the Daniudi Gas Field, Ordos Basin, China | |
| CN112526610B (en) | Three-dimensional seismic acquisition excitation well depth design method for constrained surface layer modeling | |
| Podnebesnykh et al. | New approach to the evaluation of the structure of initial reserves in Ozhginskoe gas-oil field |
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 |