CN112305585A - Three-dimensional data volume acquisition method and device - Google Patents
Three-dimensional data volume acquisition method and device Download PDFInfo
- Publication number
- CN112305585A CN112305585A CN201910680558.1A CN201910680558A CN112305585A CN 112305585 A CN112305585 A CN 112305585A CN 201910680558 A CN201910680558 A CN 201910680558A CN 112305585 A CN112305585 A CN 112305585A
- Authority
- CN
- China
- Prior art keywords
- data
- dimensional
- dimensional seismic
- seismic data
- well
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000011160 research Methods 0.000 claims description 22
- 238000010276 construction Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 16
- 230000001502 supplementing effect Effects 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 13
- 238000003780 insertion Methods 0.000 claims description 7
- 230000037431 insertion Effects 0.000 claims description 7
- 238000003860 storage Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 17
- 230000000694 effects Effects 0.000 description 10
- 238000012545 processing Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 238000012360 testing method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 239000003208 petroleum Substances 0.000 description 2
- 239000013589 supplement Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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/28—Processing seismic data, e.g. for interpretation or for event detection
-
- 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/282—Application of seismic models, synthetic seismograms
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)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a three-dimensional data volume acquisition method and a device, wherein the method comprises the following steps: acquiring a coordinate origin according to a two-dimensional seismic survey net base map of a three-dimensional data volume area to be generated and a geological task; according to the coordinate origin, taking main survey line data in the two-dimensional seismic data as a vertical axis, and taking tie line data in the two-dimensional seismic data as a horizontal axis to construct a virtual three-dimensional seismic survey network; acquiring two-dimensional seismic data of a three-dimensional data volume area to be generated, and inserting main survey line data and tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data; creating well site data according to the two-dimensional seismic data, and building a frame model according to the well site data and the sparse three-dimensional seismic data; and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model to obtain a three-dimensional data volume.
Description
Technical Field
The invention relates to the field of petroleum exploration, in particular to a method and a device for acquiring a three-dimensional data volume.
Background
The development of seismic exploration techniques is changing day by day, from two-dimensional to three-dimensional exploration in the early years. At present, three-dimensional exploration full coverage is basically realized in eastern oil fields, but two-dimensional seismic exploration does not exit from an exploration stage and is still a main means of seismic exploration on the basis of exploration cost, construction difficulty, history and other reasons. The three-dimensional seismic data has a large amount of information. The construction and reservoir prediction means for the three-dimensional data volume are far more than those of the two-dimensional data volume; at present, the two-dimensional network measurement density of a plurality of blocks reaches 500m multiplied by 500m, the data volume is very large, and when comprehensive research is carried out, the project operation time is short and the workload is very large; the three-dimensional seismic data body cannot be constructed quickly due to the complex construction and the large data volume, and workers can only analyze the massive two-dimensional data bodies one by one, so that the workload is large, and the research precision is reduced due to the fact that the data volume is increased.
Disclosure of Invention
The invention aims to provide a three-dimensional data volume acquisition method and a three-dimensional data volume acquisition device, which construct a three-dimensional seismic data volume by utilizing the existing two-dimensional seismic data.
In order to achieve the above object, the present invention provides a method for acquiring a three-dimensional data volume, the method comprising: acquiring a coordinate origin according to a two-dimensional seismic survey net base map of a three-dimensional data volume area to be generated and a geological task; according to the coordinate origin, taking main survey line data in the two-dimensional seismic data as a vertical axis, and taking tie line data in the two-dimensional seismic data as a horizontal axis to construct a virtual three-dimensional seismic survey network; acquiring two-dimensional seismic data of a three-dimensional data volume area to be generated, and inserting main survey line data and tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data; creating well site data according to the two-dimensional seismic data, and building a frame model according to the well site data and the sparse three-dimensional seismic data; and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model to obtain a three-dimensional data volume.
In the above three-dimensional data volume obtaining method, it is preferable that the inserting the main survey line data and the tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data includes: obtaining a preset formula according to the relation between the vertical axis and the horizontal axis; and obtaining data of each point where the two-dimensional seismic data pass in the three-dimensional seismic network according to main survey line data, the tie line data and the preset formula in the two-dimensional seismic data, and inserting the data of each point into the three-dimensional seismic network to obtain sparse three-dimensional seismic data.
In the above three-dimensional data volume acquiring method, preferably, the creating of well-site data from the two-dimensional seismic data includes: and supplementing the two-dimensional seismic data into a well position data table through a preset format to obtain well position data.
In the above three-dimensional data volume obtaining method, preferably, the establishing a frame model according to the well data and the sparse three-dimensional seismic data includes: loading the well location data and the sparse three-dimensional seismic data to a research work area; obtaining the false well data with complete relation of acoustic time difference, density, longitudinal wave impedance and time depth from the two-dimensional seismic data, supplementing the false well data into the research work area, and assigning values to the false well data by utilizing a false well output function to obtain a corresponding seismic curve; supplementing all the false well data in the two-dimensional seismic data into the research work area, and assigning values to the false well data in a manner of extracting a false curve so that the false well data obtain a corresponding seismic curve; and establishing a frame model according to the false well data, the well location data and the sparse three-dimensional seismic data.
In the above three-dimensional data volume obtaining method, preferably, interpolating each seismic data in the two-dimensional seismic data as well log data into the frame model to obtain a three-dimensional data volume includes: and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model by a local kriging method and a triangular weighting method to obtain a three-dimensional data volume.
In the above three-dimensional data volume obtaining method, preferably, according to the origin of coordinates, constructing a virtual three-dimensional seismic survey network with the main survey line data in the two-dimensional seismic data as a vertical axis and the tie line data in the two-dimensional seismic data as a horizontal axis includes: the main survey line and the connecting line are straight survey line data with consistent middle distance of the two-dimensional seismic data; and the distance between the main measuring line and the connecting line in the two-dimensional seismic data is the distance between the three-dimensional seismic measuring network.
The invention also provides a three-dimensional data volume acquisition device, which comprises a three-dimensional seismic network construction module, a three-dimensional seismic data insertion module, a frame model construction module and a three-dimensional data volume construction module; the three-dimensional seismic network construction module is used for obtaining a coordinate origin according to a two-dimensional seismic survey network base map of a three-dimensional data volume area to be generated and a geological task; according to the coordinate origin, taking main survey line data in the two-dimensional seismic data as a vertical axis, and taking tie line data in the two-dimensional seismic data as a horizontal axis to construct a virtual three-dimensional seismic survey network; the three-dimensional seismic data insertion module is used for acquiring two-dimensional seismic data of a three-dimensional data volume area to be generated, and inserting main survey line data and tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data; the frame model building module is used for building well bit data according to the two-dimensional seismic data and building a frame model according to the well bit data and the sparse three-dimensional seismic data; and the three-dimensional data volume construction module is used for interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model to obtain a three-dimensional data volume.
In the above three-dimensional data volume acquisition apparatus, preferably, the three-dimensional seismic data interpolation module includes: obtaining a preset formula according to the relation between the vertical axis and the horizontal axis; and obtaining data of each point where the two-dimensional seismic data pass in the three-dimensional seismic network according to main survey line data, the tie line data and the preset formula in the two-dimensional seismic data, and inserting the data of each point into the three-dimensional seismic network to obtain sparse three-dimensional seismic data.
In the above three-dimensional data volume acquiring apparatus, preferably, the frame model creating module includes: loading the well location data and the sparse three-dimensional seismic data to a research work area; obtaining the false well data with complete relation of acoustic time difference, density, longitudinal wave impedance and time depth from the two-dimensional seismic data, supplementing the false well data into the research work area, and assigning values to the false well data by utilizing a false well output function to obtain a corresponding seismic curve; supplementing all the false well data in the two-dimensional seismic data into the research work area, and assigning values to the false well data in a manner of extracting a false curve so that the false well data obtain a corresponding seismic curve; and establishing a frame model according to the false well data, the well location data and the sparse three-dimensional seismic data.
In the above three-dimensional data volume acquiring apparatus, preferably, the three-dimensional data volume constructing module includes: and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model by a local kriging method and a triangular weighting method to obtain a three-dimensional data volume.
The invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the computer program.
The present invention also provides a computer-readable storage medium storing a computer program for executing the above method.
The three-dimensional data volume acquisition method and the three-dimensional data volume acquisition device provided by the invention enrich reservoir prediction means on one hand, solve the problems of closure error correction, no seismic data beside a well and the like to a certain extent on the other hand, and obviously improve the precision of comprehensive research.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
fig. 1 is a schematic flow chart of a three-dimensional data volume acquiring method according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a two-dimensional direct survey line effect of the yellow river west bank collection provided by an embodiment of the present invention;
FIG. 3 is a schematic diagram of a three-dimensional seismic network and two-dimensional seismic lines provided in accordance with an embodiment of the present invention;
FIGS. 4A-4B are schematic diagrams of equivalent cross-sections of seismic lines provided by an embodiment of the invention;
FIG. 5 is a schematic illustration of the distribution of sparse seismic data on a plane, as provided by an embodiment of the present invention;
FIG. 6 is a schematic diagram illustrating the effect of smoothly generating a pseudo-seismic curve of a test well according to an embodiment of the present invention;
fig. 7 is a schematic diagram illustrating an interpolation effect of a local kriging interpolation effect according to an embodiment of the present invention;
FIG. 8 is a schematic diagram illustrating an interpolation effect of the triangular weighting method according to an embodiment of the present invention;
fig. 9 is a schematic diagram illustrating an interpolation effect of a local kriging interpolation effect according to an embodiment of the present invention;
FIG. 10 is a schematic diagram illustrating an interpolation effect of the triangular weighting method according to an embodiment of the present invention;
fig. 11 is a schematic structural diagram of a three-dimensional data volume acquiring apparatus according to an embodiment of the present invention;
FIG. 12 is a close-up schematic of a cross-sectional intersection of two-dimensional seismic data provided by an embodiment of the invention.
Detailed Description
The following detailed description of the embodiments of the present invention will be provided with reference to the drawings and examples, so that how to apply the technical means to solve the technical problems and achieve the technical effects can be fully understood and implemented. It should be noted that, unless otherwise specified, the embodiments and features of the embodiments of the present invention may be combined with each other, and the technical solutions formed are within the scope of the present invention.
Additionally, the steps illustrated in the flow charts of the figures may be performed in a computer system such as a set of computer-executable instructions and, although a logical order is illustrated in the flow charts, in some cases, the steps illustrated or described may be performed in an order different than here.
Referring to fig. 1, a three-dimensional data volume acquiring method provided by the present invention includes: s101, obtaining a coordinate origin according to a two-dimensional seismic survey net base map of a three-dimensional data volume area to be generated and a geological task; s102, according to the coordinate origin, taking main survey line data in two-dimensional seismic data as a vertical axis, and taking tie line data in the two-dimensional seismic data as a horizontal axis to construct a virtual three-dimensional seismic survey network; s103, acquiring two-dimensional seismic data of a three-dimensional data volume area to be generated, and inserting main survey line data and tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data; s104, well position data is created according to the two-dimensional seismic data, and a frame model is built according to the well position data and the sparse three-dimensional seismic data; s105, taking each seismic data in the two-dimensional seismic data as logging data, and interpolating the logging data into the frame model to obtain a three-dimensional data volume. In step S102, constructing a virtual three-dimensional seismic survey network by using main survey line data in the two-dimensional seismic data as a vertical axis and using tie line data in the two-dimensional seismic data as a horizontal axis according to the origin of coordinates includes: the main survey line and the connecting line are straight survey line data with consistent middle distance of the two-dimensional seismic data; and the distance between the main measuring line and the connecting line in the two-dimensional seismic data is the distance between the three-dimensional seismic measuring network. Creating well data from the two-dimensional seismic data in step S104 includes: and supplementing the two-dimensional seismic data into a well position data table through a preset format to obtain well position data. In actual work, before defining the three-dimensional seismic survey network, straight survey lines should be selected as much as possible, the main survey lines and the connecting lines are required to be uniformly distributed, the line distances are basically consistent, preferably within 500 meters (underground geologic bodies are changed, the line distances are larger, the errors are larger), and the CDP distances are also preferably consistent. Meanwhile, before defining the virtual three-dimensional measuring network, an actual two-dimensional seismic measuring network base map is opened, then a coordinate origin is selected on the base map according to geological tasks, generally speaking, the direction of a main measuring line is the Inline direction, the direction of a connecting measuring line perpendicular to the main measuring line is the Xline measuring line direction, and the CDP surface area elements of the three-dimensional measuring network are respectively equal to the CDP distance between the main measuring line and the connecting measuring line.
In the above embodiment, the two-dimensional seismic data is preferably acquired in the same year and with the same construction factor, and processed in the same flow, so that the section frequency, the signal-to-noise ratio, the amplitude variation range are consistent, the network density is uniform, and the closure difference at the intersection point of the survey lines is small, which can be referred to as fig. 12 in the specific case, and the intersection points of the three two-dimensional seismic sections in fig. 12 basically meet the requirement.
In an embodiment of the present invention, in the step S103, inserting the main survey line data and the tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network, and obtaining the sparse three-dimensional seismic data includes: obtaining a preset formula according to the relation between the vertical axis and the horizontal axis; and obtaining data of each point where the two-dimensional seismic data pass in the three-dimensional seismic network according to main survey line data, the tie line data and the preset formula in the two-dimensional seismic data, and inserting the data of each point into the three-dimensional seismic network to obtain sparse three-dimensional seismic data.
In an embodiment of the present invention, the building a frame model according to the well data and the sparse three-dimensional seismic data in step S104 includes: loading the well location data and the sparse three-dimensional seismic data to a research work area; obtaining the false well data with complete relation of acoustic time difference, density, longitudinal wave impedance and time depth from the two-dimensional seismic data, supplementing the false well data into the research work area, and assigning values to the false well data by utilizing a false well output function to obtain a corresponding seismic curve; supplementing all the false well data in the two-dimensional seismic data into the research work area, and assigning values to the false well data in a manner of extracting a false curve so that the false well data obtain a corresponding seismic curve; and establishing a frame model according to the false well data, the well location data and the sparse three-dimensional seismic data.
Interpolating each of the two-dimensional seismic data as well log data into the frame model in the step S105 to obtain a three-dimensional data volume comprises: and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model by a local kriging method and a triangular weighting method to obtain a three-dimensional data volume.
In order to make the above three-dimensional data volume acquisition method provided by the present invention more clearly understood, the above embodiments are described in the following by taking actual work as an example, and it should be understood by those skilled in the art that the example is not further limited thereto.
Depending on the geological task, several types of three-dimensional grids can be generated: 1. only the inline direction is considered; 2. only the tie line direction is considered; 3. the directions of the main measuring line and the connecting line are comprehensively considered. The first two types are relatively simple to implement: for example, the leftmost lower corner of the seismic survey line is selected as the origin of coordinates on the graph, and the CDP (continuous data processing) distance of the two-dimensional survey line is taken as the CDP surface element of the three-dimensional survey network. The third type is described in detail later.
Referring to fig. 2, fig. 2 is a two-dimensional straight survey line collected by the medium petroleum coalbed methane company 2013 on the west bank of the yellow river, and the CDP distances are all 10 m. 7 main measuring lines (in the east-west direction) and the line distance is 2 km; three measuring lines are connected, and the line distance is 4 km. In order to ensure that the virtual data can integrate information in east-west, south-north directions, data such as main and tie lines X, Y and CDP distances are analyzed on a two-dimensional measuring net, the origin of coordinates is selected to be located between DJ13-242 and DJ13-16-240 measuring lines, 600 inlines (the number of main lines is multiplied by 100 or 2 lines of tie lines are drawn by 1) and 650 Xline lines (the number of main lines is drawn by 2 lines or the number of middle tie lines is multiplied by 1 and the two lines on two sides are multiplied by 225) are defined, and the obtained measuring net and the original two-dimensional measuring net are shown in FIG. 3. After the virtual three-dimensional seismic survey network is established, reading the corresponding relation between each two-dimensional seismic survey line and the three-dimensional survey network on a base map (namely, figure 3), wherein for a main survey line, the change of a CDP number corresponds to the change of an Xline number; for the crossline, the change of CDP number corresponds to the change of Inline number; the specific correspondence between the Inline and Xline dynamic ranges and the two-dimensional seismic lines is shown in table 1.
TABLE 1
| Serial number | Line number | CDP Range of variation | Corresponding Inline number | Corresponding to the |
| 1 | DJ13-242 | 49-1348 | 50 | 1-650 |
| 2 | DJ13-244 | 149-1490 | 150 | 1-650 |
| 3 | DJ13-246 | 199-1496 | 250 | 1-650 |
| 4 | DJ13-248 | 558-1855 | 350 | 1-650 |
| 5 | DJ13-250 | 556-1853 | 450 | 1-650 |
| 6 | DJ13-252 | 554-1851 | 550 | 1-650 |
| 7 | DJ13-500 | 282-1482 | 1-600 | 120 |
| 8 | DJ13-504 | 410-1609 | 1-600 | 320 |
| 9 | DJ13-508 | 230-1430 | 1-600 | 520 |
After the three-dimensional seismic survey network is constructed through the definition, loading data of a two-dimensional seismic survey line, specifically, outputting the seismic survey lines in a work area one by using a Geoeast processing and interpretation integrated system developed by eastern geophysics corporation, and defining an Inline number and an Xline number by using a formula when outputting, for example, for a main survey line direction DJ13-242 survey line, the Inline number is defined as 50, and the Xline number is defined as (CDp-49)/2+ 1; for the tie direction DJ13-500 line measurements, Inline No. (CDP-282)/2+1, Xline No. (120) are defined. FIGS. 4A and 4B show the complete agreement between the amplitude, waveform and time of DJ13-242 (upper) and Inline50 (lower) profiles; thus, data loading on the three-dimensional seismic network is realized, but only the data is loaded, and the seismic data which are not copied are sparse, that is, in the three-dimensional seismic network, only the points which are passed by the two-dimensional seismic network contain the seismic data, and no seismic data exists (figure 5); therefore, data supplement is needed, namely, each path of seismic data is loaded as logging data and then is interpolated into a three-dimensional data body, specifically:
(1) bulk creation of well site data
Well data is first edited in the format:
Wellname1 KB1 TD1 X1 Y1；
Wellname2 KB2 TD2 X2 Y2；
wherein Wellname represents well name, KB represents core filling altitude, TD represents drilling completion depth, and X, Y represents coordinates; the X, Y coordinates are obtained from coordinates with seismic data, well names are generally given by the number of Inline and Xline, for example In001X001, In001X002, etc.
(2) Realizing a pseudo three-dimensional volume with seismic data in each channel
a. And (3) respectively loading the sparse three-dimensional seismic data obtained in the step (5.3.1) and the well data created in the step (1) to a research work area.
b. Copying a false well with complete relation of time difference, density, longitudinal wave impedance and time depth of sound waves to the work area, and assigning a value to the well by using the output function of the false well to ensure that a Seismic curve can be obtained (as shown in figure 6).
c. And (4) importing all the false wells, and copying the test well data to the false well curve by using a copy function.
d. And assigning values to all the false wells by using the Extract Pseudo wells function to ensure that each false well can obtain a Seismic curve (as shown in figure 7).
e. Establishing a frame Model, namely respectively introducing sparse seismic data, horizon data and a false well into a Model Builder (Without TDC) for each sub-work area, and establishing the frame Model, wherein the interpolation method of the well data comprises reverse distance weighting, local weighting, global kriging, local kriging, triangular weighting, natural neighbors and the like, and only the local kriging and the triangular weighting method can better supplement the seismic data of blank areas by comparison, and the interpolation effects of the two methods are shown in fig. 7 and fig. 8.
f. Generation of sub-work area interpolated seismic data
The seismic data generated by using the local kriging and triangular weighting method retains the characteristics of the original seismic data and changes naturally, but a small number of seismic traces exist, and the data change is very small due to the fact that the seismic traces are far away from a false well (fig. 9 and 10).
Referring to fig. 11, the present invention further provides a three-dimensional data volume obtaining apparatus, which includes a three-dimensional seismic network building module, a three-dimensional seismic data inserting module, a frame model building module, and a three-dimensional data volume building module; the three-dimensional seismic network construction module is used for obtaining a coordinate origin according to a two-dimensional seismic survey network base map of a three-dimensional data volume area to be generated and a geological task; according to the coordinate origin, taking main survey line data in the two-dimensional seismic data as a vertical axis, and taking tie line data in the two-dimensional seismic data as a horizontal axis to construct a virtual three-dimensional seismic survey network; the three-dimensional seismic data insertion module is used for acquiring two-dimensional seismic data of a three-dimensional data volume area to be generated, and inserting main survey line data and tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data; the frame model building module is used for building well bit data according to the two-dimensional seismic data and building a frame model according to the well bit data and the sparse three-dimensional seismic data; and the three-dimensional data volume construction module is used for interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model to obtain a three-dimensional data volume.
In the above embodiment, the three-dimensional seismic data insertion module comprises: obtaining a preset formula according to the relation between the vertical axis and the horizontal axis; and obtaining data of each point where the two-dimensional seismic data pass in the three-dimensional seismic network according to main survey line data, the tie line data and the preset formula in the two-dimensional seismic data, and inserting the data of each point into the three-dimensional seismic network to obtain sparse three-dimensional seismic data.
In the above embodiment, the framework model building module includes: loading the well location data and the sparse three-dimensional seismic data to a research work area; obtaining the false well data with complete relation of acoustic time difference, density, longitudinal wave impedance and time depth from the two-dimensional seismic data, supplementing the false well data into the research work area, and assigning values to the false well data by utilizing a false well output function to obtain a corresponding seismic curve; supplementing all the false well data in the two-dimensional seismic data into the research work area, and assigning values to the false well data in a manner of extracting a false curve so that the false well data obtain a corresponding seismic curve; and establishing a frame model according to the false well data, the well location data and the sparse three-dimensional seismic data.
In the above embodiment, the three-dimensional data volume construction module includes: and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model by a local kriging method and a triangular weighting method to obtain a three-dimensional data volume.
The invention also provides a computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method when executing the computer program.
The present invention also provides a computer-readable storage medium storing a computer program for executing the above method.
Since the prosperity of inversion techniques in the nineties of the last century, the wave impedance inversion technique has been widely used in the deldos basin and its surrounding areas. However, because the phenomenon of longitudinal wave impedance superposition of sand and mudstone in a target interval is very common, the technology has a serious multi-solution problem. If the AC curve is replaced by a characteristic curve which can identify lithology and maintain the basic skeleton of the AC curve, the reservoir prediction problem which troubles many years can be well solved. The advent of the present invention makes it practical for post-stack inversion to distinguish lithology. By using the method, the flow of reservoir prediction can be simplified, and the inversion result is more real and reliable. Characteristic curve inversion tests of the Suligonun block prove that by the method, the wave impedance inversion data volume which can visually reflect underground lithological changes can be obtained quickly. Compared with the traditional wave impedance inversion method, the method has no multi-solution and is easier to explain; compared with prestack inversion, the method is convenient and quick to operate, the required data volume is small (angle or offset stacking does not need to be carried out on prestack seismic data, a plurality of seismic wavelets do not need to be extracted, and transverse wave logging data do not need to be extracted), and therefore, the method can foresee that: the construction of the characteristic curve and the corresponding inversion technology can be used for developing the characteristics in reservoir prediction in the future, and play unique and irreplaceable roles.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (12)
1. A method for three-dimensional data volume acquisition, the method comprising:
acquiring a coordinate origin according to a two-dimensional seismic survey net base map of a three-dimensional data volume area to be generated and a geological task;
according to the coordinate origin, taking main survey line data in the two-dimensional seismic data as a vertical axis, and taking tie line data in the two-dimensional seismic data as a horizontal axis to construct a virtual three-dimensional seismic survey network;
acquiring two-dimensional seismic data of a three-dimensional data volume area to be generated, and inserting main survey line data and tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data;
creating well site data according to the two-dimensional seismic data, and building a frame model according to the well site data and the sparse three-dimensional seismic data;
and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model to obtain a three-dimensional data volume.
2. The method of claim 1, wherein the inserting inline data and inline data in the two-dimensional seismic data into the three-dimensional seismic network to obtain sparse three-dimensional seismic data comprises:
obtaining a preset formula according to the relation between the vertical axis and the horizontal axis;
and obtaining data of each point where the two-dimensional seismic data pass in the three-dimensional seismic network according to main survey line data, the tie line data and the preset formula in the two-dimensional seismic data, and inserting the data of each point into the three-dimensional seismic network to obtain sparse three-dimensional seismic data.
3. The method of claim 1, wherein creating well data from the two-dimensional seismic data comprises: and supplementing the two-dimensional seismic data into a well position data table through a preset format to obtain well position data.
4. The method of claim 1, wherein building a frame model from the well data and the sparse three-dimensional seismic data comprises:
loading the well location data and the sparse three-dimensional seismic data to a research work area;
obtaining the false well data with complete relation of acoustic time difference, density, longitudinal wave impedance and time depth from the two-dimensional seismic data, supplementing the false well data into the research work area, and assigning values to the false well data by utilizing a false well output function to obtain a corresponding seismic curve;
supplementing all the false well data in the two-dimensional seismic data into the research work area, and assigning values to the false well data in a manner of extracting a false curve so that the false well data obtain a corresponding seismic curve;
and establishing a frame model according to the false well data, the well location data and the sparse three-dimensional seismic data.
5. The method of claim 1, wherein interpolating each of the two-dimensional seismic data as well log data into the frame model to obtain a three-dimensional data volume comprises: and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model by a local kriging method and a triangular weighting method to obtain a three-dimensional data volume.
6. The method for acquiring the three-dimensional data volume according to any one of claims 1 to 5, wherein constructing the virtual three-dimensional seismic network by using the main survey line data in the two-dimensional seismic data as a vertical axis and the tie line data in the two-dimensional seismic data as a horizontal axis according to the origin of coordinates comprises: the main survey line and the connecting line are straight survey line data with consistent middle distance of the two-dimensional seismic data; and the distance between the main measuring line and the connecting line in the two-dimensional seismic data is the distance between the three-dimensional seismic measuring network.
7. The three-dimensional data volume acquisition device is characterized by comprising a three-dimensional seismic network construction module, a three-dimensional seismic data insertion module, a frame model construction module and a three-dimensional data volume construction module;
the three-dimensional seismic network construction module is used for obtaining a coordinate origin according to a two-dimensional seismic survey network base map of a three-dimensional data volume area to be generated and a geological task; according to the coordinate origin, taking main survey line data in the two-dimensional seismic data as a vertical axis, and taking tie line data in the two-dimensional seismic data as a horizontal axis to construct a virtual three-dimensional seismic survey network;
the three-dimensional seismic data insertion module is used for acquiring two-dimensional seismic data of a three-dimensional data volume area to be generated, and inserting main survey line data and tie line data in the two-dimensional seismic data into the three-dimensional seismic survey network to obtain sparse three-dimensional seismic data;
the frame model building module is used for building well bit data according to the two-dimensional seismic data and building a frame model according to the well bit data and the sparse three-dimensional seismic data;
and the three-dimensional data volume construction module is used for interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model to obtain a three-dimensional data volume.
8. The three-dimensional data volume acquisition device according to claim 7, wherein the three-dimensional seismic data insertion module comprises:
obtaining a preset formula according to the relation between the vertical axis and the horizontal axis;
and obtaining data of each point where the two-dimensional seismic data pass in the three-dimensional seismic network according to main survey line data, the tie line data and the preset formula in the two-dimensional seismic data, and inserting the data of each point into the three-dimensional seismic network to obtain sparse three-dimensional seismic data.
9. The apparatus according to claim 7, wherein the frame model building module comprises:
loading the well location data and the sparse three-dimensional seismic data to a research work area;
obtaining the false well data with complete relation of acoustic time difference, density, longitudinal wave impedance and time depth from the two-dimensional seismic data, supplementing the false well data into the research work area, and assigning values to the false well data by utilizing a false well output function to obtain a corresponding seismic curve;
supplementing all the false well data in the two-dimensional seismic data into the research work area, and assigning values to the false well data in a manner of extracting a false curve so that the false well data obtain a corresponding seismic curve;
and establishing a frame model according to the false well data, the well location data and the sparse three-dimensional seismic data.
10. The three-dimensional data volume obtaining apparatus according to claim 7, wherein the three-dimensional data volume construction module comprises: and interpolating each seismic data in the two-dimensional seismic data as logging data into the frame model by a local kriging method and a triangular weighting method to obtain a three-dimensional data volume.
11. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 6 when executing the computer program.
12. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910680558.1A CN112305585B (en) | 2019-07-26 | 2019-07-26 | Three-dimensional data volume acquisition method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910680558.1A CN112305585B (en) | 2019-07-26 | 2019-07-26 | Three-dimensional data volume acquisition method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112305585A true CN112305585A (en) | 2021-02-02 |
| CN112305585B CN112305585B (en) | 2024-08-27 |
Family
ID=74329653
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910680558.1A Active CN112305585B (en) | 2019-07-26 | 2019-07-26 | Three-dimensional data volume acquisition method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112305585B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113267815A (en) * | 2021-07-07 | 2021-08-17 | 中海油田服务股份有限公司 | Filtering method and device for repeated broken edge data |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4672545A (en) * | 1984-04-06 | 1987-06-09 | Pennzoil Company | Method and apparatus for synthesizing three dimensional seismic data |
| US20080144436A1 (en) * | 2006-12-13 | 2008-06-19 | Stephen Charles Pickering | Time-Lapse Seismic Acquisition |
| CN101663597A (en) * | 2007-02-13 | 2010-03-03 | 兰得马克制图公司 | Performing 3-D scatterer imaging from 2-D seismic data |
| CN105510967A (en) * | 2015-12-30 | 2016-04-20 | 中国石油天然气集团公司 | Processing method and device for achieving three-dimensional seismic data migration |
| CN107703542A (en) * | 2017-08-31 | 2018-02-16 | 中国石油天然气集团公司 | A kind of determination method and apparatus of seismic horizon |
| CN107783187A (en) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | It is a kind of that logging speed and seismic velocity are combined to the method for establishing 3D velocity field |
| CN108254784A (en) * | 2017-12-11 | 2018-07-06 | 中国石油天然气股份有限公司 | Fault identification method, device and system based on two-dimensional seismic data |
-
2019
- 2019-07-26 CN CN201910680558.1A patent/CN112305585B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4672545A (en) * | 1984-04-06 | 1987-06-09 | Pennzoil Company | Method and apparatus for synthesizing three dimensional seismic data |
| US20080144436A1 (en) * | 2006-12-13 | 2008-06-19 | Stephen Charles Pickering | Time-Lapse Seismic Acquisition |
| CN101663597A (en) * | 2007-02-13 | 2010-03-03 | 兰得马克制图公司 | Performing 3-D scatterer imaging from 2-D seismic data |
| CN105510967A (en) * | 2015-12-30 | 2016-04-20 | 中国石油天然气集团公司 | Processing method and device for achieving three-dimensional seismic data migration |
| CN107783187A (en) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | It is a kind of that logging speed and seismic velocity are combined to the method for establishing 3D velocity field |
| CN107703542A (en) * | 2017-08-31 | 2018-02-16 | 中国石油天然气集团公司 | A kind of determination method and apparatus of seismic horizon |
| CN108254784A (en) * | 2017-12-11 | 2018-07-06 | 中国石油天然气股份有限公司 | Fault identification method, device and system based on two-dimensional seismic data |
Non-Patent Citations (5)
| Title |
|---|
| EDUARDO TRINCHERO等: "3D seismic processing and interpretation from 2D seismic data: Application in environmentally sensitive areas of the Neuquén Basin, Argentina", THE LEADING EDGE, 31 July 2014 (2014-07-31) * |
| JASON CHIN-SEN KAO;李土杰;: "二维地震测网形成三维数据体的自动内插", 勘探地球物理进展, no. 01, 31 December 1991 (1991-12-31) * |
| 刘学清;魏福吉;王大全;秦俐;白玉春;: "二维高分辨率反演技术在苏里格气田的应用", 石油地球物理勘探, no. 01, 15 February 2013 (2013-02-15) * |
| 王鹏;钟建华;张宝权;张宏伟;安鹏;王玉珊;: "虚拟三维地震技术在绥滨拗陷油气勘探中的应用", 石油地球物理勘探, no. 1, 30 November 2011 (2011-11-30) * |
| 田连玉等: "诺木洪北虚拟三维地震技术的实现及应用", 中国石油大学胜利学院学报, 31 March 2016 (2016-03-31), pages 5 - 7 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113267815A (en) * | 2021-07-07 | 2021-08-17 | 中海油田服务股份有限公司 | Filtering method and device for repeated broken edge data |
| CN113267815B (en) * | 2021-07-07 | 2022-05-10 | 中海油田服务股份有限公司 | Method and device for filtering repeated broken edge data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112305585B (en) | 2024-08-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1397703B1 (en) | Automated system for modeling faulted multi-valued horizons | |
| US10795053B2 (en) | Systems and methods of multi-scale meshing for geologic time modeling | |
| EP3293552B1 (en) | System and method for editing geological models by switching between volume-based models and surface-based structural models augmented with stratigraphic fiber bundles | |
| US20130223187A1 (en) | Geological Structure Contour Modeling and Imaging | |
| CA2724589C (en) | Distribution of properties in a 3d volumetric model using a maximum continuity field | |
| AU2011339017A1 (en) | Constructing geologic models from geologic concepts | |
| Wuestefeld et al. | Benchmarking earthquake location algorithms: A synthetic comparison | |
| CN110501757B (en) | Point-line-surface three-dimensional point-by-point well depth design method | |
| Taranchuk | The integrated computer complex of an estimation and adapting of digital geological models | |
| CN112305585B (en) | Three-dimensional data volume acquisition method and device | |
| RU2011148308A (en) | METHOD FOR COMPREHENSIVE PROCESSING OF GEOPHYSICAL DATA AND TECHNOLOGICAL SYSTEM "LITOSCAN" FOR ITS IMPLEMENTATION | |
| WO2024064126A1 (en) | Geologic modeling framework | |
| CN107945271B (en) | Three-dimensional pressure field modeling method based on geological block tracking | |
| CN107229074A (en) | The variable speed mapping method and apparatus in low exploration activity field marine facies exploratory area | |
| CN109188522A (en) | Velocity field construction method and device | |
| Rüger et al. | Meshing for velocity modeling and ray tracing in complex velocity fields | |
| Huang et al. | Reservoir characterization by integration of seismic and dynamic data | |
| CN111815769A (en) | Modeling method, computing device and storage medium for thrust-driven tectonic belt structure | |
| Lavenu et al. | Modeling Fractures in a Heterogeneous Carbonate Reservoir Onshore UAE: A Case Study | |
| Wood et al. | Model creation based on digital borehole records and interpreted geological cross‐sections | |
| Gawith et al. | Integrating geoscience and engineering for improved field management and appraisal | |
| CN110954956B (en) | Method for evaluating acquisition trace of observation system and computer-readable storage medium | |
| Ouenes et al. | Integrated property and fracture modeling using 2D seismic data: Application to an Algerian Cambrian field | |
| Hammam et al. | Novel Dual Porosity Simulation Workflow and Solutions for the 3D Fracture Modeling | |
| CN106556864B (en) | The formation of wave equation physical efficiency amount and observation system corresponding relationship and localization method |
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 |