Construction method, device and equipment of forward modeling geophysical parameter model and medium

Abstract

The invention provides a construction method, device and equipment of a forward modeling geophysical parameter model and a medium, and belongs to the technical field of oil and gas exploration and development. According to the method, on the basis of an established complex structure model, based on the complex structure model, seismic traces and well trajectories in a model range are subjected to segmented parameterization and distributed to corresponding stratigraphic units, and a required forward geophysical parameter model is constructed based on seismic trace and well trajectory information interpolation of segmented parameterization. According to the method, the geophysical parameter distribution in the stratigraphic unit can be finely described only by explaining the main speed change interface layer without explaining a large number of horizon positions for modeling, and meanwhile, the section model does not need to be simplified, so that the construction efficiency of the forward modeling geophysical parameter model in a complex structure region is improved.

Description

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 readable storage medium. Background Art

[0002] Three-dimensional geological models are widely used in oil exploration and development. In seismic exploration, they can be used to design acquisition observation systems, test and verify new processing methods and technologies, build inversion initial models, and establish velocity models. In oil field development, they can be used for static and dynamic simulations, design development plans, and tap the potential of remaining oil.

[0003] The main purpose of the 3D seismic forward numerical simulation study for the dual-complex exploration area is to study the excitation and propagation mechanism of seismic waves in the dual-complex area and provide a data basis for the study of acquisition, processing and interpretation methods. This study mainly includes the establishment of a typical seismic geophysical parameter model and the high-precision forward geophysical parameter simulation. The geological model can be divided into acoustic wave velocity model, elastic wave velocity and density model and multi-parameter model according to different simulation equations.

[0004] The establishment of a 3D seismic forward modeling geophysical parameter model based on an undulating surface mainly includes the establishment of a near-surface geophysical parameter model and a deep underground geophysical parameter model. The near-surface geophysical parameter model can be obtained through tomographic inversion and corrected by micro-well logging. The underground geophysical parameter model usually requires the establishment of a structural model first, and then the establishment of a geophysical parameter distribution model based on the structural model.

[0005] The establishment of the seismic forward geophysical parameter model of 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 complex spatial faults is the basic guarantee for the forward geophysical parameter model. On the basis of the structural model, the velocity information such as offset velocity, VSP and logging can be used for interpolation to establish a geophysical parameter model with fine spatial changes. 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., the velocity interpolation becomes difficult.

[0006] The traditional method of establishing a forward underground geophysical parameter model for complex structures is to establish a stratigraphic block model consisting of a large number of layers and faults, and then fill in constants or vertical gradients in the fault blocks, or interpolate and fill in the lateral direction. However, this method has defects. When there are a large number of non-permeable faults that cannot be divided into blocks, the cross-section model needs to be simplified. In addition, due to the filling in the blocks, the model is rough in the vertical direction and cannot express the forward geophysical parameter model with fine spatial changes. In order to solve this problem, it is usually necessary to interpret a large number of layers for stratigraphic modeling 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. Summary of the invention

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

[0009] In order to achieve the above object, an embodiment of the present invention provides a method for constructing a forward geophysical parameter model, comprising: 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.

[0010] 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: 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.

[0011] Optionally, determine the depositional pattern of each stratum, 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.

[0012] 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 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.

[0013] 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: 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.

[0014] Optionally, the position information of each sampling point in the seismic trace corresponding to each stratum is set in sequence from top to bottom with a zero starting value, 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.

[0015] 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: 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.

[0016] Optionally, 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.

[0017] Optionally, determining the position information of the sampling point to be interpolated based on the sedimentation 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: 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.

[0018] In a second aspect of the embodiments of the present invention, a device for constructing a forward geophysical parameter model is provided, comprising: 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 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.

[0019] In a third aspect of the embodiments of the present invention, there is provided an electronic device, 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.

[0020] In a fourth aspect of the embodiments of the present invention, a computer-readable storage medium is provided, storing computer instructions, which, when executed on a computer, enable the computer to execute the above-mentioned method for constructing a forward geophysical parameter model.

[0021] 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 scope of the three-dimensional structural model are obtained; wherein each seismic trace and each well trajectory include 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; the three-dimensional structural model is divided into multiple strata, and the sedimentation mode of each stratum is determined; the multiple seismic traces and the multiple well trajectories are segmented according to the multiple strata, and each seismic trace and each well trajectory corresponding to each stratum are obtained, 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 track corresponding to each stratum to determine the overlapping sampling points as known sampling points, and the non-overlapping sampling points are determined 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 track 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 finely describe the distribution of geophysical parameters inside the stratum unit, and at the same time, the cross-section model does not need to be simplified, thereby improving the construction efficiency of the forward geophysical parameter model in complex structural areas.

