Method, device, computer equipment and storage medium for determining overburden permeability

Abstract

The present application provides a method, device, computer equipment and storage medium for determining overburden permeability, which belongs to the field of exploration and development of tight sandstone oil and gas reservoirs. The present application obtains the initial permeability of multiple samples and the overburden permeability of multiple samples under different net stress information, so as to determine the dimensionless permeability of multiple samples under different net stress information based on a first relationship calculation model, and then determine the second relationship calculation model corresponding to each sample. Based on multiple second relationship calculation models, it is possible to obtain parameter values ​​of multiple first preset parameters and parameter values ​​of multiple second preset parameters to determine a third relationship calculation model corresponding to the first preset parameters and a fourth relationship calculation model corresponding to the second preset parameters, and then determine a target relationship calculation model. Through the target relationship calculation model, based on the initial permeability and net pressure information of the sample that has not been tested for overburden permeability, the overburden permeability is determined.

Description

Technical Field

[0001] The present application relates to the field of exploration and development of tight sandstone oil and gas reservoirs, and in particular to a method, device, computer equipment and storage medium for determining overburden permeability. Background Art

[0002] The reservoir overburden permeability field is the basic parameter field for numerical simulation of oil and gas reservoirs, and is also the basic parameter in the calculation of single well productivity. The productivity calculation result is the basis for production allocation in the subsequent actual production process. Therefore, the determination of overburden permeability has important guiding significance for the subsequent oil production process.

[0003] At present, when determining the overburden permeability, the overburden permeability of a core under a set of overburden pressure conditions is mainly determined through experiments, and then the relationship calculation model is fitted with the set of overburden permeability values ​​and the corresponding normal pressure permeability, net stress or confining pressure to obtain the coefficient value in the relationship calculation model corresponding to the core participating in the experiment, and then the relationship calculation model for calculating the overburden permeability of the core is obtained, so that the overburden permeability of the core under different overburden pressure conditions can be determined through the relationship calculation model.

[0004] In the above implementation process, the parameters in the relationship calculation model corresponding to the cores with different initial permeabilities are different. Only the cores that have undergone overburden permeability experiments can calculate their overburden permeability under any other overburden conditions within the permitted range of laboratory test conditions according to their corresponding calculation models. For the cores that have not undergone overburden permeability experiments, the coefficients in the relationship calculation model between their overburden permeability values ​​and the corresponding normal pressure permeability, net stress or confining pressure fitting relationship are unknown, and therefore the relationship calculation model of the overburden permeability of this core cannot be determined, that is, the overburden permeability of the core that has not undergone overburden permeability experiments cannot be determined. Summary of the invention

[0005] The embodiment of the present application provides a method, device, computer equipment and storage medium for determining overburden permeability, which can determine the overburden permeability of a core that has not undergone an overburden permeability test. The technical solution is as follows:

[0006] In one aspect, a method for determining overburden permeability is provided, the method comprising:

[0007] Acquire a plurality of net stress information, overburden permeabilities of a plurality of samples under the plurality of net stress information, and initial permeabilities of the plurality of samples;

[0008] Determining the dimensionless permeability of the plurality of samples under the plurality of net pressure information based on the initial permeability, the overburden permeability and a first relationship calculation model, wherein the first relationship calculation model is used to indicate the correlation between the overburden permeability and the initial permeability, and the dimensionless permeability is used to indicate the ratio information of the overburden permeability to the initial permeability;

[0009] Based on the multiple net stress information and the dimensionless permeability under the multiple net stress information, determine multiple second relationship calculation models corresponding to the multiple samples, the second relationship calculation model is used to indicate the correlation between the overburden permeability and the net stress information of the corresponding sample, and the second relationship calculation model includes a first preset parameter and a second preset parameter;

[0010] Based on the multiple second relationship calculation models, determine a third relationship calculation model corresponding to the first preset parameter and a fourth relationship calculation model corresponding to the second preset parameter, the third relationship calculation model is used to indicate the correlation between the first preset parameter and the initial permeability, and the fourth relationship calculation model is used to indicate the correlation between the second preset parameter and the initial permeability;

[0011] Based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, a target relationship calculation model is determined, and the target relationship calculation model is used to determine multiple initial permeabilities and overburden permeabilities under multiple net stress conditions.

[0012] In a possible implementation, the first relationship calculation model is:

[0013]

[0014] Among them, K D is the dimensionless permeability, K is the overburden permeability, K i is the initial permeability.

[0015] In a possible implementation, determining a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information includes:

[0016] For any sample among the plurality of samples, based on the plurality of net stress information and the dimensionless permeability of the any sample under the plurality of net stress information, determining a relationship curve plate corresponding to the any sample, the relationship curve plate being used to represent a change curve of the dimensionless permeability and the net pressure information;

[0017] Based on the relationship curve chart, a second relationship calculation model corresponding to any one of the samples is determined.

[0018] In a possible implementation, the second relationship calculation model is:

[0019] K D =A(p i -p wf ) -B

[0020] Among them, KD is the dimensionless permeability, p i is the confining pressure, p ωf is the core pore pressure, A is the first preset parameter, and B is the second preset parameter.

[0021] In a possible implementation, the third relationship calculation model is:

[0022] A=aln(K i )+b

