Method and device for determining flow conductivity of crack and readable storage medium

Abstract

The embodiment of the invention provides a method and device for determining the flow conductivity of a crack and a readable storage medium, and the method comprises the steps: obtaining the laying layer number of a propping agent in the crack at a target moment according to the position information of the supporting particles in the crack at the target moment; according to the laying layer number of the propping agent in the crack at the target moment, the permeability of the crack at the target moment is obtained; according to the initial width of the crack and the embedding amount and the deformation amount of the supporting particles at the target moment, the target width of the crack at the target moment is obtained; and according to the permeability and the target width, the overall flow conductivity of the crack at the target moment is determined. According to the embodiment of the invention, the paving condition and the interaction condition of the proppant in the crack are fully considered, and the flow conductivity of the crack is accurately calculated.

Description

Technical Field

[0001] The embodiments of the present application relate to the technical field of oil and gas field development engineering, and in particular, to a method and device for determining the conductivity of a fracture, and a readable storage medium. Background Art

[0002] Hydraulic fracturing is a common technical means for exploiting low permeability oil and gas reservoirs. During hydraulic fracturing, proppants are injected to support the cracks and prevent them from closing under the action of ground stress. The placement of proppants in the cracks will directly affect the conductivity of the cracks, and thus affect the oil and gas production output of hydraulic fracturing. The placement of proppants in the cracks will directly affect the conductivity of the cracks. However, the current calculation and research on the conductivity of fractures generally directly uses two parallel rock plates to test the conductivity of a single fracture, or performs theoretical calculations based on parallel rock plate models, which ignores the placement of proppants in the fractures. Therefore, the current conductivity of fractures, whether it is actual simulation tests or theoretical simulation calculations, has particularly large limitations, and it is difficult to obtain accurate fracture conductivity. In oil and gas production operations, inaccurate calculations of fracture conductivity will have an adverse effect on a series of processes in oil and gas production operations. Therefore, how to accurately calculate the conductivity of fractures is a problem that technicians urgently need to solve. Summary of the invention

[0003] The embodiments of the present application provide a method, device and readable storage medium for determining the conductivity of a fracture, which are intended to fully consider the placement and interaction of proppants in the fracture and accurately calculate the conductivity of the fracture.

[0004] In one aspect, an embodiment of the present application provides a method for determining the conductivity of a fracture, wherein the fracture is filled with a proppant, and the proppant includes: proppant particles; the method includes:

[0005] According to the position information of the propping particles in the fracture at the target moment, obtaining the number of paving layers of the propping agent in the fracture at the target moment;

[0006] According to the number of paving layers of the proppant in the fracture at the target time, obtaining the permeability of the fracture at the target time;

[0007] Obtaining a target width of the crack at a target time according to an initial width of the crack and an embedding amount and a deformation amount of the supporting particles at a target time;

[0008] The overall conductivity of the fracture at a target moment is determined according to the permeability and the target width.

[0009] Optionally, before the step of obtaining the number of paving layers of the proppant in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment, the step further includes:

[0010] Establishing a simulation model of the proppant in the fracture;

[0011] According to the simulation model, the position information of the supporting particles in the crack at the target moment is obtained.

[0012] Optionally, the step of establishing a simulation model of the proppant in the fracture comprises:

[0013] Based on the method of coupling computational fluid dynamics and discrete element in Euler-Lagrangian model, a coupling model of the fluid flowing between the proppants is established by using preset interaction parameters as a simulation model of the proppants in the fractures;

[0014] Wherein, the preset interaction parameter includes at least one of the following:

[0015] The interaction between the proppants, the interaction between the proppants and the wall of the fracture, and the interaction between the proppants and the fluid.

[0016] Optionally, the step of obtaining the position information of the supporting particles in the crack at a target time according to the simulation model includes:

[0017] Inputting simulation condition parameters into the simulation model, and performing numerical simulation of proppant migration and placement in the fracture at different simulation times to obtain position information of the proppant particles in the fracture at a target time;

[0018] The simulation condition parameters include at least one of the following:

[0019] Proppant parameters, fluid parameters, fracturing pumping process parameters of the fluid, roughness parameters or smoothness parameters of the fracture wall and branching seam parameters.

[0020] Optionally, the step of obtaining the number of paving layers of the proppant in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment comprises:

[0021] Dividing a plurality of orthogonal grids on the fracture surface of the fracture;

[0022] Determining the number of supporting particles in each of the orthogonal grids at the target time according to the orthogonal grids and the position information of the supporting particles at the target time;

[0023] The number of layers of the proppant in the fracture at the target moment is determined according to the number of proppant particles in each of the grids at the target moment, the diameter of the proppant particles and the size of the grid.

[0024] Optionally, the step of obtaining the permeability of the fracture at the target moment according to the number of paving layers of the proppant in the fracture at the target moment comprises:

[0025] Obtaining the porosity of the proppant in the fracture at the target time according to the diameter of the proppant, the position information of the proppant in the fracture at the target time, the number of the proppant in each of the orthogonal grids at the target time, the deformation of the proppant, and the zero-pressure aperture of the fracture;

[0026] Obtaining the tortuosity of the proppant in the fracture at the target time according to the diameter of the proppant, the position information of the proppant in the fracture at the target time, the number of the proppant in each of the orthogonal grids at the target time, the deformation of the proppant, and the zero-pressure aperture of the fracture;

[0027] Obtaining the roaring diameter of the proppant in the fracture at the target time according to the diameter of the proppant particles, the position information of the proppant particles in the fracture at the target time, the number of the proppant particles in each of the orthogonal grids at the target time, the embedding amount of the proppant, the deformation amount of the proppant, the zero-pressure opening of the fracture, and the roaring diameter at each of the grids;

[0028] The permeability of the fracture is obtained according to the porosity of the proppant in the fracture, the tortuosity and the roar diameter.

[0029] Optionally, before the step of obtaining a target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the support particles at a target moment, the step includes:

[0030] Obtaining the embedding amount of the support particles at a target time according to the diameter of the support particles, the spacing coefficient of the support particles, the extrusion stress between the support particles, the Poisson's ratio of the support particles, the Poisson's ratio of the wall of the crack, the elastic model of the support particles, the elastic model of the wall of the crack, the thickness of the wall of the crack, the viscosity coefficient of the support particles, and the viscosity coefficient of the wall of the crack;