[0022] 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

[0023] 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 specific implementations, they are used to explain the embodiments of the present invention, but do not constitute a limitation on the embodiments of the present invention. In the accompanying drawings: Figure 1 It is a flow chart of a method for constructing a forward geophysical parameter model provided by an embodiment of the present invention; Figure 2 is a schematic diagram of a three-dimensional structural model provided by an embodiment of the present invention; Figure 3 is a schematic diagram of a sedimentation pattern of a stratum provided by an embodiment of the present invention; Figure 4 is a schematic diagram of segmented parameterization of seismic traces and well trajectories provided by an embodiment of the present invention; Figure 5 is a schematic diagram of searching for geophysical parameters of sampling points to be interpolated provided by an embodiment of the present invention; Figure 6 is a schematic diagram of a forward geophysical parameter model provided by an embodiment of the present invention; Figure 7 is a schematic diagram of a forward geophysical parameter model including a special geological body provided by an embodiment of the present invention; Figure 8 It is a schematic diagram of a complex structural model of a double complex exploration area in Tarim provided by an embodiment of the present invention; Fig. 9 It is the sedimentation pattern of each stratum in the complex structural model of a double complex exploration area in Tarim provided by an embodiment of the present invention; Fig.10 It is a P-wave velocity model of a complex overthrust structure area in Tarim provided by an embodiment of the present invention; Fig.11 This is the mask processing effect of a high-speed conglomerate fan in a complex overthrust nappe structure area in Tarim provided by an embodiment of the present invention; Fig.12 It is a structural schematic diagram of a device for constructing a forward geophysical parameter model provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0024] The specific implementation of the embodiment of the present invention is described in detail below in conjunction with the accompanying drawings. It should be understood that the specific implementation described here 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.

[0025] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by technicians in the technical field 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.

[0026] In the description of the embodiments of the present application, the technical terms "first", "second", etc. are only used to distinguish different objects, and cannot be understood as indicating or implying relative importance or implicitly indicating the number, specific order or primary and secondary relationship of the indicated technical features. In the description of the embodiments of the present application, the meaning of "multiple" is more than two, unless otherwise clearly and specifically defined.

[0027] Embodiment 1 Please refer to Figure 1 , Figure 1 1 is a flow chart of a method for constructing a forward geophysical parameter model provided by an embodiment of the present invention, the method comprising the following steps: S100, obtaining a three-dimensional structural model based on sections and planes, and a plurality of seismic traces and a plurality of well trajectories within the scope of the three-dimensional structural model; wherein each seismic trace and each well trajectory includes a plurality of sampling points; and each well trajectory includes a plurality of well logging curves corresponding to the position information of different sampling points and the corresponding relationship between different geophysical parameter information; A 3D structural model is a virtual model that uses computer technology to model and simulate the internal structure and geological processes of the earth. It can describe and explain the internal structure and evolution of the earth by collecting geological data, seismic data, surface morphology and other geological information, combining geological principles and mathematical models, such as Figure 2 shown.

[0028] 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.

[0029] A seismic trace is a curve chart 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. Seismic traces are usually curve charts with time as the horizontal axis and seismic wave amplitude or energy as the vertical axis. Figure 4 It can be seen that the seismic traces are perpendicular to the stratum direction, and the seismic traces can be longitudinal or transverse. It is understandable 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, one containing (nx, ny, nz) seismic traces, and each (ix, iy) position corresponds to a seismic trace containing nz points.

[0030] Well trajectory refers to the path of drilling activities in oil and gas exploration and development. In the process of oil and gas exploration and development, in order to reach the target oil and gas layer, it is necessary to obtain underground information by drilling wells. The well trajectory describes the path and direction of the drilling wellbore underground, from Figure 5 It can be seen that the well trajectory can be straight or inclined.

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

[0032] It can be understood that the well trajectory provides the location and direction information of the wellbore, including parameters such as well depth, well inclination and azimuth, while the well logging curve measures the formation in the wellbore through the logging tool. In other words, the well trajectory includes the well logging curve, and the well trajectory includes the various well logging curves mentioned above.

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

[0034] In one embodiment, the well trajectory is a well trajectory of VSP or sonic logging after normalization.