[0023] Among them, A is the first preset parameter, K i is the initial permeability, a is the first preset coefficient, and b is the second preset coefficient;

[0024] The fourth relationship calculation model is:

[0025] B=cln(K i )+d

[0026] Among them, B is the second preset parameter, K i is the initial permeability, c is the third preset coefficient, and d is the fourth preset coefficient.

[0027] In a possible implementation, the determining the target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model, and the fourth relationship calculation model includes:

[0028] Based on the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, a fifth relationship calculation model is determined, where the fifth relationship calculation model is used to indicate the correlation between the dimensionless permeability and the initial permeability and the net stress information;

[0029] Based on the first relational calculation model and the fifth relational calculation model, the target relational calculation model is determined.

[0030] In a possible implementation, the fifth relationship calculation model is:

[0031]

[0032] Among them, K D is the dimensionless permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient;

[0033] The target relationship calculation model is:

[0034]

[0035] Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

[0036] In one aspect, a device for determining overburden permeability is provided, the device comprising:

[0037] An acquisition module, used for acquiring a plurality of net stress information, overburden permeabilities of a plurality of samples under the plurality of net stress information, and initial permeabilities of the plurality of samples;

[0038] a permeability determination module, for determining the dimensionless permeability of the plurality of samples under the plurality of net pressure information based on the initial permeability, the overburden permeability and a first relationship calculation model, wherein the first relationship calculation model is used to indicate the correlation between the overburden permeability and the initial permeability, and the dimensionless permeability is used to indicate the ratio information of the overburden permeability to the initial permeability;

[0039] A model determination module, for determining a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information, the second relationship calculation model being used to indicate a correlation between the overburden permeability and the net stress information of the corresponding sample, the second relationship calculation model comprising a first preset parameter and a second preset parameter;

[0040] The determination module is further used to determine, based on the multiple second relationship calculation models, a third relationship calculation model corresponding to the first preset parameter and a fourth relationship calculation model corresponding to the second preset parameter, the third relationship calculation data being used to indicate the correlation between the first preset parameter and the initial permeability, and the fourth relationship calculation model being used to indicate the correlation between the second preset parameter and the initial permeability;

[0041] The model determination module is also used to determine a target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model. The target relationship calculation model is used to determine multiple initial permeabilities and overburden permeabilities under multiple net stress conditions.

[0042] In a possible implementation, the first relationship calculation model is:

[0043]

[0044] Among them, K D is the dimensionless permeability, K is the overburden permeability, K i is the initial permeability.

[0045] In one possible implementation, the model determination module is used to determine, for any sample among the multiple samples, a relationship curve plate corresponding to the any sample based on the multiple net stress information and the dimensionless permeability of the any sample under the multiple net stress information, wherein the relationship curve plate is used to represent a variation curve of the dimensionless permeability and the net pressure information; based on the relationship curve plate, determine a second relationship calculation model corresponding to the any sample.

[0046] In a possible implementation, the second relationship calculation model is:

[0047] K D =A( pi -p wf ) -B

[0048] Among them, K D is the dimensionless permeability, p i is the confining pressure, p wf is the core pore pressure, A is the first preset parameter, and B is the second preset parameter.

[0049] In a possible implementation, the third relationship calculation model is:

[0050] A=aln(K i )+b

[0051] Among them, A is the first preset parameter, K i is the initial permeability, a is the first preset coefficient, and b is the second preset coefficient;

[0052] The fourth relationship calculation model is:

[0053] B=cln(K i )+d

[0054] Among them, B is the second preset parameter, K i is the initial permeability, c is the third preset coefficient, and d is the fourth preset coefficient.

[0055] In one possible implementation, the model determination module is used to determine a fifth relationship calculation model based on the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, wherein the fifth relationship calculation model is used to indicate the correlation between the dimensionless permeability and the initial permeability and the net stress information; and to determine the target relationship calculation model based on the first relationship calculation model and the fifth relationship calculation model.

[0056] In a possible implementation, the fifth relationship calculation model is:

[0057]

[0058] Among them, K D is the dimensionless permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient;

[0059] The target relationship calculation model is:

[0060]

[0061] Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

[0062] On the one hand, a computer device is provided, which includes one or more processors and one or more memories, wherein the one or more memories store at least one program code, and the program code is loaded and executed by the one or more processors to implement the operations performed by the method for determining the overburden permeability.

[0063] On the one hand, a computer-readable storage medium is provided, in which at least one program code is stored, and the program code is loaded and executed by a processor to implement the operations performed by the method for determining overburden permeability.

[0064] In one aspect, a computer program product or computer program is provided, the computer program product or computer program comprising computer program code, the computer program code being stored in a computer-readable storage medium. A processor of a computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code to implement the operations performed by the method for determining overburden permeability.

[0065] The solution provided in the present application obtains the initial permeabilities of multiple samples and the overburden permeabilities of multiple samples under different net stress information, so as to determine the dimensionless permeabilities of multiple samples under different overburden pressure conditions based on a first relationship calculation model, thereby determining the second relationship calculation model corresponding to each sample based on the dimensionless permeability of each sample under different net stress information, and obtaining multiple second relationship calculation models, wherein one second relationship calculation model corresponds to a parameter value of a first preset parameter and a parameter value of a second preset parameter, so as to obtain the parameter values ​​of multiple first preset parameters and the parameter values ​​of multiple second preset parameters based on the multiple second relationship calculation models, so as to determine the third relationship calculation model corresponding to the first preset parameter and the fourth relationship calculation model corresponding to the second preset parameter, and then determine the target relationship calculation model, and through the target relationship calculation model, based on the initial permeability and net pressure information of the sample that has not been tested for overburden permeability, the overburden permeability is determined. BRIEF DESCRIPTION OF THE DRAWINGS

