Forward geophysical parameter model construction method and device, equipment and medium

Abstract

The application provides a method, device and equipment for constructing a forward geophysical parameter model and a medium, and belongs to the technical field of oil and gas exploration and development. The method is based on a complex structure model and a model range, and is used for segmentally parameterizing and distributing a seismic trace and a well track in the model range to a corresponding stratum unit, and interpolating and constructing a required forward geophysical parameter model based on the segmentally parameterized seismic trace and well track information. The application can finely describe the internal geophysical parameter distribution of the stratum unit by only interpreting a main velocity variation interface layer without interpreting a large number of horizons for modeling, and the section model does not need to be simplified, so that the construction efficiency of the forward geophysical parameter model in a complex structure area is improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of oil and gas exploration and development, 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

[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 techniques, construct inversion initial models, and establish velocity models. In oilfield development, they can be used for static and dynamic simulation, design development plans, and tap remaining oil.

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

[0004] Establishing a three-dimensional seismic forward geophysical parameter model based on a rolling surface mainly includes establishing a near-surface geophysical parameter model and a deep underground geophysical parameter model. The near-surface geophysical parameter model can be obtained by tomographic inversion and corrected by micro-logging. The underground geophysical parameter model usually needs to establish a structure model first, and then establish a geophysical parameter distribution model according to the structure model.

[0005] The establishment of the seismic forward geophysical parameter model of the underground part includes the establishment of a three-dimensional structure model and a geophysical parameter model. The structure model is the foundation, especially for the underground overthrust nappe structure of the double complex exploration area. Establishing a stratum unit that can accurately represent the spatial complex fault is the basic guarantee of the forward geophysical parameter model. On the basis of the structure model, the spatial fine change geophysical parameter model can be established by using the velocity information of migration velocity, VSP and logging. The geophysical parameter model is usually represented and stored based on the post-stack seismic data body format. Due to the complex underground structure of the double complex exploration area, the structure is crossed, overthrust nappe and reverse fault is developed, the stratum is developed with high dip angle, the deep stratum and high-speed conglomerate body and other complex conditions, and the velocity interpolation becomes difficult.

[0006] The conventional complex structure forward underground geophysical parameter model building method is to build a stratigraphic fault block model composed of a large number of layer surfaces and faults, and then fill in constants or longitudinal gradients in the fault blocks or interpolate in the horizontal direction. However, this method has defects. When a large number of non-penetrating faults cannot be blocked, the section model needs to be simplified. In addition, due to the filling in the block, the model is longitudinally rough and cannot express the spatially fine forward geophysical parameter model. In order to solve this problem, a large number of layer surfaces are usually needed to be interpreted for stratigraphic modeling to improve the longitudinal fineness of the model.

[0007] However, this will cause the modeling process to be very time-consuming and laborious, resulting in low efficiency of the forward geophysical parameter model modeling. SUMMARY

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

[0009] In order to achieve the above-mentioned purpose, the embodiments of the present application provide a forward geophysical parameter model building method, comprising:

[0010] Obtaining a three-dimensional structure model composed of sections and layer surfaces, and a plurality of seismic traces and a plurality of well tracks within the range of the three-dimensional structure model; wherein each seismic trace and each well track comprises a plurality of sampling points; each well track corresponds to a plurality of logging curves of the corresponding relationship between the position information of different sampling points and different geophysical parameter information;

[0011] Dividing the three-dimensional structure model into a plurality of strata, and determining the sedimentary mode of each stratum;

[0012] Segmenting the plurality of seismic traces and the plurality of well tracks according to the plurality of strata respectively to obtain each seismic trace and each well track corresponding to each stratum, and calculating the position information of each sampling point in each seismic trace and each well track corresponding to each stratum;

[0013] Overlapping each seismic trace corresponding to each stratum with each well track corresponding to each stratum to determine the overlapping sampling points as known sampling points and the non-overlapping sampling points as interpolation sampling points;

[0014] Matching the known sampling points with the corresponding logging curves to obtain the geophysical parameter information of each sampling point in each well track corresponding to each stratum;

[0015] Based on the sedimentary mode of the stratum where the interpolation sampling point is located, the position information of each sampling point in the seismic trace corresponding to the stratum where the interpolation sampling point is located, and the logging curves, the geophysical parameter information of the interpolation sampling point is obtained;

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

[0017] Optionally, before the step of constructing the 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 comprises:

[0018] calculating the geophysical parameter of the sampling point of the special geologic body in the three-dimensional structure model composed of the section and the layer surface;

[0019] determining the sampling points coinciding with the sampling points in the special geologic body in the known sampling points and the sampling points to be interpolated as the sampling points to be replaced;

[0020] replacing the geophysical parameter of the sampling points to be replaced with the geophysical parameter of the sampling points in the special geologic body to obtain the geophysical parameter of the sampling points to be replaced after replacement;

[0021] constructing the forward geophysical parameter model based on the geophysical parameter of the sampling points except the sampling points to be replaced and the geophysical parameter of the sampling points to be replaced after replacement.

[0022] Optionally, the deposition mode of each layer is determined, comprising:

[0023] for each layer:

[0024] if the top surface of the layer and the ground surface are in conformity, the layer is determined as the proportional mode;

[0025] if the top surface of the layer is an erosion surface, the layer is determined as the parallel bottom mode;

[0026] if the ground surface of the layer is an erosion surface, the layer is determined as the parallel top mode;

[0027] if the attribute of the layer has no obvious change rule in the spatial direction, the layer is determined as the massive mode.

[0028] Optionally, the plurality of seismic traces and the plurality of well tracks are respectively segmented according to the plurality of layers to obtain each seismic trace and each well track corresponding to each layer, and the position information of each sampling point in each seismic trace and each well track corresponding to each layer is calculated, comprising:

[0029] the plurality of seismic traces are respectively segmented according to the plurality of layers to obtain each seismic trace corresponding to each layer, and the position information of each sampling point in each seismic trace corresponding to each layer is calculated;

[0030] According to the multiple strata, the multiple well tracks are segmented to obtain each well track corresponding to each stratum, and position information of each sampling point in each well track corresponding to each stratum is calculated.

[0031] Optionally, according to the multiple strata, the multiple seismic traces are segmented 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, comprising:

[0032] For each seismic trace:

[0033] Multiple intersection points of the multiple strata and the seismic trace are determined;

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

[0035] The seismic trace is divided into multiple seismic traces based on each two sorted intersection points;

[0036] The multiple seismic traces are matched with the multiple strata to obtain seismic traces corresponding to each stratum;

[0037] From top to bottom, position information of each sampling point in the seismic trace corresponding to each stratum is sequentially set with a zero starting value; wherein distances between adjacent sampling points in the seismic trace corresponding to each stratum are the same.

[0038] Optionally, from top to bottom, position information of each sampling point in the seismic trace corresponding to each stratum is sequentially set with a zero starting value, comprising:

[0039] When a certain sampling point in the seismic trace corresponding to each stratum is an intersection point of a fault and the seismic trace, the position information of the intersection point of the fault and the seismic trace is reset based on a fault throw to obtain position information of a first intersection point of the fault and the seismic trace and position information of a second intersection point of the fault and the seismic trace, and the position information of the first intersection point of the fault and the seismic trace is taken as position information of a last sampling point of a previous stratum, and the position information of the second intersection point of the fault and the seismic trace is taken as position information of an initial sampling point of a next stratum.

[0040] Optionally, according to the multiple strata, the multiple well tracks are segmented to obtain each well track corresponding to each stratum, and position information of each sampling point in each well track corresponding to each stratum is calculated, comprising:

[0041] For each well track:

[0042] The multiple well tracks are extended to outside of the three-dimensional structure model to determine multiple intersection points of the multiple strata and the well tracks;

[0043] The multiple intersection points are sequentially sorted from top to bottom based on a quick sorting method to obtain the multiple sorted intersection points;

[0044] The seismic traces are divided into multiple well traces based on each two sorted intersection points as a division basis;

[0045] The multiple well traces are matched with the multiple strata to obtain a well trace corresponding to each stratum;

[0046] The position information of each sampling point in the well trace corresponding to each stratum is sequentially set from top to bottom based on a zero starting value; wherein the distance between adjacent sampling points in the well trace corresponding to each stratum is the same.

[0047] Optionally, the geophysical parameter information of the sampling point to be interpolated is obtained based on the sedimentary mode 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 the well curves, including:

[0048] The position information of the sampling point to be interpolated is determined based on the sedimentary mode 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;

[0049] The geophysical parameter information of the sampling point to be interpolated under the well curves is determined by matching the position information of the sampling point to be interpolated with the well curves;

[0050] The geophysical parameter information of the sampling point to be interpolated is calculated based on a geostatistical Kriging algorithm to determine the geophysical parameter information of the sampling point to be interpolated.

[0051] Optionally, the position information of the sampling point to be interpolated is determined based on the sedimentary mode 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, including:

[0052] If the stratum where the sampling point to be interpolated is located is a 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;

[0053] If the stratum where the sampling point to be interpolated is located is a parallel bottom 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 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 a 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.

[0055] In a second aspect of the embodiments of the present application, a device for constructing a forward geophysical parameter model is provided, comprising:

[0056] a data acquisition module configured to acquire a three-dimensional structure model composed of sections and layers, and a plurality of seismic traces and a plurality of well tracks within the range of the three-dimensional structure model; wherein each seismic trace and each well track comprises a plurality of sampling points; and each well track corresponds to a plurality of logging curves of the correspondence between the position information of different sampling points and different geophysical parameter information;

[0057] a model layering module configured to divide the three-dimensional structure model into a plurality of strata, and determine the sedimentary pattern of each stratum;

[0058] a model parameterization module configured to segment process the plurality of seismic traces and the plurality of well tracks according to the plurality of strata respectively, to obtain each seismic trace and each well track corresponding to each stratum, and to calculate the position information of each sampling point in each seismic trace and each well track corresponding to each stratum;

[0059] a sampling point classification module configured to superimpose the seismic trace corresponding to each stratum and the well track corresponding to each stratum, to determine the overlapping sampling points as known sampling points, and to determine the non-overlapping sampling points as interpolation sampling points;

[0060] a first parameter determination module configured to match the known sampling points with the corresponding logging curves, to obtain the geophysical parameter information of each sampling point in each well track corresponding to each stratum;

[0061] a second parameter determination module configured to obtain the geophysical parameter information of the interpolation sampling points based on the sedimentary pattern of the stratum where the interpolation sampling points are located, the position information of each sampling point in the seismic trace corresponding to the stratum where the interpolation sampling points are located, and the logging curves;

[0062] a model construction module configured to construct a forward geophysical parameter model based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the interpolation sampling points.

[0063] In a third aspect of the embodiments of the present application, an electronic device is provided, comprising a processor and a memory, wherein the memory stores machine readable instructions executable by the processor, and the machine readable instructions are executed by the processor to perform the above-mentioned method for constructing a forward geophysical parameter model.

[0064] In a fourth aspect of the embodiments of the present application, a computer readable storage medium is provided, which stores computer instructions, and when the computer instructions are executed on a computer, the computer is caused to perform the above-mentioned method for constructing a forward geophysical parameter model.

[0065] In the embodiment of the present application, a three-dimensional structure model based on a section and a layer is obtained, and a plurality of seismic traces and a plurality of well tracks in the range of the three-dimensional structure model are obtained; each seismic trace and each well track comprises a plurality of sampling points; each well track corresponds to a plurality of logging curves of the corresponding relationship between the position information of different sampling points and different geophysical parameter information; the three-dimensional structure model is divided into a plurality of strata, and the deposition mode of each stratum is determined; the plurality of seismic traces and the plurality of well tracks are segmented and processed according to the plurality of strata, to obtain each seismic trace and each well track corresponding to each stratum, and the position information of each sampling point in each seismic trace and each well track 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 as interpolation sampling points; the geophysical parameter information of each sampling point in each well track corresponding to each stratum is obtained by matching the known sampling points with the corresponding logging curves; the geophysical parameter information of the interpolation sampling points is obtained based on the deposition mode of the stratum where the interpolation sampling points are located, the position information of each sampling point in the seismic trace corresponding to the stratum where the interpolation sampling points are located, and the logging curves; and a forward geophysical parameter model is constructed based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the interpolation sampling points. The present application can finely describe the distribution of geophysical parameters in the internal stratum unit by interpreting only the main velocity variation interface layer without interpreting a large number of horizon layers for modeling, and the section model does not need to be simplified, thereby improving the construction efficiency of the forward geophysical parameter model in a complex structure region.

[0066] Other features and advantages of the embodiments of the present application will be described in detail in the following specific embodiments. BRIEF DESCRIPTION OF DRAWINGS

[0067] The accompanying drawings are included to provide a further understanding of the embodiments of the present application, and constitute a part of the specification, and are used to explain the embodiments of the present application together with the following specific embodiments, but do not constitute a limitation of the embodiments of the present application. In the drawings:

[0068] Figure 1 FIG. 1 is a flowchart of the method for constructing a forward geophysical parameter model provided by the embodiments of the present application;

[0069] Figure 2 FIG. 3 is a schematic diagram of a three-dimensional structure model provided by the embodiments of the present application;

[0070] Figure 3 FIG. 5 is a schematic diagram of a deposition mode of a stratum provided by the embodiments of the present application;

[0071] Figure 4 FIG. 7 is a schematic diagram of the segmented parameterization of a seismic trace and a well track provided by the embodiments of the present application;

[0072] Figure 5 is a schematic diagram of searching for a geophysical parameter of a to-be-interpolated sampling point provided by an embodiment of the present application;

[0073] Figure 6 is a schematic diagram of a forward geophysical parameter model provided by an embodiment of the present application;

[0074] Figure 7 is a schematic diagram of a forward geophysical parameter model containing a special geological body provided by an embodiment of the present application;

[0075] Figure 8 is a schematic diagram of a complex structure model of a certain double complex exploration area in Tarim provided by an embodiment of the present application;

[0076] Figure 9 is a schematic diagram of a sedimentary model of each stratum in the complex structure model of the certain double complex exploration area in Tarim provided by an embodiment of the present application;

[0077] Figure 10 is a P-wave velocity model of a certain complex overthrust nappe structure area in Tarim provided by an embodiment of the present application;

[0078] Figure 11 is a high-speed conglomerate fan body mask processing effect of the certain complex overthrust nappe structure area in Tarim provided by an embodiment of the present application;

[0079] Figure 12 is a structural schematic diagram of a construction device of a forward geophysical parameter model provided by an embodiment of the present application. DETAILED DESCRIPTION

[0080] The specific embodiments of the embodiments of the present application are described in detail below with reference to the accompanying drawings. It should be understood that the specific embodiments described herein are only for the purpose of describing and explaining the embodiments of the present application, and are not intended to limit the embodiments of the present application.

[0081] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this application belongs; the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the present application.

[0082] In the description of the embodiments of the present application, the technical terms "first", "second", and the like 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 "a plurality of" is two or more, unless otherwise specifically limited.

[0083] Embodiment one

[0084] Please refer to Figure 1 ,Figure 1 is a flowchart of a method for constructing a forward geophysical parameter model according to an embodiment of the present application. The method comprises the following steps:

[0085] S100, obtaining a three-dimensional structure model composed of sections and layers, and a plurality of seismic traces and a plurality of well tracks within the range of the three-dimensional structure model; wherein each seismic trace and each well track comprises a plurality of sampling points; each well track corresponds to a plurality of logging curves of the correspondence relationship between the position information of different sampling points and different geophysical parameter information;

[0086] The three-dimensional structure model refers to a virtual model of the internal structure and geological process of the earth modeled and simulated by geology using computer technology. It can describe and explain the structure and evolution of the earth's interior by collecting geological data, seismic data, surface topography and other geological information, combining with geological principles and mathematical models, such as Figure 2 .

[0087] In an embodiment, in order to improve the intersection processing of seismic traces and these spatial triangular faces, the layer and section in the three-dimensional structure model are triangular mesh element faces.

[0088] The seismic trace refers to a curve graph recording 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. The seismic trace is usually a curve 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 trace is perpendicular to the direction of the stratum, and the seismic trace can be longitudinal or transverse. Understandably, 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 kind of seismic trace containing (nx, ny, nz), each (ix, iy) position corresponds to a seismic trace containing nz samples.

[0089] The well track 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 through drilling wells. The well track describes the path and direction of the wellbore in the ground, as can be seen from Figure 5 , the well track can be straight or inclined.

[0090] The logging curve refers to the curve of the variation of the properties of the stratum with the depth of the wellbore measured by the logging tool in the wellbore, so as to reveal the lithology, oil and gas properties, etc. of the stratum, including but not limited to: such as resistivity curve, natural gamma radiation curve, P-wave velocity curve, S-wave velocity curve, density curve, VTI curve and TTI curve.

[0091] It can be understood that the well trajectory provides the position and direction information of the wellbore, including parameters such as well depth, well inclination and azimuth, and the well log is measured by the logging tool on the formation in the wellbore. That is, the well trajectory includes the well log, and the well trajectory includes the above-mentioned various well logs.

[0092] It should be noted that there are many seismic traces, but relatively few well trajectories.

[0093] In an embodiment, the well trajectory is a normalized VSP or acoustic logging well trajectory.

[0094] S200, dividing the three-dimensional structure model into multiple strata, and determining the sedimentary mode of each stratum;

[0095] The multiple strata refer to the division of the three-dimensional structure model into multiple levels. As shown in Figure 2 It can be clearly seen that there are many layers in the three-dimensional structure model, and each layer represents a stratum.

[0096] The sedimentary mode refers to the distribution and arrangement of sediments on the earth's surface or the bottom of a water body. It describes the variation of sediments in time and space, reflects the characteristics and evolution process of the sedimentary environment, including but not limited to: proportional mode, parallel top mode, parallel bottom mode and block mode, as shown in Figure 3 .

[0097] In an embodiment, the determination of the sedimentary mode of each stratum can be realized by the following way:

[0098] For each stratum:

[0099] S210, if the top surface of the stratum and the ground surface are in conformable relationship, determining that the stratum is in proportional mode;

[0100] The proportional mode refers to the deposition of sediments in a continuous and flat manner, and there is no obvious discontinuity or erosion phenomenon between the top and the underlying ground. This conformable relationship indicates that the sediments have not been subjected to significant erosion or rearrangement during the deposition process. On the contrary, a flat and continuous interface is formed between the top of the sediments and the underlying ground.

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

[0102] The parallel bottom mode refers to the existence of a clear erosion surface between the top of the stratum and the underlying ground during the deposition of sediments. This erosion surface indicates that an erosion event has occurred after the formation of the stratum, resulting in the erosion of the top of the stratum and the formation of a bottom parallel to the ground.

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

[0104] Parallel top mode refers to a clear erosion surface between the bottom of the stratum and the upper ground during the deposition of the sediments. This erosion surface indicates that an erosion event occurred before the formation of the stratum, resulting in the erosion of the bottom of the stratum and the formation of a top parallel to the ground. The erosion surface can be caused by factors such as erosion, sea level drop, etc.

[0105] S240, if the properties of the stratum have no obvious change rule in the spatial direction, determining that the stratum is in block mode.

[0106] Block mode refers to the fact that the sediments maintain relatively consistent properties in the vertical direction during the deposition process, i.e., the properties, composition, or characteristics of the stratum do not change significantly in the vertical direction. This uniformity indicates that the sediments have not been significantly affected by environmental changes or changes in deposition conditions during the deposition process.

[0107] S300, segmenting the plurality of seismic traces and the plurality of well tracks according to the plurality of strata respectively to obtain each seismic trace and each well track corresponding to each stratum, and calculating the position information of each sampling point in each seismic trace and each well track corresponding to each stratum;

[0108] It should be noted that in order to overcome the need for a large number of explicit layer modeling to refine the geophysical parameter model in the traditional method, to realize the consistency of the lateral trend change of the geophysical parameter model under the constraint of three-dimensional structure with the trend of the small layer deposition mode inside the stratum unit and the fine change of the geophysical parameter, and to determine the stratum of the deposition mode, the seismic trace corresponding to each stratum is parameterized according to the regional range of the model to determine the position information of each sampling point on the seismic trace corresponding to each stratum.

[0109] Specifically, step S300 can be implemented in the following way:

[0110] S310, segmenting the plurality of seismic traces according to the plurality of strata respectively 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;

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

[0112] For each seismic trace:

[0113] S311, determining a plurality of intersection points of the plurality of strata and the seismic trace;

[0114] As shown in the following table: Figure 4 Seismic trace 1 intersects with layer1, layer2, layer3, and layer4 at three points.

[0115] S312, sequentially sort the multiple intersection points from top to bottom based on the quick sort method to obtain the multiple intersection points after sorting processing;

[0116] Specifically, after calculating all intersection points, since the intersection point calculation process cannot determine the up-down order of all intersection points on the same seismic trace, the seismic trace cannot be divided into seismic trace segments according to the intersection points, and therefore, the intersection points on each seismic trace need to be sorted from top to bottom based on the quick sort method.

[0117] S313, divide the seismic trace into multiple seismic traces based on each two intersection points after sorting processing as a division basis;

[0118] Specifically, the seismic trace is divided into seismic trace segments according to the intersection point order, and the seismic trace segments are distributed into the corresponding strata while marking the belonging stratum number.

[0119] S314, match the multiple seismic traces with the multiple strata to obtain the seismic trace corresponding to each stratum;

[0120] Specifically, after dividing the seismic trace segment corresponding to each stratum, the midpoint coordinate of each segment is selected to determine which stratum the midpoint belongs to. The determination method is to emit a line upward and downward along the point, calculate the intersection points of the line with all stratum units, and determine which stratum the nearest intersection point belongs to by taking the belonging surface information of the nearest intersection point. In this way, all seismic trace segments can be distributed into the corresponding strata.

[0121] S315, set the position information of each sampling point in the seismic trace corresponding to each stratum from top to bottom with a zero starting value; wherein the distance between adjacent sampling points in the seismic trace corresponding to each stratum is the same.

[0122] Specifically, the top and bottom endpoint parameterization processing of the seismic trace corresponding to each stratum is performed. After determining all seismic traces corresponding to each stratum inside, the distance values from the top and bottom endpoints of these segments to the top and bottom interfaces of the stratum units also need to be determined.

[0123] In an embodiment, when a certain sampling point in the seismic trace corresponding to each stratum is the intersection point of the fault and the seismic trace, the position information of the intersection point of the fault and the seismic trace is reset based on the fault throw to obtain the position information of the first fault and seismic trace intersection point and the position information of the second fault and seismic trace intersection point, and the position information of the first fault and seismic trace intersection point is taken as the position information of the last sampling point of the previous stratum, and the position information of the second fault and seismic trace intersection point is taken as the position information of the initial sampling point of the next stratum.

[0124] Specifically, as Figure 4As shown, if the top end point is on the top interface, its distance value is set to 0.0, if the top end point is on a section, the distance value needs to be corrected by using the section interval on the section, so that the distance value is not 0.0, but a value between 0.0 and 1.0; if the bottom end point is on the bottom interface, its distance value is set to 1.0, if the top end point is on a section, the distance value needs to be corrected by using the section interval on the section, so that the distance value is not 1.0, but a value between 0.0 and 1.0.

[0125] S320, respectively, according to a plurality of strata, a plurality of well trajectories are segmented to obtain each stratum corresponding to each well trajectory, and the position information of each sampling point in each stratum corresponding to each well trajectory is calculated.

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

[0127] For each well trajectory:

[0128] S321, extending the plurality of well trajectories to the outside of the three-dimensional structure model to determine a plurality of intersection points of the plurality of strata and the well trajectories;

[0129] Specifically, in order to ensure that the well trajectory and the curve can penetrate the top and bottom of the model, the first point and the last point of the well trajectory and the curve are copied and moved up and down to the outside of the top and bottom interfaces of the model. This processing is to ensure that the well trajectory and the three-dimensional structure model are well processed, and at the same time, sufficient data can be searched in the top and bottom stratum unit part during subsequent interpolation.

[0130] Specifically, since the well trajectory includes vertical wells and inclined wells, determining the plurality of intersection points of the plurality of strata and the well trajectories includes the following two cases:

[0131] (1) If it is a vertical well, the intersection point is directly taken as the intersection of the well position and the seismic trace, because at this time the well trajectory coincides with the seismic grid section at this position, as shown in Figure 4 ;

[0132] (2) If it is an inclined well, the intersection point of the inclined well trajectory and the layer and section of the three-dimensional structure model needs to be calculated. At this time, since the well trajectory is a monotonous polyline in space, the entire intersection process is actually the intersection of the straight line segment on the well trajectory and the triangular elements of the layer and section in space. In order to speed up the intersection efficiency, a multi-level bounding box intersection detection and exclusion method is adopted, and the bounding box of the well trajectory is first used to judge the bounding box of the triangular grid surface of the section and the layer. If there is an intersection, it is judged whether the bounding box of the straight line segment of the well trajectory and the bounding box of the triangular element exist intersection, which excludes most of the invalid intersection calculation, and significantly improves the intersection efficiency.

[0133] S322, sequentially sort the multiple intersection points from top to bottom based on the quick sort method to obtain the multiple intersection points after sorting processing;

[0134] Specifically, the intersection point sorting and segmentation, after the intersection point calculation, similar to the intersection point sorting on the seismic trace, since the vertical well and the inclined well trajectory are monotonous, so based on the quick sort method from top to bottom sorting processing.

[0135] S323, divide the seismic trace into multiple well trajectories based on each two intersection points after sorting processing as the division basis;

[0136] Specifically, then the well trajectory is segmented into well trajectory segments according to the intersection point order, and each well trajectory segment belongs to a stratum, such as Figure 4 As shown.

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

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

[0139] Specifically, after the well trajectory segmentation, parameterization processing is needed, which includes two types of vertical well and inclined well processing:

[0140] (1) For the vertical well, since the well trajectory and the seismic trace at this position coincide, the position information of the seismic trace corresponding to each stratum can be used as the position information of the corresponding well trajectory;

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

[0142] S400, superimpose the seismic trace corresponding to each stratum 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 interpolation sampling points;

[0143] It can be understood that each seismic trace and each well trajectory includes multiple sampling points. Since the number of seismic traces is much larger than the number of well trajectories, some sampling points in the seismic trace will not coincide with the sampling points in the well trajectory, and these non-coincident sampling points are determined as interpolation sampling points. The sampling points in each well trajectory corresponding to each stratum are known sampling points. Since they do not have geophysical parameter information, they need to be inserted with geophysical parameter information, so as to ensure that each sampling point in the three-dimensional structure model has geophysical parameter information, thereby constructing a high-accuracy forward geophysical parameter model.

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

[0145] It can be understood that the more logging curves there are, the more geophysical parameter information of each sampling point in each well track corresponding to each stratum there is, and the more forward geophysical parameter models are constructed subsequently. For example, if the logging curves are P-wave curves and S-wave curves, then the geophysical parameter information of each sampling point in each well track corresponding to each stratum is S-wave and P-wave, and the forward geophysical parameter models subsequently constructed are a forward S-wave parameter model and a forward P-wave parameter model.

[0146] S600, obtaining geophysical parameter information of the sampling point to be interpolated based on a sedimentary mode of a stratum where the sampling point to be interpolated is located, position information of each sampling point in a seismic trace corresponding to the stratum where the sampling point to be interpolated is located, and the logging curves;

[0147] Specifically, step S600 can further include the following steps:

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

[0149] Specifically, step S610 includes:

[0150] S611, if the stratum where the sampling point to be interpolated is located is a proportional mode, obtaining position information of the sampling point to be interpolated based on position information of the sampling point to be interpolated, position information of an initial sampling point of the stratum where the sampling point to be interpolated is located, and position information of a final sampling point of the stratum where the sampling point to be interpolated is located;

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

[0152] S612, if the stratum where the sampling point to be interpolated is located is a parallel bottom mode, determining position information of the sampling point to be interpolated based on position information of the sampling point to be interpolated and position information of a final sampling point of the stratum where the sampling point to be interpolated is located;

[0153] Specifically, if the stratum is a parallel bottom mode, a distance value is obtained by subtracting the Z value of the bottom endpoint of the seismic trace from the Z value of the sampling point to be interpolated, and then adding the distance parameter value of the bottom endpoint to obtain the small layer parameter value of the sampling point to be interpolated; 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 endpoint of the seismic trace is the distance from the bottom of the bottom surface of the seismic trace to the bottom surface of the three-dimensional structural model.

[0154] In step S613, if the stratum where the sampling point to be interpolated is located is a 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.

[0155] Specifically, if the stratum is 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 endpoint of the seismic trace, and then adding the distance parameter value of the top endpoint to obtain the small layer parameter value of the Z value of the sampling point to be interpolated; 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 endpoint 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.

[0156] In step S620, the position information of the sampling point to be interpolated is matched with each well log curve to determine the geophysical parameter information of the sampling point to be interpolated under each well log curve.

[0157] Specifically, the position information of the sampling point to be interpolated is matched with all well log curves of the stratum where the sampling point to be interpolated is located to obtain the geophysical parameter information of the sampling point to be interpolated under each well log curve; wherein for each well curve segment, it is determined whether the position information of the sampling point to be interpolated is contained by using the top and bottom endpoint parameter values thereof, if not, it indicates that there is no interpolation control point on the well curve segment, and if yes, it is needed to quickly locate between which two adjacent well curve sampling points based on the dichotomy, and then linearly interpolate the geophysical parameter information of the control point.

[0158] In step S630, the geophysical parameter information of the sampling point to be interpolated is determined based on the geostatistical Kriging algorithm.

[0159] The geostatistical Kriging algorithm belongs to the prior art, and will not be described here.

[0160] In order to facilitate the understanding of step S600, the following is an example:

[0161] Referring to Figure 5 , Figure 5The sample point in the well log is a sample point to be interpolated, and then the depositional pattern of the layer where the sample point is located is determined as a proportional pattern, so as to calculate the position information of the sample point, and then the position information of the sample point is matched with the three curves of Well A, Well B and Well C, and the geophysical parameter information of the three points of A point on Well A, B point on Well B and C point on Well C is the geophysical parameter information of the sample point under each well log curve.

[0162] S700, based on the geophysical parameter information of the known sample point and the geophysical parameter information of the sample point to be interpolated, a forward geophysical parameter model is constructed.

[0163] Specifically, in order to realize parallel interpolation while reducing memory usage, multi-thread interpolation is used along the line direction to interpolate all seismic traces in one profile, and after one line profile is filled, the file is written, and the above steps are repeated to obtain the forward geophysical parameter model of the entire model area, and the constructed forward geophysical parameter model is as shown in Figure 6

[0164] In an embodiment, if special geological bodies in the three-dimensional structure model are considered, the geophysical parameter information of the known sample point and the geophysical parameter information of the sample point to be interpolated need to be replaced, and the specific process is as follows:

[0165] First step: calculate the geophysical parameters of the sample points in the special geological bodies in the three-dimensional structure model composed of cross sections and layer surfaces;

[0166] The special geological bodies include but are not limited to multi-value gypsiferous rock bodies, igneous rock bodies, high-speed conglomerate fan bodies and other velocity anomaly structures.

[0167] In this step, the method for calculating the geophysical parameters of the known sample points and the geophysical parameters of the sample points to be interpolated in the above embodiments can be used to calculate the geophysical parameters of the sample points in the special geological bodies. Of course, other methods can also be used to calculate the geophysical parameters of the sample points in the special geological bodies, and the embodiments of the present application do not make specific limitations thereto.

[0168] Second step: determine the sample points in the known sample points and the sample points to be interpolated that coincide with the sample points in the special geological bodies as sample points to be replaced;

[0169] Third step: replace the geophysical parameters of the sample points to be replaced with the geophysical parameters of the sample points in the special geological bodies to obtain the geophysical parameters of the sample points to be replaced after replacement;

[0170] ​Specifically, the sampling points in the known sampling points and the sampling points to be interpolated which coincide with the special geologic body are determined, and then the geophysical parameters of the sampling points of the special geologic body are replaced into the coincident sampling points.

[0171] The fourth step is to construct a forward geophysical parameter model based on the geophysical parameters of the sampling points except the sampling points to be replaced and the geophysical parameters of the sampling points to be replaced after replacement, and the constructed forward geophysical parameter model is as shown in FIG. 4. Figure 7

[0172] In the embodiment of the present application, a three-dimensional structure model based on a section and a layer is obtained, and a plurality of seismic traces and a plurality of well tracks in the range of the three-dimensional structure model are obtained; each seismic trace and each well track comprises a plurality of sampling points; each well track corresponds to a plurality of logging curves of the corresponding relationship between the position information of different sampling points and different geophysical parameter information; the three-dimensional structure model is divided into a plurality of strata, and the deposition mode of each stratum is determined; the plurality of seismic traces and the plurality of well tracks are processed by segments according to the plurality of strata respectively, to obtain each seismic trace and each well track corresponding to each stratum, and the position information of each sampling point in each seismic trace and each well track corresponding to each stratum is calculated; each seismic trace corresponding to each stratum is superimposed with each well track 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 logging curves to obtain the geophysical parameter information of each sampling point in each well track corresponding to each stratum; the geophysical parameter information of the sampling points to be interpolated is obtained based on the deposition mode of the stratum where the sampling points to be interpolated are located, the position information of each sampling point in the seismic trace corresponding to the stratum where the sampling points to be interpolated are located, and the logging curves; and a forward geophysical parameter model is constructed based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the sampling points to be interpolated. The present application can finely describe the distribution of geophysical parameters in the internal stratum unit without interpreting a large number of horizons for modeling and only interpreting the main velocity variation interface layer, and the section model does not need to be simplified, thereby improving the construction efficiency of the forward geophysical parameter model in a complex structure region.

[0173] Embodiment two

[0174] ​The method for constructing a forward geophysical parameter model provided by the application is applied to a double complex typical exploration area in the Tarim Basin. The main features of the area are complex surface velocity structure and highly developed underground overthrust nappe faults, which lead to complex underground velocity structure. How to establish a fine forward velocity model that can represent the double complex exploration area under the constraint of a structural framework model is one of the main technical challenges of typical forward modeling. In order to establish a fine forward model, 22 reverse faults are interpreted. Since reverse fault interpretation needs to be interpreted in upper and lower plates, a total of 85 slices of interpreted horizons are divided. At the same time, VSP logging and sonic logging in the exploration area, prestack depth migration velocity model and poststack seismic data volume information are collected for velocity model establishment. The structural model established based on the above-mentioned faults and horizon data is shown in Figure 8 . In order to establish a velocity model, the 15 VSP and sonic logging data in the area are reconstructed based on artificial intelligence curve reconstruction technology using poststack seismic data and migration velocity model before modeling to fill in the gaps. The stratigraphic unit sedimentary pattern is shown in Figure 9 . The underground velocity model established based on the above-mentioned data using the technology of the application is shown in Figure 10 and Figure 11 and

[0175] Example Three

[0176] Please refer to Figure 12 , Figure 12 is a structural schematic diagram of a device for constructing a forward geophysical parameter model provided by an embodiment of the application.

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

[0178] A data acquisition module 210 is configured to acquire a three-dimensional structural model composed of sections and horizons, and a plurality of seismic traces and a plurality of well tracks within the range of the three-dimensional structural model. Each seismic trace and each well track includes a plurality of sampling points. Each well track corresponds to a plurality of logging curves of the corresponding relationship between the position information of different sampling points and different geophysical parameter information.

[0179] A model layering module 220 is configured to divide the three-dimensional structural model into a plurality of strata and determine the sedimentary pattern of each stratum.

[0180] A model parameterization module 230 is configured to perform segmented processing on the plurality of seismic traces and the plurality of well tracks according to the plurality of strata, to obtain each seismic trace and each well track corresponding to each stratum, and to calculate the position information of each sampling point in each seismic trace and each well track corresponding to each stratum.

[0181] The sampling point classification module 240 is configured to superimpose each seismic trace corresponding to each stratum with each well track corresponding to each stratum, to determine the overlapped sampling points as known sampling points, and to determine the non-overlapped sampling points as interpolation sampling points;

[0182] The first parameter determination module 250 is configured to match the known sampling points with the corresponding logging curves, to obtain the geophysical parameter information of each sampling point in each well track corresponding to each stratum;

[0183] The second parameter determination module 260 is configured to obtain the geophysical parameter information of the interpolation sampling points based on the sedimentary mode of the stratum where the interpolation sampling points are located, the position information of each sampling point in the seismic trace corresponding to the stratum where the interpolation sampling points are located, and the logging curves.

[0184] The model construction module 270 is configured to construct a forward geophysical parameter model based on the geophysical parameter information of the known sampling points and the geophysical parameter information of the interpolation sampling points.

[0185] It should be understood that the device corresponds to the above-mentioned embodiments of the method for constructing a forward geophysical parameter model, and can perform each step involved in the above-mentioned method embodiments. The specific functions of the device can be referred to the description above, and the detailed description is appropriately omitted here to avoid repetition. The device includes at least one software function module that can be stored in the memory in the form of software or firmware or solidified in the operating system (OS) of the device.

[0186] Embodiment four

[0187] Based on the same inventive concept, the embodiments of the present application also provide an electronic device, comprising a processor and a memory, the memory storing machine readable instructions executable by the processor, and the machine readable instructions being executed by the processor to perform the above-mentioned method for constructing a forward geophysical parameter model.

[0188] The memory can include non-persistent memory in computer readable media, random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash memory (flash RAM). The memory is an example of computer readable media.

[0189] Computer-readable media includes permanent and non-permanent, movable and non-movable media, which 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 technologies, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassette, magnetic tape, magnetic disk storage or other magnetic storage devices, or any other non-transmission medium that can be used to store information accessible by a computing device. According to the definition herein, computer-readable media does not include transitory media such as modulated data signals and carriers.

[0190] Embodiment five

[0191] Based on the same inventive concept, the embodiments of the present application also provide a computer-readable storage medium, which stores computer instructions, when the computer instructions are run on a computer, the computer instructions make the computer execute the forward geophysical parameter model construction method described above.

[0192] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system or a computer program product. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0193] The present application is described with reference to flowcharts and / or block diagrams according to the method, device (system) and computer program product of the embodiments of the present application. It should be understood that each flow and / or block in the flowcharts and / or block diagrams, and the combination of the flows and / or blocks in the flowcharts and / or block diagrams can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor 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 that implements the functions specified in the flowcharts and / or block diagrams. Figure 1 The function of the device specified in one flow or multiple flows and / or blocks Figure 1 The function of the device specified in one flow or multiple flows and / or blocks