[0035] S200, dividing the three-dimensional structural model into a plurality of strata, and determining a sedimentation pattern of each stratum; Multiple strata refer to dividing the three-dimensional structural model into multiple levels. Figure 2 As shown, it can be clearly seen that there are many layers in the three-dimensional structural model, and each layer represents a stratum.

[0036] Sedimentation pattern refers to the distribution and arrangement of sediments on the earth's surface or at the bottom of a water body. It describes the changing patterns of sediments in time and space, and reflects the characteristics and evolution of the sedimentary environment, including but not limited to: proportional pattern, parallel top pattern, parallel bottom pattern and block pattern, such as Figure 3 shown.

[0037] In one embodiment, the sedimentation mode of each stratum may be determined by: For each stratum: 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; The proportional pattern means that the sediments were deposited in a continuous, flat manner, with no obvious discontinuity or erosion between the top and the ground below. This conformity indicates that the sediments did not undergo significant erosion or rearrangement during deposition. Instead, a flat, continuous interface was formed between the top of the sediment and the ground below.

[0038] S220, if the top surface of the stratum is an erosion surface, determining that the stratum is a parallel bottom mode; The parallel bottom model refers to the presence of an obvious erosion surface between the top of the stratum and the ground below during the sedimentation process. This erosion surface indicates that after the formation of the stratum, an erosion event occurred, causing the top of the stratum to be eroded, forming a bottom parallel to the ground.

[0039] S230, if the ground surface of the stratum is an erosion surface, determining that the stratum is a parallel top mode; The parallel top pattern refers to the presence of an obvious erosion surface between the bottom of the stratum and the ground above during the sedimentation process. This erosion surface indicates that an erosion event occurred before the formation of the stratum, causing the bottom of the stratum to be eroded and forming a top parallel to the ground. The erosion surface may be caused by factors such as erosion and sea level drop.

[0040] S240: If the attributes of the stratum have no obvious changing pattern in the spatial direction, it is determined that the stratum is in a block mode.

[0041] The block pattern refers to the fact that the sediment remains relatively consistent in the vertical direction during the deposition process, that is, the properties, composition or characteristics of the strata do not change significantly in the vertical direction. This uniformity indicates that the sediment was not significantly affected by environmental changes or changes in depositional conditions during the deposition process.

[0042] 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; 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 achieve the consistency of the lateral trend changes of the geophysical parameter model under three-dimensional structural constraints with the trend of the small layer sedimentation pattern within the stratigraphic unit and to represent the fine changes of geophysical parameters, it is necessary to determine the strata of the sedimentation 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.

[0043] Specifically, step S300 can be implemented in the following manner: 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; Specifically, step S310 includes the following steps: For each seismic trace: S311, determining a plurality of intersection points between a plurality of strata and seismic traces; like Figure 4 As shown, seismic trace 1 intersects with layer1, layer2, layer3 and layer4 at three points.

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

[0045] S313, dividing the seismic trace into a plurality of seismic traces based on every two intersection points after the sorting process; 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.

[0046] S314, matching the multiple seismic traces with the multiple strata to obtain a seismic trace corresponding to each stratum; Specifically, after dividing the seismic trace segments corresponding to each stratum, select the midpoint coordinates of each segment to determine which stratum the midpoint belongs to. The determination method is to emit lines upward and downward along the point, calculate the intersection of the rays with 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 corresponding strata.

[0047] S315, starting from zero, sequentially setting 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.

[0048] Specifically, the top and bottom endpoints of the seismic trace corresponding to each stratum are parameterized. After all the seismic traces corresponding to each stratum are determined, the distance values ​​from the top and bottom endpoints of these segments to the top and bottom interfaces of the stratum units need to be determined.

[0049] 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.

[0050] Specifically, Figure 4 As shown in the figure, 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.

[0051] 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.

[0052] Specifically, step S320 includes the following steps: For each well trajectory: 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; 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 upward and downward outside the top and bottom interfaces of the model. This processing is to ensure that the well trajectory and the three-dimensional structural model are well intersected, and at the same time, subsequent interpolation can ensure that sufficient data can be searched for the top and bottom stratigraphic units.

[0053] 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: (1) If it is a vertical well, then the intersection point is the point where the well position intersects with the seismic trace, because at this time the well trajectory coincides with the seismic grid trace segment at that location, such as Figure 4 As shown; (2) If it is an inclined well, it is necessary to directly calculate the intersection points of the inclined well trajectory with the planes and sections 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 with the triangular units of the spatial planes and sections. 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 further judged 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, which significantly improves the intersection efficiency.

