Method and device for determining long-term conductivity of acid-etched fractures

Abstract

The application provides a method and device for determining long-term conductivity of acid-etched fractures in carbonate rocks. The method comprises: obtaining rock mineral composition, porosity, permeability, formation Young's modulus, formation closure stress, formation Poisson's ratio, acid concentration and temperature; determining the damage factor of acid-etched fractures based on rock mineral composition, porosity, permeability, Young's modulus, acid concentration and temperature; determining the change in acid-etched fracture width after acid etching of the fracture surface based on formation Young's modulus, closure stress, Poisson's ratio and acid-rock reaction rate; determining the creep variable of acid-etched fracture conductivity based on the initial fracture width, Young's modulus, Poisson's ratio and damage factor; and finally, determining the long-term conductivity by the fracture width after acid etching and the formation creep variable. The method comprehensively considers lithology, formation physical properties, rock mechanics parameters, acid etching damage and creep effect, and can quickly and accurately determine the conductivity of acid-etched fractures.

Description

TECHNICAL FIELD

[0001] The application belongs to the technical field of oil exploitation, and particularly relates to a method and device for determining long-term conductivity of acid-etched fractures. BACKGROUND

[0002] Carbonate reservoirs have large oil and gas reserves, about 48% of global crude oil reserves and 28% of natural gas reserves, and high single-well production, and thus have attracted increasing attention in the world oil and gas industry. Acid fracturing is one of the key technologies for reservoir reconstruction. The fracture morphology, rock and mineral content, porosity and permeability conditions, and rock mechanics parameters after acid etching all affect the conductivity of acid-etched fractures, and the conductivity of acid-etched fractures has become one of the important indicators for evaluating acid fracturing.

[0003] The existing methods for characterizing the space of acid-etched fractures have defects, and more are equivalent to the characteristics of the space of fractures, such as roughness, fractal dimension, kurtosis, tortuosity, etc., by using the profile characteristics of a single side of the fracture or a single side of the fracture. Although the closure degree can represent the closure characteristics of both sides of the fracture, it cannot reflect the morphology of the fracture space. At present, the acid-etched fracture conductivity models used to determine the conductivity of acid-etched fractures are mostly empirical models or improved models, and there is a common problem that the determined results of the conductivity of acid-etched fractures are inaccurate. The space morphology of acid-etched fractures, rock mechanics parameters, closure stress, and fracture contact state can all affect the conductivity of acid-etched fractures. At present, technical solutions for accurately determining the conductivity of acid-etched fractures still need to be researched.

[0004] The acid-etched fractures generated in the acid fracturing process are affected by multiple factors such as closure stress, rock creep, and acid etching damage, and the fracture conductivity is difficult to maintain for a long time, which seriously restricts the development effect. The existing prediction methods fail to comprehensively consider the acid etching damage and creep effect, resulting in a large deviation in the evaluation of long-term conductivity. SUMMARY

[0005] The purpose of the present application is to overcome the problems existing in the prior art, and to provide an innovative long-term conductivity prediction model that comprehensively considers the acid etching damage and creep effect of carbonate rocks and quantifies the dynamic change law of fracture width with time. The model can accurately predict the long-term conductivity of fractures under different acid fracturing processes, providing a scientific basis for optimizing the development and design of carbonate reservoirs, and has important theoretical value and engineering significance.

[0006] In order to achieve the above purpose, the present application provides the following five technical solutions.

[0007] In a first aspect, the present application provides a method for determining the long-term conductivity of acid-etched fractures, which comprises:

[0008] acquire rock mineral composition, porosity, permeability, formation Young's modulus, formation closure stress, formation Poisson's ratio, acid concentration used for reconstruction, temperature of the work area;

[0009] determine acid-etched fracture damage factor of the work area based on rock mineral composition, porosity, permeability, formation Young's modulus, acid concentration, and temperature;

[0010] determine acid-etched fracture face acid-etched fracture width change of the work area based on formation Young's modulus, formation closure stress, formation Poisson's ratio, and formation acid-rock reaction rate;