[0066] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work.

[0067] Figure 1 is a flow chart of a method for determining overburden permeability provided in an embodiment of the present application;

[0068] Figure 2 is a flow chart of a method for determining overburden permeability provided in an embodiment of the present application;

[0069] Figure 3 It is a dimensionless permeability and net stress relationship curve provided in an embodiment of the present application;

[0070] Figure 4 is a relationship curve between a first preset parameter and an initial permeability provided in an embodiment of the present application;

[0071] Figure 5 is a relationship curve between a second preset parameter and an initial permeability provided in an embodiment of the present application;

[0072] Figure 6 It is a structural schematic diagram of a device for determining overburden permeability provided in an embodiment of the present application;

[0073] Figure 7 It is a structural diagram of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION

[0074] In order to make the objectives, technical solutions and advantages of the present application clearer, the implementation methods of the present application will be further described in detail below with reference to the accompanying drawings.

[0075] The embodiment of the present application provides a method for determining the overburden permeability of a tight sandstone oil and gas reservoir, which is specifically used for permeability field calculation of porous tight sandstone reservoirs and capacity calculation under stress-sensitive conditions. Relevant technicians collect multiple cores from coring wells through equipment, and drill plunger cores of a certain size on the cores as samples to be tested, and then use experimental equipment to determine the initial permeability of the samples and the corresponding overburden permeability under multiple test net stresses, and input the determined initial permeability and the corresponding overburden permeability under multiple test net stresses into a computer device, and input the initial permeability and net pressure information of samples that have not been tested for overburden permeability into a computer device, and determine the overburden permeability of samples that have not been tested for overburden permeability under different net stress information after calculation and processing by the computer device.

[0076] Figure 1 is a flow chart of a method for determining overburden permeability provided in an embodiment of the present application, see Figure 1 , the method comprising:

[0077] 101. A computer device obtains a plurality of net stress information, overburden permeabilities of a plurality of samples under the plurality of net stress information, and initial permeabilities of the plurality of samples.

[0078] 102. The computer device determines the dimensionless permeability of the multiple samples under the multiple net pressure information based on the initial permeability, the overburden permeability and the first relationship calculation model, the first relationship calculation data is used to indicate the correlation between the overburden permeability and the initial permeability, and the dimensionless permeability is used to indicate the ratio information of the overburden permeability to the initial permeability.

[0079] 103. The computer device determines a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information. The second relationship calculation model is used to indicate the correlation between the overburden permeability and the net stress information of the corresponding sample. The second relationship calculation model includes a first preset parameter and a second preset parameter.

[0080] 104. The computer device determines a third relationship calculation model corresponding to the first preset parameter and a fourth relationship calculation model corresponding to the second preset parameter based on the multiple second relationship calculation models, wherein the third relationship calculation model is used to indicate the correlation between the first preset parameter and the initial permeability, and the fourth relationship calculation model is used to indicate the correlation between the second preset parameter and the initial permeability.

[0081] 105. The computer device determines a target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model. The target relationship calculation model is used to determine multiple initial permeabilities and overburden permeabilities under multiple net stress conditions.

[0082] The solution provided by the embodiment of the present application obtains the initial permeabilities of multiple samples and the overburden permeabilities of multiple samples under different net stress information, so as to determine the dimensionless permeabilities of multiple samples under different overburden pressure conditions based on the first relationship calculation model, thereby determining the second relationship calculation model corresponding to each sample based on the dimensionless permeability of each sample under different net stress information, and obtaining multiple second relationship calculation models, wherein one second relationship calculation model corresponds to a parameter value of a first preset parameter and a parameter value of a second preset parameter, so as to obtain the parameter values ​​of multiple first preset parameters and the parameter values ​​of multiple second preset parameters based on the multiple second relationship calculation models, so as to determine the third relationship calculation model corresponding to the first preset parameter and the fourth relationship calculation model corresponding to the second preset parameter, and then determine the target relationship calculation model, and through the target relationship calculation model, based on the initial permeability and net pressure information of the sample that has not been tested for overburden permeability, the overburden permeability is determined.

[0083] In a possible implementation, the first relationship calculation model is:

[0084]

[0085] Among them, K D is the dimensionless permeability, K is the overburden permeability, K i is the initial permeability.

[0086] In a possible implementation, determining a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information includes:

[0087] For any sample among the plurality of samples, based on the plurality of net stress information and the dimensionless permeability of the any sample under the plurality of net stress information, determining a relationship curve plate corresponding to the any sample, the relationship curve plate being used to represent a change curve of the dimensionless permeability and the net pressure information;

[0088] Based on the relationship curve chart, a second relationship calculation model corresponding to any one of the samples is determined.

[0089] In a possible implementation, the second relationship calculation model is:

[0090] K D =A(p i -p wf) -B

