Fracturing filling crack scale optimization analysis method, device, medium and equipment

Abstract

The invention relates to a fracturing filling crack scale optimization analysis method and device, a medium and equipment. The method comprises the steps that the particle migration critical flow velocity of target loose sandstone reservoir rock is determined based on a flow velocity sensitivity experiment; collecting field logging interpretation data, formation pressure gradient data and production well production data; calculating the single well yield by utilizing a seepage mechanical formula; calculating the flow velocity distribution of reservoir fluid under the fracturing-free condition by utilizing a seepage mechanical formula, and determining the blockage range of the maximum particle migration of the reservoir according to the critical flow velocity of the particle migration; calculating the fluid flow velocity distribution of the reservoir under the same yield and the fracturing filling condition of the production well by utilizing a seepage mechanical formula; flow velocity distribution of fluid under different fracture half lengths is calculated, and the fracture filling fracture scale for relieving / relieving particle migration pollution is determined; and determining the scale of the fracturing filling crack considering the particle migration based on the blockage range of the maximum particle migration of the reservoir and the scale of the fracturing filling crack for retarding / relieving the particle migration pollution.

Description

Technical Field

[0001] The invention relates to a method, a device, a medium and equipment for optimizing and analyzing the scale of fractures in fracturing and filling, and belongs to the technical field of loose sandstone fracturing and filling. Background Art

[0002] The characteristics of loose sandstone are low strength and poor cementation. In addition, the mineral composition of loose sandstone has a high content of mineral components such as kaolinite and montmorillonite that are highly sensitive to speed and water. During the development process, it is easy to cause particle migration and particle swelling during water blockage pollution. In response to the above-mentioned blockage pollution problem, the on-site method generally adopts small-scale acidification to improve it. However, this method can only remove the 1-2m contamination zone near the wellbore, which is very limited. Therefore, for the blockage pollution problem with poor acidification effect, domestic researchers have proposed a deep penetration fracturing technology. After numerical simulation and field application, it has been confirmed that when the fracture passes through the contamination zone, the production capacity can be restored to achieve the effect of removing pollution. However, fracturing can not only increase production and remove blockage, but also change the flow state of the near-wellbore zone, transform the radial flow without fractures into linear flow, and change the fluid flow rate in the near-wellbore area. The reservoir rock of the shallow oilfield in Bohai Sea is loose sandstone. During the production process, it is very easy for particles to migrate, resulting in reservoir blockage and production decline. The existing practice on site is to take a series of measures to remove blockage and increase production after blockage occurs and production declines. If the appropriate fracture size is determined during the drilling and completion stage, and the production well is completed by fracturing and filling, the occurrence of particle migration blockage pollution can be greatly reduced. Therefore, it is necessary to determine a method for determining the fracture size that can remove particle migration blockage by fracturing and filling with consideration of particle migration. Summary of the invention

[0003] In view of the above technical problems, the present invention provides a method, device, medium and equipment for optimizing the analysis of the scale of fractures in fracturing and filling. The method optimizes the scale of fractures in loose sandstone reservoirs with different permeabilities to reduce or eliminate the blockage of particle migration, and provides a theoretical guidance method for on-site fracturing and filling design.

[0004] To achieve the above object, the present invention adopts the following technical solutions:

[0005] A method for optimizing and analyzing the scale of fracture filling by hydraulic fracturing comprises the following steps:

[0006] Determine the critical velocity of particle migration in the target loose sandstone reservoir rock based on velocity sensitivity experiments;

[0007] Collect on-site logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production well production data;

[0008] Based on the critical velocity of particle migration, well logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells, the single well production is calculated using the seepage mechanics formula;

[0009] The flow velocity distribution of reservoir fluid under non-fracturing conditions is calculated using the seepage mechanics formula, and the maximum particle migration blockage range of the reservoir is determined based on the critical particle migration velocity;

[0010] The fluid velocity distribution of the reservoir under the conditions of fracturing and filling of production wells at the same production rate is calculated using the seepage mechanics formula;

