阅读说明：本技术 地层参数模型建立方法及装置 (Stratum parameter model establishing method and device ) 是由 符志国 陈康 于 2020-04-09 设计创作，主要内容包括：本发明提供了一种地层参数模型建立方法及装置,其中该方法包括：获取目标地层的地震信号和目标地层的测井数据,建立目标地层的初始地层参数模型；若初始地层参数模型吻合目标地层的地震信号的变化特征,将初始地层参数模型确定为目标地层的地层参数模型；若不吻合,根据目标地层的地震信号的变化特征,调整多个测井在待求位置处的地层参数加权值,建立目标地层的地层参数模型。该方法在地层参数值求取时,需要根据每口测井在待求位置处的地层参数和地层参数加权值确定,使得求取的地层参数值与测量点间的地层参数变化相关；调整多个测井在待求位置处的地层参数加权值,使得模型贴切地反映真实地层参数分布结构特征,提高了准确性。(The invention provides a method and a device for establishing a stratum parameter model, wherein the method comprises the following steps: acquiring seismic signals of a target stratum and logging data of the target stratum, and establishing an initial stratum parameter model of the target stratum; if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as a stratum parameter model of the target stratum; if the seismic signals do not coincide with each other, adjusting stratum parameter weighted values of the multiple logging wells at the positions to be solved according to the variation characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum. When stratum parameter values are obtained, the method needs to determine according to stratum parameters of each log at a position to be obtained and stratum parameter weighted values, so that the obtained stratum parameter values are related to stratum parameter changes among measurement points; and the stratum parameter weighted values of a plurality of logging wells at the positions to be solved are adjusted, so that the model can appropriately reflect the real stratum parameter distribution structure characteristics, and the accuracy is improved.)

1. A method for establishing a formation parameter model is characterized by comprising the following steps:

acquiring seismic signals of a target stratum and logging data of the target stratum; the logging data of the target stratum comprises a plurality of stratum parameter values of the target stratum;

establishing an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum; determining the stratum parameter values of the positions to be solved in the initial stratum parameter model according to the stratum parameter values of the multiple logging wells at the positions to be solved and stratum parameter weighted values;

if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as the stratum parameter model of the target stratum;

if the initial stratum parameter model does not fit the change characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of the multiple logs at the positions to be solved according to the change characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum; and determining the formation parameter values at the positions to be solved in the formation parameter model of the target formation according to the formation parameter values of the multiple logs at the positions to be solved and the adjusted formation parameter weighted values.

2. The method of claim 1, wherein determining the formation parameter value at the desired location in the initial formation parameter model based on the formation parameter values and the formation parameter weighting values for the plurality of logs at the desired location comprises:

and taking the weighted values of the formation parameters of the multiple logging wells at the positions to be solved as the weights of the formation parameter values of the multiple logging wells at the positions to be solved, and carrying out weighted summation to obtain the formation parameter values of the positions to be solved in the initial formation parameter model.

3. The method of claim 2, wherein the formation parameter value at the desired location in the initial formation parameter model is determined based on the formation parameter values and the formation parameter weighting values for the plurality of logs at the desired location according to the formula:

M＝w₁×Pw₁+w₂×Pw₂+…+w_R×Pw_R

wherein M represents a stratum parameter value at a position to be solved in the initial stratum parameter model; w is a_iI is 1,2, …, and R represents the weighted value of the parameter of the ith logging well at the position to be solved; pw_i1,2, …, R, representing the formation parameter value of the ith log at the location to be found; r is more than 1 and represents the logging number; the position to be solved is described by a three-dimensional coordinate consisting of a ground coordinate and a stratum depth.

4. The method of claim 3 wherein the parameter weighting w for the ith well log at the desired location is determined according to the formula_i：

w_i＝a_i/S,i＝1,2,…,R

Wherein, S ═ a₁+a₂+…+a_R；f_iThe adjustment factor of the ith logging at the position to be solved is represented as a real number; d_iRepresenting the spatial horizontal distance between the position to be solved and the ith logging well; q represents a power parameter, which is a real number.

5. The method of claim 4, wherein adjusting the formation parameter weighting values of the plurality of logs at the desired location based on the varying characteristics of the seismic signals of the target formation to create a formation parameter model of the target formation comprises:

adjusting the adjustment factor f of the ith well logging at the position to be solved according to the change characteristics of the seismic signal of the target stratum_iAnd establishing a stratum parameter model of the target stratum according to the power parameter q.

6. Method according to claim 4, characterized in that the ground coordinates (X) according to the position to be found are used_n,Y_n) Determining the spatial level between the position to be determined and the ith well log according to the following formulaDistance d_i：

Wherein (X)_i,Y_i) Representing the surface coordinates of the ith log.

7. Method according to claim 6, characterized in that the ground coordinates (X) of the position to be determined are determined_n,Y_n) The method comprises the following steps:

dividing the target stratum into a plurality of seismic Trace traces₁,Trace₂,…,Trace_N(N > R); therein, Trace_iThe seismic trace at the ith logging position is taken;

determining the plurality of seismic Trace traces₁,Trace₂,…,Trace_NThe ground coordinates of (a);

determining the Trace of the seismic record of the position to be solved according to the position to be solved_n，1≤n≤N；

Tracing the seismic Trace_nIs determined as the ground coordinates (X) of the location to be determined_n,Y_n)。

8. The method of claim 7, wherein a formation parameter value Pw for the ith log at the location to be determined is determined_iThe method comprises the following steps:

obtaining an on-time window [ wt ] according to the logging data of the ith logging at the target stratum₁,wt₂]Formation parameter sampling data Pw of the ith logging well in the range at different sampling depths in the target formation_i,jI ═ 1,2, …, R; j ═ 1,2, …, L, and sampling interval ds; wherein j represents the sampling point number at different sampling depths in the target stratum; l represents the total number of sampling points;

determining the position moment t of each sampling point of the ith logging in the target stratum according to the sampling depth and the seismic signal of the target stratum_i,j；

According to whatThe formation parameter sampling data Pw of the ith well logging at different sampling depths in the target formation_i,jAnd the position moment t of each sampling point of the ith well logging in the target stratum_i,jAnd sampling time interval ds, and determining formation parameter value Pw of the ith logging well at the position to be solved by adopting linear interpolation_i。