[0031] The deformation amount of the support particles at the target moment is obtained according to the diameter of the support particles, the spacing coefficient of the support particles, the extrusion stress between the support particles, the elastic model of the support particles, the Poisson's ratio of the support particles, and the viscosity coefficient of the support particles.

[0032] Optionally, the step of determining the overall conductivity of the fracture at a target moment according to the permeability and the target width includes:

[0033] Determine the conductivity of each location in the crack at the target time according to the product of the permeability of the crack at the target time, the target width and a preset correction coefficient;

[0034] The conductivity of each location in the fracture at the target time is harmonized and averaged to determine the overall conductivity of the fracture at the target time.

[0035] In another aspect, an embodiment of the present application provides a device for determining the conductivity of a fracture, wherein the fracture is filled with a proppant, and the proppant includes: proppant particles; the device includes:

[0036] A positioning unit, used for obtaining the number of paving layers of the proppant in the fracture at the target moment according to the position information of the proppant particles in the fracture at the target moment;

[0037] a permeability calculation unit, for obtaining the permeability of the fracture at the target moment according to the number of paving layers of the proppant in the fracture at the target moment;

[0038] a width calculation unit, configured to obtain a target width of the crack at a target moment according to an initial width of the crack and an embedding amount and a deformation amount of the supporting particles at a target moment;

[0039] The conductivity calculation unit is used to determine the overall conductivity of the fracture at a target moment according to the permeability and the target width.

[0040] On the other hand, another embodiment of the present application provides a readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the steps in the method of any of the above embodiments of the present application.

[0041] On the other hand, another embodiment of the present application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, which, when executed by the processor, implements the steps in the method of any of the above embodiments of the present application.

[0042] On the other hand, another embodiment of the present application provides a computer program product comprising instructions, and when the computer program product is run on a signal detection device, the signal detection device executes the steps of the signal detection method in any of the above embodiments.

[0043] On the other hand, another embodiment of the present application provides a chip, the chip includes a processor and a communication interface, the communication interface and the processor are coupled, and the processor is used to run a computer program or instructions to implement the steps of the signal detection method in any of the above embodiments.

[0044] Optionally, the chip provided in the embodiment of the present application also includes a memory for storing computer programs or instructions.

[0045] Compared with the prior art, the advantages of the embodiments of the present application are:

[0046] Considering that the position distribution and laying concentration of proppant in the crack are closely related to the conductivity of the crack, in the embodiment of the present application, when the crack is filled with proppant, and the proppant includes proppant particles, the number of laying layers of the proppant in the crack at the target time is obtained according to the position information of the proppant particles in the crack at the target time; the permeability of the crack at the target time is obtained according to the number of laying layers of the proppant in the crack at the target time; the target width of the crack at the target time is obtained according to the initial width of the crack and the embedding amount and deformation amount of the proppant particles at the target time; the overall conductivity of the crack at the target time is determined according to the permeability and the target width. In the case of fully considering the laying conditions and interaction of the proppant in the crack, the embedding and deformation of the proppant particles are calculated by simulating the change of the laying conditions of the proppant over time, so as to obtain the change of the width of the crack, and the conductivity of the crack is calculated in combination with the permeability of the crack, which can effectively improve the accuracy of the calculation of the conductivity of the crack, and then effectively improve the work efficiency of hydraulic fracturing oil and gas production operations. BRIEF DESCRIPTION OF THE DRAWINGS

[0047] The accompanying drawings are for reference and illustration only and are not intended to limit the scope of protection of this application. The technical solutions in the embodiments of this application will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of this application. Obviously, the described embodiments are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative work are within the scope of protection of this application.

[0048] Figure 1 A flowchart showing a method for determining the conductivity of a fracture in an embodiment provided by the present application is shown;

[0049] Figure 2 A flowchart showing a method for determining the conductivity of a fracture in another embodiment provided by the present application is shown;

[0050] Figure 3A flowchart showing a method for determining the conductivity of a fracture in another embodiment provided by the present application is shown;

[0051] Figure 4 A schematic diagram of orthogonal grid division in another embodiment provided by the present application is shown;

[0052] Figure 5 A flowchart showing a method for determining the conductivity of a fracture in another embodiment provided by the present application is shown;

[0053] Figure 6 A flowchart showing a method for determining the conductivity of a fracture in another embodiment provided by the present application is shown;

[0054] Figure 7 A flowchart showing a method for determining the conductivity of a fracture in another embodiment provided by the present application is shown;

[0055] Figure 8 A flowchart of determining the conductivity of a fracture in another embodiment provided by the present application is shown;

[0056] Fig. 9 A structural block diagram of a device for determining the conductivity of a fracture in an embodiment provided by the present application is shown. DETAILED DESCRIPTION

[0057] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0058] In the process of exploiting low permeability oil and gas reservoirs by hydraulic fracturing, proppants are often pre-injected to support the fractures. However, the placement of proppants at different locations in the fracture is often non-uniform. In some locations, the number of proppants is large, resulting in a high concentration of sand, while in some locations, the number of proppants is small, resulting in a low concentration of sand, and there is even no proppants in some locations. Therefore, ignoring the placement of proppants in the fracture, a data result obtained by simply testing with parallel rock plates or theoretically calculating based on the parallel rock plate model cannot characterize the conductivity of the entire fracture, but can only characterize the fracture conductivity under the sand concentration. In other words, the current theoretical calculation of conductivity has limitations, and directly ignores the placement of proppants in the fracture. In actual situations, the fracture conductivity is closely related to the placement of proppants in the fracture, including factors such as the distribution of proppant positions and the concentration of proppants. The placement of proppants in the fracture should be fully considered when calculating the fracture conductivity.

[0059] In view of this, the embodiments of the present application provide a method, device and readable storage medium for determining the conductivity of a fracture. When the fracture is filled with proppant, and the proppant includes proppant particles, the number of proppant layers in the fracture at the target moment is obtained according to the position information of the proppant particles in the fracture at the target moment; the permeability of the fracture at the target moment is obtained according to the number of proppant layers in the fracture at the target moment; the target width of the fracture at the target moment is obtained according to the initial width of the fracture and the embedding amount and deformation amount of the proppant particles at the target moment; the overall conductivity of the fracture at the target moment is determined according to the permeability and the target width. Under the condition of fully considering the laying condition and interaction of the proppant in the fracture, the embedding and deformation of the proppant particles are calculated by simulating the change of the laying condition of the proppant over time, so as to obtain the change of the width of the fracture, and the conductivity of the fracture is calculated in combination with the permeability of the fracture, which can effectively improve the accuracy of the calculation of the conductivity of the fracture, and then effectively improve the work efficiency of hydraulic fracturing oil and gas production operations.