[0011] determine acid-etched fracture width creep amount of the work area based on acid-etched fracture width after acid etching, formation Young's modulus, Poisson's ratio, and acid-etched fracture damage factor of the work area;

[0012] determine acid-etched fracture conductivity of the work area based on acid-etched fracture width after acid etching and fracture width after creep of the work area;

[0013] In a second aspect, the present application provides a long-term acid-etched fracture conductivity determination device, which comprises:

[0014] The data acquisition module is configured to acquire rock mineral composition, porosity, permeability, formation Young's modulus, formation closure stress, formation Poisson's ratio, acid concentration used for reconstruction, and temperature of the work area.

[0015] The acid-etched damage determination module is configured to determine acid-etched fracture damage factor of the work area based on rock mineral composition, porosity, permeability, mechanical properties, acid concentration used for reconstruction, and temperature of the work area.

[0016] The acid-etched fracture width determination module is configured to determine acid-etched fracture width of the work area based on formation Young's modulus, formation closure stress, formation Poisson's ratio, and formation acid-rock reaction rate.

[0017] The fracture width after creep determination module is configured to determine acid-etched fracture width after creep of the work area based on acid-etched fracture damage factor and acid-etched fracture width.

[0018] The conductivity determination module is configured to determine long-term acid-etched fracture conductivity of the work area based on acid-etched fracture width after acid etching and fracture width after creep.

[0019] In a third aspect, the present application provides a computer device, which comprises a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the acid-etched fracture long-term conductivity determination method provided in the first aspect when executing the computer program.

[0020] In a fourth aspect, the present application provides a computer readable storage medium, which stores a computer program, and the computer program, when executed by a processor, implements the acid-etched fracture long-term conductivity determination method provided in the first aspect.

[0021] In a fifth aspect, the present application provides a computer program product, which comprises a computer program, and the computer program, when executed by a processor, implements the acid-etched fracture long-term conductivity determination method provided in the first aspect.

[0022] Advantages and beneficial effects of the present application:

[0023] The present application proposes an innovative long-term conductivity prediction model, which comprehensively considers the carbonate rock acid-etching damage and creep effect, and quantifies the dynamic change law of fracture width with time. The model can accurately predict the long-term conductivity of fractures under different acid fracturing processes, provides a scientific basis for optimizing the development and design of carbonate reservoirs, and has important theoretical value and engineering significance. BRIEF DESCRIPTION OF DRAWINGS

[0024] Figure 1 is a fracture width deformation-closure stress relationship diagram;

[0025] Figure 2 is a fracture width deformation-reaction time relationship diagram;

[0026] Figure 3 is a fracture width deformation-time relationship diagram;

[0027] Figure 4 is a conductivity-time relationship diagram;

[0028] Figure 5 is an acid-etched fracture long-term conductivity determination device structure diagram;

[0029] Figure 6 is an acid-etching damage determination module structure diagram;

[0030] Figure 7 is an acid-etched fracture width determination module structure diagram. DETAILED DESCRIPTION

[0031] In order to make the purpose, technical scheme and advantages of the embodiments of the present application clearer, further detailed description of the embodiments of the present application will be given below with reference to the accompanying drawings. Herein, the illustrative embodiments of the present application and their descriptions are used to explain the present application, but not as a limitation of the present application.

[0032] In the description of the present specification, "include", "includes", "have", "has", and the like are open terms, that is, mean including but not limited to. The description referring to the terms "one embodiment", "one specific embodiment", "some embodiments", "for example", and the like means that the specific features, structures, or characteristics described in connection with the embodiment or example are included in at least one embodiment or example of the present application. In the present specification, the illustrative description of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, or characteristics described can be combined in any one or more embodiments or examples in a suitable manner. The order of the steps involved in each embodiment is used to illustrate the implementation of the present application, and the order of the steps is not limited and can be appropriately adjusted as needed.

[0033] The embodiment of the present application provides a method for determining long-term conductivity of acid-etched cracks, which comprises the following steps:

[0034] Step 101, obtaining rock mineral composition, porosity, permeability, formation Young's modulus, formation closure stress, formation Poisson's ratio, acid concentration used for reconstruction, and temperature of a work area;

[0035] Step 102, determining an acid-etched crack damage factor of the work area based on the rock mineral composition, the porosity, the permeability, the formation Young's modulus, the acid concentration, and the temperature;

[0036] Step 103, determining a change in crack width after acid etching of a crack surface of the acid-etched crack of the work area based on the formation Young's modulus, the formation closure stress, the formation Poisson's ratio, and an acid-rock reaction rate of the formation;

[0037] Step 104, determining a crack width creep amount of the acid-etched crack of the work area based on the crack width after acid etching of the acid-etched crack of the work area, the formation Young's modulus, the Poisson's ratio, and the acid-etched crack damage factor;

[0038] Step 105, determining a conductivity of the acid-etched crack of the work area based on the crack width after acid etching of the acid-etched crack of the work area and the crack width after creep.

[0039] In one embodiment, in step 102, the acid-etched crack damage factor of the work area is determined based on the rock mineral composition, the porosity, the permeability, the formation Young's modulus, the acid concentration, and the temperature, and comprises the following steps:

[0040] Step 1021, determining an acid-rock reaction rate of the work area by using an acid-rock reaction kinetics model based on the acid concentration used for reconstruction and the temperature of the work area;

[0041] Step 1022, determining a porosity change of the work area by using an acid-rock reaction rate and porosity change relationship model based on the acid-rock reaction rate of the work area;

[0042] Step 1023, based on the acid-rock reaction rate of the work area, the acid-rock reaction rate and the Young's modulus change relationship model is used to determine the Young's modulus change of the work area.

[0043] Step 1024, based on the acid-rock reaction rate of the work area, the porosity change, the Young's modulus change, the acid-rock reaction rate, the porosity change, the Young's modulus change relationship and the acid etching damage relationship model are used to determine the acid etching crack damage factor of the work area.

[0044] In one embodiment, in step 1021, the acid-rock reaction kinetics model is:

[0045]

[0046] In the formula, T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid concentration, mol / L; A, b, m are constants, dimensionless;

[0047] For example, the relationship model between the acid-rock reaction rate and the acid concentration and temperature used in the reconstruction is:

[0048]

[0049] In the formula, T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid concentration, mol / L.

[0050] In one embodiment, in step 1022, the acid-rock reaction rate and porosity change relationship model is:

[0051]

[0052] In the formula, The porosity change is %, T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid concentration, mol / L; t is the reaction time, H; S is the reaction area, m 2 ; A, b, m are constants, dimensionless;

[0053] For example, the acid-rock reaction rate and porosity change relationship model is:

[0054]

[0055] In the formula, The porosity change is %, T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid concentration, mol / L; t is the reaction time, H; S is the reaction area, m 2 ;

[0056] In one embodiment, in step 1023, the acid-rock reaction rate and Young's modulus change relationship model is:

[0057]

[0058] wherein E is the Young's modulus after acid etching, GPa; E0 is the initial Young's modulus, GPa; T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid liquid concentration, mol / L; t is the reaction time, H; S is the reaction area, m c 2 ; A, b, m, n are constants, dimensionless;

[0059] For example, the acid rock reaction rate and Young's modulus change relationship model is:

[0060]

[0061] wherein E is the Young's modulus after acid etching, GPa; E0 is the initial Young's modulus, GPa; T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid liquid concentration, mol / L; t is the reaction time, H; S is the reaction area, m c 2 ;

[0062] In one embodiment, in step 1024, the acid rock reaction rate, porosity change, Young's modulus change relationship and acid etching damage relationship model is:

[0063]

[0064] wherein T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid liquid concentration, mol / L; t is the reaction time, H; S is the reaction area, m 2 ; A, b, m, n are constants, dimensionless;

[0065] For example, the acid rock reaction rate, porosity change, Young's modulus change relationship and acid etching damage relationship model is:

[0066]

[0067] wherein T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid liquid concentration, mol / L; t is the reaction time, H; S is the reaction area, m 2 ;

[0068] In one embodiment, in step 103, the acid etching fracture surface acid etching fracture width change in the working area is determined based on the formation Young's modulus, the formation closure stress, the formation Poisson's ratio and the formation acid rock reaction rate, which includes:

[0069] Step 1031, based on the formation Young's modulus, the formation closure stress and the formation Poisson's ratio in the working area, the closure stress fracture width in the working area is determined.

[0070] ​​In step 1032, based on the acid-rock reaction rate of the work area, the acid-rock reaction crack width of the work area is determined by using the acid-rock reaction rate and the crack width relationship model.

[0071] In step 1033, based on the closure stress crack width, the acid-rock reaction crack width and the acid-etched crack width relationship model, the acid-etched crack width is determined.

[0072] In one embodiment, in step 1031, the relationship model between the closure stress crack width and the formation Young's modulus, the formation closure stress and the formation Poisson's ratio is:

[0073]

[0074] In the formula, N is the number of micro-convex units per unit area of acid-etched crack, dimensionless; w is the initial crack width, m; R is the micro-convex radius, m; E is the Young's modulus, GPa; v is the Poisson's ratio, dimensionless; D is the acid-etched crack damage factor, dimensionless; σ1 is the closure stress, MPa; μ is the average value of micro-convex height distribution, m; σ is the standard deviation of micro-convex height distribution, m; f is the friction coefficient, dimensionless; is the cumulative distribution value of the micro-convex distribution function; dimensionless.

[0075] In one embodiment, in step 1032, the relationship model between the acid-rock reaction crack width and the acid-rock reaction rate is:

[0076]

[0077] In the formula, J is the acid-rock reaction rate, g / (cm 2 ·s); t is the reaction time, h; ρ is the rock density, kg / m 3 ; S is the reaction area; m 2 .

[0078] In one embodiment, in step 1033, the relationship model between the acid-etched crack width and the closure stress crack width and the acid-rock reaction crack width is:

[0079] w = w0-w(σ)-w(J)

[0080] In the formula, w is the acid-etched crack width, m; w0 is the initial crack width, m; w(σ) is the closure stress crack width, m; w(J) is the acid-rock reaction crack width, m.

[0081] In one embodiment, in step 104, the creep crack width determination module determines the post-creep crack width of the work area by the following formula:

[0082]

[0083] wherein w0 is the initial fracture width, m; σ1 is the closure stress, MPa; σ3 is the confining pressure, MPa; K is the bulk modulus, GPa; G is the shear modulus, GPa; B is the bulk parameter, dimensionless; β is the material homogeneity parameter, dimensionless; η is the viscosity parameter, dimensionless; α is the loading level strength parameter, dimensionless; t σ is the loading time, H; D is the acid-etched fracture damage factor, dimensionless; t B is the accelerated creep time, H; t f is the collapse time, H.

[0084] In one embodiment, in step 105, the fracture conductivity determination module determines the acid-etched fracture conductivity of the work area by the following formula:

[0085]

[0086] wherein α is the acid-etched fracture roughness correction coefficient, dimensionless; w is the fracture width after acid etching, m; W kf is the acid-etched fracture conductivity, μm 2 ·cm.

[0087] In one embodiment, the rock mineral composition of the work area is determined by an X-ray diffraction experiment.

[0088] In one embodiment, the porosity and permeability of the work area are determined by a porosity and permeability determination experiment.

[0089] In one embodiment, the formation Young's modulus, the formation closure stress, and the formation Poisson's ratio are determined by a conventional Young's modulus, closure stress, and Poisson's ratio test, which is not limited herein.

[0090] The embodiment of the present application also provides an acid-etched fracture long-term conductivity determination device, as described in the following embodiment. Since the principle of solving the problem of the device is similar to that of the acid-etched fracture long-term conductivity determination method, the implementation of the device can be referred to the implementation of the acid-etched fracture conductivity determination method, and the repeated parts will not be described herein again.