[0194] These computer program instructions can also be stored in a computer readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer readable memory produce an article of manufacture including instructions which implement the Figure 1 function specified in the flow or flows and / or blocks Figure 1 of the block or blocks.

[0195] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions that are executed on the computer or other programmable apparatus provide steps for implementing the functions specified in the flow or flows and / or blocks Figure 1 function specified in the flow or flows and / or blocks Figure 1 of the block or blocks.

[0196] It should also be noted that each of the various illustrative technical features described in the above specific embodiments can be implemented in any suitable manner without departing from the disclosure. To avoid unnecessary repetition, the various possible combinations of the technical features are not described separately in the embodiments of the disclosure.

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

[0198] It should also be noted that the terms "comprising", "including", or any other variant thereof are intended to cover non-exclusive inclusions, such that processes, methods, articles, or devices that include a series of elements not only include those elements, but also include other elements not explicitly listed, or other elements inherent in such processes, methods, articles, or devices. Without more limitations, the element defined by the phrase "comprising a" does not exclude the presence of additional identical elements in the process, method, article, or device that includes the element.

[0199] The above is only an embodiment of the present application and is not intended to limit the present application. The present application can have various modifications and changes for those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present application shall be included in the scope of 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, 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; Divide the 3D structural model into multiple strata and determine the sedimentary pattern of each stratum; 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; 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; 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 the sampling points to be interpolated, a forward geophysical parameter model is constructed; 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 point to be interpolated under each logging curve is calculated to determine the geophysical parameter information of the sampling point to be interpolated.

