Method and apparatus for constructing forward geophysical parameter model, and device and medium

Abstract

A method for constructing a forward geophysical parameter model, which method belongs to the technical field of oil and gas exploration and development. On the basis of an established complex structural model, seismic traces and well trajectories within the range of the model are segmented and parameterized on the basis of the complex structural model, and are distributed into corresponding stratum units, interpolation is performed on the basis of segmented and parameterized seismic trace and well trajectory information, and a required forward geophysical parameter model is constructed. A geophysical parameter distribution inside a stratum unit can be finely described simply by means of explaining main speed change interface layers instead of explaining a large number of horizons to perform modeling, and it is not necessary to simplify a section model, thereby improving the construction efficiency of a forward geophysical parameter model for a complex structural region. Further provided are an apparatus for constructing a forward geophysical parameter model, and a device and a medium.

Description

Method, device, equipment and medium for constructing forward geophysical parameter model Technical Field

[0001] The present invention relates to the technical field of oil and gas exploration and development, and in particular to a method for constructing a forward geophysical parameter model, a device for constructing a forward geophysical parameter model, an electronic device, and a computer-readable storage medium. Background Art

[0002] 3D geological models are widely used in oil exploration and development. In seismic exploration, they can be used to design acquisition and observation systems, test and validate new processing methods and technologies, construct initial inversion models, and establish velocity models. In oilfield development, they can be used for static and dynamic simulations, development plan design, and remaining oil potential.

[0003] The primary purpose of 3D seismic forward numerical simulation research in dual-complexity exploration areas is to investigate the excitation and propagation mechanisms of seismic waves in these areas, providing a data foundation for research into acquisition, processing, and interpretation methods. This research primarily involves establishing typical seismic geophysical parameter models and performing high-precision forward geophysical parameter simulations. Geological models can be categorized into acoustic wave velocity models, elastic wave velocity and density models, and multi-parameter models, depending on the simulation equations.

[0004] Establishing a 3D seismic forward geophysical parameter model based on an undulating surface primarily involves building a near-surface geophysical parameter model and a deep subsurface geophysical parameter model. The near-surface geophysical parameter model can be obtained through tomographic inversion and calibrated using micro-well logging. However, developing a subsurface geophysical parameter model typically requires first establishing a structural model, and then building a geophysical parameter distribution model based on the structural model.

[0005] The establishment of the seismic forward geophysical parameter model for the underground part includes the establishment of a three-dimensional structural model and a geophysical parameter model. The structural model is the foundation, especially for the underground overthrust and thrust structures in the double-complex exploration area. The establishment of stratigraphic units that can accurately represent the complex spatial faults is the basic guarantee for the forward geophysical parameter model. On the basis of the structural model, the offset velocity, VSP and well logging velocity information can be used for interpolation to establish a geophysical parameter model with fine spatial variations. The geophysical parameter model is usually represented and stored in the format of post-stack seismic data. Due to the very complex underground structure of the double-complex exploration area, the structural intersection, the development of overthrust and thrust faults, the high dip angle of the strata, the deep strata and the high-speed conglomerate bodies, etc., velocity interpolation becomes difficult.

[0006] The traditional method for constructing forward geophysical parameter models for complex structures involves creating a stratigraphic block model consisting of a large number of planes and faults. Constants or vertical gradients are then filled in within the blocks, or interpolated values ​​are used laterally. However, this approach has drawbacks. When a large number of impermeable faults exist that cannot be divided into blocks, the cross-sectional model must be simplified. Furthermore, the in-block filling results in a coarse vertical model, making it impossible to express the spatially detailed forward geophysical parameter model. To address this issue, stratigraphic modeling typically requires interpreting a large number of planes to improve the vertical refinement of the model.

[0007] However, doing so will make the modeling process very time-consuming and labor-intensive, resulting in low efficiency in forward geophysical parameter modeling.

[0008] Summary of the Invention

[0009] The purpose of the embodiments of the present invention is to provide a method, apparatus, device and medium for constructing a forward geophysical parameter model to solve the problem of low efficiency in forward geophysical parameter modeling.

[0010] To achieve the above objectives, an embodiment of the present invention provides a method for constructing a forward geophysical parameter model, comprising:

[0011] Acquire a three-dimensional structural model based on sections and planes, as well as multiple seismic traces and multiple well trajectories within the three-dimensional structural model; wherein each seismic trace and each well trajectory includes multiple sampling points; and each well trajectory includes multiple well logging curves corresponding to the position information of different sampling points and different geophysical parameter information;

[0012] Divide the 3D structural model into multiple strata and determine the sedimentary pattern of each stratum;

[0013] Segmenting multiple seismic traces and multiple well trajectories according to multiple strata to obtain each seismic trace and each well trajectory corresponding to each stratum, and calculating the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum;

[0014] Overlapping the seismic trace corresponding to each stratum with the well trajectory corresponding to each stratum, so as to determine the overlapping sampling points as known sampling points and determine the non-overlapping sampling points as sampling points to be interpolated;

[0015] Match the known sampling points with the corresponding well logging curves to obtain the geophysical parameter information of each sampling point in each well trajectory corresponding to each formation;

[0016] Based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each logging curve, the geophysical parameter information of the sampling point to be interpolated is obtained;

[0017] Based on the geophysical parameter information of known sampling points and the geophysical parameter information of the sampling points to be interpolated, a forward geophysical parameter model is constructed.

[0018] Optionally, before the step of constructing a forward geophysical parameter model based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated, the method further includes:

[0019] Calculate geophysical parameters of sampling points of specific geological bodies in a three-dimensional structural model composed of sections and layers;

[0020] The sampling points among the known sampling points and the sampling points to be interpolated that coincide with the sampling points in the special geological body are determined as the sampling points to be replaced;

[0021] Replacing the geophysical parameters of the sampling point to be replaced with the geophysical parameters of the sampling point in the special geological body to obtain the replaced geophysical parameters of the sampling point to be replaced;

[0022] A forward geophysical parameter model is constructed based on the geophysical parameters of the sampling points other than the sampling point to be replaced and the geophysical parameters of the replaced sampling point.

[0023] Optionally, determine the depositional pattern of each formation, including:

[0024] For each stratum:

[0025] If the top surface of the stratum and the ground surface are in an integrated relationship, the stratum is determined to be in a proportional mode;

[0026] If the top surface of the stratum is an erosion surface, the stratum is determined to be in a parallel bottom mode;

[0027] If the ground surface of the formation is an erosion surface, the formation is determined to be a parallel top model;

[0028] If the properties of the stratum have no obvious variation pattern in the spatial direction, the stratum is determined to be a blocky pattern.