[0091] Among them, K D is the dimensionless permeability, p i is the confining pressure, p wf is the core pore pressure, A is the first preset parameter, and B is the second preset parameter.

[0092] In a possible implementation, the third relationship calculation model is:

[0093] A=aln(K i )+b

[0094] Among them, A is the first preset parameter, K i is the initial permeability, a is the first preset coefficient, and b is the second preset coefficient;

[0095] The fourth relationship calculation model is:

[0096] B=cln(K i )+d

[0097] Among them, B is the second preset parameter, K i is the initial permeability, c is the third preset coefficient, and d is the fourth preset coefficient.

[0098] In a possible implementation, the determining the target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model, and the fourth relationship calculation model includes:

[0099] Based on the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, a fifth relationship calculation model is determined, where the fifth relationship calculation model is used to indicate the correlation between the dimensionless permeability and the initial permeability and the net stress information;

[0100] Based on the first relationship calculation data and the fifth relationship calculation model, the target relationship calculation model is determined.

[0101] In a possible implementation, the fifth relationship calculation model is:

[0102]

[0103] Among them, K D is the dimensionless permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient;

[0104] The target relationship calculation model is:

[0105]

[0106] Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

[0107] Figure 2 is a flow chart of a method for determining overburden permeability provided in an embodiment of the present application, see Figure 1 , the method comprising:

[0108] 201. A computer device obtains a plurality of net stress information, overburden permeabilities of a plurality of samples under the plurality of net stress information, and initial permeabilities of the plurality of samples.

[0109] In one possible implementation, relevant technical personnel prepare multiple cores as the multiple samples in accordance with industry standards, and measure the initial permeability of each core, that is, each sample, and then measure the overburden permeability corresponding to each sample under different net stress information in accordance with industry standards, and then input the measured initial permeability of each sample and the overburden permeability corresponding to the sample, as well as the net stress information used to measure the overburden permeability into a computer device, and the computer device obtains the multiple net stress information input by the relevant technical personnel, the overburden permeabilities of the multiple samples under the multiple net stress information, and the initial permeabilities of the multiple samples.

[0110] It should be noted that when determining the initial permeability of the rock and the overburden permeability under different net stress information, the determination is carried out by experimental means, and the specific method of the experiment is not limited in the embodiment of the present application. It should be emphasized that the net stress values ​​corresponding to the multiple net stress information gradually increase from the required minimum net stress value to the upper limit net stress value, and the embodiment of the present application does not limit the values ​​of the minimum net stress value and the upper limit net stress value.

[0111] For example, the initial permeability of the core prepared according to industry standards, that is, the sample, is shown in Table 1:

[0112] Table 1

[0113] Core number Initial permeability (mD) 9 0.309 23 0.417 57 0.447 37 0.524 49 0.529 67 0.757 91 1.07 73 1.42 104 1.43 113 1.79 122 1.93 81 2.11 107 2.13

[0114] The cores, i.e. samples, prepared according to industry standards, have overburden permeabilities measured under different net stress information, see Table 2:

[0115] Table 2

[0116]

[0117] 202. The computer device determines the dimensionless permeability of the multiple samples under the multiple net pressure information based on the initial permeability, the overburden permeability and a first relationship calculation model, wherein the first relationship calculation model is used to indicate the correlation between the overburden permeability and the initial permeability, and the dimensionless permeability is used to indicate the ratio information of the overburden permeability to the initial permeability.

[0118] It should be noted that the first relationship calculation model refers to formula (1):

[0119]

[0120] Among them, K D is the dimensionless permeability, K is the overburden permeability, K i is the initial permeability.

[0121] In a possible implementation, the computer device inputs the acquired initial permeabilities of the multiple samples and the overburden permeabilities of the multiple samples into a first relationship calculation model, and determines the dimensionless permeabilities corresponding to the initial permeability of each sample and the multiple overburden permeabilities of the sample based on the first relationship calculation model. For example, for any sample among the multiple samples, the computer device inputs the multiple overburden permeabilities of the sample under different net stress information and the initial permeability of the sample into the first relationship calculation model, and determines multiple ratio information of the multiple overburden permeabilities of the any sample to the initial permeability based on the first relationship calculation model as the dimensionless permeability of the sample under different net stress information.

[0122] By taking the samples corresponding to the initial permeability and the overburden permeability as a group and inputting the initial permeability and the overburden permeability of the same sample into the first relationship calculation model, the data of each sample can be sorted, and then the dimensionless permeability of each sample under different overburden pressure conditions can be obtained.

[0123] For example, the dimensionless permeability of each sample in Table 2 under different net pressure information is shown in Table 3 below:

[0124] Table 3

[0125]

[0126] 203. The computer device determines a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information. The second relationship calculation model is used to indicate the correlation between the overburden permeability and the net stress information of the corresponding sample. The second relationship calculation model includes a first preset parameter and a second preset parameter.

[0127] In a possible implementation, for any sample among the multiple samples, based on the multiple net stress information and the dimensionless permeability of the any sample under the multiple net stress information, a relationship curve plate corresponding to the any sample is determined, and based on the relationship curve plate, a second relationship calculation model corresponding to the any sample is determined, and the relationship curve plate is used to represent the change curve of the dimensionless permeability and the net pressure information. That is, for each sample, based on the dimensionless permeability under different net stress information and the net stress value indicated by the corresponding net stress information, a relationship curve plate between the dimensionless permeability and the net stress is determined, and then a power function relationship calculation model between the overburden permeability and the net stress is determined, that is, the second relationship calculation model, so as to obtain the second relationship calculation model corresponding to each sample.

