Acid rock reaction rate determination method and device

Abstract

The invention provides an acid rock reaction rate determination method and device. The method comprises the following steps: determining the acid etching crack width of a work area based on rock mineral composition, porosity and permeability of the work area; correcting the surface area of the rock initially participating in the reaction by using the average pore diameter of the work area, the number of pores in unit area and the wall surface roughness coefficient of the crack formed by the reaction of the acid liquor and the rock in the work area to obtain the surface area of the rock participating in the reaction in the work area; determining the reaction volume of the acid liquor in the work area based on the filtration coefficient of the acid liquor in the reaction process of the acid liquor and the rock in the work area, the reaction time of the acid liquor and the rock in the work area and the volume of the acid liquor injected into the work area; determining the H < + > diffusion coefficient of the work area based on the reaction temperature of the acid liquor and the rock in the work area and the speed and concentration of the acid liquor injected into the work area; and determining the reaction rate of the acid rock in the work area based on the width of the acid etching crack in the work area, the surface area of the rock participating in the reaction, the reaction volume of the acid liquor, the H < + > diffusion coefficient, the injection speed and viscosity of the acid liquor and the concentration gradient of the diffusion boundary layer.

Description

Technical Field

[0001] The invention belongs to the technical field of petroleum extraction, and in particular relates to a method and a device for determining an acid-rock reaction rate. Background Art

[0002] Conventional acid-rock reaction rate determination technology does not consider the impact of lithology and physical properties on acid-rock reaction rate; when the mineral composition in carbonate rock lithology is very different, the acid-rock reaction rate will be very different. For example, the reaction rate is fast when the calcite content is high, and the reaction rate is slow when the content of mud, organic matter, quartz and other minerals is high. Conventional acid-rock reaction rate determination technology does not consider the change of acid-rock reaction area. For example, the reaction rate of cracks and caves in carbonate rock physical properties is faster than that of pores. During the acid-rock reaction, as the pores and cracks are dissolved, the specific surface area increases, resulting in an accelerated reaction rate.

[0003] When the lithology and physical properties of the reservoir vary greatly, the results of the acid-rock reaction rate in the core sampling test will vary greatly; if a unified conventional acid-rock reaction kinetic equation is established, it will lead to inaccurate predictions of the acid-rock reaction rate and the effective action distance of the acid etching, which will in turn lead to deviations in the design ideas of the acidizing process.

[0004] In summary, there is still a need to study technical solutions that can accurately determine the acid-rock reaction rate. Summary of the invention

[0005] The purpose of the present invention is to provide a technical solution that can accurately determine the acid-rock reaction rate. In order to achieve the above purpose, the present invention provides the following five technical solutions.

[0006] In a first aspect, the present invention provides a method for determining an acid-rock reaction rate, the method comprising:

[0007] Determine the width of acid-etched cracks in the work area based on the rock mineral composition, porosity, and permeability of the work area;

[0008] The rock surface area initially involved in the reaction is corrected using the average pore diameter and the number of pores per unit area in the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid solution and the rock in the work area to obtain the corrected rock surface area involved in the reaction;

[0009] The acid reaction volume of the work area is determined based on the acid loss coefficient during the reaction between the acid and the rock in the work area, the time for the acid to react with the rock in the work area, and the volume of the acid injected into the work area;

[0010] Based on the temperature of the acid and the rock in the work area, the speed and concentration of the acid injected into the work area, the H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid determines the H of the work area. + Diffusion coefficient;

[0011] Based on the width of the acid-etched cracks in the work area, the surface area of ​​the rock involved in the reaction in the work area, the acid reaction volume in the work area, and the H + The diffusion coefficient, the velocity and viscosity of the acid injected into the work area, the concentration gradient of the diffusion boundary layer in the work area, and the acid-rock reaction rate and fracture width, the rock surface area involved in the reaction, the acid reaction volume, H + The relationship model between the diffusion coefficient, the velocity and viscosity of the acid fluid, and the concentration gradient of the diffusion boundary layer determines the acid-rock reaction rate in the work area.

[0012] In a second aspect, the present invention provides a device for determining an acid-rock reaction rate, the device comprising:

[0013] Fracture width determination module: used to determine the width of acid-etched cracks in the work area based on the rock mineral composition, porosity, and permeability of the work area;

[0014] The module for determining the rock surface area participating in the reaction is used to correct the initial rock surface area participating in the reaction by using the average pore diameter and the number of pores per unit area of ​​the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid solution and the rock in the work area. The corrected rock surface area participating in the reaction is the rock surface area participating in the reaction in the work area.

[0015] Acid reaction volume determination module: used to determine the acid reaction volume of the work area based on the acid loss coefficient during the reaction between the acid and the rock in the work area, the time for the acid to react with the rock in the work area, and the volume of the acid injected into the work area;

[0016] H + Diffusion coefficient determination module: It is used to determine the diffusion coefficient based on the temperature of the acid and the rock in the work area, the speed and concentration of the acid injected into the work area, and the H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid determines the H of the work area. + Diffusion coefficient;

[0017] Acid rock reaction rate determination module: used to determine the acid-rock reaction rate based on the width of the acid-etched cracks in the work area, the surface area of ​​the rock involved in the reaction in the work area, the acid reaction volume in the work area, and the H + The diffusion coefficient, the velocity and viscosity of the acid injected into the work area, the concentration gradient of the diffusion boundary layer in the work area, and the acid-rock reaction rate and fracture width, the rock surface area involved in the reaction, the acid reaction volume, H + The relationship model between the diffusion coefficient, the velocity and viscosity of the acid fluid, and the concentration gradient of the diffusion boundary layer determines the acid-rock reaction rate in the work area.

[0018] In a third aspect, the present invention provides a computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the acid-rock reaction rate determination method provided in the first aspect when executing the computer program.

[0019] In a fourth aspect, the present invention provides a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the acid-rock reaction rate determination method provided in the first aspect.