[0029] Optionally, multiple seismic traces and multiple well trajectories are segmented according to multiple strata to obtain each seismic trace and each well trajectory corresponding to each stratum, and position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum is calculated, including:

[0030] Segmenting multiple seismic traces according to multiple strata to obtain each seismic trace corresponding to each stratum, and calculating the position information of each sampling point in each seismic trace corresponding to each stratum;

[0031] The multiple well trajectories are segmented according to the multiple strata to obtain each well trajectory corresponding to each stratum, and the position information of each sampling point in each well trajectory corresponding to each stratum is calculated.

[0032] Optionally, multiple seismic traces are segmented according to multiple strata to obtain each seismic trace corresponding to each stratum, and position information of each sampling point in each seismic trace corresponding to each stratum is calculated, including:

[0033] For each seismic trace:

[0034] Determine multiple intersection points of multiple strata and seismic traces;

[0035] Based on the quick sorting method, multiple intersection points are sorted from top to bottom to obtain multiple sorted intersection points;

[0036] The seismic trace is divided into multiple seismic traces based on every two intersection points after sorting;

[0037] Match multiple seismic traces with multiple strata to obtain the seismic trace corresponding to each stratum;

[0038] Starting from zero, the position information of each sampling point in the seismic trace corresponding to each stratum is set sequentially from top to bottom; wherein, the distance between adjacent sampling points in the seismic trace corresponding to each stratum is the same.

[0039] Optionally, starting from zero, the position information of each sampling point in the seismic trace corresponding to each stratum is sequentially set from top to bottom, including:

[0040] When a certain sampling point in the seismic trace corresponding to each stratum is the intersection of the fault and the seismic trace, the position information of the intersection of the fault and the seismic trace is reset based on the fault distance to obtain the position information of the intersection of the first fault and the seismic trace and the position information of the intersection of the second fault and the seismic trace. The position information of the intersection of the first fault and the seismic trace is used as the position information of the final sampling point of the previous stratum, and the position information of the intersection of the second fault and the seismic trace is used as the position information of the initial sampling point of the next stratum.

[0041] Optionally, multiple well trajectories are segmented according to multiple strata to obtain each well trajectory corresponding to each stratum, and the position information of each sampling point in each well trajectory corresponding to each stratum is calculated, including:

[0042] For each well trajectory:

[0043] Extending multiple well trajectories outside the three-dimensional structural model to determine multiple intersection points of multiple formations with the well trajectories;

[0044] Based on the quick sorting method, multiple intersection points are sorted from top to bottom to obtain multiple sorted intersection points;

[0045] The seismic trace is divided into multiple well trajectories based on every two intersection points after sorting;

[0046] Match multiple well trajectories with multiple strata to obtain the well trajectory corresponding to each stratum;

[0047] Starting from zero, the position information of each sampling point in the well trajectory corresponding to each stratum is sequentially set from top to bottom; wherein, the distance between adjacent sampling points in the well trajectory corresponding to each stratum is the same.

[0048] Optionally, geophysical parameter information of the sampling point to be interpolated is obtained based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each well logging curve, including:

[0049] Determine the position information of the sampling point to be interpolated based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located and the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located;

[0050] Matching the position information of the sampling points to be interpolated with each well logging curve to determine the geophysical parameter information of the sampling points to be interpolated under each well logging curve;

[0051] Based on the geostatistical Kriging algorithm, the geophysical parameter information of the sampling point to be interpolated under each logging curve is calculated to determine the geophysical parameter information of the sampling point to be interpolated.

[0052] Optionally, determining the position information of the sampling point to be interpolated based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located and the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located includes:

[0053] If the stratum where the sampling point to be interpolated is located is in proportional mode, the position information of the sampling point to be interpolated is obtained based on the position information of the sampling point to be interpolated, the position information of the initial sampling point of the stratum where the sampling point to be interpolated is located, and the position information of the final sampling point of the stratum where the sampling point to be interpolated is located;

[0054] If the stratum where the sampling point to be interpolated is located is in a parallel bottom mode, determining the position information of the sampling point to be interpolated based on the position information of the sampling point to be interpolated and the position information of the last-level sampling point of the stratum where the sampling point to be interpolated is located;

[0055] If the stratum where the sampling point to be interpolated is located is in parallel top mode, the position information of the sampling point to be interpolated is determined based on the position information of the sampling point to be interpolated and the position information of the initial sampling point of the stratum where the sampling point to be interpolated is located.

[0056] In a second aspect of an embodiment of the present invention, a device for constructing a forward geophysical parameter model is provided, comprising:

[0057] A data acquisition module is used to acquire a three-dimensional structural model based on sections and planes, as well as multiple seismic traces and multiple well trajectories within the scope of the three-dimensional structural model; wherein each seismic trace and each well trajectory includes multiple sampling points; each well trajectory corresponds to multiple well logging curves corresponding to the position information of different sampling points and different geophysical parameter information;

[0058] Model layering module, used to divide the 3D structural model into multiple strata and determine the sedimentary pattern of each stratum;

[0059] The model parameterization module is used to segment multiple seismic traces and multiple well trajectories according to multiple strata, obtain each seismic trace and each well trajectory corresponding to each stratum, and calculate the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum;

[0060] A sampling point classification module is used to overlap the seismic trace corresponding to each stratum with the well trajectory corresponding to each stratum, so as to determine the overlapping sampling points as known sampling points and determine the non-overlapping sampling points as sampling points to be interpolated;

[0061] The first parameter determination module is used to match the known sampling points with the corresponding well logging curves to obtain the geophysical parameter information of each sampling point in each well trajectory corresponding to each formation;

[0062] The second parameter determination module is used to obtain geophysical parameter information of the sampling point to be interpolated based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each well logging curve;

[0063] The model building module is used to build a forward geophysical parameter model based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated.

[0064] In a third aspect of the embodiments of the present invention, an electronic device is provided, comprising: a processor and a memory, wherein the memory stores machine-readable instructions executable by the processor, and when the machine-readable instructions are executed by the processor, the method for constructing a forward geophysical parameter model is performed.

[0065] In a fourth aspect of the embodiments of the present invention, a computer-readable storage medium is provided, storing computer instructions. When the computer instructions are executed on a computer, the computer executes the above-mentioned method for constructing a forward geophysical parameter model.