[0128] It should be noted that the second relationship calculation model is shown in formula (2):

[0129] K D =A(p i -p wf ) -B (2)

[0130] Among them, K D is the dimensionless permeability, p i is the confining pressure, p wf is the core pore pressure, A is the first preset parameter, and B is the second preset parameter.

[0131] For example, see Figure 3 , Figure 3 It is a relationship curve chart of dimensionless permeability and net stress provided in an embodiment of the present application. The relationship curve chart includes relationship curves of dimensionless permeability and net stress of samples with initial permeabilities of 0.524 millidarcy (mD), 0.757mD, 0.309mD, 1.42mD, 1.43mD, 1.79mD, 1.93mD, and 2.13mD. Through these relationship curves, the second relationship calculation model of samples with initial permeabilities of 0.524 millidarcy (mD), 0.757mD, 0.309mD, 1.42mD, 1.43mD, 1.79mD, 1.93mD, and 2.13mD can be determined.

[0132] 204. The computer device determines a third relationship calculation model corresponding to the first preset parameter and a fourth relationship calculation model corresponding to the second preset parameter based on the multiple second relationship calculation models, wherein the third relationship calculation model is used to indicate the correlation between the first preset parameter and the initial permeability, and the fourth relationship calculation model is used to indicate the correlation between the second preset parameter and the initial permeability.

[0133] It should be noted that one second relationship calculation model corresponds to a parameter value of a first preset parameter and a parameter value of a second preset parameter, and thus the plurality of second relationship calculation models correspond to the parameter values ​​of the plurality of first preset parameters, and the correlation between the first preset parameter and the initial permeability can be determined by the parameter values ​​of the plurality of first preset parameters. The plurality of second relationship calculation models also correspond to the parameter values ​​of the plurality of second preset parameters, and the correlation between the second preset parameter and the initial permeability can be determined by the parameter values ​​of the plurality of second preset parameters.

[0134] For example, Figure 3 Taking the second relational calculation model in as an example, Figure 3 The parameter values ​​of the first preset parameters and the second preset parameters corresponding to each second relationship calculation model in are shown in Table 4 below:

[0135] Table 4

[0136] Core number Initial permeability (mD) First preset parameter Second preset parameter 122 1.93 1.3235 0.375 113 1.79 1.4962 0.448 107 2.13 1.332 0.431 104 1.43 1.4388 0.504 73 1.42 1.739 0.574 67 0.757 1.647 0.716 37 0.524 2.4325 0.95 9 0.309 2.5262 0.961

[0137] In one possible implementation, the computer device determines the third relationship calculation model based on the parameter values ​​of the multiple first preset parameters corresponding to the multiple second relationship calculation models, and determines the fourth relationship calculation model based on the parameter values ​​of the multiple second preset parameters corresponding to the multiple second relationship calculation models.

[0138] It should be noted that the third relationship calculation model refers to formula (3), and the fourth relationship calculation model refers to formula (4):

[0139] A=aln(K i )+b (3)

[0140] Among them, A is the first preset parameter, K i is the initial permeability, a is the first preset coefficient, and b is the second preset coefficient;

[0141] B=cln(K i )+d (4)

[0142] Among them, B is the second preset parameter, K i is the initial permeability, c is the third preset coefficient, and d is the fourth preset coefficient.

[0143] It should be noted that when determining the third relationship calculation model according to the parameter values ​​of the multiple first preset parameters corresponding to the multiple second relationship calculation models, the computer device determines the values ​​of the first preset coefficient and the second preset coefficient in the third relationship calculation model based on the parameter values ​​of the multiple first preset parameters and the initial permeability of each sample. Still taking the parameter values ​​of the multiple first preset parameters in Table 4 as an example, based on the parameter values ​​of the multiple first preset parameters and the initial permeability of each sample, the relationship curve plate of the third relationship calculation model is determined, see Figure 4 , Figure 4 is a relationship curve between a first preset parameter and an initial permeability provided in an embodiment of the present application. Through the relationship curve, it can be determined that the value of the first preset coefficient is -0.637, and the value of the second preset coefficient is 1.7894. The third relationship calculation model is shown in formula (5):

[0144] A=-0.637ln(K i )+1.7894 (5)

[0145] Among them, A is the first preset parameter, K i is the initial permeability.

[0146] Taking the parameter values ​​of the second preset parameters in Table 4 as an example, based on the parameter values ​​of the second preset parameters and the initial permeability of each sample, the relationship curve chart of the fourth relationship calculation model is determined, see Figure 5 , Figure 5 is a relationship curve between a second preset parameter and an initial permeability provided in an embodiment of the present application. Through the relationship curve, it can be determined that the value of the third preset coefficient is -0.332, and the value of the fourth preset coefficient is 0.654. ​​The fourth relationship calculation model is shown in formula (6):

[0147] B=-0.322ln(K i )+0.654 (6)

[0148] Among them, B is the second preset parameter, K i is the initial permeability.