[0011] Calculate the velocity distribution of fluids at different fracture half-lengths, draw contour maps of velocity, and determine the scale of fracture filling that can slow down / eliminate particle migration contamination;

[0012] Based on the blockage range of the maximum particle migration in the reservoir and the scale of the fracturing and filling fractures that can slow down / eliminate the pollution caused by particle migration, the scale of the fracturing and filling fractures that takes particle migration into consideration is determined.

[0013] The fracturing filling fracture scale optimization analysis method is preferably based on the flow rate sensitivity experiment to determine the critical flow rate of particle migration of the target loose sandstone reservoir rock, and the calculation formula of the rock sample permeability change rate corresponding to different flow rates is as follows:

[0014]

[0015] Where D is the change rate of rock permeability corresponding to different flow rates; K n is the permeability of rock samples at different flow rates; K1 is the initial permeability.

[0016] The method for optimizing the scale of fracture filling is preferably to calculate the single well production using a seepage mechanics formula, and the specific calculation method is as follows:

[0017]

[0018] Where: q(t) is the production of the production well; K is the reservoir permeability; h is the reservoir thickness; p e is the initial formation pressure; p w is the bottom hole flowing pressure; r e is the supply radius; r w The diameter of the production well.

[0019] The method for optimizing the analysis of the scale of fracture filling is preferably to use a seepage mechanics formula to calculate the flow velocity distribution of the reservoir fluid under non-fracture conditions, and the calculation formula is as follows:

[0020]

[0021] Where: V(z) is the flow velocity of the reservoir fluid; z is the particle position in the rectangular coordinate system; z c is the wellbore position in the rectangular coordinate system.

[0022] The method for optimizing the analysis of the scale of fracture filling is preferably to determine the maximum particle migration blockage range of the reservoir according to the critical velocity of particle migration, as follows:

[0023] Draw a contour map of the critical velocity and record the wellbore z c The distance r to the critical velocity boundary is the blockage range of the maximum particle migration in the reservoir.

[0024] The method for optimizing the analysis of fracture scale of fracturing filling preferably uses a seepage mechanics formula to calculate the fluid velocity distribution of the reservoir under the conditions of fracturing filling of the production well at the same production rate. The specific calculation formula is as follows:

[0025]

[0026] Where: β is the angle between the crack and the horizontal coordinate axis; L is the crack length; i is an imaginary number.

[0027] The method for optimizing and analyzing the scale of fracture filling by fracturing preferably calculates the flow velocity distribution of the fluid at different fracture half-lengths, and draws contour maps of the flow velocity to determine the scale of fracture filling by fracturing that slows down / eliminates particle migration pollution, as follows:

[0028] The velocity distribution of fluids at different fracture half-lengths is calculated by programming, and contour maps of the velocity are drawn. The fracture half-length when the critical velocity boundary of the reservoir / fracture half-length is less than 0.1 is recorded as the optimized fracturing filling fracture for removing particle migration pollution. When the calculated fracture length reaches the maximum fracture length that can be achieved in on-site construction and the critical velocity boundary / fracture half-length is still not less than 0.1, the optimization calculation is stopped. At this time, the fracture half-length is the optimal fracture half-length for slowing down particle migration blockage.

[0029] A second aspect of the present invention provides a device for optimizing and analyzing the scale of fracture filling, comprising:

[0030] The first processing unit is used to determine the critical velocity of particle migration of the target loose sandstone reservoir rock based on the velocity sensitivity experiment;

[0031] The second processing unit is used to collect on-site logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells;

[0032] The third processing unit is used to calculate the single well production using the seepage mechanics formula based on the critical velocity of particle migration, well logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells;

[0033] The fourth processing unit is used to calculate the flow velocity distribution of the reservoir fluid under the condition of no fracturing by using the seepage mechanics formula, and determine the blockage range of the maximum particle migration of the reservoir according to the critical flow velocity of particle migration;

[0034] The fifth processing unit is used to calculate the fluid velocity distribution of the reservoir under the condition of fracturing and filling of the production well under the same production rate by using the seepage mechanics formula;