[0066] In an embodiment of the present invention, a three-dimensional structural model based on sections and planes, and multiple seismic traces and multiple well trajectories within the three-dimensional structural model are obtained; each seismic trace and each well trajectory includes multiple sampling points; and each well trajectory corresponds to multiple well logging curves representing the correspondence between position information of different sampling points and different geophysical parameter information. The three-dimensional structural model is divided into multiple strata, and the sedimentary pattern of each stratum is determined. The multiple seismic traces and multiple well trajectories are segmented according to the multiple strata to obtain each seismic trace and each well trajectory corresponding to each stratum, and the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum is calculated. The seismic trace corresponding to each stratum is superimposed with the well trajectory corresponding to each stratum to determine the overlapping sampling points as known sampling points and the non-overlapping sampling points as sampling points to be interpolated; the known sampling points are matched with the corresponding well logging curves to obtain the geophysical parameter information of each sampling point in each well trajectory corresponding to each stratum; based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each well logging curve, the geophysical parameter information of the sampling point to be interpolated is obtained; based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated, a forward geophysical parameter model is constructed. The present invention does not need to explain a large number of horizons for modeling, but only needs to explain the main velocity change interface layer to accurately describe the distribution of geophysical parameters within the stratum unit, and at the same time, the cross-section model does not need to be simplified, thereby improving the efficiency of constructing the forward geophysical parameter model in complex tectonic areas.

[0067] Other features and advantages of the embodiments of the present invention will be described in detail in the subsequent detailed description. BRIEF DESCRIPTION OF THE DRAWINGS

[0068] The accompanying drawings are used to provide a further understanding of the embodiments of the present invention and constitute a part of the specification. Together with the following detailed description, they are used to explain the embodiments of the present invention but do not constitute a limitation of the embodiments of the present invention. In the accompanying drawings:

[0069] FIG1 is a schematic flow chart of a method for constructing a forward geophysical parameter model provided by an embodiment of the present invention;

[0070] FIG2 is a schematic diagram of a three-dimensional structural model provided by an embodiment of the present invention;

[0071] FIG3 is a schematic diagram of a sedimentation pattern of a stratum provided by an embodiment of the present invention;

[0072] FIG4 is a schematic diagram of segmented parameterization of seismic traces and well trajectories provided by an embodiment of the present invention;

[0073] FIG5 is a schematic diagram of searching for geophysical parameters of sampling points to be interpolated according to an embodiment of the present invention;

[0074] FIG6 is a schematic diagram of a forward geophysical parameter model provided by an embodiment of the present invention;

[0075] FIG7 is a schematic diagram of a forward geophysical parameter model including a special geological body provided by an embodiment of the present invention;

[0076] FIG8 is a schematic diagram of a complex structural model of a dual-complex exploration area in Tarim provided by an embodiment of the present invention;

[0077] FIG9 is a sedimentary pattern of each stratum in a complex structural model of a dual-complex exploration area in Tarim provided by an embodiment of the present invention;

[0078] FIG10 is a P-wave velocity model of a complex overthrust structure area in Tarim provided by an embodiment of the present invention;

[0079] FIG11 is a mask processing effect of a high-speed conglomerate fan in a complex overthrust structure area in Tarim provided by an embodiment of the present invention;

[0080] FIG12 is a schematic structural diagram of a device for constructing a forward geophysical parameter model provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0081] The following describes the specific implementation of the embodiment of the present invention in detail with reference to the accompanying drawings. It should be understood that the specific implementation described herein is only used to illustrate and explain the embodiment of the present invention and is not used to limit the embodiment of the present invention.

[0082] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by those skilled in the art to which this application belongs. The terms used herein are only for the purpose of describing specific embodiments and are not intended to limit this application.

[0083] In the description of the embodiments of this application, the technical terms "first," "second," etc. are used only to distinguish different objects and should not be understood to indicate or imply relative importance or to implicitly indicate the quantity, specific order, or primary and secondary relationship of the indicated technical features. In the description of the embodiments of this application, the meaning of "plurality" is more than two, unless otherwise specifically defined.

[0084] Example 1

[0085] Please refer to FIG1 , which is a flow chart of a method for constructing a forward geophysical parameter model provided by an embodiment of the present invention. The method includes the following steps:

[0086] S100, obtaining a three-dimensional structural model based on sections and planes, and multiple seismic traces and multiple well trajectories within the three-dimensional structural model; wherein each seismic trace and each well trajectory includes multiple sampling points; and each well trajectory includes multiple well logging curves corresponding to position information of different sampling points and different geophysical parameter information;

[0087] A 3D tectonic model is a virtual model used by geologists to model and simulate the Earth's internal structure and geological processes using computer technology. It uses geological data, seismic data, surface topography, and other geological information, combined with geological principles and mathematical models, to describe and explain the Earth's internal structure and evolutionary processes, as shown in Figure 2.

[0088] In one embodiment, in order to improve the intersection processing between the seismic traces and these spatial triangular surfaces, the planes and sections in the three-dimensional structural model are all triangular mesh unit surfaces.

[0089] A seismic trace refers to a graph that records the propagation process of seismic waves in seismic exploration. In seismic exploration, seismic energy propagates through underground media, and seismic instruments convert the received seismic signals into seismic traces. A seismic trace is usually a graph with time as the horizontal axis and seismic wave amplitude or energy as the vertical axis. As can be seen from Figure 4, the seismic traces are all perpendicular to the direction of the stratum, and the seismic traces can be longitudinal or transverse. It can be understood that the seismic trace is a three-dimensional forward model data body, that is, a set of regular seismic grids are defined according to the grid scale of dx, dy, dz, a seismic trace containing (nx, ny, nz), and each (ix, iy) position corresponds to a seismic trace containing various points of nz.

[0090] A well trajectory refers to the path of drilling activities in oil and gas exploration and development. In order to reach the target oil and gas formations, drilling wells is necessary to obtain subsurface information. The well trajectory describes the path and direction of the wellbore underground. As shown in Figure 5, a well trajectory can be either straight or inclined.

[0091] Well logging curves refer to curves obtained by measuring formation properties with well depth in the wellbore using logging tools, thereby revealing the lithology, oil and gas properties, etc. of the formation, including but not limited to: resistivity curves, natural gamma radiation curves, compressional wave velocity curves, shear wave velocity curves, density curves, VTI curves and TTI curves.

[0092] It can be understood that the well trajectory provides information on the location and orientation of the wellbore, including parameters such as well depth, well inclination, and azimuth, while the well log curve measures the formation within the wellbore using logging tools. In other words, the well trajectory includes the well log curve, and when the well trajectory includes the various well log curves mentioned above.

[0093] It should be noted that there can be many seismic traces, but relatively fewer well trajectories.

[0094] In one embodiment, the well trajectory is a normalized well trajectory of VSP or sonic logging.

[0095] S200, dividing the three-dimensional structural model into multiple strata and determining the sedimentary pattern of each stratum;