[0091] One embodiment of the present application provides an acid-etched fracture long-term conductivity determination device, which specifically can include as shown in Figure 5

[0092] The data acquisition module 201 is configured to acquire the rock mineral composition, the porosity, the permeability, the formation Young's modulus, the formation closure stress, the formation Poisson's ratio, the acid concentration used for reconstruction, and the temperature of the work area.

[0093] The acid-etched damage determination module 202 is configured to determine the acid-etched fracture damage factor of the work area based on the rock mineral composition, the porosity, the permeability, the mechanical property, the acid concentration used for reconstruction, and the temperature of the work area. ​

[0094] Module 203 for determining fracture width after acid etching: This module is used to determine the fracture width after acid etching in the work area based on Young's modulus of the strata, closure stress of the strata, Poisson's ratio of the strata, and acid-rock reaction rate of the strata.

[0095] Creep-after crack width determination module 204: Used to determine the creep-after crack width of acid-etched cracks in the work area based on the acid-etched crack damage factor and the creep-after crack width.

[0096] Flow capacity determination module 205: Used to determine the long-term flow capacity of acid-etched cracks in the work area based on the crack width after acid etching and the crack width after creep.

[0097] In one embodiment, the acid etching damage determination module 202 includes, as follows: Figure 6 As shown:

[0098] Acid-rock reaction rate determination submodule 2021: Based on the acid concentration and temperature used in the modification of the work area, the acid-rock reaction rate of the work area is determined using the acid-rock reaction kinetic model;

[0099] Porosity Change Determination Submodule 2022: Used to determine the porosity change in the work area based on the acid-rock reaction rate and the relationship model between the acid-rock reaction rate and porosity change.

[0100] Young's modulus variation determination submodule 2023: Used to determine the Young's modulus variation in the work area based on the acid-rock reaction rate and the relationship model between the acid-rock reaction rate and the Young's modulus variation.

[0101] The acid etching damage determination submodule 2024 is used to determine the acid etching crack damage factor in the work area based on the acid rock reaction rate, porosity change, Young's modulus change, and acid etching damage relationship model.

[0102] In one embodiment, the etched seam width determination module 203 includes, as follows: Figure 7 As shown:

[0103] Closed stress joint width determination submodule 2031: Based on the Young's modulus, formation closure stress, and formation Poisson's ratio of the work area, determine the width of the closed stress joint in the work area;

[0104] Submodule 2032 for determining the width of acid-rock reaction joints in the work area: Based on the acid-rock reaction rate in the work area, the width of acid-rock reaction joints in the work area is determined using a model that relates the acid-rock reaction rate and the joint width.

[0105] Submodule 2033 for determining the width of the fracture after acid etching: This module is used to determine the width of the fracture after acid etching based on the relationship model between the width of the closed stress fracture in the work area, the width of the acid-rock reaction fracture, and the width of the fracture after acid etching.

[0106] The embodiment of the present application further provides a computer device, which comprises a memory, a processor and a computer program stored in the memory and executable on the processor, and the processor implements the acid-etched fracture conductivity determination method when executing the computer program.

[0107] The embodiment of the present application further provides a computer readable storage medium, which stores a computer program, and the computer program implements the acid-etched fracture conductivity determination method when executed by a processor.

[0108] The embodiment of the present application further provides a computer program product, which comprises a computer program, and the computer program implements the acid-etched fracture conductivity determination method when executed by a processor.

[0109] Embodiment 1:

[0110] The embodiment provides an acid-etched fracture long-term conductivity determination method, which comprises the following steps:

[0111] 1. Obtain the rock mineral composition, porosity, permeability, formation Young's modulus, formation closure stress, formation Poisson's ratio, acid concentration used for reconstruction and temperature of a work area.

[0112] In the method, the rock mineral composition is determined by X-ray diffraction experiment, and the porosity and permeability are determined by porosity and permeability determination experiment.