[0020] In a fifth aspect, the present invention provides a computer program product, which includes a computer program, and when the computer program is executed by a processor, it implements the acid-rock reaction rate determination method provided in the first aspect.

[0021] The technical solution provided by the present invention takes into account the influence of multiple factors such as crack wall roughness, acid liquid properties, rock mineral composition, reaction temperature, etc. on the acid-rock reaction rate, and can accurately determine the acid-rock reaction rate, providing a basis for accurately analyzing the effective action distance of acid etching and optimizing the liquid parameters of the acid fracturing process. BRIEF DESCRIPTION OF THE DRAWINGS

[0022] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the prior art descriptions. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without creative work. In the drawings:

[0023] Figure 1 This is the acid-rock reaction rate-temperature relationship diagram.

[0024] Figure 2 This is the relationship diagram between acid-rock reaction rate and acid concentration.

[0025] Figure 3 This is the relationship diagram between acid rock reaction rate and mineral content.

[0026] Figure 4 This is the relationship diagram between acid-rock reaction rate and acid viscosity.

[0027] Figure 5 This is the acid-rock reaction rate-pressure difference relationship diagram.

[0028] Figure 6 This is the relationship diagram between acid-rock reaction rate and acid flow rate. DETAILED DESCRIPTION

[0029] To make the purpose, technical solution and advantages of the embodiments of the present invention more clear, the embodiments of the present invention are further described in detail below in conjunction with the accompanying drawings. Here, the exemplary embodiments of the present invention and their descriptions are used to explain the present invention, but are not intended to limit the present invention.

[0030] In the description of this specification, the terms "include", "including", "have", "contain", etc. are all open terms, which mean including but not limited to. The descriptions with reference to the terms "one embodiment", "a specific embodiment", "some embodiments", "for example", etc. mean that the specific features, structures or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present application. In this specification, the schematic representations of the above terms do not necessarily refer to the same embodiment or example. Moreover, the specific features, structures or characteristics described may be combined in a suitable manner in any one or more embodiments or examples. The order of steps involved in each embodiment is used to schematically illustrate the implementation of the present application, and the order of steps is not limited and can be appropriately adjusted as needed.

[0031] In order to construct a shale composite model that is well representative of actual shale, the embodiments of the present invention propose a shale composite pore and fracture medium molecular model construction method, a shale composite pore and fracture medium molecular model construction device, a computer device, a computer-readable storage medium, and a computer program product.

[0032] An embodiment of the present invention provides a method for determining an acid-rock reaction rate, the method comprising the following steps:

[0033] Step 101, determining the width of the acid-etched cracks in the work area based on the rock mineral composition, porosity, and permeability of the work area;

[0034] Step 102, using the average pore diameter and the number of pores per unit area of ​​the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid solution and the rock in the work area, the initial rock surface area participating in the reaction (i.e., the projected area of ​​the rock wall participating in the reaction) is corrected to obtain the corrected rock surface area participating in the reaction, which is the rock surface area participating in the reaction in the work area;

[0035] Step 103, determining the reaction volume of the acid solution in the work area based on the loss coefficient of the acid solution during the reaction between the acid solution and the rock in the work area, the reaction time between the acid solution and the rock in the work area, and the volume of the acid solution injected into the work area;

[0036] Step 104, based on the temperature of the reaction between the acid and the rock in the work area, the speed and concentration of the acid injected into the work area, the relationship model between the H+ diffusion coefficient and the temperature of the reaction between the acid and the rock, the speed and concentration of the acid is used to determine the H+ in the work area. + Diffusion coefficient;

[0037] Step 105, based on the width of the acid-etched cracks in the work area, the surface area of ​​the rock involved in the reaction in the work area, the acid reaction volume in the work area, and the H +The acid-rock reaction rate in the work area is determined by the diffusion coefficient, the velocity and viscosity of the acid injected into the work area, and the concentration gradient of the diffusion boundary layer in the work area, and the relationship model between the acid-rock reaction rate and the fracture width, the corrected rock surface area participating in the reaction, the acid reaction volume, the H+ diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer.

[0038] The existing methods for determining the acid-rock reaction rate do not take into account the influence of lithology and physical properties on the acid-rock reaction rate, and fail to reflect the problem that as the holes and cracks are dissolved during the acid-rock reaction, the specific surface area increases, resulting in an accelerated reaction rate. The above-mentioned technical scheme provided by the present invention takes into account the influence of multiple factors such as the roughness of the crack wall, the physical properties of the acid liquid, the mineral composition of the rock, and the reaction temperature on the acid-rock reaction rate. It can accurately determine the acid-rock reaction rate, and provide a basis for accurately analyzing the effective action distance of the acid etching and optimizing the liquid parameters of the acid fracturing process.

[0039] In one embodiment, in step 101, determining the width of the acid-etched cracks in the work area based on the rock mineral composition, porosity, and permeability of the work area includes:

[0040] Step 1011, based on the rock mineral composition, porosity, and permeability of the work area, the fractal dimension of the fracture wall of the work area is determined by using a relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity, and permeability;

[0041] Step 1012, based on the fractal dimension of the crack wall in the work area, the average curvature of the acid-etched crack space in the work area is determined using a relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall;

[0042] Step 1013, based on the average curvature of the acid-etched crack space in the work area, the width of the acid-etched crack in the work area is determined using a relationship model between the width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area.

[0043] In one embodiment, in step 1011, the relationship model between the fracture wall fractal dimension and the rock mineral composition, porosity, and permeability is:

[0044]

[0045] Where, X i is the content of the ith mineral, dimensionless; is porosity, dimensionless; K is permeability, mD; a i , b1, b2 are coefficients; D is the fractal dimension of the crack wall, dimensionless;

[0046] Furthermore, the rock mineral composition includes N-1 carbonate minerals, X1 to X N-1 are the contents of these N-1 carbonate minerals, X Nis the non-carbonate mineral content;

[0047] For example, the relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity, and permeability is:

[0048]