[0060] The embodiments of the present application are described below with reference to the accompanying drawings.

[0061] Embodiment 1

[0062] Reference Figure 1 , Figure 1 FIG. 1 is a flowchart showing a method for determining the conductivity of a fracture in an embodiment of the present application. Figure 1 As shown, an embodiment of the present application provides a method for determining the conductivity of a fracture, which can be applied to fractures, specifically fractures in low permeability oil and gas reservoirs that are exploited by hydraulic fracturing in oil and gas production operations.

[0063] Specifically, the fracture is filled with proppant, and the proppant includes: proppant particles.

[0064] The types of proppants and supporting particles therein can be selected according to the material properties required for oil and gas production, and will not be elaborated here.

[0065] The methods specifically include:

[0066] Step S210, obtaining the number of proppant layers in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment.

[0067] Step S220, obtaining the permeability of the fracture at the target moment according to the number of proppant layers laid in the fracture at the target moment.

[0068] Step S230 , obtaining a target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target moment.

[0069] Step S240, determining the overall conductivity of the fracture at the target time according to the permeability and the target width.

[0070] Through the above embodiments, considering that the position distribution and laying concentration of proppant in the crack are closely related to the conductivity of the crack, the present application obtains the number of laying layers of proppant in the crack at the target time according to the position information of the supporting particles in the crack at the target time when the crack is filled with proppant, and the proppant includes supporting particles; obtains the permeability of the crack at the target time according to the number of laying layers of proppant in the crack at the target time; obtains the target width of the crack at the target time according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target time; determines the overall conductivity of the crack at the target time according to the permeability and target width. In the case of fully considering the laying conditions and interaction of the proppant in the crack, the embedding and deformation of the supporting particles are calculated by simulating the change of the laying conditions of the proppant over time, so as to obtain the change of the width of the crack, and the conductivity of the crack is calculated in combination with the permeability of the crack, which can effectively improve the accuracy of the calculation of the conductivity of the crack, and then can effectively improve the work efficiency of hydraulic fracturing oil and gas production operations.

[0071] Embodiment 2

[0072] Reference Figure 2 , Figure 2 FIG. 2 is a flowchart showing a method for determining the conductivity of a fracture in another embodiment of the present application. Figure 2 As shown, the embodiment of the present application provides another method for determining the conductivity of a fracture, which can be applied to fractures, specifically fractures in low permeability oil and gas reservoirs that are exploited by hydraulic fracturing in oil and gas production operations.

[0073] Specifically, the fracture is filled with proppant, and the proppant includes: proppant particles.

[0074] The types of proppants and supporting particles therein can be selected according to the material properties required for oil and gas production, and will not be elaborated here.

[0075] The methods specifically include:

[0076] Step S210, obtaining the number of proppant layers in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment.

[0077] In the embodiment of the present application, the following steps may be further included before step S210:

[0078] Step S110, establishing a simulation model of proppant in the fracture.

[0079] Step S120, obtaining position information of the supporting particles in the crack at the target moment according to the simulation model.

[0080] Furthermore, in the embodiment of the present application, step S110 may specifically include:

[0081] Step S111, based on the computational fluid dynamics and discrete element coupling method in the Euler-Lagrangian model, using preset interaction parameters, a coupling model of fluid flow between proppants is established as a simulation model of proppants in fractures.

[0082] The preset interaction parameters include at least one of the following:

[0083] The interaction between proppants and proppants, the interaction between proppants and the wall of the fracture, and the interaction between proppants and fluids.

[0084] In the embodiments of the present application, the simulation model of proppant in the fracture can be a discrete element coupled CFD-DEM model. The present application can establish a coupled model of proppant flow in the fracturing fluid by considering the interaction between proppant and proppant, proppant and fracture wall, and proppant and fracturing fluid, and simulate the proppant migration and laying morphology in the fracture, so as to derive the position information of each proppant particle.

[0085] Step S112, inputting simulation condition parameters into the simulation model, and performing numerical simulation of proppant migration and placement in the fracture at different simulation times to obtain position information of proppant particles in the fracture at the target time.

[0086] The simulation condition parameters include at least one of the following:

[0087] Proppant parameters, fluid parameters, fluid fracturing pumping process parameters, fracture wall roughness parameters or smoothness parameters and branch fracture parameters.

[0088] In an embodiment of the present application, the position information of each supporting particle may change with time. By setting different simulation times for the simulation model, the position coordinates of the supporting particles at each time node can be derived under the above-mentioned simulation condition parameters, which can then be used to characterize the position information of the supporting particles in the crack at the target moment.

[0089] Step S202, obtaining the permeability of the fracture at the target moment according to the number of proppant layers laid in the fracture at the target moment.

[0090] Step S203, obtaining a target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target moment.

[0091] Step S204, determining the overall conductivity of the fracture at the target time according to the permeability and the target width.

[0092] Through the above embodiments, considering that the position distribution and laying concentration of proppant in the crack are closely related to the conductivity of the crack, the present application establishes a discrete element coupling CFD-DEM model when the crack is filled with proppant, and the proppant includes proppant particles, fully considers the interaction between proppant and proppant, proppant and crack wall, and proppant and fracturing fluid, establishes a coupling model of proppant flow in fracturing fluid, simulates the proppant migration and laying morphology in the crack, and derives the position information of each supporting particle. Then, according to the position information of the supporting particles in the crack at the target time, the number of laying layers of the proppant in the crack at the target time is obtained; according to the number of laying layers of the proppant in the crack at the target time, the permeability of the crack at the target time is obtained; according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target time, the target width of the crack at the target time is obtained; according to the permeability and target width, the overall conductivity of the crack at the target time is determined. Under full consideration of the placement and interaction of proppants in the fractures, the embedding and deformation of proppant particles are calculated by simulating the changes in the placement of proppants over time, thereby obtaining the change in the width of the fractures. Combined with the permeability of the fractures, the conductivity of the fractures is calculated, which can effectively improve the accuracy of the calculation of the conductivity of the fractures, and then effectively improve the work efficiency of hydraulic fracturing oil and gas production operations.