[0149] 205. The computer device determines a fifth relationship calculation model based on the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model. The fifth relationship calculation model is used to indicate the correlation between the dimensionless permeability and the initial permeability and the net stress information.

[0150] In a possible implementation, the computer device inputs the third relational calculation model and the fourth relational calculation model into the second relational calculation model, and outputs the fifth relational calculation model.

[0151] It should be noted that the fifth relationship calculation model is shown in formula (7):

[0152]

[0153] Among them, K D is the dimensionless permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, p i -p wf is the net stress, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

[0154] 206. The computer device determines a target relationship calculation model based on the first relationship calculation model and the fifth relationship calculation model, where the target relationship calculation model is used to determine multiple initial permeabilities and overburden permeabilities under multiple net stress conditions.

[0155] In a possible implementation, the computer device inputs the first relational calculation model into the fifth relational calculation model, and outputs the target relational calculation model.

[0156] It should be noted that the target relationship calculation model is shown in formula (8):

[0157]

[0158] Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, p i -p wf is the net stress, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

[0159] Taking the third relationship calculation model as shown in formula (5) and the fourth relationship calculation model as shown in formula (6) as an example, the target relationship calculation model is shown in formula (9):

[0160]

[0161] Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, p i -p wf is the net stress.

[0162] It should be noted that the relevant technicians input the initial permeability of the sample that has not been tested for overburden permeability and the net stress value corresponding to the net stress information that needs to be tested into the computer device, and the computer device obtains the initial permeability and net stress value input by the relevant technicians, inputs the obtained initial permeability and net stress value into the target relationship calculation model, and outputs the overburden permeability determined by the target relationship calculation model, so as to obtain the overburden permeability of the sample under the corresponding net pressure information. For example, when it is necessary to determine the overburden permeability of any sample under any net stress information, the initial permeability of any sample and the net stress value corresponding to the any net stress information are input into the target relationship calculation model, and the overburden permeability of any sample under the any net stress information can be determined.

[0163] It should be emphasized that the input initial permeability and net stress values ​​must be within the permissible range of laboratory test conditions to ensure the accuracy of the determined overburden permeability.

[0164] For example, taking the target relationship calculation model as shown in formula (9) as an example, the overburden permeability of each sample that has not been tested for overburden permeability under different net stress information is shown in Table 5 below:

[0165] Table 5

[0166]

[0167] The solution provided by the embodiment of the present application obtains the initial permeability of multiple samples and the overburden permeability of multiple samples under different net stress information, so as to determine the dimensionless permeability of multiple samples under different overburden conditions based on the first relationship calculation model, thereby determining the second relationship calculation model corresponding to each sample based on the dimensionless permeability of each sample under different net stress information, and obtaining multiple second relationship calculation models, one second relationship calculation model corresponds to a parameter value of a first preset parameter and a parameter value of a second preset parameter, so as to obtain the parameter values ​​of multiple first preset parameters and the parameter values ​​of multiple second preset parameters based on multiple second relationship calculation models, so as to determine the third relationship calculation model corresponding to the first preset parameter and the fourth relationship calculation model corresponding to the second preset parameter, and then determine the target relationship calculation model, and through the target relationship calculation model, based on the initial permeability and net pressure information of the sample without overburden permeability test, the overburden permeability is determined. According to the method combining physical experiments, gas reservoir engineering and higher mathematics, the overburden permeability is determined by computer equipment, and the manual processing process is reduced, thereby improving the efficiency of overburden permeability determination.

[0168] All the above optional technical solutions can be arbitrarily combined to form optional embodiments of the present application, which will not be described one by one here.

[0169] Figure 6 is a schematic diagram of a device for determining overburden permeability provided in an embodiment of the present application, see Figure 6 , the device comprises:

[0170] An acquisition module 601 is used to acquire a plurality of net stress information, overburden permeabilities of a plurality of samples under the plurality of net stress information, and initial permeabilities of the plurality of samples;

[0171] A permeability determination module 602 is used to determine the dimensionless permeability of the multiple samples under the multiple net pressure information based on the initial permeability, the overburden permeability and a first relationship calculation model, wherein the first relationship calculation model is used to indicate the correlation between the overburden permeability and the initial permeability, and the dimensionless permeability is used to indicate the ratio information of the overburden permeability to the initial permeability;

[0172] A model determination module 603 is used to determine a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information, wherein the second relationship calculation model is used to indicate the correlation between the overburden permeability and the net stress information of the corresponding sample, and the second relationship calculation model includes a first preset parameter and a second preset parameter;

[0173] The model determination module 603 is further used to determine, based on the multiple second relationship calculation models, a third relationship calculation model corresponding to the first preset parameter and a fourth relationship calculation model corresponding to the second preset parameter, wherein the third relationship calculation model is used to indicate the correlation between the first preset parameter and the initial permeability, and the fourth relationship calculation model is used to indicate the correlation between the second preset parameter and the initial permeability;

[0174] The model determination module 603 is also used to determine a target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model. The target relationship calculation model is used to determine multiple initial permeabilities and overburden permeabilities under multiple net stress conditions.