2. The method for constructing a forward geophysical parameter model according to claim 1, wherein: 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 specific geological bodies in a three-dimensional structural model 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.

3. The method for constructing a forward geophysical parameter model according to claim 1, wherein: 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 a parallel bottom mode; If the ground surface of the formation is an erosion surface, the formation is determined to be a parallel top model; If the properties of the stratum have no obvious variation pattern in the spatial direction, the stratum is determined to be a blocky pattern.

4. The method for constructing a forward geophysical parameter model according to claim 1, wherein: 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: 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; 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, wherein: Multiple seismic traces are segmented according to 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, 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 sorted intersection points; 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 distance between adjacent sampling points in the seismic trace corresponding to each stratum is the same.

6. The method for constructing a forward geophysical parameter model according to claim 5, wherein: Starting from zero, set the position information of each sampling point in the seismic trace corresponding to each stratum from top to bottom, including: 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.

7. The method for constructing a forward geophysical parameter model according to claim 4, wherein: 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 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 sorted intersection points; The seismic trace is divided into multiple well trajectories based on every two intersection points after sorting; Match multiple well trajectories with multiple strata to obtain the 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 sequentially set from top to bottom; wherein, the distance between adjacent sampling points in the well trajectory corresponding to each stratum is the same.

8. The method for constructing a forward geophysical parameter model according to claim 1, wherein: 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: 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 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; 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.

9. A device for constructing a forward geophysical parameter model, characterized in that: include: The data acquisition module is used to obtain a three-dimensional structural model composed of sections and layers, 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 logging curves of the corresponding relationship between the position information of different sampling points and different geophysical parameter information Model layering module, used to divide the 3D structural model into multiple strata and determine the sedimentary pattern of each stratum; 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; 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 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 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; A model building module is used to build a forward geophysical parameter model based on the geophysical parameter information of known sampling points and the geophysical parameter information of the sampling points to be interpolated; The second parameter determination module is specifically used for: 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 point to be interpolated under each logging curve is calculated to determine the geophysical parameter information of the sampling point to be interpolated.

10. 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 according to any one of claims 1 to 8 is executed.

11. 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 8.

Images