[0093] Embodiment 3

[0094] Reference Figure 3 , Figure 3 FIG. 2 is a flowchart showing a method for determining the conductivity of a fracture in another embodiment of the present application. Figure 3 As shown, the embodiment of the present application provides another method for determining the conductivity of a fracture, which can be applied to fractures, specifically fractures in low permeability oil and gas reservoirs that are exploited by hydraulic fracturing in oil and gas production operations.

[0095] Specifically, the fracture is filled with proppant, and the proppant includes: proppant particles.

[0096] The types of proppants and supporting particles therein can be selected according to the material properties required for oil and gas production, and will not be elaborated here.

[0097] The methods specifically include:

[0098] Step S210, obtaining the number of proppant layers in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment.

[0099] Reference Figure 4 , Figure 4 FIG. 2 shows a schematic diagram of orthogonal grid division in another embodiment provided by the present application. Figure 4As shown, in the embodiment of the present application, step S210 may specifically include the following steps:

[0100] Step S211, dividing a plurality of orthogonal grids on the crack surface of the crack.

[0101] The crack surface may refer to a virtual plane obtained by overlapping the crack opening with the plane where the opening is located, and the boundary of the crack surface is the opening boundary of the crack.

[0102] Among them, orthogonal grids can be divided on the crack surface, and the grid size of each orthogonal grid is set to dx×dy.

[0103] Step S212, determining the number of supporting particles in each orthogonal grid at the target moment according to the orthogonal grid and the position information of the supporting particles at the target moment.

[0104] like Figure 4 As shown in the figure, assuming that the proppant in each grid is evenly distributed and stacked in a simple cubic manner, the number of proppant layers can be simply estimated based on the number of proppant particles in each grid. Furthermore, after obtaining the number of proppant particles in each grid at the target time, the number of proppant layers in the fracture can be more accurately obtained.

[0105] Step S213, determining the number of proppant layers to be laid in the fracture at the target time according to the number of proppant particles in each grid at the target time, the diameter of the proppant particles and the size of the grid.

[0106] In an optional example, the number of proppant layers in the fracture at the target time can be obtained according to formula (1):

[0107]

[0108] Where t is the target time; i represents the i-th grid in the X direction; j represents the j-th grid in the Z direction; () int represents rounding; dx is the grid length in the X direction, and the unit can be uniformly set to meters; dy is the grid length in the Y direction, and the unit can be uniformly set to meters; n ij (t) is the number of proppant layers at the grid of row i and column j at time t; N ij (t) is the number of supporting particles at the grid of row i and column j at time t; D p To support particle diameter, the unit can be uniformly set to meters.

[0109] Step S202, obtaining the permeability of the fracture at the target moment according to the number of proppant layers laid in the fracture at the target moment.

[0110] Step S203, obtaining a target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target moment.

[0111] Step S204, determining the overall conductivity of the fracture at the target time according to the permeability and the target width.

[0112] Through the above embodiments, the present application divides the orthogonal grids, analyzes and counts the number of supporting particles falling into each grid, and combines the size of the grid and the diameter of the supporting particles to obtain the number of proppant laying layers at different times in the cracks; then obtain the permeability of the cracks at the target time according to the number of proppant laying layers in the cracks at the target time; obtain the target width of the cracks at the target time according to the initial width of the cracks and the embedding amount and deformation amount of the supporting particles at the target time; determine the overall conductivity of the cracks at the target time according to the permeability and target width. In the case of fully considering the laying conditions and interactions of the proppant in the cracks, the embedding and deformation of the supporting particles are calculated by simulating the changes in the laying conditions of the proppant over time, thereby obtaining the change in the width of the cracks, and combining the permeability of the cracks to calculate the conductivity of the cracks, which can effectively improve the accuracy of the calculation of the conductivity of the cracks, and then effectively improve the work efficiency of hydraulic fracturing oil and gas production operations.

[0113] Embodiment 4

[0114] Reference Figure 5 , Figure 5 FIG. 2 is a flowchart showing a method for determining the conductivity of a fracture in another embodiment of the present application. Figure 5 As shown, the embodiment of the present application provides another method for determining the conductivity of a fracture, which can be applied to fractures, specifically fractures in low permeability oil and gas reservoirs that are exploited by hydraulic fracturing in oil and gas production operations.

[0115] Specifically, the fracture is filled with proppant, and the proppant includes: proppant particles.

[0116] The types of proppants and supporting particles therein can be selected according to the material properties required for oil and gas production, and will not be elaborated here.

[0117] The methods specifically include:

[0118] Step S210, obtaining the number of proppant layers in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment.

[0119] Step S220, obtaining the permeability of the fracture at the target moment according to the number of proppant layers laid in the fracture at the target moment.

[0120] Assuming that the proppant is evenly laid in the fracture, the fracture permeability is a function of the fracture porosity, the roaring channel diameter and the tortuosity. Considering the embedding and deformation of the proppant, the porosity, the tortuosity and the roaring channel diameter must be dynamically corrected. To this end, in the embodiment of the present application, step S220 may specifically include the following steps:

[0121] Step S221, obtaining the porosity of the proppant in the fracture at the target time according to the diameter of the proppant, the position information of the proppant in the fracture at the target time, the number of the proppant in each orthogonal grid at the target time, the deformation of the proppant and the zero-pressure aperture of the fracture.

[0122] In an optional example, the porosity of the proppant in the fracture at the target time can be calculated according to formula (2):

[0123]

[0124] Where t is the target time; i represents the i-th grid in the X direction; j represents the j-th grid in the Z direction; dx is the grid length in the X direction, and the unit can be uniformly set to meters; dy is the grid length in the Y direction, and the unit can be uniformly set to meters; φ ij (t) is the porosity at the grid of row i and column j at time t, dimensionless; D is the zero-pressure aperture of the fracture, that is, the aperture of the fracture when the closing pressure is 0 and the fracture is just closed and pressurized onto the proppant, and the unit can be uniformly set to meter; D p To support particle diameter, the unit can be uniformly set to meter; L d (t) is the deformation of the proppant at time t; N ij (t) is the number of supporting particles at the grid of row i and column j at time t.