[0096] Multiple strata refer to dividing the 3D structural model into multiple levels. As shown in Figure 2, it can be clearly seen that there are many layers in the 3D structural model, and each layer represents a stratum.

[0097] Sedimentation patterns refer to the distribution and arrangement of sediments on the Earth's surface or at the bottom of a body of water. They describe how sediments change over time and space, reflecting the characteristics and evolution of the sedimentary environment. These patterns include, but are not limited to, proportional patterns, parallel-top patterns, parallel-bottom patterns, and massive patterns, as shown in Figure 3.

[0098] In one embodiment, the sedimentation pattern of each stratum may be determined by:

[0099] For each stratum:

[0100] S210, if the top surface of the stratum and the ground surface are in an integrated relationship, determining that the stratum is in a proportional mode;

[0101] The proportional pattern describes sediment deposition in a continuous, flat pattern, with no apparent discontinuity or erosion between the top and the underlying ground. This conformity suggests that the sediments did not undergo significant erosion or rearrangement during deposition. Instead, a flat, continuous interface forms between the top of the sediment and the underlying ground.

[0102] S220, if the top surface of the stratum is an erosion surface, determining that the stratum is in a parallel bottom mode;

[0103] The parallel-bottom pattern refers to a situation where a distinct erosion surface exists between the top of a stratum and the underlying ground during sediment deposition. This erosion surface indicates that an erosion event occurred after the formation of the stratum, causing the top of the stratum to be eroded, forming a bottom parallel to the ground.

[0104] S230, if the ground surface of the stratum is an erosion surface, determining that the stratum is a parallel top mode;

[0105] A parallel-top pattern occurs when a distinct erosion surface exists between the bottom of a formation and the surface above it during sediment deposition. This surface suggests that an erosion event occurred before the formation, eroding the bottom of the formation and creating a top parallel to the surface. This erosion surface may be caused by factors such as erosion and sea level drop.

[0106] S240: If the attributes of the stratum have no obvious change pattern in the spatial direction, the stratum is determined to be a block pattern.

[0107] Massive patterns refer to the relatively consistent vertical properties of sediments during their deposition. This means that the properties, composition, or characteristics of the formations do not change significantly in the vertical direction. This uniformity indicates that the sediments were not significantly affected by environmental changes or changes in depositional conditions during their deposition.

[0108] S300, segmenting the plurality of seismic traces and the plurality of well trajectories according to the plurality of strata to obtain each seismic trace and each well trajectory corresponding to each stratum, and calculating the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum;

[0109] It should be noted that in order to overcome the need for a large amount of explicit layer modeling to refine the geophysical parameter model in traditional methods, and to ensure that the lateral trend changes of the geophysical parameter model under three-dimensional structural constraints are consistent with the sedimentary pattern trends of small layers within the stratigraphic unit and can represent the fine changes in geophysical parameters, it is necessary to determine the strata of the sedimentary pattern, and parameterize the seismic traces corresponding to each stratum according to the regional scope of the model to determine the position information of each sample point on the seismic trace corresponding to each stratum.

[0110] Specifically, step S300 can be implemented in the following manner:

[0111] S310, segmenting the plurality of seismic traces according to the plurality of strata to obtain each seismic trace corresponding to each stratum, and calculating the position information of each sampling point in each seismic trace corresponding to each stratum;

[0112] Specifically, step S310 includes the following steps:

[0113] For each seismic trace:

[0114] S311, determining a plurality of intersection points between a plurality of strata and seismic traces;

[0115] As shown in Figure 4, seismic trace 1 intersects with layer 1, layer 2, layer 3, and layer 4 at three points.

[0116] S312, sorting the plurality of intersection points from top to bottom based on a quick sorting method to obtain a plurality of sorted intersection points;

[0117] Specifically, after all the intersections are calculated, the intersection calculation process cannot determine the upper and lower order of all the intersections on the same seismic trace, so the seismic trace cannot be divided into seismic trace segments according to the intersections. Therefore, it is necessary to sort the intersections on each seismic trace from top to bottom based on the quick sort method.

[0118] S313, dividing the seismic trace into multiple seismic traces based on every two intersection points after the sorting process;

[0119] Specifically, the seismic traces are divided into seismic trace segments according to the order of intersection points, and the seismic trace segments are assigned to corresponding strata, and the stratum numbers to which they belong are marked.

[0120] S314, matching the multiple seismic traces with the multiple strata to obtain a seismic trace corresponding to each stratum;

[0121] Specifically, after dividing the seismic trace segments corresponding to each stratum, the midpoint coordinates of each segment are selected to determine which stratum the midpoint belongs to. The judgment method is to emit lines upward and downward along the point, calculate the intersection of the rays and all stratigraphic units, and take out the surface information of the nearest intersection to determine which stratum it belongs to. In this way, all seismic trace segments can be assigned to the corresponding stratum.

[0122] S315 , starting from zero, sequentially set the position information of each sampling point in the seismic trace corresponding to each stratum from top to bottom; wherein the distances between adjacent sampling points in the seismic trace corresponding to each stratum are the same.

[0123] Specifically, the top and bottom endpoints of the seismic trace corresponding to each stratum are parameterized. After determining all the seismic traces corresponding to each stratum, it is also necessary to determine the distance values ​​from the top and bottom endpoints of these segments to the top and bottom interfaces of the stratum units.

[0124] In one embodiment, when a sampling point in the seismic trace corresponding to each stratum is the intersection of a fault and the seismic trace, the position information of the intersection of the fault and the seismic trace is reset based on the fault distance to obtain the position information of the intersection of the first fault and the seismic trace and the position information of the intersection of the second fault and the seismic trace, and the position information of the intersection of the first fault and the seismic trace is used as the position information of the final sampling point of the previous stratum, and the position information of the intersection of the second fault and the seismic trace is used as the position information of the initial sampling point of the next stratum.

[0125] Specifically, as shown in Figure 4, if the top point is on the top interface, then its distance value is set to 0.0. If the top point is on a certain section, it is necessary to use the fault distance on the section to correct the distance value, so that the distance value is no longer 0.0, but a value between 0.0-1.0; if the bottom point is on the bottom interface, then its distance value is set to 1.0. If the top point is on a certain section, it is necessary to use the fault distance on the section to correct the distance value, so that the distance value is no longer 1.0, but a value between 0.0-1.0.

[0126] S320 , segmenting the plurality of well trajectories according to the plurality of strata to obtain each well trajectory corresponding to each stratum, and calculating the position information of each sampling point in each well trajectory corresponding to each stratum.

[0127] Specifically, step S320 includes the following steps:

[0128] For each well trajectory:

[0129] S321, extending the multiple well trajectories outside the three-dimensional structural model to determine multiple intersection points between the multiple formations and the well trajectories;