[0113] The results are shown in Table 1.

[0114] Table 1

[0115] Physical parameters Unit symbol Value Formation closure stress (MPa) P 5-50 Permeability (um 3 )]]> K 1 Porosity (%) Φ 0.7 Formation Young's modulus (GPa) E 50 Formation Poisson's ratio ν 0.25 Acid concentration C 1.5-6 Acid temperature T 25-100 Dolomite content ％ 90-100

[0116] 2. Determine the acid-rock reaction kinetics model, determine the relationship model between acid-rock reaction rate and porosity change, determine the relationship model between acid-rock reaction rate and Young's modulus change, and determine the relationship model among acid-rock reaction rate, porosity change, Young's modulus change and acid-etched damage; specifically comprising:

[0117] Twelve cores numbered 1-12 are obtained from the work area, and the porosity difference before and after acid etching, permeability, Young's modulus ratio, Poisson's ratio and acid-rock reaction rate of the 12 cores are determined by indoor experiment, and the results are shown in Table 2.

[0118] Table 2

[0119]

[0120]

[0121] 3、Based on the concentration and temperature of acid used in the reconstruction of the work area, the acid-rock reaction kinetics model is used to determine the acid-rock reaction rate and acid-etched fracture damage factor of the work area, as shown in Table 3.

[0122] Table 3

[0123]

[0124]

[0125] 4、Through laboratory experiments, the number of micro-convex bodies per unit area of acid-etched fractures, initial fracture width, micro-convex body radius, acid-etched fracture damage factor, closure stress, average value of micro-convex body height distribution, standard deviation of micro-convex body height distribution, and friction coefficient are determined, as shown in Table 4.

[0126] Table 4

[0127] Physical parameters Unit symbol Value Acid-etched fracture asperity number per unit area N 1000000 Initial fracture width (m) w 0.1 Aperity radius (m) R 0.000001 Acid-etched fracture damage factor D 0.3 Closure stress (MPa) ​ 0-100 Aperity height distribution average μ 0.3 Aperity height distribution standard deviation σ 0.3 Friction coefficient f 0.6

[0128] 5、Based on the formation Young's modulus, formation closure stress, and formation Poisson's ratio, the fracture width after closure in the work area is determined by formula, as shown in Table 5 and Figure 1 .

[0129] Table 5

[0130]

[0131] 5、Based on the formation Young's modulus, formation closure stress, and formation Poisson's ratio, the fracture width after acid etching in the work area is determined by formula, as shown in Table 6 and Figure 2 .

[0132] Table 6

[0133] t / s 10 50 100 150 200 250 300 W / mm 0.04 0.21 0.42 0.62 0.83 1.04 1.25

[0134] 6、Based on the formation Young's modulus, formation closure stress, and formation Poisson's ratio, the fracture width after acid etching in the work area is determined by formula, as shown in Table 7 and Figure 3 .

[0135] Table 7

[0136] t / s 10 1000 2000 3000 4000 5000 10000 75000 W / mm 1.54 3.05 2.70 2.44 2.44 2.44 2.44 5.28

[0137] 7、Based on the fracture width after acid etching and the fracture width after creep, the long-term conductivity of acid-etched fractures in the work area is determined by formula, as shown in Table 8 and Figure 4 .

[0138] Table 8

[0139]

[0140] The above merely describes the preferred embodiments of the present application, and it should be noted that, for those skilled in the art, several modifications and improvements can be made without departing from the inventive concept, and these all belong to the protection scope of the present application.

Claims