[0175] The device provided in the embodiment of the present application obtains the initial permeabilities of multiple samples and the overburden permeabilities of multiple samples under different net stress information, so as to determine the dimensionless permeabilities of multiple samples under different overburden pressure conditions based on the first relationship calculation model, thereby determining the second relationship calculation model corresponding to each sample based on the dimensionless permeability of each sample under different net stress information, and obtaining multiple second relationship calculation models, wherein one second relationship calculation model corresponds to a parameter value of a first preset parameter and a parameter value of a second preset parameter, so as to obtain the parameter values ​​of multiple first preset parameters and the parameter values ​​of multiple second preset parameters based on the multiple second relationship calculation models, so as to determine the third relationship calculation model corresponding to the first preset parameter and the fourth relationship calculation model corresponding to the second preset parameter, and then determine the target relationship calculation model, and through the target relationship calculation model, based on the initial permeability and net pressure information of the sample that has not been tested for overburden permeability, the overburden permeability is determined.

[0176] In a possible implementation, the first relationship calculation model is:

[0177]

[0178] Among them, K D is the dimensionless permeability, K is the overburden permeability, K i is the initial permeability.

[0179] In a possible implementation, the model determination module 603 is used to determine, for any sample among the multiple samples, a relationship curve plate corresponding to the any sample based on the multiple net stress information and the dimensionless permeability of the any sample under the multiple net stress information, wherein the relationship curve plate is used to represent a curve of changes in the dimensionless permeability and the net pressure information; based on the relationship curve plate, determine a second relationship calculation model corresponding to the any sample.

[0180] In a possible implementation, the second relationship calculation model is:

[0181] K D =A(p i -p wf ) -B

[0182] Among them, K D is the dimensionless permeability, p i is the confining pressure, p wf is the core pore pressure, A is the first preset parameter, and B is the second preset parameter.

[0183] In a possible implementation, the third relationship calculation model is:

[0184] A=aln(Ki )+b

[0185] Among them, A is the first preset parameter, K i is the initial permeability, a is the first preset coefficient, and b is the second preset coefficient;

[0186] The fourth relationship calculation model is:

[0187] B=cln(K i )+d

[0188] Among them, B is the second preset parameter, K i is the initial permeability, c is the third preset coefficient, and d is the fourth preset coefficient.

[0189] In one possible implementation, the model determination module 603 is used to determine a fifth relationship calculation model based on the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, wherein the fifth relationship calculation model is used to indicate the correlation between the dimensionless permeability and the initial permeability and the net stress information; and to determine the target relationship calculation model based on the first relationship calculation model and the fifth relationship calculation model.

[0190] In a possible implementation, the fifth relationship calculation model is:

[0191]

[0192] Among them, K D is the dimensionless permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient;

[0193] The target relationship calculation model is:

[0194]

[0195] Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

[0196] It should be noted that: the overburden permeability determination device provided in the above embodiment only uses the division of the above functional modules as an example to illustrate when determining the overburden permeability of the reservoir. In actual applications, the above functions can be assigned to different functional modules as needed, that is, the internal structure of the computer device is divided into different functional modules to complete all or part of the functions described above. In addition, the overburden permeability determination device provided in the above embodiment and the overburden permeability determination method embodiment belong to the same concept, and the specific implementation process is detailed in the method embodiment, which will not be repeated here.

[0197] Figure 7 It is a structural diagram of a computer device provided in an embodiment of the present application. The computer device 700 may have relatively large differences due to different configurations or performances, and may include one or more processors (Central Processing Units, CPU) 701 and one or more memories 702, wherein the one or more memories 702 store at least one program code, and the at least one program code is loaded and executed by the one or more processors 701 to implement the methods provided in the above-mentioned various method embodiments. Of course, the computer device 700 may also have components such as a wired or wireless network interface, a keyboard, and an input and output interface for input and output. The computer device 700 may also include other components for implementing device functions, which will not be described in detail here.

[0198] In an exemplary embodiment, a computer-readable storage medium is also provided, such as a memory including a program code, and the program code can be executed by a processor to complete the method for determining the overburden permeability in the above embodiment. For example, the computer-readable storage medium can be a read-only memory (ROM), a random access memory (RAM), a compact disc read-only memory (CD-ROM), a magnetic tape, a floppy disk, an optical data storage device, etc.

[0199] In an exemplary embodiment, a computer program product or a computer program is also provided, which includes a computer program code, which is stored in a computer-readable storage medium. A processor of a computer device reads the computer program code from the computer-readable storage medium, and the processor executes the computer program code, so that the computer device executes the method steps of the method for determining the overburden permeability provided in the above-mentioned embodiment.

[0200] A person skilled in the art will appreciate that all or part of the steps for implementing the above embodiments may be accomplished by hardware or by hardware associated with a program code, and the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a disk or an optical disk, etc.

[0201] The above are only optional embodiments of the present application and are not intended to limit the present application. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.

Claims