[0035] The sixth processing unit is used to calculate the flow velocity distribution of the fluid at different fracture half-lengths, draw contour maps of the flow velocity, and determine the scale of the fracture filling for slowing down / relieving particle migration pollution;

[0036] The seventh processing unit is used to determine the scale of the fracturing filling cracks considering particle migration based on the maximum particle migration blockage range of the reservoir and the scale of the fracturing filling cracks that slow down / eliminate particle migration pollution.

[0037] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the steps of any one of the above-mentioned methods for optimizing the analysis of the scale of fracture filling and fracturing are implemented.

[0038] A fourth aspect of 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 when the processor executes the computer program, the steps of any one of the above-mentioned methods for optimizing the analysis of the scale of fracture filling and fracturing are implemented.

[0039] The present invention adopts the above technical solution, which has the following advantages:

[0040] The method of the present invention optimizes the scale of cracks in loose sandstone reservoirs with different permeabilities to reduce or relieve particle migration blockage, and provides a theoretical guidance method for on-site fracturing and filling design. The present invention can better optimize the scale of fracturing and filling cracks when considering relieving or slowing down particle migration blockage pollution. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Figure 1 It is a flow chart of a method for optimizing the analysis of fracture scale of hydraulic fracturing and filling considering particle migration provided by an embodiment of the present invention;

[0042] Figure 2 Schematic diagram of the crack not breaking through the particle migration blocking boundary and the crack breaking through the particle migration boundary provided by the embodiment of the present invention, wherein Figure a shows no breaking through and Figure b shows breaking through;

[0043] Figure 3 This is a flow half-length / fracture half-length curve diagram of the critical flow velocity of particle migration under different permeabilities and different half-lengths provided by this embodiment of the present invention. DETAILED DESCRIPTION

[0044] In order to make the purpose, technical solution and advantages of the present invention clearer, the technical solution of the present invention is described clearly and completely below. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in the art without creative work are within the scope of protection of the present invention.

[0045] The reservoir rock of the shallow oilfield in Bohai Sea is loose sandstone, which is very easy to cause particle migration during the production process, resulting in reservoir blockage and production decline. The existing practice on site is to take a series of measures to remove blockage and increase production after blockage occurs and production declines. If the appropriate fracture size is determined during the drilling and completion stage, and the production well is completed by fracturing and filling, the occurrence of particle migration blockage pollution can be greatly reduced. Therefore, it is necessary to determine a method for determining the fracture size that takes particle migration into consideration to remove particle migration blockage by fracturing and filling.

[0046] In view of the above technical problems, the present invention provides a method, device, medium and equipment for optimizing the analysis of the scale of fractures in fracturing and filling. The method optimizes the scale of fractures in loose sandstone reservoirs with different permeabilities to reduce or eliminate the blockage of particle migration, and provides a theoretical guidance method for on-site fracturing and filling design.

[0047] like Figure 1 As shown, the method for optimizing the analysis of the scale of fracture filling provided by the present invention is as follows:

[0048] Step 1: Velocity sensitivity experiment to determine the critical velocity of reservoir rock

[0049] The experimental method of flow velocity sensitivity experiment refers to "Reservoir Sensitivity Flow Experiment Evaluation Method" SY / T5358-2010.

[0050] Through the standard experimental core (diameter 2.5cm, length 5cm) according to the liquid flow rate 0.1cm 3 / min, 0.25cm 3 / min, 0.50cm 3 / min, 0.1cm 3 / min, 0.75cm 3 / min, 1.0cm 3 / min, 1.5cm 3 / min, 2.0cm 3 / min, 3.0cm 3 / min, 4.0cm 3 / min, 5.0cm 3 / min and 6.0cm 3 / min liquid flow rate, and the rock permeability is measured in turn. During the test, the initial confining pressure of the rock is 2.0MPa, and the confining pressure value is always kept greater than the rock inlet pressure by 1.5MPa~2.0MPa. According to formula (1), the permeability change rate of the rock at different flow rates is calculated, and the flow rate corresponding to the previous point when the rock permeability change rate D is greater than 20% is recorded as the critical flow rate:

[0051]

[0052] Where: D: the change rate of rock sample permeability at different flow rates; K n : Rock sample permeability (corresponding to different flow rates in the experiment), unit 10 -3 μm 2 ; K1: initial permeability (the permeability of the rock sample corresponding to the minimum flow rate in the experiment), unit 10 -3 μm 2 .

[0053] Step 2: Basic information collection

[0054] Collect the field logging interpretation data (as shown in Table 1), formation pressure gradient data, reservoir fluid viscosity and production data of production wells (wellbore radius, bottom hole flowing pressure).

[0055] Table 1 Schematic diagram of field logging interpretation results

[0056]

[0057]

[0058] Step 3: Calculate single well production using the seepage mechanics formula

[0059] The basic data collected in step 1 and step 2 are used to calculate the production q(t) of the production well under the basic data conditions according to formula (2):

[0060]

[0061] Where: q(t): production well output, unit: m 3 / s; K: reservoir permeability, unit: μm 2 (In this formula, the unit is m 2 ); h: reservoir thickness, unit: m; p e : initial formation pressure (supply boundary pressure), unit: Pa; p w : Bottom hole pressure, unit: Pa; r e : Supply radius, unit: m; r w : Diameter of production well, unit: m.

[0062] Formula 3 is used to calculate the reservoir fluid velocity distribution under non-fracturing conditions:

[0063]

[0064] Where: V(z): reservoir fluid velocity, unit: m / s; z: particle position in rectangular coordinate system; z c : Wellbore position in rectangular coordinate system.

[0065] Draw a contour map of the critical velocity and record the wellbore z c The distance r to the critical velocity boundary is the blockage range of the maximum particle migration in the reservoir (particle migration boundary).

[0066] Formula (4) is used to calculate the fluid velocity distribution of the reservoir under the conditions of fracturing and filling of production wells at the same production rate:

[0067]

[0068] Where: β: the angle between the crack and the horizontal coordinate axis, unit: rad; L: crack length, unit: m; i: imaginary number.