9. The method of claim 8, further comprising:

determining a top interface seismic horizon and a bottom interface seismic horizon of the target stratum according to ground seismic exploration;

determining a plurality of seismic Trace traces according to the seismic signals of the target stratum, and the top interface seismic horizon and the bottom interface seismic horizon of the target stratum₁,Trace₂,…,Trace_N(N > R) Top interface time Top (Trace)_K) K1, 2, …, N and bottom interface time bot (Trace)_K),K＝1,2,…,N；

According to the following relation, according to the multiple seismic record Trace traces₁,Trace₂,…,Trace_N(N > R) Top interface time Top (Trace)_K) Bottom interface time bot (Trace)_K) Determining the time window [ wt ]₁,wt₂]：

wt₁＜min(top(Trace_K))＜max(bot(Trace_K))＜wt₂

Wherein, min (Trace)_K) Represents the top interface time top (Trace) of multiple seismic traces_K) K is the minimum of 1,2, …, N; max (Trace)_K) Represents the bottom boundary time bot (Trace) of a plurality of seismic traces_K) K is the maximum of 1,2, …, N.

10. The method of claim 9, wherein the formation parameter sample data Pw is at different sample depths in the formation of interest based on the ith log_i,jAnd the position moment t of each sampling point of the ith well logging in the target stratum_i,jAnd sampling time interval ds, and determining the stratum of the ith logging at the position to be solved by adopting linear interpolationParameter value Pw_iThe method comprises the following steps:

determining the position moment t of the position to be solved according to the stratum depth of the position to be solved;

according to the position moment t and the time window [ wt ] of the position to be solved₁,wt₂]The Trace of the seismic record of the position to be solved_nTop interface time top (Trace)_n) And the seismic Trace of the position to be solved_nBottom interface time bot (Trace)_n) Determining the equal proportional position time T of the position to be obtained at the ith well logging position_i；

According to the formation parameter sampling data Pw of the ith logging at different sampling depths in the target formation_i,jAnd the equal proportional position time T of the position to be solved at the ith well logging position_iThe position time t of each sampling point in the target stratum of the ith well logging_i,jAnd the sampling time interval ds, according to the following formula, determining the formation parameter value Pw of the ith well logging at the position to be solved_i：

Pw_i＝Pw_i,j+(Pw_i,j+1-Pw_i,j)×(T_i-t_i,j)/ds

Wherein, wt₁≤t_i,j≤T_i≤t_i,j+1≤wt₂；j＝1,2,…,L-1。

11. Method according to claim 10, characterized in that the position time t and the time window [ wt ] are dependent on the position of the position to be determined₁,wt₂]The Trace of the seismic record of the position to be solved_nTop interface time top (Trace)_n) And the seismic Trace of the position to be solved_nBottom interface time bot (Trace)_n) Determining the equal proportional position time T of the position to be obtained at the ith well logging position_iThe method comprises the following steps:

if wt₁≤t＜top(Trace_n) And then:

T_i＝wt₁+(top(Trace_i)-wt₁)×(t-wt₁)/(top(Trace_n)-wt₁)；

if top (Trace)_n)≤t＜bot(Trace_n) And then:

if bot (Trace)_n)≤t≤wt₂And then:

T_i＝bot(Trace_n)+(wt₂-bot(Trace_i))×(t-bot(Trace_n))/(wt₂-bot(Trace_n))

wherein, top (Trace)_i) Representing the top interface time of the seismic trace at the ith log; bot (Trace)_i) Representing the bottom interface time of the seismic trace at the ith log.

12. A formation parameter model building apparatus, comprising:

the data acquisition module is used for acquiring seismic signals of a target stratum and logging data of the target stratum; the logging data of the target stratum comprises a plurality of stratum parameter values of the target stratum;

the initial model building module is used for building an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum; determining the stratum parameter values of the positions to be solved in the initial stratum parameter model according to the stratum parameter values of the multiple logging wells at the positions to be solved and stratum parameter weighted values;

the first model establishing module is used for determining the initial stratum parameter model as the stratum parameter model of the target stratum if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum;

the second model establishing module is used for adjusting stratum parameter weighted values of the multiple logs at positions to be solved according to the change characteristics of the seismic signals of the target stratum and establishing a stratum parameter model of the target stratum if the initial stratum parameter model is not matched with the change characteristics of the seismic signals of the target stratum; and determining the formation parameter values at the positions to be solved in the formation parameter model of the target formation according to the formation parameter values of the multiple logs at the positions to be solved and the adjusted formation parameter weighted values.

13. The apparatus of claim 12, wherein the initial model building module comprises:

and the stratum parameter value calculation unit is used for taking the stratum parameter weighted values of the multiple logging wells at the positions to be solved as the weights of the stratum parameter values of the multiple logging wells at the positions to be solved, and carrying out weighted summation to obtain the stratum parameter values of the positions to be solved in the initial stratum parameter model.

14. The apparatus of claim 13, wherein the formation parameter value calculation unit is specifically configured to:

determining the formation parameter value at the position to be solved in the initial formation parameter model according to the formation parameter values and the formation parameter weighted values of the plurality of logs at the position to be solved according to the following formula:

M＝w₁×Pw₁+w₂×Pw₂+…+w_R×Pw_R

wherein M represents a stratum parameter value at a position to be solved in the initial stratum parameter model; w is a_iI is 1,2, …, and R represents the weighted value of the parameter of the ith logging well at the position to be solved; pw_i1,2, …, R, representing the formation parameter value of the ith log at the location to be found; r is more than 1 and represents the logging number; the position to be solved is described by a three-dimensional coordinate consisting of a ground coordinate and a stratum depth.

15. The apparatus of claim 14, wherein the initial model building module comprises: a weight value calculation unit for:

determining the weighted value w of the parameter of the ith well logging at the position to be determined according to the following formula_i：

w_i＝a_i/S,i＝1,2,…,R

Wherein, S ═ a₁+a₂+…+a_R；f_iThe adjustment factor of the ith logging at the position to be solved is represented as a real number; d_iRepresenting the spatial horizontal distance between the position to be solved and the ith logging well; q represents a power parameter, which is a real number.

16. The apparatus of claim 15, wherein the second model building module is specifically configured to:

adjusting the adjustment factor f of the ith well logging at the position to be solved according to the change characteristics of the seismic signal of the target stratum_iAnd establishing a stratum parameter model of the target stratum according to the power parameter q.

17. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 11 when executing the computer program.

18. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 11.

Technical Field

The invention relates to the technical field of geophysical exploration, in particular to a method and a device for establishing a stratum parameter model.

Background

One of the major efforts in petroleum seismic exploration is to predict the spatial distribution of underground petroleum and gas reservoirs. In this work, model-based seismic data inversion methods are widely used. The implementation process of the inversion method needs to use a stratum parameter model when being initialized. Many studies have shown that the processing results of model-based seismic data inversion methods depend on this formation parameter model. The structural characteristics of the formation parameter model determine the structural characteristics of the inversion result. Therefore, the production of a stratum parameter model is a key link for predicting the petroleum and natural gas reservoir.

In the prior art, a stratum parameter model is established by using an inverse distance weighted interpolation method, but the result of the method only depends on the distance between an interpolation point and a measuring point, but is irrelevant to the stratum parameter change between the measuring points, the real stratum parameter distribution structure characteristics of a target stratum cannot be reflected in a proper way, and compared with the reality, the error is large.

Disclosure of Invention

The embodiment of the invention provides a stratum parameter model establishing method, which is used for establishing a stratum parameter model according to stratum parameter change among measuring points, appropriately reflecting real stratum parameter distribution structure characteristics of a target stratum and improving the accuracy of the established stratum parameter model, and comprises the following steps:

acquiring seismic signals of a target stratum and logging data of the target stratum; the logging data of the target stratum comprises a plurality of stratum parameter values of the target stratum;

establishing an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum; determining the stratum parameter values of the positions to be solved in the initial stratum parameter model according to the stratum parameter values of the multiple logging wells at the positions to be solved and stratum parameter weighted values;

if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as the stratum parameter model of the target stratum;

if the initial stratum parameter model does not fit the change characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of the multiple logs at the positions to be solved according to the change characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum; and determining the formation parameter values at the positions to be solved in the formation parameter model of the target formation according to the formation parameter values of the multiple logs at the positions to be solved and the adjusted formation parameter weighted values.

The embodiment of the invention also provides a device for establishing a stratum parameter model, which is used for establishing the stratum parameter model according to the stratum parameter change among the measuring points, appropriately reflecting the real stratum parameter distribution structure characteristics of the target stratum and improving the accuracy of the established stratum parameter model, and comprises the following steps:

the data acquisition module is used for acquiring seismic signals of a target stratum and logging data of the target stratum; the logging data of the target stratum comprises a plurality of stratum parameter values of the target stratum;

the initial model building module is used for building an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum; determining the stratum parameter values of the positions to be solved in the initial stratum parameter model according to the stratum parameter values of the multiple logging wells at the positions to be solved and stratum parameter weighted values;

the first model establishing module is used for determining the initial stratum parameter model as the stratum parameter model of the target stratum if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum;

the second model establishing module is used for adjusting stratum parameter weighted values of the multiple logs at positions to be solved according to the change characteristics of the seismic signals of the target stratum and establishing a stratum parameter model of the target stratum if the initial stratum parameter model is not matched with the change characteristics of the seismic signals of the target stratum; and determining the formation parameter values at the positions to be solved in the formation parameter model of the target formation according to the formation parameter values of the multiple logs at the positions to be solved and the adjusted formation parameter weighted values.

The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the formation parameter model building method when executing the computer program.

Embodiments of the present invention also provide a computer-readable storage medium, which stores a computer program for executing the above-mentioned formation parameter model building method.

In the embodiment of the invention, seismic signals of a target stratum and logging data of the target stratum are obtained; establishing an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum; if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as a stratum parameter model of the target stratum; if the initial stratum parameter model does not fit the change characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of a plurality of logs at positions to be solved according to the change characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum; when the formation parameter value at each position to be solved is solved, the formation parameter value at each position to be solved is determined according to the formation parameter of each logging at the position to be solved and the formation parameter weighted value, so that the formation parameter value at each position to be solved is related to the formation parameter change between the measurement points; compared with the prior art of establishing a stratum parameter model by an inverse distance weighted interpolation method, the established stratum parameter model can more closely reflect the real stratum parameter distribution structure characteristics of the target stratum by adjusting the stratum parameter weighted values of a plurality of logs at the positions to be solved, the error between the actual stratum parameters and the established stratum parameter model is small, and the accuracy of the established stratum parameter model is improved.

Drawings

In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.

FIG. 1 is a schematic diagram of a method for establishing a formation parameter model according to an embodiment of the present invention.

FIG. 2 is a diagram illustrating an exemplary embodiment of the present invention for determining a formation parameter value Pw for an ith log at a desired location_iThe method of (1) is a schematic flow diagram.

FIG. 3 is a schematic diagram of a multi-well logging tool connection in accordance with an embodiment of the present invention.

FIG. 4 is a schematic diagram of A, B, C three log locations and line locations in another embodiment of the present invention.

FIG. 5 is a block diagram of a device for establishing a formation parameter model according to an embodiment of the present invention.

Detailed Description

The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

The embodiment of the invention provides a stratum parameter model establishing method, which is used for establishing a stratum parameter model according to stratum parameter change among measuring points, appropriately reflecting real stratum parameter distribution structure characteristics of a target stratum and improving the accuracy of the established stratum parameter model, and as shown in figure 1, the method comprises the following steps:

step 101: acquiring seismic signals of a target stratum and logging data of the target stratum; the logging data of the target stratum comprise a plurality of stratum parameter values of the logging in the target stratum;

step 102: establishing an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum;

the stratum parameter value of the position to be solved in the initial stratum parameter model is determined according to the stratum parameter values of the plurality of logging wells at the position to be solved and the stratum parameter weighted value;

step 103: if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as a stratum parameter model of the target stratum;

step 104: if the initial stratum parameter model does not fit the change characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of a plurality of logs at positions to be solved according to the change characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum;

and determining the formation parameter value at the position to be solved in the formation parameter model of the target formation according to the formation parameter values of the multiple logging wells at the position to be solved and the adjusted formation parameter weighted value.