[0130] Specifically, in order to ensure that the well trajectory and curve can penetrate the top and bottom of the model, it is necessary to copy the first and last points of the well trajectory and curve and move them up and down outside the top and bottom interfaces of the model. This is done to ensure that the well trajectory and the three-dimensional structural model are well intersected, and that sufficient data can be searched for the top and bottom stratigraphic units during subsequent interpolation.

[0131] Specifically, since there are vertical wells and inclined wells in the well trajectory, the determination of multiple intersection points between multiple formations and the well trajectory includes the following two situations:

[0132] (1) If it is a vertical well, then the intersection point is directly taken as the point where the well position intersects the seismic trace, because the well trajectory coincides with the seismic grid trace segment at that location, as shown in Figure 4;

[0133] (2) If it is an inclined well, it is necessary to directly calculate the intersection of the inclined well trajectory with the plane and section of the three-dimensional structural model. At this time, since the well trajectory is a monotonous broken line in space, the entire intersection process is actually the intersection of the straight line segment on the well trajectory and the triangular unit of the spatial plane and section. In order to speed up the intersection efficiency, a multi-level bounding box detection and exclusion method is adopted, and the bounding box of the well trajectory is first used to make a judgment with the bounding box of the triangular mesh surface of the section and layer. If there is an intersection, it is continued to judge whether the bounding box of the straight line segment of the well trajectory and the bounding box of the triangular unit intersect. In this way, most of the invalid intersection calculations are eliminated, and the intersection efficiency is significantly improved.

[0134] S322, sorting the plurality of intersection points from top to bottom based on a quick sorting method to obtain a plurality of sorted intersection points;

[0135] Specifically, after the intersection points are calculated, the intersection point sorting and segmentation are similar to the intersection point sorting on the seismic trace. Since the vertical well and inclined well trajectories are monotonous, they are sorted from top to bottom based on the quick sorting method.

[0136] S323, dividing the seismic trace into a plurality of well trajectories based on every two intersection points after the sorting process;

[0137] Specifically, the well trajectory is then divided into well trajectory segments according to the order of intersection points, and each well trajectory segment belongs to a formation, as shown in FIG4 .

[0138] S324, matching the multiple well trajectories with the multiple strata to obtain the well trajectory corresponding to each stratum;

[0139] S325 , starting from zero, sequentially set the position information of each sampling point in the well trajectory corresponding to each stratum from top to bottom; wherein, the distance between adjacent sampling points in the well trajectory corresponding to each stratum is the same.

[0140] Specifically, after the well trajectory is segmented, parameter processing is required, which is divided into two types: vertical well processing and inclined well processing:

[0141] (1) For vertical wells, since the well trajectory and the seismic trace at that location coincide, the location information of the seismic trace corresponding to each formation can be used as the location information of the corresponding well trajectory;

[0142] (2) For inclined wells, the position information of the sampling points in the seismic trace cannot be used, and a parameterization processing method similar to that of the seismic trace segment is required to parameterize the position information of each sampling point in the well trajectory.

[0143] S400, superimposing the seismic trace corresponding to each stratum with the well trajectory corresponding to each stratum, so as to determine the overlapping sampling points as known sampling points, and determine the non-overlapping sampling points as sampling points to be interpolated;

[0144] It is understood that each seismic trace and each well trajectory includes multiple sampling points. Since the number of seismic traces far exceeds the number of well trajectories, some sampling points in the seismic traces may not overlap with sampling points in the well trajectories. These non-overlapping sampling points are identified as sampling points to be interpolated. The sampling points in each well trajectory corresponding to each stratum are known sampling points. Since these points do not have geophysical parameter information, geophysical parameter information must be interpolated to ensure that each sampling point in the 3D structural model has geophysical parameter information, thereby constructing a highly accurate forward geophysical parameter model.

[0145] S500, matching the known sampling points with the corresponding well logging curves to obtain geophysical parameter information of each sampling point in each well trajectory corresponding to each formation;

[0146] It's understandable that for every well logging curve, there are as many different geophysical parameter information types for each sampling point in each well trajectory corresponding to each formation, and as many forward geophysical parameter models subsequently constructed. For example, if the well logging curves are P-wave and S-wave curves, then for every sampling point in each well trajectory corresponding to each formation, there are S-wave and P-wave geophysical parameter information types, and the forward geophysical parameter models subsequently constructed are a forward S-wave parameter model and a forward P-wave parameter model.

[0147] S600, obtaining geophysical parameter information of the sampling point to be interpolated based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each well logging curve;

[0148] Specifically, step S600 may further include the following steps:

[0149] S610, determining the position information of the sampling point to be interpolated based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located and the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located;

[0150] Specifically, step S610 includes:

[0151] S611: If the stratum where the sampling point to be interpolated is located is in proportional mode, obtain the position information of the sampling point to be interpolated based on the position information of the sampling point to be interpolated, the position information of the initial sampling point of the stratum where the sampling point to be interpolated is located, and the position information of the final sampling point of the stratum where the sampling point to be interpolated is located;

[0152] Specifically, if the stratum is in proportional mode, the position parameter of the sampling point to be interpolated can be obtained by using the sum of the distance parameter value between the seismic trace and the top and the distance parameter value between the top and the bottom corresponding to the stratum where the sampling point to be interpolated is located;

[0153] S612: If the stratum where the sampling point to be interpolated is located is in a parallel bottom mode, determine the position information of the sampling point to be interpolated based on the position information of the sampling point to be interpolated and the position information of the last-level sampling point of the stratum where the sampling point to be interpolated is located;

[0154] Specifically, if the stratum is in a parallel bottom mode, the Z value of the sampling point to be interpolated is subtracted from the Z value of the bottom end point of the seismic trace to obtain a distance value, and then the distance parameter value of the bottom end point is added to obtain the small layer parameter value of the sampling point to be interpolated; among them, the Z value of the sampling point to be interpolated is the distance from the sampling point to be interpolated to the bottom surface of the three-dimensional structural model, and the Z value of the bottom end point of the seismic trace is the distance from the bottom surface of the bottom end of the seismic trace to the bottom surface of the three-dimensional structural model.

[0155] S613: If the stratum where the sampling point to be interpolated is located is in parallel top mode, determine the position information of the sampling point to be interpolated based on the position information of the sampling point to be interpolated and the position information of the initial sampling point of the stratum where the sampling point to be interpolated is located.