[0069] The flow velocity distribution of the fluid under different fracture half-lengths (0.5L) is calculated by programming (the fracture half-length starts from 10m and the calculation interval is 5m) and the contour map of the flow velocity is drawn. The fracture half-length ( Figure 2 The figure shows the fracture filling crack optimized by the method to relieve particle migration contamination. When the calculated fracture length reaches the maximum fracture length that can be achieved in field construction and the critical velocity boundary / fracture half-length is still not less than 0.1, the optimization calculation is stopped. At this time, the fracture half-length is the optimal fracture half-length for slowing down particle migration blockage.

[0070] The technical solution of the present invention is described in detail below in conjunction with specific application examples.

[0071] Table 2 shows the logging interpretation results of a loose sandstone reservoir in Bohai:

[0072] Table 2 Well logging interpretation results of an oil field in Bohai Sea

[0073]

[0074]

[0075] The production well and reservoir conditions are shown in Table 3:

[0076] Table 3 Key parameters used in calculation

[0077] <![CDATA[Critical flow velocity (v r : 10 -5 m / s)]]> 8.19 <![CDATA[Initial formation pressure (p e : 10 6 Pa)]]> 13 <![CDATA[Bottom-hole flowing pressure (p w : 10 6 Pa)]]> 8 Fluid viscosity (μ: Pa·s) 0.001 <![CDATA[Hole diameter (r w : m)]]> 0.1 <![CDATA[Supply radius (r e : m)]]> 250

[0078] According to formulas 2 and 3, the maximum range of particle migration blockage in the reservoir under the conditions in Table 3 can be calculated. The results are shown in Table 4:

[0079] Table 4 Well logging interpretation results of an oil field in Bohai Sea

[0080]

[0081]

[0082] Taking reservoirs 7, 8, 12, 13 and 14 as examples, the optimal fracture scale for removing / slowing down particle migration is optimized and calculated using formula 2 and formula 4. The calculation results are plotted with "critical flow half length of particle migration / fracture half length" as the vertical axis and "fracture half length" as the horizontal axis. It is stipulated that when the particle migration boundary > fracture half length, the critical flow velocity half length of particle migration / fracture half length is 1. The calculation results are shown in Figure 3 shown.

[0083] The calculation optimization results show that the optimal half-length of reservoir 7 is 15m, the optimal half-length of reservoirs 8 and 12 is 30m, the optimal half-length of reservoir 13 is 20m, and the optimal half-length of reservoir 14 is 50m.

[0084] Formulas 2 to 3 involved in the optimization calculation in this method require that all parameter units involved in the calculation are in SI units.

[0085] A second aspect of the present invention provides a device for optimizing and analyzing the scale of fracture filling, comprising:

[0086] The first processing unit is used to determine the critical velocity of particle migration of the target loose sandstone reservoir rock based on the velocity sensitivity experiment;

[0087] The second processing unit is used to collect on-site logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells;

[0088] The third processing unit is used to calculate the single well production using the seepage mechanics formula based on the critical velocity of particle migration, well logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells;

[0089] The fourth processing unit is used to calculate the flow velocity distribution of the reservoir fluid under the condition of no fracturing by using the seepage mechanics formula, and determine the blockage range of the maximum particle migration of the reservoir according to the critical flow velocity of particle migration;

[0090] The fifth processing unit is used to calculate the fluid velocity distribution of the reservoir under the condition of fracturing and filling of the production well under the same production rate by using the seepage mechanics formula;

[0091] The sixth processing unit is used to calculate the flow velocity distribution of the fluid at different fracture half-lengths, draw contour maps of the flow velocity, and determine the scale of the fracture filling for slowing down / relieving particle migration pollution;

[0092] The seventh processing unit is used to determine the scale of the fracturing filling cracks considering particle migration based on the maximum particle migration blockage range of the reservoir and the scale of the fracturing filling cracks that slow down / eliminate particle migration pollution.

[0093] A third aspect of the present invention provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the steps of any one of the above-mentioned methods for optimizing the analysis of the scale of fracture filling and fracturing are implemented.

[0094] A fourth aspect of 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 when the processor executes the computer program, the steps of any one of the above-mentioned methods for optimizing the analysis of the scale of fracture filling and fracturing are implemented.

[0095] The present invention is described in terms of flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to specific embodiments. 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 generate 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.

[0096] 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 1 A function specified in one or more boxes.

[0097] 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 1The steps for the functions specified in one or more boxes.

[0098] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit it. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they can still modify the technical solutions described in the aforementioned embodiments, or make equivalent replacements for some of the technical features therein. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention.

Claims

1. A method for optimizing the scale of fracture filling, characterized in that: The steps include: Determine the critical velocity of particle migration in the target loose sandstone reservoir rock based on velocity sensitivity experiments; Collect on-site logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production well production data; Based on the critical velocity of particle migration, well logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells, the single well production is calculated using the seepage mechanics formula; The flow velocity distribution of reservoir fluid under non-fracturing conditions is calculated using the seepage mechanics formula, and the maximum particle migration blockage range of the reservoir is determined based on the critical particle migration velocity; The fluid velocity distribution of the reservoir under the conditions of fracturing and filling of production wells at the same production rate is calculated using the seepage mechanics formula; Calculate the velocity distribution of fluids at different fracture half-lengths, draw contour maps of velocity, and determine the scale of fracture filling that can slow down / eliminate particle migration contamination; Based on the blockage range of the maximum particle migration in the reservoir and the scale of the fracturing and filling fractures that can slow down / eliminate the pollution caused by particle migration, the scale of the fracturing and filling fractures that takes particle migration into consideration is determined.

2. The method for optimizing the scale of fracture filling according to claim 1, characterized in that: Based on the velocity sensitivity experiment, the critical velocity of particle migration in the target loose sandstone reservoir rock is determined. The calculation formula of the critical velocity of particle migration is as follows: Where D is the change rate of rock permeability corresponding to different flow rates; K n is the rock sample permeability; K1 is the initial permeability.

3. The method for optimizing the scale of fracture filling according to claim 1, characterized in that: The single well production is calculated using the seepage mechanics formula. The specific calculation method is as follows: Where: q(t) is the production of the production well; K is the reservoir permeability; h is the reservoir thickness; p e is the initial formation pressure; p w is the bottom hole flowing pressure; r e is the supply radius; r w The diameter of the production well.

4. The method for optimizing the scale of fracture filling according to claim 1, characterized in that: The flow velocity distribution of reservoir fluid under non-fracturing conditions is calculated using the seepage mechanics formula. The calculation formula is as follows: Where: V(z) is the flow velocity of the reservoir fluid; z is the particle position in the rectangular coordinate system; z c is the wellbore position in the rectangular coordinate system.

5. The method for optimizing the scale of fracture filling according to claim 4, characterized in that: The maximum particle migration blockage range of the reservoir is determined based on the critical velocity of particle migration, as follows: Draw a contour map of the critical velocity and record the wellbore z c The distance r to the critical velocity boundary is the blockage range of the maximum particle migration in the reservoir.

6. The method for optimizing the scale of fracture filling according to claim 5, characterized in that: The fluid velocity distribution of the reservoir under the conditions of fracturing and filling of production wells at the same production rate is calculated using the seepage mechanics formula. The specific calculation formula is as follows: Where: β is the angle between the crack and the horizontal coordinate axis; L is the crack length; i is an imaginary number.

7. The method for optimizing the scale of fracture filling according to claim 6, characterized in that: Calculate the velocity distribution of fluids at different fracture half-lengths, draw contour maps of the velocity, and determine the scale of the fracture filling that can slow down / eliminate particle migration pollution, as follows: The velocity distribution of the fluid under different fracture half-lengths is calculated, and the contour map of the velocity is drawn. The fracture half-length when the critical velocity boundary of the reservoir / fracture half-length is less than 0.1 is recorded as the optimized fracturing filling fracture for removing particle migration pollution. When the calculated fracture length reaches the maximum fracture length that can be achieved in on-site construction, the optimization calculation is stopped when the critical velocity boundary / fracture half-length is still not less than 0.

1. At this time, the fracture half-length is the optimal fracture half-length for slowing down particle migration blockage.

8. A device for optimizing and analyzing the scale of fracture filling, characterized in that: include: The first processing unit is used to determine the critical velocity of particle migration of the target loose sandstone reservoir rock based on the velocity sensitivity experiment; The second processing unit is used to collect on-site logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells; The third processing unit is used to calculate the single well production using the seepage mechanics formula based on the critical velocity of particle migration, well logging interpretation data, formation pressure gradient data, reservoir fluid viscosity and production data of production wells; The fourth processing unit is used to calculate the flow velocity distribution of the reservoir fluid under the condition of no fracturing by using the seepage mechanics formula, and determine the blockage range of the maximum particle migration of the reservoir according to the critical flow velocity of particle migration; The fifth processing unit is used to calculate the fluid velocity distribution of the reservoir under the condition of fracturing and filling of the production well under the same production rate by using the seepage mechanics formula; The sixth processing unit is used to calculate the flow velocity distribution of the fluid at different fracture half-lengths, draw contour maps of the flow velocity, and determine the scale of the fracture filling for slowing down / relieving particle migration pollution; The seventh processing unit is used to determine the scale of the fracturing filling cracks considering particle migration based on the maximum particle migration blockage range of the reservoir and the scale of the fracturing filling cracks that slow down / eliminate particle migration pollution.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for optimizing the analysis of the scale of fracture filling and fracturing are implemented.

10. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that: When the processor executes the computer program, the steps of the method for optimizing the analysis of the scale of fracture filling and fracturing are implemented as described in any one of claims 1 to 7.

Images