As can be seen from FIG. 1, in the embodiment of the invention, the seismic signals of the target stratum and the logging data of the target stratum are acquired; establishing an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum; if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as a stratum parameter model of the target stratum; if the initial stratum parameter model does not fit the change characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of a plurality of logs at positions to be solved according to the change characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum; when the formation parameter value at each position to be solved is solved, the formation parameter value at each position to be solved is determined according to the formation parameter of each logging at the position to be solved and the formation parameter weighted value, so that the formation parameter value at each position to be solved is related to the formation parameter change between the measurement points; compared with the prior art of establishing a stratum parameter model by an inverse distance weighted interpolation method, the established stratum parameter model can more closely reflect the real stratum parameter distribution structure characteristics of the target stratum by adjusting the stratum parameter weighted values of a plurality of logs at the positions to be solved, the error between the actual stratum parameters and the established stratum parameter model is small, and the accuracy of the established stratum parameter model is improved.

In specific implementation, firstly, a seismic signal of a target stratum and logging data of the target stratum are obtained, wherein the logging data of the target stratum comprise stratum parameter values of a plurality of logs in the target stratum, and the stratum parameters can be basic data provided for reservoir description, such as parameter data including lithology, shale content, porosity, permeability and the like. The seismic signal is a signal obtained by transmitting a group of continuously vibrating elastic wave signals (also called scanning signals) to the underground by a controllable seismic source in a reaction mode, processing and identifying reflected wave signals received by the ground, and is used for explaining the structural form and the occurrence of an underground geological target.

After the seismic signal of the target stratum and the logging data of the target stratum are obtained, an initial stratum parameter model of the target stratum is established according to the seismic signal of the target stratum and the logging data of the target stratum. In a specific embodiment, the position to be solved in the initial formation parameter model is described by a three-dimensional coordinate composed of a surface coordinate and a formation depth, weighted values of formation parameters of a plurality of logs at the position to be solved are taken as weights of formation parameter values of the plurality of logs at the position to be solved, weighted summation is performed to obtain the formation parameter value of the position to be solved in the initial formation parameter model, and the formation parameter value is determined according to the following formula:

M＝w₁×Pw₁+w₂×Pw₂+…+w_R×Pw_R (1)

wherein M represents a stratum parameter value at a position to be solved in the initial stratum parameter model; w is a_iI is 1,2, …, and R represents the weighted value of the parameter of the ith logging well at the position to be solved; pw_i1,2, …, R, representing the formation parameter value of the ith log at the location to be found; r > 1, representing the number of logs.

The weighted value w of the parameter of the ith well logging in the formula (1) at the position to be solved_iCalculated according to the following formula:

w_i＝a_i/S,i＝1,2,…,R (2)

wherein, S ═ a₁+a₂+…+a_R；f_iThe adjustment factor of the ith logging at the position to be solved is represented as a real number; d_iRepresenting the spatial horizontal distance between the position to be solved and the ith logging well; q represents a power parameter, which is a real number.

Dividing target stratum into a plurality of seismic recording traces₁,Trace₂,…,Trace_N(N > R); therein, Trace_iThe seismic trace at the ith logging position is taken; determining multiple seismic Trace traces₁,Trace₂,…,Trace_NThe ground coordinates of (a); determining the Trace of the seismic record of the position to be solved according to the position to be solved_nN is more than or equal to 1 and less than or equal to N; trace of seismic Trace_nIs determined as the ground coordinates (X) of the location to be determined_n,Y_n)。

Ground coordinates (X) according to the position to be determined_n,Y_n) And the surface coordinates (X) of the ith well log_i,Y_i) D is determined according to the following formula (3)_i：

D is determined according to the formula (3)_iAfter a, a_iIs then determined, then w_iAnd (4) determining.

The formation parameter Pw of the ith well logging in the above formula (1) at the position to be found_iAs shown in fig. 2, the determination process includes:

step 201: obtaining an on-time window [ wt ] according to the logging data of the ith logging at the target stratum₁,wt₂]Formation parameter sampling data Pw of the ith logging well in the range at different sampling depths in the target formation_i,jI ═ 1,2, …, R; j ═ 1,2, …, L, and sampling interval ds; wherein j represents the sampling point number at different sampling depths in the target stratum; l meterShowing the total number of sampling points;

step 202: determining the position moment t of each sampling point of the ith logging in the target stratum according to the sampling depth and the seismic signal of the target stratum_i,j；

Step 203: according to formation parameter sampling data Pw of the ith logging at different sampling depths in the target formation_i,jAnd the position moment t of each sampling point of the ith well logging in the target stratum_i,jAnd sampling time interval ds, and determining formation parameter value Pw of the ith logging well at the position to be solved by adopting linear interpolation_i。

In a specific implementation process, on the basis of the method flow shown in fig. 2, the method further includes:

determining a top interface seismic horizon and a bottom interface seismic horizon of a target stratum according to ground seismic exploration;

determining the Trace of a plurality of seismic record tracks according to the seismic signals of the target stratum, the top interface seismic horizon and the bottom interface seismic horizon of the target stratum₁,Trace₂,…,Trace_N(N > R) Top interface time Top (Trace)_K) K1, 2, …, N and bottom interface time bot (Trace)_K),K＝1,2,…,N；

According to the following relation, according to the multiple seismic record Trace traces₁,Trace₂,…,Trace_N(N > R) Top interface time Top (Trace)_K) Bottom interface time bot (Trace)_K) Determining the time window [ wt ]₁,wt₂]：

wt₁＜min(top(Trace_K))＜max(bot(Trace_K))＜wt₂ (4)

Wherein, min (Trace)_K) Represents the top interface time top (Trace) of multiple seismic traces_K) K is the minimum of 1,2, …, N; max (Trace)_K) Represents the bottom boundary time bot (Trace) of a plurality of seismic traces_K) K is the maximum of 1,2, …, N.

The time of the top interface and the time of the bottom interface of the plurality of seismic traces are the time of the top interface and the time of the bottom interface of the plurality of seismic traces, the stratum depth and the bottom interface stratum depth of the top interface of each seismic trace are determined on a seismic section according to the top interface and the bottom interface of each seismic trace, and the time value corresponding to the stratum depth of the top interface of each seismic trace and the time value corresponding to the bottom interface stratum depth are found according to the corresponding relation between the stratum depth and the seismic signal time axis, namely the time of the top interface and the time of the bottom interface of the plurality of seismic traces.