[0156] Specifically, if the stratum is in a parallel-top mode, the Z value of the top point of the seismic trace minus the Z value of the sampling point to be interpolated is a distance value, and then the distance parameter value of the top point is added to obtain the small layer parameter value of the Z value of the sampling point to be interpolated; among them, the Z value of the sampling point to be interpolated is the distance from the sampling point to be interpolated to the bottom surface of the three-dimensional structural model, and the Z value of the top point of the seismic trace is the distance from the top surface of the seismic trace to the bottom surface of the three-dimensional structural model.

[0157] S620, matching the position information of the sampling point to be interpolated with each well logging curve to determine the geophysical parameter information of the sampling point to be interpolated under each well logging curve;

[0158] Specifically, the location information of the sampling point to be interpolated is matched with all the well logging curves of the formation where the sampling point to be interpolated is located, and the geophysical parameter information of the sampling point to be interpolated under each well logging curve is obtained. For each well curve segment, the top and bottom endpoint parameter values ​​are used to determine whether it contains the location information of the sampling point to be interpolated. If not, it means that there is no interpolation control point on the well curve segment. If it is included, it is necessary to quickly locate the sampling points between two adjacent well curves based on the bisection method, and then linearly interpolate the geophysical parameter information of the control point.

[0159] S630 , calculating the geophysical parameter information of the sampling point to be interpolated under each well logging curve based on the geostatistical kriging algorithm to determine the geophysical parameter information of the sampling point to be interpolated.

[0160] The geostatistical kriging algorithm belongs to the existing technology and will not be described in detail here.

[0161] To facilitate understanding of step S600, an example is given below:

[0162] Referring to Figure 5, the sample point in Figure 5 is the sampling point to be interpolated. Then, the sedimentary mode of the stratum where the sample point is located is determined to be the proportional mode, so as to calculate the position information of the sample point. Then, the position information of the sample point is matched with the three curves of WellA, WellB and WellC. The geophysical parameter information of point A on WellA, point B on WellB and point C on WellC is the geophysical parameter information of the sample point under each logging curve.

[0163] S700: Construct a forward geophysical parameter model based on geophysical parameter information of known sampling points and geophysical parameter information of sampling points to be interpolated.

[0164] Specifically, to achieve parallel interpolation and reduce memory usage, multi-threaded interpolation is used to interpolate all seismic traces within a profile along the line direction. Each time a line profile is fully interpolated, it is written to a file. By repeating the above steps, the forward geophysical parameter model of the entire model area can be obtained. The constructed forward geophysical parameter model is shown in Figure 6.

[0165] In one embodiment, if a special geological body is considered in the three-dimensional structural model, it is necessary to replace the geophysical parameter information of the sampling points to be interpolated with the geophysical parameter information of the known sampling points. The specific process is as follows:

[0166] Step 1: Calculate the geophysical parameters of sampling points of special geological bodies in the three-dimensional structural model composed of sections and layers;

[0167] Special geological bodies include but are not limited to: multi-value gypsum bodies, igneous rock bodies, high-speed conglomerate fans and other velocity abnormal structures.

[0168] In this step, the method for calculating the geophysical parameters of known sampling points and the geophysical parameters of the sampling points to be interpolated in the above-described embodiment can be used to calculate the geophysical parameters of the sampling points of the specific geological body. Of course, other methods can also be used to calculate the geophysical parameters of the sampling points of the specific geological body, and this embodiment of the present invention is not specifically limited to this method.

[0169] Step 2: The sampling points among the known sampling points and the sampling points to be interpolated that coincide with the sampling points in the special geological body are determined as the sampling points to be replaced;

[0170] Step 3: replacing the geophysical parameters of the sampling point to be replaced with the geophysical parameters of the sampling point in the special geological body to obtain the geophysical parameters of the sampling point to be replaced after replacement;

[0171] Specifically, the sampling points among the known sampling points and the sampling points to be interpolated that coincide with the special geological body are determined, and then the geophysical parameters of the sampling points of the special geological body are replaced with the coincident sampling points.

[0172] The fourth step is to construct a forward geophysical parameter model based on the geophysical parameters of the sampling points other than the sampling point to be replaced and the geophysical parameters of the replaced sampling point. The constructed forward geophysical parameter model is shown in FIG7 .

[0173] In an embodiment of the present invention, a three-dimensional structural model based on sections and planes, and multiple seismic traces and multiple well trajectories within the three-dimensional structural model are obtained; each seismic trace and each well trajectory includes multiple sampling points; and each well trajectory corresponds to multiple well logging curves representing the correspondence between position information of different sampling points and different geophysical parameter information. The three-dimensional structural model is divided into multiple strata, and the sedimentary pattern of each stratum is determined. The multiple seismic traces and multiple well trajectories are segmented according to the multiple strata to obtain each seismic trace and each well trajectory corresponding to each stratum, and the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum is calculated. The seismic trace corresponding to each stratum is superimposed with the well trajectory corresponding to each stratum to determine the overlapping sampling points as known sampling points and the non-overlapping sampling points as sampling points to be interpolated; the known sampling points are matched with the corresponding well logging curves to obtain the geophysical parameter information of each sampling point in each well trajectory corresponding to each stratum; based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each well logging curve, the geophysical parameter information of the sampling point to be interpolated is obtained; based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated, a forward geophysical parameter model is constructed. The present invention does not need to explain a large number of horizons for modeling, but only needs to explain the main velocity change interface layer to accurately describe the distribution of geophysical parameters within the stratum unit, and at the same time, the cross-section model does not need to be simplified, thereby improving the efficiency of constructing the forward geophysical parameter model in complex tectonic areas.

[0174] Example 2

[0175] The forward geophysical parameter model construction method provided by the present invention was applied to a typical dual-complex exploration area in the Tarim Basin. This area is characterized by a complex surface velocity structure and a highly developed subsurface thrust-nappe fault system, resulting in a complex subsurface velocity structure. One of the main technical challenges of typical forward modeling is how to establish a detailed subsurface forward velocity model that represents this dual-complex exploration area within the constraints of a tectonic framework model. To establish a detailed forward model, 22 reverse faults were interpreted. Because reverse fault interpretation requires interpretation of both the upper and lower walls, a total of 85 horizons were interpreted in slices. VSP and sonic logging data from the exploration area, along with prestack depth migration velocity models and poststack seismic data, were collected for velocity model construction. The structural model established based on these fault and horizon data is shown in Figure 8. To establish the velocity model, VSP and sonic logging data from 15 wells in the area were reconstructed using poststack seismic data and migration velocity models using artificial intelligence curve reconstruction technology before modeling. The stratigraphic unit sedimentary patterns are shown in Figure 9. Using these data, the subsurface velocity model established using the present invention is shown in Figures 10 and 11.

[0176] Example 3

[0177] Please refer to FIG. 12 , which is a schematic diagram of the structure of a device for constructing a forward geophysical parameter model provided in an embodiment of the present application.