[0125] Step S222, obtaining the tortuosity of the proppant in the fracture at the target time according to the diameter of the proppant, the position information of the proppant in the fracture at the target time, the number of the proppant in each orthogonal grid at the target time, the deformation of the proppant and the zero-pressure aperture of the fracture.

[0126] In an optional example, the tortuosity of the proppant in the fracture at the target time can be calculated according to formula (3) and formula (4):

[0127]

[0128]

[0129] Where t is the target time; i represents the i-th grid in the X direction; j represents the j-th grid in the Z direction; dx is the grid length in the X direction, and the unit can be uniformly set to meters; dy is the grid length in the Y direction, and the unit can be uniformly set to meters; τ ij(t) is the tortuosity at the grid of row i and column j at time t, dimensionless; D is the zero-pressure opening of the fracture, that is, the opening of the fracture when the closing pressure is 0 and the fracture is just closed and pressurized onto the proppant, and the unit can be uniformly set to meter; L d (t) is the deformation of the proppant at time t; N ij (t) is the number of supporting particles at the grid of row i and column j at time t; D p To support particle diameter, the unit can be uniformly set to meters.

[0130] Step S223, according to the diameter of the supporting particles, the position information of the supporting particles in the fracture at the target time, the number of supporting particles in each orthogonal grid at the target time, the embedding amount of the proppant, the deformation amount of the proppant, the zero-pressure opening of the fracture and the number of laying layers of the proppant at each orthogonal grid at the target time, the roar diameter of the proppant in the fracture at the target time is obtained.

[0131] In an optional example, the roar diameter of the proppant in the fracture at the target time can be calculated according to formula (5):

[0132]

[0133] Wherein, t is the target time; i represents the i-th grid in the X direction; j represents the j-th grid in the Z direction; dx is the grid length in the X direction, and the unit can be uniformly set to meters; dy is the grid length in the Y direction, and the unit can be uniformly set to meters; D is the zero-pressure opening of the fracture, that is, the opening of the fracture when the closing pressure is 0 and the fracture is just closed and pressurized onto the proppant, and the unit can be uniformly set to meters; r ij (t) is the diameter of the roar channel at the grid of row i and column j at time t, and the unit can be uniformly set to meter; N in is the number of supporting particles embedded in the crack wall at the i-row and j-column grid; C B is a coefficient related to the number of proppant layers n.

[0134] In an alternative example, C B There may be the following relationship between and the number of proppant layers n:

[0135] If n=3m-2, then B=2(m-1); if n=3m-1, then B=2(m-1); if n=3m, B=2m-1.

[0136] Step S224, obtaining the permeability of the fracture according to the porosity, tortuosity and roar diameter of the proppant in the fracture.

[0137] In an optional example, the permeability of the fracture can be calculated according to formula (6):

[0138]

[0139] Where t is the target time; i represents the i-th grid in the X direction; j represents the j-th grid in the Z direction; φ ij (t) is the porosity of the grid at row i and column j at time t, dimensionless; τ ij (t) is the tortuosity of the grid at row i and column j at time t, dimensionless; r ij (t) is the diameter of the roar channel at the grid of row i and column j at time t, and the unit can be uniformly set to meter; k ij (t) is the permeability of the crack at the grid in row i and column j at time t.

[0140] Step S230 , obtaining a target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target moment.

[0141] Step S240, determining the overall conductivity of the fracture at the target time according to the permeability and the target width.

[0142] Through the above embodiments, considering that the position distribution and laying concentration of proppant in the crack are closely related to the conductivity of the crack, the present application obtains the number of laying layers of proppant in the crack at the target time according to the position information of the supporting particles in the crack at the target time when the crack is filled with proppant, and the proppant includes supporting particles; obtains the permeability of the crack at the target time according to the number of laying layers of proppant in the crack at the target time; obtains the target width of the crack at the target time according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target time; determines the overall conductivity of the crack at the target time according to the permeability and target width. In the case of fully considering the laying conditions and interaction of the proppant in the crack, the embedding and deformation of the supporting particles are calculated by simulating the change of the laying conditions of the proppant over time, so as to obtain the change of the width of the crack, and the conductivity of the crack is calculated in combination with the permeability of the crack, which can effectively improve the accuracy of the calculation of the conductivity of the crack, and then can effectively improve the work efficiency of hydraulic fracturing oil and gas production operations.

[0143] Embodiment 5

[0144] Reference Figure 6 , Figure 6 FIG. 2 is a flowchart showing a method for determining the conductivity of a fracture in another embodiment of the present application. Figure 6 As shown, an embodiment of the present application provides a method for determining the conductivity of a fracture, which can be applied to fractures, specifically fractures in low permeability oil and gas reservoirs that are exploited by hydraulic fracturing in oil and gas production operations.

[0145] Specifically, the fracture is filled with proppant, and the proppant includes: proppant particles.

[0146] The types of proppants and supporting particles therein can be selected according to the material properties required for oil and gas production, and will not be elaborated here.

[0147] The methods specifically include:

[0148] Step S210, obtaining the number of proppant layers in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment.

[0149] Step S220, obtaining the permeability of the fracture at the target moment according to the number of proppant layers laid in the fracture at the target moment.

[0150] Under the crack closure pressure, the supporting particles and the supporting particles and the crack wall will deform and embed, thereby causing the crack width to change. To this end, in the embodiment of the present application, the target width of the crack at each time node can be more accurately calculated based on the deformation and embedding amount of the supporting particles.

[0151] In some optional embodiments, the deformation amount and embedding amount of the support particles may be calculated by the following steps, specifically including:

[0152] Step 2301, based on the diameter of the supporting particles, the spacing coefficient of the supporting particles, the extrusion stress between the supporting particles, the Poisson's ratio of the supporting particles, the Poisson's ratio of the crack wall, the elastic model of the supporting particles, the elastic model of the crack wall, the thickness of the crack wall, the viscosity coefficient of the supporting particles and the viscosity coefficient of the crack wall, obtain the embedded amount of the supporting particles at the target moment.

[0153] In an optional example, the embedding amount of the support particles at the target time can be calculated according to formula (7):

[0154]