Determining the time window [ wt ] according to the relation (4)₁,wt₂]Determining a good time window [ wt ]₁,wt₂]Then, the time window [ wt ] is obtained₁,wt₂]Pw within the range_i,jAnd ds, determining the position moment t of each sampling point of the ith well in the target stratum on the seismic signal of the target stratum according to the depth of each sampling point of the ith well in the target stratum_i,jObtaining Pw_i,jDs and t_i,jThen, linear interpolation is adopted to determine the formation parameter value Pw of the ith well logging at the position to be solved_iThe specific process comprises the following steps:

determining the position moment t of the position to be solved in the seismic signal of the target stratum according to the stratum depth of the position to be solved;

according to the position moment t and the time window [ wt ] of the position to be obtained₁,wt₂]Seismic Trace to be located_nTop interface time top (Trace)_n) And the Trace of the seismic record of the position to be solved_nBottom interface time bot (Trace)_n) Determining the equal proportional position time T of the position to be obtained at the ith well logging position_i；

According to formation parameter sampling data Pw of the ith logging at different sampling depths in the target formation_i,jAnd the equal proportional position time T of the position to be solved at the ith well logging position_iThe position time t of each sampling point of the ith well logging in the target stratum_i,jAnd the sampling time interval ds, according to the following formula, determining the formation parameter value Pw of the ith well logging at the position to be solved_i：

Pw_i＝Pw_i,j+(Pw_i,j+1-Pw_i,j)×(T_i-t_i,j)/ds (5)

Wherein, wt₁≤t_i,j≤T_i≤t_i,j+1≤wt₂；j＝1,2,…,L-1。

T in the above formula (5)_iThe determination process of (1), comprising:

if wt₁≤t＜top(Trace_n) And then:

T_i＝wt₁+(top(Trace_i)-wt₁)×(t-wt₁)/(top(Trace_n)-wt₁) (6)

if top (Trace)_n)≤t＜bot(Trace_n) And then:

if bot (Trace)_n)≤t≤wt₂And then:

T_i＝bot(Trace_n)+(wt₂-bot(Trace_i))×(t-bot(Trace_n))/(wt₂-bot(Trace_n)) (8)

wherein, top (Trace)_i) Representing the top interface time of the seismic trace at the ith log; bot (Trace)_i) Representing the bottom interface time of the seismic trace at the ith log.

According to the position moment t and the time window [ wt ] of the position to be obtained₁,wt₂]Seismic Trace to be located_nTop interface time top (Trace)_n) And the Trace of the seismic record of the position to be solved_nBottom interface time bot (Trace)_n) Is different, T is determined by the above formula (6), formula (7) or formula (8)_i，T_iAfter the determination, the formation parameter value Pw of the ith well logging at the position to be solved can be determined according to the formula (5)_i。

w_iAnd Pw_iAnd (3) after all the stratum parameters are determined, determining the stratum parameter values of the positions to be solved in the initial stratum parameter model according to the formula (1), and after the stratum parameters of all the positions are determined according to the process, establishing the initial stratum parameter model.

And if the established initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as the stratum parameter model of the target stratum.

And if the established initial stratum parameter model does not conform to the variation characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of the multiple logs at the positions to be solved according to the variation characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum. Wherein, the formation parameter value M' at the position to be solved in the formation parameter model of the target formation is obtained according to the formation parameter values Pw of a plurality of logs at the position to be solved_iAnd adjusted formation parameter weighting value w_i' determination. The specific calculation formula is as follows:

M'＝w₁'×Pw₁+w₂'×Pw₂+…+w_R'×Pw_R (9)

wherein M' represents a stratum parameter value at a position to be solved in a stratum parameter model of the target stratum; w is a_i', i-1, 2, …, wherein R represents the weighted value of the parameter of the adjusted ith logging well at the position to be found.

In a particular embodiment of the present invention,

w_i'＝a_i'/S',i＝1,2,…,R (10)

wherein, S' ═ a₁'+a₂'+…+a_R'；f_i' represents the adjusting factor of the adjusted ith logging in the position to be solved, and is a real number; q' represents the adjusted power parameter, which is a real number.

Adjusting the adjustment factor f of the ith well logging at the position to be solved according to the change characteristics of the seismic signal of the target stratum_iAnd the power parameter q is adjusted according to the adjustment factor f of the ith well logging at the position to be solved_i'establishing a stratum parameter model of the target stratum with the adjusted power parameter q', so that the established stratum parameter model of the target stratum is matched with the change characteristics of the seismic signals of the target stratum. In specific implementation, the parameter f may need to be adjusted for multiple times_iAnd q, up to the ground of the target stratumThe layer parameter model is matched with the variation characteristics of the seismic signals of the target stratum.

A specific example is given below to illustrate how embodiments of the present invention model formation parameters. This example applies to the zone Station1, where the zone has two logs A and B.

Firstly, acquiring seismic signals of a target stratum and logging data of the target stratum; the logging data of the target stratum comprise a plurality of stratum parameter values of the logging in the target stratum;

and secondly, establishing an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum. The specific process is as follows:

according to the surface coordinates (X) of A and B logs_A,Y_A) And (X)_B,Y_B) Determining a straight line section AB of a continuous well, extracting a seismic section of a target stratum according to the straight line section AB, and determining top and bottom interface seismic horizons of the target stratum, wherein the seismic section of the stratum is obtained by ground seismic exploration through data acquisition, processing and interpretation. And extracting 50 seismic traces from the seismic section: trace₁,Trace₂,…,Trace₅₀Wherein, the seismic record Trace at the two well logging positions of A and B is Trace₁And Trace₂. The ground coordinate of each seismic trace is (X)_n,Y_n) N is 1. ltoreq. n.ltoreq.50, and (X)₁,Y₁)＝(X_A,Y_A)，(X₂,Y₂)＝(X_B,Y_B)。

Determining the Trace of the seismic record of the position to be solved according to the position to be solved_nN is more than or equal to 1 and less than or equal to 50; trace of seismic Trace_nIs determined as the ground coordinates (X) of the location to be determined_n,Y_n). The three-dimensional coordinate representation of the position to be found in the initial formation parameter model is (X)_n,Y_nH), where H represents the depth of layer at the location of the location.