[0049] Where, X1 is the calcite content, dimensionless; X2 is the dolomite content, dimensionless; X3 is the non-carbonate mineral content, dimensionless; is the porosity, dimensionless; K is the permeability, mD; D is the fractal dimension of the fracture wall, dimensionless.

[0050] In one embodiment, in step 1012, the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall is:

[0051] K a =mD+n

[0052] In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; m and n are coefficients; D is the fractal dimension of the crack wall, dimensionless;

[0053] For example, the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall is:

[0054] K a =0.4524D-0.3937

[0055] In the formula, K a is the spatial average curvature of the acid-etched crack, dimensionless; D is the fractal dimension of the crack wall, dimensionless.

[0056] In one embodiment, in step 1013, the relationship model between the acid-etched crack width and the average curvature of the acid-etched crack space in the work area is:

[0057] w=pK a +q

[0058] In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; p and q are coefficients; w is the width of the acid-etched crack, dimensionless;

[0059] For example, the relationship model between the width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area is:

[0060] w=8.0667K a -0.5974

[0061] In the formula, K ais the spatial average curvature of the acid-etched crack, dimensionless; D is the fractal dimension of the crack wall, dimensionless; w is the width of the acid-etched crack, unit: mm.

[0062] In one embodiment, in step 102, the surface area of ​​the rock initially involved in the reaction is corrected using the average pore diameter and the number of pores per unit area of ​​the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid solution with the rock in the work area, using the following formula:

[0063]

[0064] Where s′ is the corrected rock surface area participating in the reaction, cm 2 ; s is the initial rock surface area involved in the reaction, cm 2 ;n i is the number of pores per unit area in the work area, dimensionless; U is the wall roughness coefficient of the cracks formed by the reaction of acid and rock in the work area, dimensionless; d is the average pore diameter in the work area, micrometers.

[0065] In one embodiment, in step 103, the acid solution reaction volume in the working area is determined by the following formula:

[0066]

[0067] Where V' is the acid reaction volume, L; V is the initial acid volume, L; C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, t is the time for the acid to react with the rocks in the work area, s.

[0068] In one embodiment, in step 103, the loss coefficient of the acid solution during the reaction between the acid solution and the rock in the work area is determined based on the porosity and permeability of the work area and the pore volume (i.e., the ratio of the pore volume to the volume of mineral particles), the viscosity of the acid solution injected into the work area, and the pressure difference between the fracture in the work area and the formation (i.e., the difference between the bottom hole pressure and the formation pressure during construction);

[0069] Furthermore, the filtration coefficient of the acid during the reaction between the acid and the rock in the work area is calculated by the following formula:

[0070]

[0071] Where C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, is the porosity, dimensionless; K is the permeability, mD; μ is the acid viscosity, mPa·s; Q i is the pore volume of the work area, dimensionless; Δp is the pressure difference between the fracture in the work area and the formation, MPa.

[0072] In one embodiment, in step 104, H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid is:

[0073]

[0074] In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid injected into the work area, mol / L; u is the velocity of the acid, m / s; m1, m2, n1, n2 are coefficients;

[0075] Furthermore, H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid is:

[0076]

[0077] In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; α, β, γ are correction coefficients;

[0078] This preferred technical solution is commonly used in H + The diffusion coefficient is determined based on the model for H + The diffusion coefficient was corrected, and the H obtained by using this preferred technical solution + The acid-rock reaction rate is determined by the diffusion coefficient, and the obtained acid-rock reaction rate is more accurate;

[0079] Furthermore, α = 0.4, γ=10.

[0080] In one embodiment, in step 105, the acid-rock reaction rate is related to the fracture width, the surface area of ​​the rock involved in the reaction, the acid reaction volume, H + The relationship model between the diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer is:

[0081]

[0082] In the formula, H + Diffusion coefficient, m 2 / s; w is the crack width, cm; u is the velocity of the acid, m / s; μ is the viscosity of the acid, mPa·s; s′ is the surface area of ​​the rock involved in the reaction, cm 2 ; V′ is the volume of acid solution reaction, L; is the concentration gradient of the diffusion boundary layer, mol / (L·cm); K is the acid-rock reaction rate, mol / (cm 2 ·s).

[0083] In one embodiment, the method further comprises a data acquisition step, wherein the data acquisition step comprises:

[0084] Obtain the mineral composition, porosity, permeability, average pore diameter, and number of pores per unit area of ​​the rock in the work area; obtain the volume, concentration, viscosity, and velocity of the acid injected into the work area; obtain the time and temperature of the reaction between the acid and the rock in the work area, the wall roughness coefficient of the formed cracks, and the initial rock surface area participating in the reaction;

[0085] Furthermore, the data acquisition step includes: obtaining the rock mineral composition, porosity, permeability, average pore diameter, pore volume, and number of pores per unit area of ​​the work area; obtaining the volume, concentration, viscosity, and velocity of the acid injected into the work area; obtaining the time and temperature of the reaction between the acid and the rock in the work area, the wall roughness coefficient of the formed cracks, and the rock surface area initially involved in the reaction; and obtaining the pressure difference between the cracks in the work area and the formation.

[0086] The present invention also provides an acid-rock reaction rate determination device, as described in the following embodiments. Since the principle of the device to solve the problem is similar to the acid-rock reaction rate determination method, the implementation of the device can refer to the implementation of the acid-rock reaction rate determination method, and the repeated parts will not be repeated.

[0087] An embodiment of the present invention provides a device for determining an acid-rock reaction rate, which may specifically include:

[0088] The crack width determination module 201 is used to determine the crack width of the acid-etched cracks in the work area based on the rock mineral composition, porosity and permeability of the work area;

[0089] The rock surface area determination module 202 is used to correct the initial rock surface area participating in the reaction by using the average pore diameter and the number of pores per unit area of ​​the work area and the wall roughness coefficient of the cracks formed by the reaction of the acid solution and the rock in the work area to obtain the corrected rock surface area participating in the reaction, which is the rock surface area participating in the reaction in the work area;