[0178] Based on the same inventive concept, an embodiment of the present invention further provides a device 200 for constructing a forward geophysical parameter model, comprising:

[0179] The data acquisition module 210 is used to acquire a three-dimensional structural model based on sections and planes, as well as multiple seismic traces and multiple well trajectories within the three-dimensional structural model; each seismic trace and each well trajectory includes multiple sampling points; each well trajectory includes multiple well logging curves corresponding to the position information of different sampling points and different geophysical parameter information;

[0180] The model layering module 220 is used to divide the three-dimensional structural model into multiple strata and determine the sedimentation pattern of each stratum;

[0181] The model parameterization module 230 is used to segment the multiple seismic traces and the multiple well trajectories according to the multiple strata, obtain each seismic trace and each well trajectory corresponding to each stratum, and calculate the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum;

[0182] The sampling point classification module 240 is used to overlap the seismic trace corresponding to each stratum with the well trajectory corresponding to each stratum, so as to determine the overlapping sampling points as known sampling points and determine the non-overlapping sampling points as sampling points to be interpolated;

[0183] The first parameter determination module 250 is used to match the known sampling points with the corresponding well logging curves to obtain the geophysical parameter information of each sampling point in each well trajectory corresponding to each formation;

[0184] The second parameter determination module 260 is configured to obtain geophysical parameter information of the sampling point to be interpolated based on the sedimentary pattern of the formation where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the formation where the sampling point to be interpolated is located, and each well logging curve;

[0185] The model building module 270 is used to build a forward geophysical parameter model based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated.

[0186] It should be understood that the device corresponds to the aforementioned embodiment of the method for constructing a forward geophysical parameter model and is capable of executing each of the steps involved in the aforementioned method embodiment. The specific functions of the device can be found in the description above, and a detailed description is omitted here to avoid repetition. The device includes at least one software functional module that can be stored in a memory in the form of software or firmware or embedded in the device's operating system (OS).

[0187] Example 4

[0188] Based on the same inventive concept, an embodiment of the present invention also provides an electronic device, including: a processor and a memory, the memory storing machine-readable instructions executable by the processor, and the machine-readable instructions executing the above-mentioned forward geophysical parameter model construction method when executed by the processor.

[0189] The memory may include non-permanent memory in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. The memory is an example of a computer-readable medium.

[0190] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory computer-readable media (transitory media), such as modulated data signals and carrier waves.

[0191] Example 5

[0192] Based on the same inventive concept, an embodiment of the present invention further provides a computer-readable storage medium storing computer instructions. When the computer instructions are executed on a computer, the computer executes the above-mentioned method for constructing a forward geophysical parameter model.

[0193] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to magnetic disk storage, CD-ROM, optical storage, etc.) that contain computer-usable program code.

[0194] The present application is described with reference to the flow chart and / or block diagram of the method, device (system), and computer program product according to the embodiment of the present application. It should be understood that each flow process and / or box in the flow chart and / or block diagram and the combination of the flow process and / or box in the flow chart and / or block diagram can be realized by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processing machine or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device produce a device for realizing the function specified in one flow chart flow or multiple flows and / or one box or multiple boxes of the block diagram.

[0195] These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce a product including an instruction device that implements the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

[0196] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, so that the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more processes in the flowchart and / or one or more boxes in the block diagram.

[0197] It should also be noted that the various specific technical features described in the above specific embodiments can be combined in any appropriate manner without contradiction. To avoid unnecessary repetition, the embodiments of the present invention will not further describe various possible combinations.

[0198] In addition, the functional modules in each embodiment of the present application can be integrated together to form an independent part, or each module can exist independently, or two or more modules can be integrated to form an independent part.

[0199] It should also be noted that the terms "comprises," "includes," or any other variations thereof are intended to encompass non-exclusive inclusion, such that a process, method, commodity, or apparatus that includes a series of elements includes not only those elements but also other elements not explicitly listed, or includes elements inherent to such process, method, commodity, or apparatus. In the absence of further limitations, an element defined by the phrase "comprises a ..." does not exclude the presence of other identical elements in the process, method, commodity, or apparatus that includes the element.

[0200] The above are merely embodiments of the present application and are not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application should all be included within the scope of the claims of the present application.

Claims

1. A method for constructing a forward geophysical parameter model, characterized in that: include: Acquire a three-dimensional structural model based on sections and planes, and multiple seismic traces and multiple well trajectories within the scope of the three-dimensional structural model; wherein each seismic trace and each well trajectory includes multiple sampling points; and each well trajectory includes multiple well logging curves corresponding to the position information of different sampling points and the corresponding relationship between different geophysical parameter information; Divide the three-dimensional structural model into multiple strata and determine the sedimentary pattern of each stratum; According to the multiple strata, the multiple seismic traces and the multiple well trajectories are segmented to obtain each seismic trace and each well trajectory corresponding to each stratum, and the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum is calculated; Superimposing the seismic trace corresponding to each stratum with the well trajectory corresponding to each stratum, so as to determine the overlapping sampling points as known sampling points, and determine the non-overlapping sampling points as sampling points to be interpolated; Match the known sampling points with the corresponding logging curves to obtain the geophysical parameter information of each sampling point in each well trajectory corresponding to each formation; Based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each logging curve, the geophysical parameter information of the sampling point to be interpolated is obtained; Based on the geophysical parameter information of known sampling points and the geophysical parameter information of sampling points to be interpolated, a forward geophysical parameter model is constructed.

2. The method for constructing a forward geophysical parameter model according to claim 1, characterized in that: Before the step of constructing a forward geophysical parameter model based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated, the method further includes: Calculate geophysical parameters of sampling points of special geological bodies in three-dimensional structural models composed of sections and layers; The sampling points among the known sampling points and the sampling points to be interpolated that coincide with the sampling points in the special geological body are determined as the sampling points to be replaced; Replacing the geophysical parameters of the sampling point to be replaced with the geophysical parameters of the sampling point in the special geological body to obtain the replaced geophysical parameters of the sampling point to be replaced; A forward geophysical parameter model is constructed based on the geophysical parameters of the sampling points other than the sampling point to be replaced and the geophysical parameters of the replaced sampling point to be replaced.

3. The method for constructing a forward geophysical parameter model according to claim 1, characterized in that: Determine the depositional pattern of each formation, including: For each stratum: If the top surface of the stratum and the ground surface are in an integrated relationship, the stratum is determined to be in a proportional mode; If the top surface of the stratum is an erosion surface, the stratum is determined to be in the parallel bottom mode; If the ground surface of the formation is an erosion surface, the formation is determined to be in parallel top mode; If the properties of the stratum have no obvious change pattern in the spatial direction, the stratum is determined to be a block pattern.