Ground coordinates (X) of the location to be determined_n,Y_n) After the determination, the spatial horizontal distance between the position to be solved and the A and B logs can be determined:

wherein d is₁Representing the spatial horizontal distance between the position to be solved and the A well logging; d₂Representing the spatial horizontal distance between the location to be sought and the B-log.

Given f₁、f₂And the value of q, and f₁Representing the adjustment factor of the A log at the position to be determined, as a real number, f₂Indicating the adjustment factor of the B-log at the desired position, as a real number, e.g. f, may be retrieved for the first time₁＝1，f₂When q is 1 and q is 2, then:

w₁＝a₁/S，w₂＝a₂/S

obtaining the weighted value w of the parameter of the A logging at the position to be solved₁And B logging parameter weighted value w at the position to be determined₂。

Determining 50 seismic Trace traces according to the seismic signals of the target stratum, the top interface seismic horizon and the bottom interface seismic horizon of the target stratum₁,Trace₂,…,Trace₅₀Top interface time top (Trace)_K) K ═ 1,2, …,50 and bottom interface time bot (Trace)_K),K＝1,2,…,50；

Determining a time window [ wt ] according to the following relation₁,wt₂]The range of (A):

wt₁＜min(top(Trace_K))＜max(bot(Trace_K))＜wt₂

wherein, min (Trace)_K) Represents the minimum in time of the top interface of 50 seismic traces; max (Trace)_K) Represents 50 seismic tracesMaximum in bottom interface time.

Obtaining the time window [ wt ] of the A logging according to the logging data of the A logging and the B logging at the target stratum₁,wt₂]Formation parameter sampling data Pw at different sampling depths in a target formation within a range_1,jJ-1, 2, …,100, B log in time window wt₁,wt₂]Formation parameter sampling data Pw at different sampling depths in a target formation within a range_2,jJ equals 1,2, …,100 and the sampling interval ds equals 1 second; where j represents the sample point number at different sample depths in the destination formation. Determining the position moment t of each sampling point in the target stratum of the A logging according to the sampling depth of each sampling point and the seismic signal of the target stratum_1,jAnd B logging the position time t of each sampling point in the target stratum_2,jIn specific implementation, the stratum depth and the time axis of the seismic signal for converting the data have a one-to-one correspondence, and only the value at the time axis of the seismic signal corresponding to the sampling depth of each sampling point needs to be found, namely the position and the time of each sampling point.

And determining the position moment t of the position to be solved in the seismic signal of the target stratum according to the stratum depth H of the position to be solved. According to the position moment t and the time window [ wt ] of the position to be obtained₁,wt₂]Seismic Trace to be located_nTop interface time top (Trace)_n) And the Trace of the seismic record of the position to be solved_nBottom interface time bot (Trace)_n) Determining the proportional position time T of the position to be obtained at the logging position A₁：

If wt₁≤t＜top(Trace_n) And then:

T₁＝wt₁+(top(Trace₁)-wt₁)×(t-wt₁)/(top(Trace_n)-wt₁)；

if top (Trace)_n)≤t＜bot(Trace_n) And then:

if bot (Trace)_n)≤t≤wt₂And then:

T₁＝bot(Trace_n)+(wt₂-bot(Trace₁))×(t-bot(Trace_n))/(wt₂-bot(Trace_n))

wherein, top (Trace)₁) Representing the top interface time of the seismic trace at the well A logging position; bot (Trace)₁) Representing the bottom interface time of the seismic trace at the a-log.

The same method can determine the equal proportion position time T of the position to be solved at the B logging position₂：

If wt₁≤t＜top(Trace_n) And then:

T₂＝wt₁+(top(Trace₂)-wt₁)×(t-wt₁)/(top(Trace_n)-wt₁)；

if top (Trace)_n)≤t＜bot(Trace_n) And then:

if bot (Trace)_n)≤t≤wt₂And then:

T₂＝bot(Trace_n)+(wt₂-bot(Trace₂))×(t-bot(Trace_n))/(wt₂-bot(Trace_n))

wherein, top (Trace)₂) Representing the top interface time of the seismic trace at the B log; bot (Trace)₂) Representing the bottom interface time of the seismic trace at the B log.

Pw obtained by the above calculation_1,j、ds、t_1,jAnd T₁Calculating to obtain the formation parameter value Pw of the A logging at the position to be obtained₁：

Pw₁＝Pw_1,j+(Pw_1,j+1-Pw_1,j)×(T₁-t_1,j)/ds

According to wt₁≤t_1,j≤T₁≤t_1,j+1≤wt₂And determining the value of j by using the relational expression, wherein in order to ensure that j +1 is less than or equal to L, the maximum value of j can only be L-1.

The formation parameter value Pw of the B logging at the position to be solved can be solved in the same way₂：

Pw₂＝Pw_2,j+(Pw_2,j+1-Pw_2,j)×(T₂-t_2,j)/ds

Combining the obtained parameter weighted value w of the A logging at the position to be determined₁And B logging parameter weighted value w at the position to be determined₂Then, the formation parameter value at the position to be solved can be obtained:

M(X_n,Y_n,H)＝w₁×Pw₁+w₂×Pw₂

after the stratum parameter values of all the position points are obtained according to the steps, an initial stratum parameter model of the region Station1 can be established.

And if the established initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum in the region Station1, determining the initial stratum parameter model as the stratum parameter model of the target stratum.

If the established initial stratum parameter model does not conform to the variation characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of the A logging and the B logging at the positions to be solved according to the variation characteristics of the seismic signals of the target stratum, and adjusting an adjusting factor f during specific implementation₁、f₂And the value of the power parameter q according to the adjusted regulating factor f₁'、f₂'establishing a stratum parameter model of the target stratum with the adjusted power parameter q', so that the established stratum parameter model of the target stratum is matched with the change characteristics of the seismic signals of the target stratum.

Another specific application implementation is given below to illustrate how embodiments of the present invention model formation parameters. In the previous example, one more log was added, i.e., the zone had three logs A, B and C.