[0155] Where t is time, the unit can be set to seconds; L in (t) is the embedding amount of proppant at time t, and the unit can be uniformly set to meter; D p To support particle diameter, the unit can be uniformly set to meter; C k is the particle spacing coefficient, dimensionless; σ f is the extrusion stress, the unit can be set to MPa; E f is the elastic model of the supporting particles, the unit can be set to MPa; E p is the elastic model of the crack wall, the unit can be set to MPa; v p is the Poisson's ratio of the supporting particles; v f is the Poisson's ratio of the crack wall; D f is the crack wall thickness, the unit can be uniformly set to meter; S pis the viscosity coefficient of the supporting particles, dimensionless; S f is the viscosity coefficient of the crack wall and is dimensionless.

[0156] Step 2302, obtaining the deformation of the support particles at the target moment according to the diameter of the support particles, the spacing coefficient of the support particles, the extrusion stress between the support particles, the elastic model of the support particles, the Poisson's ratio of the support particles, and the viscosity coefficient of the support particles.

[0157] In an optional example, the deformation of the supporting particles at the target time can be calculated according to formula (8):

[0158]

[0159] Where t is time, the unit can be set to seconds; L d (t) is the deformation of the proppant at time t, and the unit can be uniformly set to meter; D p is the diameter of the supporting particles, the unit can be uniformly set to meter; C k is the particle spacing coefficient, dimensionless; σ f is the extrusion stress, the unit can be set to MPa; E p is the elastic model of the crack wall, the unit can be set to MPa; v p is the Poisson's ratio of the supporting particles; S p is the viscosity coefficient of the supporting particles and is dimensionless.

[0160] Step S230 , obtaining a target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target moment.

[0161] In combination with the above implementation example, in an optional example, the target width of the crack at the target time can be calculated according to formula (9):

[0162] w(t)=w ini -2(L in (t)+L d (t)) (9)

[0163] Where t is time, the unit can be set to seconds; L in (t) is the embedding amount of proppant at time t, and the unit can be uniformly set to meter; L d (t) is the deformation of the proppant at time t, and the unit can be uniformly set to meter; w ini is the initial width of the crack, and the unit can be uniformly set to meter; w(t) is the target width of the crack at the target time, that is, the width of the crack at time t, and the unit can be uniformly set to meter.

[0164] Step S240, determining the overall conductivity of the fracture at the target time according to the permeability and the target width.

[0165] Through the above embodiments, considering that the position distribution and laying concentration of proppant in the crack are closely related to the conductivity of the crack, the present application obtains the number of laying layers of proppant in the crack at the target time according to the position information of the supporting particles in the crack at the target time when the crack is filled with proppant, and the proppant includes supporting particles; obtains the permeability of the crack at the target time according to the number of laying layers of proppant in the crack at the target time; obtains the target width of the crack at the target time according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target time; determines the overall conductivity of the crack at the target time according to the permeability and target width. In the case of fully considering the laying conditions and interaction of the proppant in the crack, the embedding and deformation of the supporting particles are calculated by simulating the change of the laying conditions of the proppant over time, so as to obtain the change of the width of the crack, and the conductivity of the crack is calculated in combination with the permeability of the crack, which can effectively improve the accuracy of the calculation of the conductivity of the crack, and then can effectively improve the work efficiency of hydraulic fracturing oil and gas production operations.

[0166] Embodiment 6

[0167] Reference Figure 7 , Figure 7 FIG. 2 is a flowchart showing a method for determining the conductivity of a fracture in another embodiment of the present application. Figure 7 As shown, an embodiment of the present application provides a method for determining the conductivity of a fracture, which can be applied to fractures, specifically fractures in low permeability oil and gas reservoirs that are exploited by hydraulic fracturing in oil and gas production operations.

[0168] Specifically, the fracture is filled with proppant, and the proppant includes: proppant particles.

[0169] The types of proppants and supporting particles therein can be selected according to the material properties required for oil and gas production, and will not be elaborated here.

[0170] The methods specifically include:

[0171] Step S210, obtaining the number of proppant layers in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment.

[0172] Step S220, obtaining the permeability of the fracture at the target moment according to the number of proppant layers laid in the fracture at the target moment.

[0173] Step S230 , obtaining a target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target moment.

[0174] Step S240, determining the overall conductivity of the fracture at the target time according to the permeability and the target width.

[0175] In some optional embodiments, step S240 may specifically include:

[0176] Step S241, determining the conductivity of each location in the crack at the target time according to the product of the crack permeability at the target time, the target width and a preset correction coefficient.

[0177] Fracture conductivity is used to evaluate the flow capacity of fluids. Fracture conductivity can be defined as the product of fracture permeability and fracture width.

[0178] In an optional example, the conductivity of each location in the crack at the target time may be equal to the product of the permeability of the crack at the target time, the target width, and a preset correction coefficient.

[0179] In another optional example, referring to formula (6) and formula (9), the conductivity of the fracture can be obtained according to the initial width of the fracture, the deformation and shallow penetration of the proppant at time t, the porosity at the grid of row i and column j at time t, the tortuosity at the grid of row i and column j at time t, and the roar channel diameter at the grid of row i and column j at time t. To this end, the conductivity of each point in the fracture at the target time can be calculated according to formula (10):

[0180]

[0181] Where t is the time, the unit can be set to seconds; i represents the i-th grid in the X direction; j represents the j-th grid in the Z direction; L in (t) is the embedding amount of proppant at time t, and the unit can be uniformly set to meter; L d (t) is the deformation of the proppant at time t, and the unit can be uniformly set to meter; w ini is the initial width of the crack, and the unit can be uniformly set to meters; w(t) is the target width of the crack at the target time, that is, the width of the crack at time t, and the unit can be uniformly set to meters; FCD ij (t) is the fracture conductivity of the grid at row i and column j at time t, and the unit can be set to μm 2 cm; c0 is the correction factor, dimensionless; φ ij (t) is the porosity at the grid of row i and column j at time t, dimensionless; D is the zero-pressure aperture of the fracture, that is, the aperture of the fracture when the closing pressure is 0 and the fracture is just closed and pressurized onto the proppant, and the unit can be uniformly set to meter; τ ij (t) is the tortuosity of the grid at row i and column j at time t, dimensionless; r ij (t) is the diameter of the roar channel at the grid of row i and column j at time t, and the unit can be uniformly set to meter.

[0182] Step S242, harmonizing and averaging the conductivity of each location in the crack at the target time to determine the overall conductivity of the crack at the target time.