1. A method for determining long-term conductivity of acid-etched fractures, the method comprising: obtaining, at step 101, rock mineral composition, porosity, permeability, formation Young's modulus, formation closure stress, formation Poisson's ratio, acid concentration used in the treatment, and temperature of the treatment area; determining, at step 102, acid-etched fracture damage factor of the treatment area based on the rock mineral composition, porosity, permeability, formation Young's modulus, acid concentration, and temperature; determining, at step 103, acid-etched fracture face acid-etched fracture width change of the treatment area based on the formation Young's modulus, formation closure stress, formation Poisson's ratio, and acid-rock reaction rate; determining, at step 104, acid-etched fracture width creep of the treatment area based on the acid-etched fracture width, formation Young's modulus, Poisson's ratio, and acid-etched fracture damage factor; determining, at step 105, acid-etched fracture conductivity of the treatment area based on the acid-etched fracture width and acid-etched fracture width after creep; the acid-etched fracture damage factor of the treatment area determined at step 102 comprises: determining, at step 1021, acid-rock reaction rate of the treatment area based on the acid concentration used in the treatment and temperature; determining, at step 1022, porosity change of the treatment area based on the acid-rock reaction rate and a model of relationship between acid-rock reaction rate and porosity change; determining, at step 1023, formation Young's modulus change of the treatment area based on the acid-rock reaction rate and a model of relationship between acid-rock reaction rate and formation Young's modulus change; determining, at step 1024, acid-etched fracture damage factor of the treatment area based on the acid-rock reaction rate, porosity change, and formation Young's modulus change; the model of relationship between acid-rock reaction rate and porosity change at step 1022 is as follows: the model of relationship between acid-rock reaction rate and formation Young's modulus change at step 1023 is as follows: the acid-etched fracture damage model at step 1024 is as follows: wherein, Dc is acid-etched fracture damage factor; the acid-etched fracture face acid-etched fracture width change of the treatment area determined at step 103 comprises: determining, at step 1031, closure stress fracture width of the treatment area based on the formation Young's modulus, formation closure stress, and formation Poisson's ratio; determining, at step 1032, acid-rock reaction fracture width of the treatment area based on the acid-rock reaction rate and a model of relationship between acid-rock reaction rate and fracture width; determining, at step 1033, acid-etched fracture width after acid-etching based on the closure stress fracture width, acid-rock reaction fracture width, and a model of relationship between closure stress fracture width, acid-rock reaction fracture width, and acid-etched fracture width; the model of relationship between closure stress fracture width and the formation Young's modulus, formation closure stress, and formation Poisson's ratio at step 1031 is as follows: the model of relationship between acid-rock reaction fracture width and acid-rock reaction rate at step 1032 is as follows: the model of relationship between acid-etched fracture width and the closure stress fracture width and acid-rock reaction fracture width at step 1033 is as follows: determining, at step 104, acid-etched fracture width after creep of the treatment area based on the acid-etched fracture width, formation Young's modulus, Poisson's ratio, and acid-etched fracture damage factor; the model of relationship between formation creep and acid-etched fracture width at step 104 is as follows: wherein, Dc is acid-etched fracture damage factor. 3.The method of claim 2, wherein: the model of relationship between closure stress fracture width and the formation Young's modulus, formation closure stress, and formation Poisson's ratio at step 1031 is as follows: the model of relationship between acid-rock reaction fracture width and acid-rock reaction rate at step 1032 is as follows: the model of relationship between acid-etched fracture width and the closure stress fracture width and acid-rock reaction fracture width at step 1033 is as follows: determining, at step 104, acid-etched fracture width after creep of the treatment area based on the acid-etched fracture width, formation Young's modulus, Poisson's ratio, and acid-etched fracture damage factor; the model of relationship between formation creep and acid-etched fracture width at step 104 is as follows: wherein, Dc is acid-etched fracture damage factor. ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ; In the formula, wherein, is the change of porosity; T0 is the temperature at normal temperature, K; T is the reaction temperature, K; C is the acid concentration, mol / L; t is the reaction time, h; S is the reaction area, m 2 ; A, b, m are constants, dimensionless; ​ ; wherein E c is the acid-etched Young's modulus, GPa; E0is the initial Young's modulus, GPa; n is a constant, dimensionless; ​ ; ​ 2. The method of claim 1, wherein, ​ ​ ​ ​ ​ ​ ; where N is the number of asperities per unit area of etched crack, dimensionless; w is the initial crack width, m; R is the asperity radius, m; E is the Young's modulus, GPa; v is the Poisson's ratio, dimensionless; and D is the etched crack damage factor, dimensionless. is the closure stress, MPa; and μ is the average of the asperity height distribution, m. is the standard deviation of the asperity height distribution, m. f is the friction coefficient, dimensionless; is the cumulative distribution value of the micro-asperity distribution function; dimensionless; ​ ; where J is the acid-rock reaction rate, g / (cm 2 ·s); t is the reaction time, h; p is the rock density, kg / m 3 ; S is the reaction area; m 2 ; ​ ; where w is the acid-etched seam width, m; w0is the initial seam width, m; is the closed stress seam width, m; is the acid-rock reaction seam width, m.