[0054] S322, based on a quick sorting method, sorting the multiple intersection points from top to bottom to obtain a plurality of sorted intersection points; Specifically, the intersection sorting and segmentation, after the intersections are calculated, is similar to the intersection sorting on the seismic trace. Since the vertical and inclined well trajectories are monotonous, they are sorted from top to bottom based on the quick sort method.

[0055] S323, dividing the seismic trace into a plurality of well trajectories based on every two intersection points after the sorting process; 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, such as Figure 4 shown.

[0056] S324, matching the multiple well trajectories with the multiple strata to obtain the well trajectory corresponding to each stratum; S325, starting from zero, sequentially setting the position information of each sampling point in the well trajectory corresponding to each stratum from top to bottom; wherein the distances between adjacent sampling points in the well trajectory corresponding to each stratum are the same.

[0057] Specifically, after the well trajectory is segmented, parameter processing is required, which is divided into two types: vertical well processing and inclined well processing: (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 stratum can be used as the location information of the corresponding well trajectory; (2) For inclined wells, the location information of the sampling points in the seismic trace cannot be used, and a parameterization method similar to that of the seismic trace segment is required to parameterize the location information of each sampling point in the well trajectory.

[0058] 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; It is understandable that each seismic trace and each well trajectory includes multiple sampling points. Among them, since the number of seismic traces is much greater than the number of well trajectories, some sampling points in the seismic traces will not overlap with the sampling points in the well trajectory, and these non-overlapping sampling points are determined as sampling points to be interpolated, and the sampling points in each well trajectory corresponding to each stratum are known sampling points. Because they do not have geophysical parameter information, they need to be inserted with geophysical parameter information to ensure that each sampling point in the three-dimensional structural model has geophysical parameter information, so that a forward geophysical parameter model with higher accuracy can be constructed.

[0059] 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; It is understandable that the number of well logging curves is the number of geophysical parameter information for each sampling point in each well trajectory corresponding to each stratum, and the number of forward geophysical parameter models constructed subsequently is the number of forward geophysical parameter models constructed subsequently. For example: the well logging curve is a longitudinal wave curve and a shear wave curve, then the geophysical parameter information for each sampling point in each well trajectory corresponding to each stratum is a shear wave and a longitudinal wave, and the forward geophysical parameter models constructed subsequently are a forward shear wave parameter model and a forward longitudinal wave parameter model.

[0060] S600, obtaining 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; Specifically, step S600 may further include the following steps: S610, determining the position information of the sampling point to be interpolated based on the sedimentation 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; Specifically, step S610 includes: S611, if the stratum where the sampling point to be interpolated is located is in proportional mode, 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, obtain the position information of the sampling point to be interpolated; 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 values ​​of the seismic trace and the top corresponding to the stratum where the sampling point to be interpolated is located and the distance parameter values ​​of the top and bottom; S612, if the stratum where the sampling point to be interpolated is located is in 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; Specifically, if the stratum is in a parallel bottom mode, a distance value is obtained by subtracting the Z value of the bottom end point of the seismic trace from the Z value of the sampling point to be interpolated, and then the small layer parameter value of the sampling point to be interpolated is obtained by adding the distance parameter value of the bottom end point; wherein, 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.

[0061] 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.

[0062] Specifically, if the stratum is in a parallel-top mode, a distance value is obtained by subtracting the Z value of the sampling point to be interpolated from the Z value of the top point of the seismic trace, and then the small layer parameter value of the Z value of the sampling point to be interpolated is obtained by adding the distance parameter value of the top point; wherein 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.

[0063] S620, 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; Specifically, the position 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; wherein, for each well curve segment, its top and bottom endpoint parameter values ​​are used to determine whether it contains the position 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 between the sampling points of two adjacent well curves based on the dichotomy method, and then linearly interpolate the geophysical parameter information of the control point; S630, based on the geostatistical Kriging algorithm, the geophysical parameter information of the sampling point to be interpolated under each well logging curve is calculated to determine the geophysical parameter information of the sampling point to be interpolated.

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

[0065] To facilitate understanding of step S600, an example is given below: Reference Figure 5 , Figure 5 The sample point in is the sampling point to be interpolated, and 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, and then match the position information of the sample point with the three curves of WellA, WellB and WellC, and 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.

[0066] 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.

[0067] Specifically, in order to achieve parallel interpolation and reduce memory usage, multi-threaded interpolation is used to interpolate all seismic traces in 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 as follows: Figure 6 shown.

[0068] In one embodiment, if a special geological body exists 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: Step 1: Calculate the geophysical parameters of sampling points of special geological bodies in the three-dimensional structural model composed of sections and layers; Special geological bodies include but are not limited to: multi-value gypsum bodies, igneous rock bodies, high-speed conglomerate fans and other velocity anomaly structures.

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

[0070] 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; Step 3: Replace 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; 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.

[0071] 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 sampling point to be replaced after replacement. The constructed forward geophysical parameter model is as follows: Figure 7 shown.

[0072] 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 scope of the three-dimensional structural model are obtained; wherein each seismic trace and each well trajectory include 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; the three-dimensional structural model is divided into multiple strata, and the sedimentation mode of each stratum is determined; the multiple seismic traces and the multiple well trajectories are segmented according to the multiple strata, and each seismic trace and each well trajectory corresponding to each stratum are obtained, 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 track corresponding to each stratum to determine the overlapping sampling points as known sampling points, and the non-overlapping sampling points are determined 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 track 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 finely describe the distribution of geophysical parameters inside the stratum unit, and at the same time, the cross-section model does not need to be simplified, thereby improving the construction efficiency of the forward geophysical parameter model in complex structural areas.

[0073] Embodiment 2 The method for constructing the forward geophysical parameter model provided by the present invention is applied to a typical double-complex exploration area in the Tarim Basin. The main feature of this area is the complex surface velocity structure and the highly developed underground reverse thrust faults, which leads to a complex underground velocity structure. How to establish a fine underground forward velocity model that can represent the double-complex exploration area under the constraints of the structural framework model is one of the main technical challenges of typical forward simulation. In order to establish a fine forward model, a total of 22 reverse faults were interpreted. Since the interpretation of reverse faults needs to be divided into upper and lower plate interpretations, a total of 85 strata were interpreted in pieces. At the same time, VSP logging and acoustic logging, pre-stack depth migration velocity model and post-stack seismic data volume information of the exploration area were collected for velocity model establishment. The structural model established based on the above-mentioned fault and stratum data is as follows: Figure 8 In order to establish the velocity model, the VSP and acoustic logging data of 15 wells in the area were reconstructed and processed using post-stack seismic data and migration velocity models based on artificial intelligence curve reconstruction technology before modeling. The sedimentary pattern of the stratigraphic unit is shown in Fig. 9 The underground velocity model established based on the above data and the technology of the present invention is shown as follows. Fig.10 and Fig.11 and as shown.

[0074] Embodiment 3 Please refer to Fig.12 , Fig.12 It is a structural schematic diagram of a device for constructing a forward geophysical parameter model provided in an embodiment of the present application.

[0075] 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: The data acquisition module 210 is used to acquire a three-dimensional structural model composed of 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; 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; A model stratification module 220 is used to divide the three-dimensional structural model into multiple strata and determine the sedimentation pattern of each stratum; The model parameterization module 230 is used to perform segmentation processing on 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; 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; 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; The second parameter determination module 260 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 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.

[0076] It should be understood that the device corresponds to the above-mentioned forward geophysical parameter model construction method embodiment, and can execute the various steps involved in the above-mentioned method embodiment. The specific functions of the device can be found in the above description. To avoid repetition, the detailed description is appropriately omitted here. The device includes at least one software function module that can be stored in a memory in the form of software or firmware or solidified in the operating system (OS) of the device.

[0077] Embodiment 4 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, when executed by the processor, execute the above-mentioned method for constructing a forward geophysical parameter model.

[0078] 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.

[0079] Computer readable media include permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. 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 disk read-only memory (CD-ROM), digital versatile disk (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 temporary computer readable media (transitory media), such as modulated data signals and carrier waves.

[0080] Embodiment 5 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.

[0081] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may 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 may adopt the form of a computer program product implemented in one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) that include computer-usable program code.

[0082] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented 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 processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0083] These computer program instructions may also be stored in a computer-readable memory capable of directing 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 an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0084] 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, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

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

[0086] 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 separately, or two or more modules can be integrated to form an independent part.

[0087] It should also be noted that the terms "include", "comprises" or any other variations thereof are intended to cover non-exclusive inclusion, so that a process, method, commodity or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, commodity or device. In the absence of more restrictions, the elements defined by the sentence "comprises a ..." do not exclude the existence of other identical elements in the process, method, commodity or device including the elements.

[0088] The above are only 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 modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application should 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