[0183] In an optional example, the overall conductivity of the fracture at the target time can be calculated by referring to formula (11):

[0184]

[0185] Where t is the time, the unit can be set to seconds; i represents the i-th grid in the X direction; j represents the j-th grid in the Z direction; FCD ij (t) is the fracture conductivity of the grid at row i and column j at time t, and the unit can be set to μm 2 ·cm; FCD(t) is the overall conductivity of the crack at time t.

[0186] Considering that during the oil and gas production operation, the conductivity of each fracture can also be intuitively displayed, which can be displayed by drawing a distribution cloud map. To this end, in some optional embodiments, the following step S243 can also be used to draw a distribution cloud map of the conductivity of the fracture.

[0187] Step S243, performing linear interpolation according to the conductivity of each location in the crack at the target time, and drawing a distribution cloud diagram of the conductivity of the crack based on the crack surface of the crack.

[0188] Through this embodiment, the conductivity of each part of the crack at the target moment can be intuitively displayed.

[0189] Through the above embodiments, considering that the position distribution and laying concentration of proppant in the crack are closely related to the conductivity of the crack, the present application obtains the number of laying layers of proppant in the crack at the target time according to the position information of the supporting particles in the crack at the target time when the crack is filled with proppant, and the proppant includes supporting particles; obtains the permeability of the crack at the target time according to the number of laying layers of proppant in the crack at the target time; obtains the target width of the crack at the target time according to the initial width of the crack and the embedding amount and deformation amount of the supporting particles at the target time; determines the overall conductivity of the crack at the target time according to the permeability and target width. In the case of fully considering the laying conditions and interaction of the proppant in the crack, the embedding and deformation of the supporting particles are calculated by simulating the change of the laying conditions of the proppant over time, so as to obtain the change of the width of the crack, and the conductivity of the crack is calculated in combination with the permeability of the crack, which can effectively improve the accuracy of the calculation of the conductivity of the crack, and then can effectively improve the work efficiency of hydraulic fracturing oil and gas production operations.

[0190] Through the above embodiments, the problem that the current theoretical calculation of fracture conductivity ignores the placement of proppants in fractures can be solved, and the fracture conductivity when the proppants are placed at different positions in the fracture can be more intuitively and truly reflected.

[0191] Reference Figure 8 , Figure 8 A flowchart of determining the conductivity of a fracture in another embodiment of the present application is shown. Referring to Table 1, Table 1 is a basic parameter table for calculating the conductivity of a fracture provided in an embodiment of the present application.

[0192] Table 1

[0193]

[0194]

[0195] like Figure 8 As shown in Table 1, the present application also provides an example, which can use the parameters in Table 1 in combination with the above embodiments to determine the long-term conductivity of the fracture.

[0196] Step 401, simulating the proppant migration and placement morphology in the fracture.

[0197] Step 402, deriving proppant position coordinate information.

[0198] Step 403: Count the number of proppant particles in each grid, estimate the number of proppant layers, calculate the proppant filling porosity, tortuosity and roar diameter at each grid, and calculate the permeability.

[0199] Step 404, calculate the deformation of the proppant embedding and calculate the change of the seam width.

[0200] The execution order of step 403 and step 403 may not be in particular order or may be performed simultaneously.

[0201] Step 405, calculate the conductivity at each grid, obtain the conductivity at any position on the fracture surface based on the linear difference, and draw a conductivity distribution cloud map.

[0202] Step 406, calculating the harmonic mean of the conductivity at all grids to obtain the overall conductivity of the fracture.

[0203] Embodiment 7

[0204] Reference Fig. 9 , Fig. 9 FIG. 1 is a block diagram showing a device for determining the conductivity of a fracture in an embodiment of the present application. Fig. 9 As shown, combined Figure 3As shown, based on the same inventive concept, the embodiment of the present application further provides a device 600 for determining the conductivity of a fracture, wherein the fracture is filled with a proppant, and the proppant includes: supporting particles; the device includes:

[0205] A positioning unit 601 is used to obtain the number of proppant layers in the fracture at the target moment according to the position information of the proppant particles in the fracture at the target moment;

[0206] The permeability calculation unit 602 is used to obtain the permeability of the fracture at the target time according to the number of proppant layers in the fracture at the target time;

[0207] A width calculation unit 603 is used to obtain a target width of the crack at a target time according to an initial width of the crack and an embedding amount and a deformation amount of supporting particles at a target time;

[0208] The conductivity calculation unit 604 is used to determine the overall conductivity of the fracture at a target time according to the permeability and the target width.

[0209] Based on the same inventive concept, an embodiment of the present application provides a computer program product comprising instructions. When the computer program product is run on a computer, the computer is enabled to execute the perforation signal detection method in the above method embodiment.

[0210] Based on the same inventive concept, an embodiment of the present application provides a readable storage medium, in which instructions are stored, and is characterized in that when a computer executes the instructions, the computer executes the steps of the signal detection method in any of the above embodiments.

[0211] Among them, the readable storage medium is a computer-readable storage medium.

[0212] Among them, the computer readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, device or device, or any combination of the above. More specific examples of computer readable storage media (a non-exhaustive list) include: an electrical connection with one or more wires, a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM), a register, a hard disk, an optical fiber, a portable compact disk read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above, or any other form of computer readable storage medium known in the art. An exemplary storage medium is coupled to a processor so that the processor can read information from the storage medium and write information to the storage medium. Of course, the storage medium can also be a component of the processor. The processor and the storage medium can be located in an application-specific integrated circuit (ASIC). In the embodiments of the present application, a computer-readable storage medium may be any tangible medium that contains or stores a program, which may be used by or in conjunction with an instruction execution system, apparatus, or device.

[0213] Based on the same inventive concept, another embodiment of the present application provides an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, which, when executed by the processor, implements the steps in the method of any of the above embodiments of the present application.

[0214] As for the device embodiments or product embodiments, since they are basically similar to the method embodiments, the description is relatively simple, and the relevant parts can be referred to the partial description of the method embodiments.

[0215] The various embodiments in this specification are described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the various embodiments can be referenced to each other.

[0216] Those skilled in the art will appreciate that the embodiments of the present application can be provided as methods, devices, or computer program products. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment in combination with software and hardware. Moreover, the present application can adopt 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.