1. A method for determining overburden permeability, characterized in that: The method comprises: Acquire a plurality of net stress information, overburden permeabilities of a plurality of samples under the plurality of net stress information, and initial permeabilities of the plurality of samples; Determining the dimensionless permeability of the plurality of samples under the plurality of net stress information based on the initial permeability, the overburden permeability and a first relationship calculation model, wherein the first relationship calculation model is used to indicate the correlation between the overburden permeability and the initial permeability, and the dimensionless permeability is used to indicate the ratio information of the overburden permeability to the initial permeability; Based on the multiple net stress information and the dimensionless permeability under the multiple net stress information, determining multiple second relationship calculation models corresponding to the multiple samples, the second relationship calculation model is used to indicate the correlation calculation between the overburden permeability and the net stress information of the corresponding sample, and the second relationship calculation model includes a first preset parameter and a second preset parameter; Based on the multiple second relationship calculation models, determine a third relationship calculation model corresponding to the first preset parameter and a fourth relationship calculation model corresponding to the second preset parameter, the third relationship calculation model is used to indicate the correlation between the first preset parameter and the initial permeability, and the fourth relationship calculation model is used to indicate the correlation between the second preset parameter and the initial permeability; Determine a target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, wherein the target relationship calculation model is used to determine multiple initial permeabilities and overburden permeabilities under multiple net stress conditions; Wherein, the first relationship calculation model is: Among them, K D is the dimensionless permeability, K is the overburden permeability, K i is the initial permeability; the second relationship calculation model is: K D =A(p i -p wf ) -B Among them, K D is the dimensionless permeability, p i is the confining pressure, p wf is the core pore pressure, A is the first preset parameter, and B is the second preset parameter; The third relationship calculation model is: A=a ln(K i )+b Among them, A is the first preset parameter, K i is the initial permeability, a is the first preset coefficient, and b is the second preset coefficient; The fourth relationship calculation model is: B=c ln(K i )+d Among them, B is the second preset parameter, K i is the initial permeability, c is the third preset coefficient, and d is the fourth preset coefficient; The target relationship calculation model is: Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

2. The method according to claim 1, characterized in that: Determining a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information includes: For any sample among the plurality of samples, based on the plurality of net stress information and the dimensionless permeability of the any sample under the plurality of net stress information, determining a relationship curve plate corresponding to the any sample, wherein the relationship curve plate is used to represent a change curve of the dimensionless permeability and the net stress information; Based on the relationship curve plate, a second relationship calculation model corresponding to any one of the samples is determined.

3. The method according to claim 1, characterized in that The determining of the target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model, and the fourth relationship calculation model includes: Based on the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, determining a fifth relationship calculation model, wherein the fifth relationship calculation model is used to indicate the correlation between the dimensionless permeability and the initial permeability and the net stress information; Based on the first relationship calculation model and the fifth relationship calculation model, the target relationship calculation model is determined.

4. The method according to claim 3, characterized in that The fifth relationship calculation model is: Among them, K D is the dimensionless permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

5. A device for determining overburden permeability, characterized in that: The device comprises: An acquisition module, used for acquiring a plurality of net stress information, overburden permeabilities of a plurality of samples under the plurality of net stress information, and initial permeabilities of the plurality of samples; a permeability determination module, configured to determine the dimensionless permeabilities of the plurality of samples under the plurality of net stress information based on the initial permeability, the overburden permeability, and a first relationship calculation model, wherein the first relationship calculation model is configured to indicate a correlation between the overburden permeability and the initial permeability, and the dimensionless permeability is configured to indicate information on a ratio of the overburden permeability to the initial permeability; A model determination module, used to determine a plurality of second relationship calculation models corresponding to the plurality of samples based on the plurality of net stress information and the dimensionless permeability under the plurality of net stress information, wherein the second relationship calculation model is used to indicate the calculation of the correlation relationship between the overburden permeability and the net stress information of the corresponding sample, and the second relationship calculation model includes a first preset parameter and a second preset parameter; The model determination module is further used to determine, based on the multiple second relationship calculation models, a third relationship calculation model corresponding to the first preset parameter and a fourth relationship calculation model corresponding to the second preset parameter, wherein the third relationship calculation model is used to indicate the correlation between the first preset parameter and the initial permeability, and the fourth relationship calculation model is used to indicate the correlation between the second preset parameter and the initial permeability; The model determination module is further used to determine a target relationship calculation model based on the first relationship calculation model, the second relationship calculation model, the third relationship calculation model and the fourth relationship calculation model, wherein the target relationship calculation model is used to determine multiple initial permeabilities and overburden permeabilities under multiple net stress conditions; The second relationship calculation model is: K D =A(p i -p wf ) -B Among them, K D is the dimensionless permeability, p i is the confining pressure, p wf is the core pore pressure, A is the first preset parameter, and B is the second preset parameter; The third relationship calculation model is: A=a ln(K i )+b Among them, A is the first preset parameter, K i is the initial permeability, a is the first preset coefficient, and b is the second preset coefficient; The fourth relationship calculation model is: B=c ln(K i )+d Among them, B is the second preset parameter, K i is the initial permeability, c is the third preset coefficient, and d is the fourth preset coefficient; The target relationship calculation model is: Where K is the overburden permeability, K i is the initial permeability, p i is the confining pressure, p wf is the core pore pressure, a is the first preset coefficient, b is the second preset coefficient, c is the third preset coefficient, and d is the fourth preset coefficient.

6. A computer device, characterized in that: The computer device includes one or more processors and one or more memories, wherein at least one program code is stored in the one or more memories, and the program code is loaded and executed by the one or more processors to implement the operations performed by the method for determining the overburden permeability as described in any one of claims 1 to 4.

7. A computer-readable storage medium, characterized in that: At least one program code is stored in the computer-readable storage medium, and the program code is loaded and executed by the processor to implement the operations performed by the method for determining overburden permeability according to any one of claims 1 to 4.

Images