4. The method of claim 3, wherein, ​ ​ ; where w0 is the initial slit width, m; is the closure stress, MPa; is the confining pressure, MPa; K is the bulk modulus, GPa; G is the shear modulus, GPa; B is the bulk parameter, dimensionless; β is the material homogeneity parameter, dimensionless; η is the viscosity parameter, dimensionless; a is the load level strength parameter, dimensionless; is the loading time, H; D is the acid-etch crack damage factor, dimensionless; t B To accelerate the creep time, H; t f To collapse time, H.

5. The method of claim 4, wherein, In step 105, the relationship model between the acid-etched fracture width and the flow conductivity is determined by the following formula: 。 6. An acid-etched fracture flow conductivity determination device according to the method of any one of claims 1-5, comprising: a data acquisition module for acquiring the rock mineral composition, porosity, permeability, formation Young's modulus, formation closure stress, formation Poisson's ratio, acid concentration used for reconstruction, temperature of the work area; an acid-etched damage determination module for determining the acid-etched fracture damage factor of the work area based on the rock mineral composition, porosity, permeability, mechanical properties, acid concentration used for reconstruction, temperature of the work area; an acid-etched fracture width determination module for determining the acid-etched fracture width of the work area based on the formation Young's modulus, formation closure stress, formation Poisson's ratio, and acid-rock reaction rate of the formation; a post-creep fracture width determination module for determining the post-creep fracture width of the acid-etched fracture of the work area based on the acid-etched fracture damage factor and the acid-etched fracture width; a flow conductivity determination module for determining the long-term flow conductivity of the acid-etched fracture of the work area based on the acid-etched fracture width and the post-creep fracture width.

7. The device of claim 6, wherein: the acid-etched damage determination module comprises: an acid-rock reaction rate determination submodule for determining the acid-rock reaction rate of the work area based on the acid concentration used for reconstruction, temperature of the work area, and using an acid-rock reaction kinetics model; a porosity change determination submodule for determining the porosity change of the work area based on the acid-rock reaction rate of the work area, and using an acid-rock reaction rate and porosity change relationship model; a Young's modulus change determination submodule for determining the Young's modulus change of the work area based on the acid-rock reaction rate of the work area, and using an acid-rock reaction rate and Young's modulus change relationship model; an acid-etched damage determination submodule for determining the acid-etched fracture damage factor of the work area based on the acid-rock reaction rate, porosity change, and Young's modulus change of the work area, and using an acid-rock reaction rate, porosity change, and Young's modulus change relationship and acid-etched damage relationship model; the acid-etched fracture width determination module comprises: a closure stress fracture width determination submodule for determining the closure stress fracture width of the work area based on the formation Young's modulus, formation closure stress, and formation Poisson's ratio of the work area; an acid-rock reaction fracture width determination submodule for determining the acid-rock reaction fracture width of the work area based on the acid-rock reaction rate of the work area, and using an acid-rock reaction rate and fracture width relationship model; an acid-etched fracture width determination submodule for determining the acid-etched fracture width of the work area based on the closure stress fracture width, acid-rock reaction fracture width, and acid-etched fracture width relationship model.

8. A computer-readable storage medium storing a computer program, wherein the computer program is executed by a processor to implement the acid-etched fracture flow conductivity determination method of any one of claims 1-5.

Images