4. The method for constructing a forward geophysical parameter model according to claim 1, characterized in that: According to multiple strata, multiple seismic traces and multiple well trajectories are segmented to obtain each seismic trace and each well trajectory corresponding to each stratum, and the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum is calculated, including: The multiple seismic traces are segmented according to the multiple strata to obtain each seismic trace corresponding to each stratum, and the position information of each sampling point in each seismic trace corresponding to each stratum is calculated; The multiple well trajectories are segmented according to the multiple strata to obtain each well trajectory corresponding to each stratum, and the position information of each sampling point in each well trajectory corresponding to each stratum is calculated.

5. The method for constructing a forward geophysical parameter model according to claim 4, characterized in that: According to multiple strata, multiple seismic traces are segmented to obtain each seismic trace corresponding to each stratum, and the position information of each sampling point in each seismic trace corresponding to each stratum is calculated, including: For each seismic trace: Determine multiple intersection points of multiple strata and seismic traces; Based on the quick sorting method, multiple intersection points are sorted from top to bottom to obtain multiple intersection points after sorting; The seismic trace is divided into multiple seismic traces based on every two intersection points after sorting; Match multiple seismic traces with multiple strata to obtain the seismic trace corresponding to each stratum; Starting from zero, the position information of each sampling point in the seismic trace corresponding to each stratum is set sequentially from top to bottom; wherein the distances between adjacent sampling points in the seismic trace corresponding to each stratum are the same.

6. The method for constructing a forward geophysical parameter model according to claim 5, characterized in that: Starting from zero, set the location information of each sampling point in the seismic trace corresponding to each stratum from top to bottom, including: When a sampling point in the seismic trace corresponding to each stratum is the intersection of the fault and the seismic trace, the position information of the intersection of the fault and the seismic trace is reset based on the fault distance to obtain the position information of the intersection of the first fault and the seismic trace and the position information of the intersection of the second fault and the seismic trace, and the position information of the intersection of the first fault and the seismic trace is used as the position information of the final sampling point of the previous stratum, and the position information of the intersection of the second fault and the seismic trace is used as the position information of the initial sampling point of the next stratum.

7. The method for constructing a forward geophysical parameter model according to claim 4, characterized in that: The multiple well trajectories are segmented according to the multiple strata to obtain each well trajectory corresponding to each stratum, and the position information of each sampling point in each well trajectory corresponding to each stratum is calculated, including: For each well trajectory: Extending the multiple well trajectories outside the three-dimensional structural model to determine multiple intersection points of multiple formations with the well trajectories; Based on the quick sorting method, multiple intersection points are sorted from top to bottom to obtain multiple intersection points after sorting; The seismic trace is divided into multiple well tracks based on every two intersection points after sorting; Matching multiple well trajectories with multiple strata to obtain a well trajectory corresponding to each stratum; Starting from zero, the position information of each sampling point in the well trajectory corresponding to each stratum is set in sequence from top to bottom; wherein the distances between adjacent sampling points in the well trajectory corresponding to each stratum are the same.

8. The method for constructing a forward geophysical parameter model according to claim 1, characterized in that: Based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and each logging curve, the geophysical parameter information of the sampling point to be interpolated is obtained, including: Determine the position information of the sampling point to be interpolated based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located and the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located; Matching the position information of the sampling points to be interpolated with each well logging curve to determine the geophysical parameter information of the sampling points to be interpolated under each well logging curve; Based on the geostatistical Kriging algorithm, the geophysical parameter information of the sampling points to be interpolated under each well logging curve is calculated to determine the geophysical parameter information of the sampling points to be interpolated.

9. The method for constructing a forward geophysical parameter model according to claim 8, characterized in that: Based on the sedimentary pattern of the stratum where the sampling point to be interpolated is located and the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling point to be interpolated is located, the position information of the sampling point to be interpolated is determined, including: If the stratum where the sampling point to be interpolated is located is in proportional mode, the position information of the sampling point to be interpolated is obtained based on the position information of the sampling point to be interpolated, the position information of the initial sampling point of the stratum where the sampling point to be interpolated is located, and the position information of the final sampling point of the stratum where the sampling point to be interpolated is located; If the stratum where the sampling point to be interpolated is located is in parallel bottom mode, based on the position information of the sampling point to be interpolated and the position information of the last-level sampling point of the stratum where the sampling point to be interpolated is located, determine the position information of the sampling point to be interpolated; If the stratum where the sampling point to be interpolated is located is in parallel top mode, the position information of the sampling point to be interpolated is determined based on the position information of the sampling point to be interpolated and the position information of the initial sampling point of the stratum where the sampling point to be interpolated is located.

10. A device for constructing a forward geophysical parameter model, characterized in that: include: The data acquisition module is used to acquire a three-dimensional structural model based on sections and planes, as well as multiple seismic traces and multiple well trajectories within the scope of the three-dimensional structural model; wherein each seismic trace and each well trajectory includes multiple sampling points; each well trajectory corresponds to multiple well logging curves of the corresponding relationship between the position information of different sampling points and different geophysical parameter information Model stratification module, used to divide the 3D structural model into multiple strata and determine the sedimentation pattern of each stratum; The model parameterization module is used to process multiple seismic traces and multiple well trajectories in sections according to multiple strata, obtain each seismic trace and each well trajectory corresponding to each stratum, and calculate the position information of each sampling point in each seismic trace and each well trajectory corresponding to each stratum; A sampling point classification module is used to overlap the seismic trace corresponding to each stratum with the well trajectory corresponding to each stratum, so as to determine the overlapping sampling points as known sampling points and determine the non-overlapping sampling points as sampling points to be interpolated; The first parameter determination module is used to match the known sampling points with the corresponding logging curves to obtain the geophysical parameter information of each sampling point in each well trajectory corresponding to each formation; The second parameter determination module is used to obtain geophysical parameter information of the sampling point to be interpolated based on the sedimentary pattern of the formation where the sampling point to be interpolated is located, the position information of each sampling point in the seismic trace corresponding to the formation where the sampling point to be interpolated is located, and each logging curve; The model building module is used to build a forward geophysical parameter model based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated.

11. An electronic device, characterized in that: include: A processor and a memory, wherein the memory stores machine-readable instructions executable by the processor, and when the machine-readable instructions are executed by the processor, the method for constructing a forward geophysical parameter model described in any one of claims 1 to 9 is executed.

12. A computer-readable storage medium storing computer instructions, characterized in that: When the computer instructions are executed on a computer, the computer is enabled to execute the method for constructing a forward geophysical parameter model according to any one of claims 1 to 9.

Images