Under the condition that more than two logging wells exist, the logging wells can be connected according to the order with the prior logging line arrangement direction, and the ground of the well connecting line segment can be extracted according to the orderThe general situation of the seismic profile is shown in fig. 3, the specific situation schematic diagram in this example is shown in fig. 4, the example extracts the seismic profile of the target stratum according to the broken line segment a-C-B, determines the top and bottom interface seismic horizons of the target stratum, and extracts 120 seismic traces in the seismic profile: trace₁,Trace₂,…,Trace₁₂₀Wherein, the seismic record channels at the A, C and B three well logging positions are Trace respectively₁、Trace₂And Trace₃. The ground coordinate of each seismic trace is (X)_n,Y_n) N is 1. ltoreq. n.ltoreq.120, and (X)₁,Y₁)＝(X_A,Y_A)，(X₂,Y₂)＝(X_C,Y_C)，(X₃,Y₃)＝(X_B,Y_B)。

Determining the Trace of the seismic record of the position to be solved according to the position to be solved_nN is more than or equal to 1 and less than or equal to 120; trace of seismic Trace_nIs determined as the ground coordinates (X) of the location to be determined_n,Y_n). The three-dimensional coordinate representation of the position to be found in the initial formation parameter model is (X)_n,Y_nH), where H represents the depth of layer at the location of the location.

Ground coordinates (X) of the location to be determined_n,Y_n) After the determination, the spatial horizontal distance between the position to be sought and the A, C and B logs can be determined:

wherein d is₁Representing the spatial horizontal distance between the position to be solved and the A well logging; d₂Representing the spatial horizontal distance between the position to be solved and the well C; d₃Representing the spatial horizontal distance between the location to be sought and the B-log.

f₁Representing the adjustment factor of the A log at the position to be determined, as a real number, f₂Representing the adjustment factor of the C log at the position to be found, as a real number, f₃Representing the adjustment factor of the B log at the position to be found, as a real number, given f₁、f₂、f₃And the value of q, e.g. f, may be taken for the first time₁＝1，f₂＝1，f₃When q is 1 and q is 2, then:

w₁＝a₁/S，w₂＝a₂/S，w₃＝a₃/S

obtaining the weighted value w of the parameter of the A logging at the position to be solved₁C, weighing the parameter w of the well logging at the position to be solved₂And B logging parameter weighted value w at the position to be determined₃。

Determining 120 seismic Trace traces according to the seismic signals of the target stratum, the top interface seismic horizon and the bottom interface seismic horizon of the target stratum₁,Trace₂,…,Trace₁₂₀Top interface time top (Trace)_K) K ═ 1,2, …,120 and bottom interface time bot (Trace)_K),K＝1,2,…,120；

Determining a time window [ wt ] according to the following relation₁,wt₂]The range of (A):

wt₁＜min(top(Trace_K))＜max(bot(Trace_K))＜wt₂

wherein, min (Trace)_K) Represents the minimum value in time for the top interface of 120 seismic traces; max (Trace)_K) Represents the maximum of the bottom interface times for 120 seismic traces.

Obtaining the time window [ wt ] of the A logging according to the logging data of the A logging, the C logging and the B logging at the target stratum₁,wt₂]Within range at different sampling depths in the formation of interestFormation parameter sampling data Pw_1,jJ-1, 2, …,100, C log in time window wt₁,wt₂]Formation parameter sampling data Pw at different sampling depths in a target formation within a range_2,jJ-1, 2, …,100 and B logs in time window wt₁,wt₂]Formation parameter sampling data Pw at different sampling depths in a target formation within a range_3,jJ equals 1,2, …,100 and the sampling interval ds equals 1 second; where j represents the sample point number at different sample depths in the destination formation. Determining the position moment t of each sampling point in the target stratum of the A logging according to the sampling depth of each sampling point and the seismic signal of the target stratum_1,jC logging position moment t of each sampling point in target stratum_2,jAnd B logging the position time t of each sampling point in the target stratum_3,j。

And determining the position moment t of the position to be solved in the seismic signal of the target stratum according to the stratum depth H of the position to be solved. According to the position moment t and the time window [ wt ] of the position to be obtained₁,wt₂]Seismic Trace to be located_nTop interface time top (Trace)_n) And the Trace of the seismic record of the position to be solved_nBottom interface time bot (Trace)_n) Determining the proportional position time T of the position to be obtained at the logging position A₁：

If wt₁≤t＜top(Trace_n) And then:

T₁＝wt₁+(top(Trace₁)-wt₁)×(t-wt₁)/(top(Trace_n)-wt₁)；

if top (Trace)_n)≤t＜bot(Trace_n) And then:

if bot (Trace)_n)≤t≤wt₂And then:

T₁＝bot(Trace_n)+(wt₂-bot(Trace₁))×(t-bot(Trace_n))/(wt₂-bot(Trace_n))

wherein, top (Trace)₁) Representing the top interface time of the seismic trace at the well A logging position; bot (Trace)₁) Representing the bottom interface time of the seismic trace at the a-log.

Determining the equal proportion position moment T of the position to be solved at the well logging position C₂：

If wt₁≤t＜top(Trace_n) And then:

T₂＝wt₁+(top(Trace₂)-wt₁)×(t-wt₁)/(top(Trace_n)-wt₁)；

if top (Trace)_n)≤t＜bot(Trace_n) And then:

if bot (Trace)_n)≤t≤wt₂And then:

T₂＝bot(Trace_n)+(wt₂-bot(Trace₂))×(t-bot(Trace_n))/(wt₂-bot(Trace_n))

wherein, top (Trace)₂) Representing the top interface time of the seismic trace at the C log; bot (Trace)₂) Representing the bottom interface time of the seismic trace at the C log.

The same method can determine the equal proportion position time T of the position to be solved at the B logging position₃And will not be described herein.

Pw obtained by the above calculation_1,j、ds、t_1,jAnd T₁Calculating to obtain the formation parameter value Pw of the A logging at the position to be obtained₁：

Pw₁＝Pw_1,j+(Pw_1,j+1-Pw_1,j)×(T₁-t_1,j)/ds

According to wt₁≤t_1,j≤T₁≤t_1,j+1≤wt₂Determining the value of j by using a relational expression, wherein in order to ensure that j +1 is less than or equal to L, the maximum value of j can only be obtainedTaking L-1.

In the same way, the formation parameter value Pw of the C logging at the position to be solved can be solved₂：