[0090] Acid liquid reaction volume determination module 203: used to determine the acid liquid reaction volume of the work area based on the acid liquid loss coefficient during the reaction process between the acid liquid and the work area rock, the time for the acid liquid to react with the work area rock, and the volume of the acid liquid injected into the work area;

[0091] H + Diffusion coefficient determination module: It is used to determine the diffusion coefficient based on the temperature of the acid and the rock in the work area, the speed and concentration of the acid injected into the work area, and the H+ The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid determines the H of the work area. + Diffusion coefficient;

[0092] Acid rock reaction rate determination module 204: used to determine the acid-rock reaction rate based on the width of the acid-etched cracks in the work area, the surface area of ​​the rock involved in the reaction in the work area, the acid reaction volume in the work area, and the H + The acid-rock reaction rate in the work area is determined by the diffusion coefficient, the velocity and viscosity of the acid injected into the work area, and the concentration gradient of the diffusion boundary layer in the work area, and the relationship model between the acid-rock reaction rate and the fracture width, the corrected rock surface area participating in the reaction, the acid reaction volume, the H+ diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer.

[0093] In one embodiment, the seam width determination module 201 includes:

[0094] Fractal dimension determination submodule 2011: used to determine the fractal dimension of the fracture wall in the work area based on the rock mineral composition, porosity and permeability of the work area and using the relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity and permeability;

[0095] The average curvature determination submodule 2012 is used to determine the average curvature of the acid-etched crack space in the work area based on the fractal dimension of the crack wall in the work area and using the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall;

[0096] The crack width determination submodule 2013 is used to determine the crack width of the acid-etched cracks in the work area based on the average curvature of the acid-etched crack space in the work area and using the relationship model between the crack width of the acid-etched cracks and the average curvature of the acid-etched crack space in the work area.

[0097] In one embodiment, the relationship model between the fracture wall fractal dimension and the rock mineral composition, porosity, and permeability is:

[0098]

[0099] In the formula, X i is the content of the ith mineral, dimensionless; is porosity, dimensionless; K is permeability, mD; a i , b1, b2 are coefficients; D is the fractal dimension of the crack wall, dimensionless;

[0100] Furthermore, the rock mineral composition includes N-1 carbonate minerals, X1 to X N-1 are the contents of these N-1 carbonate minerals, X N is the non-carbonate mineral content;

[0101] For example, the relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity, and permeability is:

[0102]

[0103] Where, X1 is the calcite content, dimensionless; X2 is the dolomite content, dimensionless; X3 is the non-carbonate mineral content, dimensionless; is the porosity, dimensionless; K is the permeability, mD; D is the fractal dimension of the fracture wall, dimensionless.

[0104] In one embodiment, the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall is:

[0105] K a =mD+n

[0106] In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; m and n are coefficients; D is the fractal dimension of the crack wall, dimensionless;

[0107] For example, the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall is:

[0108] K a =0.4524D-0.3937

[0109] In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; D is the fractal dimension of the crack wall, dimensionless. In one embodiment, the relationship model between the width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area is:

[0110] w=pK a +q

[0111] In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; p and q are coefficients; w is the width of the acid-etched crack, dimensionless;

[0112] For example, the relationship model between the width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area is:

[0113] w=8.0667K a -0.5974

[0114] In the formula, K a is the spatial average curvature of the acid-etched crack, dimensionless; D is the fractal dimension of the crack wall, dimensionless; w is the width of the acid-etched crack, unit: mm.

[0115] In one embodiment, the rock surface area determination module 202 performs initial rock surface area correction for the reaction by using the following formula:

[0116]

[0117] Where s′ is the corrected rock surface area participating in the reaction, cm 2 ; s is the initial rock surface area involved in the reaction, cm 2 ;n i is the number of pores per unit area in the work area, dimensionless; U is the wall roughness coefficient of the cracks formed by the reaction of acid and rock in the work area, dimensionless; d is the average pore diameter in the work area, micrometers.

[0118] In one embodiment, the acid liquid reaction volume determination module 203 determines the acid liquid reaction volume of the work area by the following formula:

[0119]

[0120] Where V' is the acid reaction volume, L; V is the initial acid volume, L; C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, t is the time for the acid to react with the rocks in the work area, s.

[0121] In one embodiment, the acid reaction volume determination module 203 is used to determine the loss coefficient of the acid during the reaction between the acid and the rock in the work area based on the porosity and permeability and pore volume of the work area, the viscosity of the acid injected into the work area, and the pressure difference between the fracture in the work area and the formation;

[0122] Furthermore, the acid liquid reaction volume determination module 203 determines the filtration coefficient of the acid liquid during the reaction between the acid liquid and the work area rock by the following formula:

[0123]

[0124] Where C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, is the porosity, dimensionless; K is the permeability, mD; μ is the acid viscosity, mPa·s; Q i is the pore volume of the work area, dimensionless; Δp is the pressure difference between the fracture in the work area and the formation, MPa.

[0125] In one embodiment, H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid is:

[0126]

[0127] In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; m1, m2, n1, n2 are coefficients;

[0128] Furthermore, H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid is:

[0129]

[0130] In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; α, β, γ are correction coefficients;

[0131] This preferred technical solution is commonly used in H + The diffusion coefficient is determined based on the model for H + The diffusion coefficient was corrected, and the H obtained by using this preferred technical solution + The acid-rock reaction rate is determined by the diffusion coefficient, and the obtained acid-rock reaction rate is more accurate;

[0132] Furthermore, α = 0.4, γ=10.

[0133] In one embodiment, the acid-rock reaction rate is related to the fracture width, the rock surface area involved in the reaction, the acid reaction volume, H + The relationship model between the diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer is:

[0134]