[0217] The embodiments of the present application are described with reference to the flowcharts and / or block diagrams of the methods, terminal devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes 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 terminal device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing terminal 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.

[0218] Although the preferred embodiments of the present application have been described, those skilled in the art may make additional changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications that fall within the scope of the embodiments of the present application.

[0219] Finally, it should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or terminal device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or terminal device. In the absence of further restrictions, the elements defined by the sentence "comprise a ..." do not exclude the existence of other identical elements in the process, method, article or terminal device including the elements.

Claims

1. A method for determining the conductivity of a crack, characterized in that: The fracture is filled with a proppant, the proppant comprising: propping particles; the method comprising: According to the position information of the propping particles in the fracture at the target moment, obtaining the number of paving layers of the propping agent in the fracture at the target moment; According to the number of paving layers of the proppant in the fracture at the target time, obtaining the permeability of the fracture at the target time; Obtaining a target width of the crack at a target time according to an initial width of the crack and an embedding amount and a deformation amount of the supporting particles at a target time; The overall conductivity of the fracture at a target moment is determined according to the permeability and the target width.

2. The method for determining the conductivity of a crack according to claim 1, characterized in that: Before the step of obtaining the number of paving layers of the proppant in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment, the method further includes: Establishing a simulation model of the proppant in the fracture; According to the simulation model, the position information of the supporting particles in the crack at the target moment is obtained.

3. The method for determining the conductivity of a crack according to claim 2, characterized in that: The step of establishing a simulation model of the proppant in the fracture comprises: Based on the method of coupling computational fluid dynamics and discrete element in Euler-Lagrangian model, a coupling model of the fluid flowing between the proppants is established by using preset interaction parameters as a simulation model of the proppants in the fractures; Wherein, the preset interaction parameter includes at least one of the following: The interaction between the proppants, the interaction between the proppants and the wall of the fracture, and the interaction between the proppants and the fluid.

4. The method for determining the conductivity of a crack according to claim 2, characterized in that: The step of obtaining the position information of the supporting particles in the crack at a target time according to the simulation model comprises: Inputting simulation condition parameters into the simulation model, and performing numerical simulation of proppant migration and placement in the fracture at different simulation times to obtain position information of the proppant particles in the fracture at a target time; The simulation condition parameters include at least one of the following: Proppant parameters, fluid parameters, fracturing pumping process parameters of the fluid, roughness parameters or smoothness parameters of the fracture wall and branching seam parameters.

5. The method for determining the conductivity of a crack according to claim 1, characterized in that: The step of obtaining the number of paving layers of the proppant in the fracture at the target moment according to the position information of the propping particles in the fracture at the target moment comprises: Dividing a plurality of orthogonal grids on the fracture surface of the fracture; Determining the number of supporting particles in each of the orthogonal grids at the target time according to the orthogonal grids and the position information of the supporting particles at the target time; The number of layers of the proppant in the fracture at the target moment is determined according to the number of proppant particles in each of the grids at the target moment, the diameter of the proppant particles and the size of the grid.

6. The method for determining the conductivity of a crack according to claim 5, characterized in that: The step of obtaining the permeability of the fracture at the target moment according to the number of paving layers of the proppant in the fracture at the target moment comprises: Obtaining the porosity of the proppant in the fracture at the target time according to the diameter of the proppant, the position information of the proppant in the fracture at the target time, the number of proppant in each orthogonal grid at the target time, the deformation of the proppant, and the zero-pressure aperture of the fracture; Obtaining the tortuosity of the proppant in the fracture at the target time according to the diameter of the proppant, the position information of the proppant in the fracture at the target time, the number of the proppant in each of the orthogonal grids at the target time, the deformation of the proppant, and the zero-pressure aperture of the fracture; Obtaining the roaring diameter of the proppant in the fracture at the target time according to the diameter of the proppant particles, the position information of the proppant particles in the fracture at the target time, the number of the proppant particles in each of the orthogonal grids at the target time, the embedding amount of the proppant, the deformation amount of the proppant, the zero-pressure opening of the fracture, and the roaring diameter at each of the grids; The permeability of the fracture is obtained according to the porosity of the proppant in the fracture, the tortuosity and the roar diameter.

7. The method for determining the conductivity of a crack according to claim 1, characterized in that: Before the step of obtaining the target width of the crack at a target moment according to the initial width of the crack and the embedding amount and deformation amount of the support particles at a target moment, the method includes: Obtaining the embedding amount of the support particles at a target time according to the diameter of the support particles, the spacing coefficient of the support particles, the extrusion stress between the support particles, the Poisson's ratio of the support particles, the Poisson's ratio of the wall of the crack, the elastic model of the support particles, the elastic model of the wall of the crack, the thickness of the wall of the crack, the viscosity coefficient of the support particles, and the viscosity coefficient of the wall of the crack; The deformation amount of the support particles at the target moment is obtained according to the diameter of the support particles, the spacing coefficient of the support particles, the extrusion stress between the support particles, the elastic model of the support particles, the Poisson's ratio of the support particles, and the viscosity coefficient of the support particles.

8. The method for determining the conductivity of a crack according to claim 1, characterized in that: The step of determining the overall conductivity of the fracture at a target time according to the permeability and the target width comprises: Determine the conductivity of each location in the crack at the target time according to the product of the permeability of the crack at the target time, the target width and a preset correction coefficient; The conductivity of each location in the fracture at the target time is harmonized and averaged to determine the overall conductivity of the fracture at the target time.

9. A device for determining the conductivity of a crack, characterized in that: The fracture is filled with a proppant, the proppant comprising: support particles; the device comprises: A positioning unit, used for obtaining the number of paving layers of the proppant in the fracture at the target moment according to the position information of the proppant particles in the fracture at the target moment; a permeability calculation unit, for obtaining the permeability of the fracture at the target moment according to the number of paving layers of the proppant in the fracture at the target moment; a width calculation unit, configured to obtain a target width of the crack at a target moment according to an initial width of the crack and an embedding amount and a deformation amount of the supporting particles at a target moment; The conductivity calculation unit is used to determine the overall conductivity of the fracture at a target moment according to the permeability and the target width.

10. A readable storage medium, characterized in that: A computer program is stored thereon, and when the program is executed by a processor, the steps in the method according to any one of claims 1 to 8 are implemented.

Images