Pw₂＝Pw_2,j+(Pw_2,j+1-Pw_2,j)×(T₂-t_2,j)/ds

B logging formation parameter value Pw at position to be solved₃：

Pw₃＝Pw_3,j+(Pw_3,j+1-Pw_3,j)×(T₃-t_3,j)/ds

Combining the obtained parameter weighted value w of the A logging at the position to be determined₁C, weighing the parameter w of the well logging at the position to be solved₂And B logging parameter weighted value w at the position to be determined₃Then, the formation parameter value at the position to be solved can be obtained:

M(X_n,Y_n,H)＝w₁×Pw₁+w₂×Pw₂+w₃×Pw₃

and after the stratum parameter values of all the position points are obtained according to the steps, an initial stratum parameter model of the area can be established.

And if the established initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as the stratum parameter model of the target stratum.

If the established initial stratum parameter model does not conform to the variation characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of the A well logging, the C well logging and the B well logging at the positions to be solved according to the variation characteristics of the seismic signals of the target stratum, and adjusting an adjusting factor f during specific implementation₁、f₂、f₃And the value of the power parameter q according to the adjusted regulating factor f₁'、f₂'、f₃'establishing a stratum parameter model of the target stratum with the adjusted power parameter q', so that the established stratum parameter model of the target stratum is matched with the change characteristics of the seismic signals of the target stratum.

Based on the same inventive concept, embodiments of the present invention further provide a formation parameter model building apparatus, and since the principle of the problem solved by the formation parameter model building apparatus is similar to that of the formation parameter model building method, the implementation of the formation parameter model building apparatus may refer to the implementation of the formation parameter model building method, and the repeated parts are not repeated, and the specific structure is as shown in fig. 5:

the data acquisition module 501 is configured to acquire a seismic signal of a target stratum and logging data of the target stratum; the logging data of the target stratum comprises a plurality of stratum parameter values of the target stratum;

an initial model building module 502, configured to build an initial formation parameter model of the target formation according to the seismic signal of the target formation and the logging data of the target formation; the stratum parameter value of the position to be solved in the initial stratum parameter model is determined according to the stratum parameter values of the plurality of logging wells at the position to be solved and the stratum parameter weighted value;

a first model establishing module 503, configured to determine the initial formation parameter model as a formation parameter model of the target formation if the initial formation parameter model matches the change characteristic of the seismic signal of the target formation;

a second model building module 504, configured to adjust stratum parameter weighted values of the multiple logs at the positions to be solved according to the variation characteristic of the seismic signal of the target stratum if the initial stratum parameter model does not match the variation characteristic of the seismic signal of the target stratum, and build a stratum parameter model of the target stratum; and determining the formation parameter value at the position to be solved in the formation parameter model of the target formation according to the formation parameter values of the multiple logging wells at the position to be solved and the adjusted formation parameter weighted value.

In specific implementation, the initial model building module 502 includes: and the stratum parameter value calculation unit is used for taking the stratum parameter weighted values of the multiple logging wells at the positions to be solved as the weights of the stratum parameter values of the multiple logging wells at the positions to be solved, and carrying out weighted summation to obtain the stratum parameter values of the positions to be solved in the initial stratum parameter model.

In a specific embodiment, the formation parameter value calculation unit is specifically configured to:

determining the stratum parameter value at the position to be solved in the initial stratum parameter model according to the stratum parameter values and the stratum parameter weighted values of the positions to be solved of the multiple logs according to the following formula:

M＝w₁×Pw₁+w₂×Pw₂+…+w_R×Pw_R

wherein M represents a stratum parameter value at a position to be solved in the initial stratum parameter model; w is a_iI is 1,2, …, and R represents the weighted value of the parameter of the ith logging well at the position to be solved; pw_i1,2, …, R, representing the formation parameter value of the ith log at the location to be found; r is more than 1 and represents the logging number; the position to be solved is described by a three-dimensional coordinate consisting of a ground coordinate and a stratum depth.

In a specific embodiment, the initial model building module 502 includes: a weight value calculation unit for:

determining the weighted value w of the parameter of the ith well logging at the position to be determined according to the following formula_i：

w_i＝a_i/S,i＝1,2,…,R

Wherein, S ═ a₁+a₂+…+a_R；f_iThe adjustment factor of the ith logging at the position to be solved is represented as a real number; d_iRepresenting the spatial horizontal distance between the position to be solved and the ith logging well; q represents a power parameter, which is a real number. In specific implementation, the second model establishing module 504 is specifically configured to:

adjusting the adjustment factor f of the ith well logging at the position to be solved according to the change characteristics of the seismic signal of the target stratum_iAnd establishing a stratum parameter model of the target stratum according to the power parameter q.

The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the formation parameter model building method when executing the computer program.

An embodiment of the present invention further provides a computer-readable storage medium, in which a computer program for executing the above-mentioned formation parameter model building method is stored.

In summary, the method and the device for establishing the formation parameter model provided in the embodiments of the present invention have the following advantages:

acquiring seismic signals of a target stratum and logging data of the target stratum; establishing an initial stratum parameter model of the target stratum according to the seismic signals of the target stratum and the logging data of the target stratum; if the initial stratum parameter model is matched with the change characteristics of the seismic signals of the target stratum, determining the initial stratum parameter model as a stratum parameter model of the target stratum; if the initial stratum parameter model does not fit the change characteristics of the seismic signals of the target stratum, adjusting stratum parameter weighted values of a plurality of logs at positions to be solved according to the change characteristics of the seismic signals of the target stratum, and establishing a stratum parameter model of the target stratum; when the formation parameter value at each position to be solved is solved, the formation parameter value at each position to be solved is determined according to the formation parameter of each logging at the position to be solved and the formation parameter weighted value, so that the formation parameter value at each position to be solved is related to the formation parameter change between the measurement points; compared with the prior art of establishing a stratum parameter model by an inverse distance weighted interpolation method, the established stratum parameter model can more closely reflect the real stratum parameter distribution structural characteristics of the target stratum by adjusting the stratum parameter weighted values of a plurality of logs at the positions to be solved, the error between the established stratum parameter model and the actual stratum parameters is small, the accuracy of the established stratum parameter model is improved, and the accuracy of the spatial distribution prediction of the underground petroleum and natural gas reservoir stratum is improved.

As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, apparatus, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.

The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.