[0135] In the formula, H + Diffusion coefficient, m 2 / s; w is the crack width, cm; u is the velocity of the acid, m / s; μ is the viscosity of the acid, mPa·s; s′ is the surface area of ​​the rock involved in the reaction, cm 2 ; V′ is the volume of acid solution reaction, L; is the concentration gradient of the diffusion boundary layer, mol / (L·cm); K is the acid-rock reaction rate, mol / (cm 2 ·s).

[0136] In one embodiment, the device further includes a data acquisition module, which is used to obtain the mineral composition, porosity, permeability, average pore diameter, and number of pores per unit area of ​​the rock in the work area; obtain the volume, concentration, viscosity, and speed of the acid solution injected into the work area; obtain the time and temperature of the reaction between the acid solution and the rock in the work area, the wall roughness coefficient of the formed cracks, and the surface area of ​​the rock initially participating in the reaction;

[0137] Furthermore, the data acquisition module is used to obtain the rock mineral composition, porosity, permeability, average pore diameter, pore volume, and number of pores per unit area of ​​the work area; obtain the volume, concentration, viscosity, and velocity of the acid injected into the work area; obtain the time and temperature of the reaction between the acid and the rock in the work area, the wall roughness coefficient of the formed cracks, and the rock surface area initially involved in the reaction; and obtain the pressure difference between the cracks in the work area and the formation.

[0138] An embodiment of the present invention also provides a computer device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, and the processor implements the above-mentioned acid-rock reaction rate determination method when executing the computer program.

[0139] An embodiment of the present invention further provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned acid-rock reaction rate determination method is implemented.

[0140] An embodiment of the present invention further provides a computer program product, which includes a computer program. When the computer program is executed by a processor, the above-mentioned acid-rock reaction rate determination method is implemented.

[0141] Embodiment 1:

[0142] This embodiment provides a method for determining an acid-rock reaction rate, the method comprising:

[0143] 1. Obtain the mineral composition, porosity, permeability, average pore diameter, pore volume, and number of pores per unit area of ​​the rock in the work area; obtain the volume, concentration, viscosity, and velocity of the acid injected into the work area; obtain the time and temperature of the reaction between the acid and the rock in the work area, the wall roughness coefficient of the formed cracks, and the initial rock surface area involved in the reaction; obtain the pressure difference between the cracks in the work area and the formation;

[0144] The results are shown in Table 1:

[0145] Table 1

[0146]

[0147]

[0148] 2. Determine the width of the acid-etched cracks in the work area based on the rock mineral composition, porosity and permeability of the work area; specifically include:

[0149] Based on the rock mineral composition, porosity and permeability of the work area, the fractal dimension of the fracture wall in the work area is determined by using the relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity and permeability; among them, the relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity and permeability is:

[0150]

[0151] Where, X1 is the calcite content, dimensionless; X2 is the dolomite content, dimensionless; X3 is the non-carbonate mineral content, dimensionless; is the porosity, dimensionless; K is the permeability, mD; D is the fractal dimension of the fracture wall, dimensionless;

[0152] Based on the fractal dimension of the crack wall in the work area, the average curvature of the acid-etched crack space in the work area is determined by using the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall; among them, the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall is:

[0153] K a =0.4524D-0.3937

[0154] In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; D is the fractal dimension of the crack wall, dimensionless;

[0155] Based on the average curvature of the acid-etched crack space in the work area, the relationship model between the acid-etched crack width and the average curvature of the acid-etched crack space in the work area is used to determine the acid-etched crack width in the work area; among them, the relationship model between the acid-etched crack width and the average curvature of the acid-etched crack space in the work area is:

[0156] w=8.0667K a -0.5974

[0157] In the formula, K a is the spatial average curvature of the acid-etched crack, dimensionless; D is the fractal dimension of the crack wall, dimensionless; w is the width of the acid-etched crack, unit: mm.

[0158] 3. The rock surface area initially participating in the reaction is corrected by using the average pore diameter and the number of pores per unit area in the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid solution and the rock in the work area to obtain the corrected rock surface area participating in the reaction; wherein the rock surface area initially participating in the reaction is corrected by the following formula:

[0159]

[0160] Where s′ is the corrected rock surface area participating in the reaction, cm 2 ; s is the initial surface area of ​​the rock involved in the reaction, cm 2 ;n i is the number of pores per unit area in the work area, dimensionless; U is the wall roughness coefficient of the cracks formed by the reaction of acid and rock in the work area, dimensionless; d is the average pore diameter in the work area, micrometers.

[0161] 4. Based on the porosity, permeability and pore volume of the work area, the viscosity of the acid injected into the work area, and the gap between the work area fractures and the formation, the following formula is used to determine the acid loss coefficient during the reaction between the acid and the work area rock:

[0162]

[0163] Where C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, is the porosity, dimensionless; K is the permeability, mD; μ is the acid viscosity, mPa·s; Q i is the pore volume of the work area, dimensionless; Δp is the pressure difference between the fracture in the work area and the formation, MPa;

[0164] Based on the acid loss coefficient during the reaction between the acid and the rock in the work area, the time for the acid to react with the rock in the work area, and the volume of the acid injected into the work area, the following formula is used to determine the acid reaction volume in the work area:

[0165]

[0166] Where V' is the acid reaction volume, L; V is the initial acid volume, L; C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, t is the time for the acid to react with the rocks in the work area, s.

[0167] 5. Based on the temperature of the acid and the rock in the work area, the speed and concentration of the acid injected into the work area, use H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid determines the H of the work area. + Diffusion coefficient; where H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid is:

[0168]

[0169] In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; α, β, γ are correction coefficients, α = 0.4, γ=10.

[0170] 6. Based on the width of the acid-etched cracks in the work area, the surface area of ​​the rock involved in the reaction in the work area, the acid reaction volume in the work area, and the H + The diffusion coefficient, the velocity and viscosity of the acid injected into the work area, the concentration gradient of the diffusion boundary layer in the work area, and the acid-rock reaction rate and fracture width, the rock surface area involved in the reaction, the acid reaction volume, H + The acid-rock reaction rate in the work area is determined by the relationship model between the diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer; the acid-rock reaction rate is related to the fracture width, the surface area of ​​the rock involved in the reaction, the acid reaction volume, H + The relationship model between the diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer is:

[0171]

[0172] In the formula, H + Diffusion coefficient, m 2 / s; w is the crack width, cm; u is the velocity of the acid, m / s; μ is the viscosity of the acid, mPa·s; s′ is the surface area of ​​the rock involved in the reaction, cm 2 ; V′ is the volume of acid solution reaction, L; is the concentration gradient of the diffusion boundary layer, mol / (L·cm); K is the acid-rock reaction rate, mol / (cm 2 ·s).

[0173] The results of the acid-rock reaction rate at different acid and rock reaction temperatures are as follows: Figure 1 shown.

[0174] In order to verify the accuracy of the results, the inventors conducted indoor experiments to determine the acid-rock reaction rate at different acid and rock reaction temperatures in the work area. The results are as follows: Figure 1 shown.

[0175] like Figure 1 As shown, the acid-rock reaction rate determined by using the technical solution provided by the present invention is very close to the actual experimental data, which verifies the accuracy of the technical solution provided by the present invention.

[0176] Example 2

[0177] This embodiment provides a method for determining an acid-rock reaction rate. The difference between this method and Embodiment 1 is that the values ​​of the parameters obtained in step 1 are different. The results are shown in Table 2:

[0178] Table 2

[0179]

[0180]

[0181] The final results of the acid-rock reaction rate under different acid concentrations in the injection area are as follows: Figure 2 shown.

[0182] In order to verify the accuracy of the results, the inventors conducted indoor experiments to determine the acid-rock reaction rate under different acid concentrations in the injection area. The results are as follows: Figure 2 shown.

[0183] like Figure 2 As shown, the acid-rock reaction rate determined by using the technical solution provided by the present invention is very close to the actual experimental data, which verifies the accuracy of the technical solution provided by the present invention.

[0184] Example 3

[0185] This embodiment provides a method for determining an acid-rock reaction rate. The difference between this method and Embodiment 1 is that the values ​​of the parameters obtained in step 1 are different. The results are shown in Table 3:

[0186] Table 3

[0187] Physical parameters Unit symbol Numeric Temperature of acid liquid reacting with rock in the work area (℃) T 150 <![CDATA[Permeability (um 3 )]]> K 7.14 Porosity(%) Φ 0.01 Number of pores per unit area <![CDATA[n i ]]> 207.5 The time for the acid to react with the rocks in the work area (s) t 300 Volume of acid injected into the work area (L) V 0.5 Velocity of acid injected into the work area (m / s) μ 0.25 Viscosity of the acid solution injected into the work area (mPa·s) v 30 Pore ​​volume Qi 5.554346 Pressure difference between the fracture and the formation in the work area (MPa) △P 7 Concentration of acid injected into the work area (mol / L) C 5.91 Calcite content (decimal) <![CDATA[X1]]> 0.1-0.8 Dolomite content (decimal) <![CDATA[X2]]> 0.1-0.8 Non-carbonate mineral content (decimal) <![CDATA[X3]]> 0.1

[0188] The final results of the acid-rock reaction rate under different mineral compositions are as follows: Figure 3 shown.

[0189] Example 4

[0190] This embodiment provides a method for determining an acid-rock reaction rate. The difference between this method and Embodiment 1 is that the values ​​of the parameters obtained in step 1 are different. The results are shown in Table 4:

[0191] Table 4

[0192] Physical parameters Unit symbol Numeric Temperature of acid liquid reacting with rock in the work area (℃) T 150 <![CDATA[Permeability (um 3 )]]> K 7.14 Porosity(%) Φ 0.01 Number of pores per unit area <![CDATA[n i ]]> 207.5 The time for the acid to react with the rocks in the work area (s) t 300 Volume of acid injected into the work area (L) V 0.5 Velocity of acid injected into the work area (m / s) μ 0.25 Viscosity of the acid solution injected into the work area (mPa·s) v 10-100 Pore ​​volume Qi 5.554346 Pressure difference between the fracture and the formation in the work area (MPa) △P 7 Concentration of acid injected into the work area (mol / L) C 5.91 Calcite content (decimal) <![CDATA[X1]]> 0.7 Dolomite content (decimal) <![CDATA[X2]]> 0.2 Non-carbonate mineral content (decimal) <![CDATA[X3]]> 0.1

[0193] The final results of the acid-rock reaction rate under different acid viscosity in different injection areas are as follows: Figure 4 shown.

[0194] Example 5

[0195] This embodiment provides a method for determining an acid-rock reaction rate. The difference between this method and Embodiment 1 is that the values ​​of the parameters obtained in step 1 are different. The results are shown in Table 5:

[0196] Table 5

[0197]

[0198]

[0199] The final results of the acid-rock reaction rate under the pressure difference between the fracture and the formation in different working areas are as follows: Figure 5 shown.

[0200] Example 5

[0201] This embodiment provides a method for determining an acid-rock reaction rate. The difference between this method and Embodiment 1 is that the values ​​of the parameters obtained in step 1 are different. The results are shown in Table 6:

[0202] Table 6

[0203] Physical parameters Unit symbol Numeric Temperature of acid liquid reacting with rock in the work area (℃) T 150 <![CDATA[Permeability (um 3 )]]> K 7.14 Porosity(%) Φ 0.01 Number of pores per unit area <![CDATA[n i ]]> 207.5 The time for the acid to react with the rocks in the work area (s) t 300 Volume of acid injected into the work area (L) V 0.5 Velocity of acid injected into the work area (m / s) μ 0.05-0.5 Viscosity of the acid solution injected into the work area (mPa·s) v 30 Pore ​​volume Qi 5.554346 Pressure difference between the fracture and the formation in the work area (MPa) △P 7 Concentration of acid injected into the work area (mol / L) C 5.91 Calcite content (decimal) <![CDATA[X1]]> 0.7 Dolomite content (decimal) <![CDATA[X2]]> 0.2 Non-carbonate mineral content (decimal) <![CDATA[X3]]> 0.1

[0204] The final results of the acid-rock reaction rate under different acid injection rates are as follows: Figure 6 shown.

[0205] Those skilled in the art will appreciate that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0206] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0207] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1A function specified in one or more boxes.

[0208] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0209] The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the scope of protection of the present invention.

Claims

1. A method for determining an acid-rock reaction rate, the method comprising: Determine the width of acid-etched cracks in the work area based on the rock mineral composition, porosity, and permeability of the work area; The rock surface area initially participating in the reaction is corrected using the average pore diameter and the number of pores per unit area in the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid solution and the rock in the work area. The corrected rock surface area participating in the reaction is the rock surface area participating in the reaction in the work area. The acid reaction volume of the work area is determined based on the acid loss coefficient during the reaction between the acid and the rock in the work area, the time for the acid to react with the rock in the work area, and the volume of the acid injected into the work area; Based on the temperature of the acid and the rock in the work area, the speed and concentration of the acid injected into the work area, the H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid determines the H of the work area. + Diffusion coefficient; Based on the width of the acid-etched cracks in the work area, the surface area of ​​the rock involved in the reaction in the work area, the acid reaction volume in the work area, and the H + The acid-rock reaction rate in the work area is determined by the diffusion coefficient, the velocity and viscosity of the acid injected into the work area, and the concentration gradient of the diffusion boundary layer in the work area, and the relationship model between the acid-rock reaction rate and the fracture width, the surface area of ​​the rock involved in the reaction, the acid reaction volume, the H+ diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer.

2. The method according to claim 1, wherein determining the width of the acid-etched cracks in the work area based on the rock mineral composition, porosity, and permeability of the work area comprises: Based on the rock mineral composition, porosity and permeability of the work area, the fractal dimension of the fracture wall in the work area is determined by using the relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity and permeability. Based on the fractal dimension of the crack wall in the work area, the average curvature of the acid-etched crack space in the work area is determined using the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall. Based on the average curvature of the acid-etched crack space in the work area, the width of the acid-etched crack in the work area is determined by using the relationship model between the width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area. Preferably, the relationship model between the fracture wall fractal dimension and the rock mineral composition, porosity and permeability is: Where, X i is the content of the ith mineral, dimensionless; is porosity, dimensionless; K is permeability, mD; a i , b1, b2 are coefficients; D is the fractal dimension of the crack wall, dimensionless; More preferably, the rock mineral composition includes N-1 carbonate minerals, X1 to X N-1 are the contents of these N-1 carbonate minerals, X N is the non-carbonate mineral content; Preferably, the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall is: K a =mD+n In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; m and n are coefficients; D is the fractal dimension of the crack wall, dimensionless; Preferably, the width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area are: w=pK a +q In the formula, K a is the spatial average curvature of the acid-etched crack, dimensionless; p and q are coefficients; w is the width of the acid-etched crack, dimensionless.

3. The method according to claim 1, wherein: The average pore diameter and the number of pores per unit area in the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid and the rock in the work area are used to correct the initial rock surface area involved in the reaction through the following formula: Where s′ is the corrected rock surface area participating in the reaction, cm 2 ; s is the initial rock surface area involved in the reaction, cm 2 ;n i is the number of pores per unit area in the work area, dimensionless; U is the wall roughness coefficient of the cracks formed by the reaction of acid and rock in the work area, dimensionless; d is the average pore diameter in the work area, micrometers.

4. The method according to claim 1, wherein: The acid reaction volume in the work area is determined by the following formula: Where V' is the acid reaction volume, L; V is the initial acid volume, L; C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, t is the time for the acid to react with the rocks in the work area, s.

5. The method according to claim 1, wherein: The acid loss coefficient during the reaction between the acid and the rock in the work area is determined based on the porosity, permeability and pore volume of the work area, the viscosity of the acid injected into the work area, and the pressure difference between the fractures in the work area and the formation; Preferably, the filtration coefficient of the acid solution during the reaction between the acid solution and the rock in the work area is calculated by the following formula: Where C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, is the porosity, dimensionless; K is the permeability, mD; μ is the acid viscosity, mPa·s; Q i is the pore volume of the work area, dimensionless; Δp is the pressure difference between the fracture in the work area and the formation, MPa.

6. The method according to claim 1, wherein: H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid is: In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; m1, m2, n1, n2 are coefficients; Preferably, H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid is: In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; α, β, γ are correction coefficients; More preferably, α=0.4, γ=10.

7. The method according to claim 1, wherein: The acid-rock reaction rate is related to the fracture width, the rock surface area involved in the reaction, the acid reaction volume, H + The relationship model between the diffusion coefficient, the injection velocity and viscosity of the acid solution, and the concentration gradient of the diffusion boundary layer is: In the formula, H + Diffusion coefficient, m 2 / s; w is the slit width, cm; u is the speed of the acid, m / s; μ is the viscosity of the acid solution, mPa·s; s′ is the surface area of ​​the rock involved in the reaction, cm 2 ; V′ is the volume of acid solution reaction, L; is the concentration gradient of the diffusion boundary layer, mol / (L·cm); K is the acid-rock reaction rate, mol / (cm 2 ·s).

8. A device for determining an acid-rock reaction rate, the device specifically comprising: Fracture width determination module: used to determine the width of acid-etched cracks in the work area based on the rock mineral composition, porosity, and permeability of the work area; The module for determining the rock surface area participating in the reaction is used to correct the initial rock surface area participating in the reaction by using the average pore diameter and the number of pores per unit area of ​​the work area, and the wall roughness coefficient of the cracks formed by the reaction of the acid solution and the rock in the work area. The corrected rock surface area participating in the reaction is the rock surface area participating in the reaction in the work area. Acid reaction volume determination module: used to determine the acid reaction volume of the work area based on the acid loss coefficient during the reaction between the acid and the rock in the work area, the time for the acid to react with the rock in the work area, and the volume of the acid injected into the work area; H + Diffusion coefficient determination module: It is used to determine the diffusion coefficient based on the temperature of the acid and the rock in the work area, the speed and concentration of the acid injected into the work area, and the H + The relationship model between the diffusion coefficient and the temperature of the acid-rock reaction, the velocity and concentration of the acid determines the H of the work area. + Diffusion coefficient; Acid rock reaction rate determination module: used to determine the acid-rock reaction rate based on the width of the acid-etched cracks in the work area, the surface area of ​​the rock involved in the reaction in the work area, the acid reaction volume in the work area, and the H + The acid-rock reaction rate in the work area is determined by the diffusion coefficient, the velocity and viscosity of the acid injected into the work area, and the concentration gradient of the diffusion boundary layer in the work area, and the relationship model between the acid-rock reaction rate and the fracture width, the surface area of ​​the rock involved in the reaction, the acid reaction volume, the H+ diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer.

9. The device according to claim 8, wherein: The seam width determination module includes: Fractal dimension determination submodule: It is used to determine the fractal dimension of the fracture wall in the work area based on the rock mineral composition, porosity and permeability of the work area, and by using the relationship model between the fractal dimension of the fracture wall and the rock mineral composition, porosity and permeability; Average curvature determination submodule: used to determine the average curvature of the acid-etched crack space in the work area based on the fractal dimension of the crack wall in the work area and using the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall; The crack width determination submodule is used to determine the crack width of the acid-etched crack in the work area based on the average curvature of the acid-etched crack space in the work area and using the relationship model between the crack width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area; Preferably, the relationship model between fracture width and rock mineral composition, porosity, and permeability is: In the formula, Xi is the content of the i-th mineral, dimensionless; is porosity, dimensionless; K is permeability, mD; a i , b1, b2 are coefficients; w is the seam width, cm; More preferably, the rock mineral composition includes N-1 carbonate minerals, X1 to X N-1 are the contents of these N-1 carbonate minerals, X N is the non-carbonate mineral content; Preferably, the relationship model between the average curvature of the acid-etched crack space and the fractal dimension of the crack wall is: K a =mD+n In the formula, K a is the average curvature of the acid-etched crack space, dimensionless; m and n are coefficients; D is the fractal dimension of the crack wall, dimensionless; Preferably, the width of the acid-etched crack and the average curvature of the acid-etched crack space in the work area are: w=pK a +q In the formula, K a is the spatial average curvature of the acid-etched crack, dimensionless; p and q are coefficients; w is the width of the acid-etched crack, dimensionless.

10. The device according to claim 8, wherein: The rock surface area determination module involved in the reaction performs initial correction of the rock surface area involved in the reaction through the following formula: Where s′ is the corrected rock surface area participating in the reaction, cm 2 ; s is the initial rock surface area involved in the reaction, cm 2 ;n i is the number of pores per unit area in the work area, dimensionless; U is the wall roughness coefficient of the cracks formed by the reaction of acid and rock in the work area, dimensionless; d is the average pore diameter in the work area, micrometers.

11. The device according to claim 8, wherein: The acid liquid reaction volume determination module determines the acid liquid reaction volume in the work area through the following formula: Where V' is the acid reaction volume, L; V is the initial acid volume, L; C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, t is the time for the acid to react with the rocks in the work area, s.

12. The device according to claim 8, wherein: The acid reaction volume determination module is used to determine the acid loss coefficient during the reaction between the acid and the rock in the work area based on the porosity, permeability and pore volume of the work area, the viscosity of the acid injected into the work area, and the pressure difference between the fractures in the work area and the formation; Preferably, the acid liquid reaction volume determination module determines the filtration coefficient of the acid liquid during the reaction between the acid liquid and the work area rock by the following formula: Where C1 is the acid loss coefficient during the reaction between the acid and the rock in the work area, is the porosity, dimensionless; K is the permeability, mD; μ is the acid viscosity, mPa·s; Q i is the pore volume of the work area, dimensionless; Δp is the pressure difference between the fracture in the work area and the formation, MPa.

13. The device according to claim 8, wherein: The relationship model between the H+ diffusion coefficient and the temperature of the acid-rock reaction, the acid velocity and concentration is: In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; m1, m2, n1, n2 are coefficients; Preferably, the relationship model between the H+ diffusion coefficient and the temperature of the reaction between the acid and the rock, the velocity and concentration of the acid is: In the formula, H + Diffusion coefficient, m 2 / s; T is the temperature of the reaction between the acid and the rock, °C; C is the concentration of the acid, mol / L; u is the velocity of the acid, m / s; α, β, γ are correction coefficients; More preferably, α=0.4, γ=10.

14. The device according to claim 8, wherein: The relationship model between the acid-rock reaction rate and the fracture width, the rock surface area involved in the reaction, the acid reaction volume, the H+ diffusion coefficient, the velocity and viscosity of the acid, and the concentration gradient of the diffusion boundary layer is: In the formula, H + Diffusion coefficient, m 2 / s; w is the slit width, cm; u is the speed of the acid, m / s; μ is the viscosity of the acid solution, mPa·s; s′ is the surface area of ​​the rock involved in the reaction, cm 2 ; V′ is the volume of acid solution reaction, L; is the concentration gradient of the diffusion boundary layer, mol / (L·cm); K is the acid-rock reaction rate, mol / (cm 2 ·s).

15. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: When the processor executes the computer program, the method for determining the acid-rock reaction rate according to any one of claims 1 to 7 is implemented.

16. A computer-readable storage medium storing a computer program, wherein the computer program, when executed by a processor, implements the method for determining the acid-rock reaction rate according to any one of claims 1 to 7.

17. A computer program product, comprising a computer program, wherein when the computer program is executed by a processor, the method for determining the acid-rock reaction rate according to any one of claims 1 to 7 is implemented.

Images