A method for determining development plans based on the seepage resistance of two-phase fluids in shale nanopores

Abstract

The embodiment of this specification provides a method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores. The two-phase fluid includes water-phase fluid and oil-phase fluid; the method includes: collecting fluid parameters corresponding to at least one group of experimental samples; the experimental samples have different test factor attribute values ​​corresponding to the test factors; the fluid parameters include contact angle parameters between the water-phase fluid and the oil-phase fluid; screening resistance control factors from the test factors based on the fluid parameters; calculating the two-phase fluid seepage resistance of the target formation based on the resistance control factors; and determining a development plan using the two-phase fluid seepage resistance. The above method provides guidance for the determination of the development plan, improves the effectiveness of the plan, and ensures the normal progress of production and development.

Description

Technical Field

[0001] The embodiments of this specification relate to the technical field of formation exploration and development, and in particular to a method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores. Background Art

[0002] Before exploiting oil resources in a reservoir, it is often necessary to first obtain the geological parameters of the reservoir and, based on these parameters, simulate the parameter changes during the development process to determine the appropriate exploitation plan. During the exploitation of oil resources, both the starting pressure at which crude oil begins to flow from a static state and the flow resistance encountered by crude oil during its flow through the formation have a certain impact on the exploitation process. This is especially true when exploiting shale oil, especially in the micro- and nano-pores of shale. Because the capillary pressure of shale oil reservoirs is much greater than that of conventional oil reservoirs, and due to the influence of capillary pressure and the Jamin effect, the starting pressure of the oil-water two-phase system is much greater than that of the single-phase system, the calculation of the starting pressure and flow resistance of the two-phase fluid is of great significance and is one of the decisive factors in the overall exploitation process.

[0003] However, the current mining process can only adapt to the effects of the seepage resistance of two-phase fluids through experimental simulations and adjustments during production. There is no solution that can qualitatively calculate the seepage resistance of two-phase fluids during development, which interferes with actual production results and is not conducive to the formulation of mining plans. Therefore, there is an urgent need for a method that can effectively calculate the seepage resistance of two-phase fluids in the formation and then formulate development plans. Summary of the Invention

[0004] The purpose of the embodiments of this specification is to provide a method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores, so as to solve the problem of how to quantitatively analyze and calculate the resistance of two-phase fluid.

[0005] In order to solve the above technical problems, the embodiments of this specification propose a method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores, wherein the two-phase fluid includes an aqueous fluid and an oil-phase fluid; the method includes: collecting fluid parameters corresponding to at least one group of experimental samples; the experimental samples have different test factor attribute values ​​corresponding to the test factors; the fluid parameters include contact angle parameters between the aqueous fluid and the oil-phase fluid; screening resistance control factors from the test factors according to the fluid parameters; calculating the two-phase fluid seepage resistance of the target formation based on the resistance control factors; and determining a development plan using the two-phase fluid seepage resistance.

[0006] The embodiments of this specification also propose a development plan determination device based on the seepage resistance of two-phase fluid in shale nanopores, wherein the two-phase fluid includes an aqueous fluid and an oil-phase fluid; the device includes: a parameter acquisition module for acquiring fluid parameters corresponding to at least one group of experimental samples; the experimental samples have different test factor attribute values ​​corresponding to the test factors; the fluid parameters include contact angle parameters between the aqueous fluid and the oil-phase fluid; a resistance control factor screening module for screening resistance control factors from the test factors according to the fluid parameters; a two-phase fluid seepage resistance calculation module for calculating the two-phase fluid seepage resistance of the target formation based on the resistance control factors; and a development plan determination module for determining a development plan using the two-phase fluid seepage resistance.

[0007] It can be seen from the technical solutions provided in the above embodiments of this specification that the embodiments of this specification first screen the resistance control factors from the test factors by analyzing the physical parameters and fluid parameters of the experimental samples, especially the contact angle parameters among the fluid parameters, and then calculate the two-phase fluid seepage resistance in the target formation based on the resistance control factors obtained by screening, and finally use the two-phase fluid seepage resistance to determine the final development plan. Through the above method, when determining the development plan based on the resistance encountered by the oil phase liquid and the water phase liquid during flow in the formation, it is possible to predetermine the control factors that have a strong influence on the resistance, and then complete the calculation of the seepage resistance in the target formation by combining the control factors and the contact angle between the fluids, thereby providing guidance for the determination of the development plan, improving the effectiveness of the plan, and ensuring the normal progress of production and development. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0009] Figure 1 This is a flow chart of a method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores according to an embodiment of this specification;

[0010] Figure 2 This is a schematic diagram of the shape of a liquid droplet when it is about to move under the influence of static wetting hysteresis in an embodiment of this specification;

[0011] Figure 3 This is a schematic diagram of the change of dynamic contact angle with speed under a dynamic wetting condition in an embodiment of this specification;

[0012] Figure 4AThis is a schematic diagram of a curve showing the change in the starting pressure of an alkane-water system with pore size under different nanopore sizes in an embodiment of this specification;

[0013] Figure 4B This is a schematic diagram of a curve showing the change in the starting pressure of an alkane-water system with pore size under different nanopore sizes in an embodiment of this specification;

[0014] Figure 4C This is a schematic diagram of a curve showing the change in the starting pressure of an alkane-water system with pore size under different nanopore sizes in an embodiment of this specification;

[0015] Figure 4D This is a schematic diagram of a curve showing the change in the starting pressure of an alkane-water system with pore size under different nanopore sizes in an embodiment of this specification;

[0016] Figure 4E This is a schematic diagram of a curve showing the change in the starting pressure of an alkane-water system with pore size under different nanopore sizes in an embodiment of this specification;

[0017] Figure 4F This is a schematic diagram of a curve showing the change in the starting pressure of an alkane-water system with pore size under different nanopore sizes in an embodiment of this specification;

[0018] Figure 5A This is a schematic diagram of a curve showing the change of flow resistance with flow velocity for different fluid systems according to the embodiment of this specification;

[0019] Figure 5B This is a schematic diagram of a curve showing the change of flow resistance with flow velocity for different fluid systems according to the embodiment of this specification;

[0020] Figure 5C This is a schematic diagram of a curve showing the change of flow resistance with flow velocity for different fluid systems according to the embodiment of this specification;

[0021] Figure 5D This is a schematic diagram of a curve showing the change of flow resistance with flow velocity for different fluid systems according to the embodiment of this specification;

[0022] Figure 5E This is a schematic diagram of a curve showing the change of flow resistance with flow velocity for different fluid systems according to the embodiment of this specification;

[0023] Figure 6A This is a schematic diagram of flow resistance caused by dynamic wetting hysteresis under pores of different scales according to an embodiment of this specification;

[0024] Figure 6B This is a schematic diagram of flow resistance caused by dynamic wetting hysteresis under pores of different scales according to an embodiment of this specification;

[0025] Figure 7 This is a module diagram of a device for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores according to an embodiment of this specification. DETAILED DESCRIPTION

[0026] The following will be combined with the drawings in the embodiments of this specification to clearly and completely describe the technical solutions in the embodiments of this specification. Obviously, the embodiments described are only part of the embodiments of this specification, not all of the embodiments. Based on the embodiments in this specification, all other embodiments obtained by ordinary technicians in this field without making creative efforts should fall within the scope of protection of this specification.

[0027] In order to solve the above technical problems, the embodiments of this specification propose a method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores. The execution subjects of the method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores include but are not limited to servers, industrial computers, PCs, etc. Figure 1 As shown, the method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores may specifically include the following steps.

[0028] S110: Collecting fluid parameters corresponding to at least one group of experimental samples; the experimental samples have different test factor attribute values ​​corresponding to the test factors; the fluid parameters include contact angle parameters between the water phase fluid and the oil phase fluid.

[0029] Two-phase fluids can include aqueous and oil-phase fluids. In addition to crude oil, aqueous fluids may also be present in formations, either naturally present or artificially injected during extraction. Because aqueous and oil-phase fluids are immiscible and have certain differences in physical properties, they interact with each other during fluid flow. Therefore, when calculating the seepage resistance of two-phase fluids, the impact of these various factors must be comprehensively considered.

[0030] The experimental samples can be simulated core samples prepared in order to determine the factors that are more obvious for the seepage resistance of two-phase fluids. Specifically, at least one group of experimental samples can be prepared, and each group of experimental samples corresponds to a different test factor. For example, group A experimental samples correspond to temperature, group B experimental samples correspond to pore size, and group C experimental samples correspond to permeability. Accordingly, each experimental sample in each group of experimental samples is controlled to have a different test factor attribute value corresponding to the test factor, and the values ​​of other factors except the test factor are maintained the same, thereby playing a control role. For example, the temperature of each experimental sample in group A experimental samples is controlled to be 0℃, 10℃, 20℃, 30℃, etc. The specific selection of test factors and the size of the test factor attribute values ​​can be set separately based on actual test requirements, and are not limited to the above examples, and will not be repeated here.

[0031] After completing the setup for each set of experimental samples, the fluid parameters of each set of experimental samples can be collected. Fluid parameters can be used to reflect the properties of the two-phase fluid in the experimental samples. Specifically, they can be parameters that guide the calculation of the seepage resistance of the two-phase fluid, such as the fluid's surface tension parameter, molecular average transition frequency, and molecular average transition distance.

[0032] Specifically, surface tension parameters can be obtained using the pendant drop method, while the average molecular transition frequency and average molecular transition distance can also be obtained using corresponding test methods. These parameters can also be directly obtained by calling existing test data. In practical applications, there is no restriction on the acquisition method.

[0033] In some embodiments, the fluid parameters may include contact angle parameters between the aqueous and oil phases. In a two-phase flow system, due to the presence of wetting hysteresis, two capillary pressures exist at the two interfaces between the oil and water phases. The difference between the two capillary forces acts as an obstacle to fluid flow. For systems with static wetting hysteresis, the difference between the two capillary forces is represented by the starting pressure, i.e., the maximum resistance that the fluid must overcome to begin moving from a static state. For systems with dynamic wetting hysteresis, the difference between the two capillary pressures is represented by the flow resistance, i.e., the resistance that the fluid must overcome during movement.

[0034] The aforementioned starting pressure and flow resistance have a significant impact on the movement of two-phase fluid in the formation, and therefore need to be calculated separately. Therefore, the contact angle parameter among the fluid parameters can be used for specific calculations.

[0035] The contact angle parameters include advancing contact angle and receding contact angle; the advancing contact angle includes the contact angle formed when water displaces oil and water occupies part of the original space of oil, and the receding contact angle includes the contact angle formed when oil displaces water and oil occupies part of the original space of water. Figure 2The figure shows the contact angle of the two-phase fluid. The oil droplet moves in the left direction. The contact angle between the oil phase fluid and the water phase fluid is θ. A , the receding contact angle is θ R , where θ A >θ o >θ R ,θ o is the contact angle at rest. The difference in capillary pressure between the two ends of the oil droplet on the oil-water meniscus creates additional resistance in the opposite direction of the oil's movement. To move the oil droplet, the external force must be greater than this additional resistance, i.e., the starting pressure. When the oil droplet begins to move in the capillary, both the advancing and receding contact angles change with the speed of the oil-water meniscus, resulting in dynamic wetting hysteresis.

[0036] The contact angle reflects the resistance currently experienced by a two-phase fluid. A larger contact angle indicates greater resistance. Therefore, when analyzing the seepage resistance of a two-phase fluid, the contact angle can be used to determine the relative resistance.

[0037] It should be noted that when the executor of this method is a computing device, the direct acquisition of physical parameters and fluid parameters may be performed by other detection devices, and the executor may simply accept the parameter values ​​measured by other detection devices for subsequent calculations.

[0038] In order to obtain better test results, when setting experimental samples for starting pressure and flow resistance respectively, the specific types of experimental samples can be adjusted according to corresponding needs.

[0039] Specifically, when testing for starting pressure, the experimental sample includes at least one of an octane-water system, a dodecane-water system, and a hexadecane-water system on a fluoropolymer coating on a dry mica wafer. When testing for flow resistance, the experimental sample includes at least one of a hexane-water system and a squalane-water system on a single crystal silicon wafer, and a C8-water system, a C12-water system, and a C13-water system on apatite.

[0040] In practical applications, the type of the experimental sample can be adjusted according to specific needs and is not limited to the above examples.

[0041] S120: Filtering resistance control factors from test factors according to the fluid parameters.

[0042] After obtaining the fluid parameters, the resistance control factors can be screened from the test factors based on the contact angle of the fluid parameters. Since the contact angle can be used to reflect the resistance of the two-phase fluid, the resistance control factors can be screened by analyzing the contact angles of different groups of experimental samples.

[0043] Specifically, a curve graph of the contact angle parameters corresponding to the test factors can be drawn. By analyzing the rate of change in the curve graph, if the rate of change is greater than the change threshold, it means that the test factor corresponding to the current experimental sample has a greater effect on the contact angle, and the test factor can be used as a resistance control factor.

[0044] In some embodiments, the resistance control factor determined by the above method may include at least one of temperature and pore radius. In practical applications, the resistance control factors screened out may not be limited to the above examples and will not be described in detail here.

[0045] S130: Calculating the two-phase fluid seepage resistance of the target formation based on the resistance control factor.

[0046] After determining the resistance control factors, a corresponding calculation model can be constructed based on the resistance control factors to specifically calculate the two-phase fluid seepage resistance. Since the two-phase fluid seepage resistance includes the starting pressure and the flow resistance, corresponding formulas can be constructed to calculate the above two resistances respectively.

[0047] Specifically, for the starting pressure of the two-phase fluid in a static state, combined with the capillary pressure equation, the starting pressure equation corresponding to the static wetting hysteresis can be constructed as follows: Where, P f1 is the starting pressure, γ 12 is the interfacial tension between the wetting and non-wetting phases, r is the pore radius, θ R is the receding contact angle, θ A is the advancing contact angle. By substituting the corresponding parameter values ​​into the above formula, the calculation of the starting pressure in the static state can be completed.

[0048] Regarding the flow resistance of two-phase fluid in motion, since there is a linear relationship between the moving speed of the oil-water curved interface and the dynamic contact angle, when the moving speed increases, the advancing contact angle increases and the receding contact angle decreases, that is, the difference in dynamic angles increases with the increase in moving speed, such as Figure 3 As shown in the figure, U is the moving speed of the three contact lines between the rock surface, water phase and oil phase, θ D is the dynamic contact angle of the three-phase contact line.

[0049] Based on the changes between the moving speed and contact angle of the three-phase contact line, the mathematical relationship between the dynamic contact angle and the moving speed of the three-phase contact line can be constructed as follows: Where, U is the moving speed of the three-phase contact line, m / s; K 0 is the average transition frequency of the molecule, 1 / s; λ is the average transition distance of the molecule, m; k B is the Boltzmann constant, m 2 kg / s 2 K; T is temperature, K; θ o is the equilibrium or static contact angle, deg; θ D is the dynamic contact angle, deg; sinhx is the hyperbolic sine function, which can be defined with the help of the exponential function, that is,

[0050] U>0 corresponds to the advancing contact angle θ A , U<0 corresponds to the receding contact angle θ R , and the corresponding formulas are and After subtracting the above two equations, we get Combining this formula with the above starting pressure equation, the formula for the change of flow resistance caused by dynamic wetting hysteresis with moving speed is obtained as follows: Where, P f2 is the flow resistance, k B is the Boltzmann constant, T is the temperature, r is the pore radius, λ is the average molecular transition distance, U is the moving speed of the three-phase contact line, K 0 is the average molecular transition frequency.

[0051] The data in the above formula can be obtained from published experimental data, resulting in parameters including the average molecular transition frequency and average molecular transition distance for the hexane-water and squalane-water systems on single-crystal silicon wafers, and the C8-water, C12-water, and C13-water systems on apatite. The specific data can be seen in Table 1 below.

[0052] Table 1

[0053] solid surface liquid <![CDATA[θ0 / deg]]> <![CDATA[K 0 / 10 5 Hz]]> λ / nm T / K Monocrystalline silicon wafers <![CDATA[C6H 14 ]]> 32.7 2.68 1.91 293.15 Monocrystalline silicon wafers <![CDATA[C 30 H 62 ]]> 28.0 4.15 1.11 293.15 apatite <![CDATA[C8]]> 50.0 2.93 0.245 293.15 apatite <![CDATA[C 12 ]]> 37.0 5.823 0.278 293.15 apatite <![CDATA[C 13 ]]> 21.0 9.604 0.237 293.15

[0054] Based on the formula derived above, combined with the geological parameters collected from the target formation, the two-phase fluid seepage resistance in the target formation can be calculated, which can then be used to formulate development plans in subsequent steps.

[0055] S140: Determine a development plan using the two-phase fluid seepage resistance.

[0056] After obtaining the two-phase fluid seepage resistance, the development plan can be directly formulated based on the calculated two-phase fluid seepage resistance. For example, during the specification of the development plan, if there is a need to directly use the two-phase fluid seepage resistance to perform corresponding calculations, the two-phase fluid seepage resistance can be directly brought in to obtain the corresponding calculation results.

[0057] In some embodiments, based on the derived two-phase resistance calculation formula, the factors affecting the starting pressure and flow resistance can be further analyzed. For example, different fluid systems and nanopore sizes may also affect the size of the two-phase fluid seepage resistance. The experimental data obtained from the measurement can further more accurately verify the influence of other factors.

[0058] The following uses the oil-water two-phase flow resistance calculated from each group of actual samples to illustrate the corresponding technical effects. Figure 4A 、 Figure 4B 、 Figure 4C 、 Figure 4D 、 Figure 4E and Figure 4F The figure below shows a graph of the starting pressure of an alkane-water system as a function of pore size for different nanopore sizes. As can be seen from the figure, the starting pressure exhibits different trends with pore size for different nanopore sizes, demonstrating the impact of nanopore size on starting pressure.

[0059] like Figure 5A 、 Figure 5B 、 Figure 5C 、 Figure 5D and Figure 5E The figure below shows the curves of flow resistance changing with fluidity under different fluid systems. It can be seen that different fluid systems have a significant impact on the change of flow resistance, which also proves the effect of fluid system on flow resistance.

[0060] In some embodiments, after obtaining accurate two-phase fluid flow resistance, these data can be quantitatively analyzed to obtain corresponding analysis conclusions based on actual changes and comparisons of the data.

[0061] Figure 6A and Figure 6B Schematic diagrams showing the starting pressure and flow resistance caused by dynamic wetting hysteresis at different pore sizes are shown. It can be seen intuitively from the figure that the seepage resistance of the two-phase fluid increases with the decrease of pore size, and the flow resistance is greater than the starting pressure. Therefore, the flow resistance in the nano-scale pores of shale has a greater impact on the ultimate recovery rate, and this conclusion can be used to guide production development.

[0062] Based on the above implementation, the calculated two-phase fluid seepage resistance can be used to guide the formulation of development plans from multiple levels, thereby more effectively determining the development plan and effectively improving the recovery rate.

[0063] Through the introduction of the above embodiments and scenario examples, it can be seen that the method first analyzes the physical parameters and fluid parameters of the experimental samples, especially the contact angle parameters among the fluid parameters, and screens the resistance control factors from the test factors, and then calculates the two-phase fluid seepage resistance in the target formation based on the resistance control factors obtained by screening, and finally uses the two-phase fluid seepage resistance to determine the final development plan. Through the above method, when determining the development plan based on the resistance encountered by the oil phase liquid and the water phase liquid in the formation when flowing, it is possible to predetermine the control factors that have a strong influence on the resistance, and then complete the calculation of the seepage resistance in the target formation by combining the control factors and the contact angle between the fluids, thereby providing guidance for the determination of the development plan, improving the effectiveness of the plan, and ensuring the normal progress of production and development.

[0064] Based on the above-mentioned method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores, this specification also proposes an embodiment of a device for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores. Figure 7 As shown, the device for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores specifically includes the following modules.

[0065] The parameter acquisition module 710 is used to acquire fluid parameters corresponding to at least one group of experimental samples; the experimental samples have different test factor attribute values ​​corresponding to the test factors; the fluid parameters include the contact angle parameters between the water phase fluid and the oil phase fluid.

[0066] The resistance control factor screening module 720 is used to screen the resistance control factor from the test factors according to the fluid parameters.

[0067] The two-phase fluid seepage resistance calculation module 730 is used to calculate the two-phase fluid seepage resistance of the target formation based on the resistance control factor.

[0068] The development plan determination module 740 is used to determine a development plan using the two-phase fluid seepage resistance.

[0069] Although the process flows described above include multiple operations occurring in a particular order, it should be understood that these processes may include more or fewer operations, which may be performed sequentially or in parallel (eg, using parallel processors or a multi-threaded environment).

[0070] Although the process flows described above include multiple operations occurring in a particular order, it should be understood that these processes may include more or fewer operations, which may be performed sequentially or in parallel (eg, using parallel processors or a multi-threaded environment).

[0071] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present specification. 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 device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 a process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0072] These computer program instructions may also be stored in a computer readable memory that can direct a computer or other programmable data processing device to work 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 The function specified in one or more boxes.

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

[0074] In a typical configuration, a computing device includes one or more processors (CPUs), input / output interfaces, network interfaces, and memory.

[0075] Memory may include non-permanent storage in a computer-readable medium, random access memory (RAM) and / or non-volatile memory in the form of read-only memory (ROM) or flash RAM. Memory is an example of a computer-readable medium.

[0076] Computer-readable media includes permanent and non-permanent, removable and non-removable media that can be implemented by any method or technology to store information. The information can be computer-readable instructions, data structures, program modules or other data. Examples of computer storage media include, but are not limited to, phase change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technology, compact disc read-only memory (CD-ROM), digital versatile disc (DVD) or other optical storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices or any other non-transmission media that can be used to store information that can be accessed by a computing device. As defined herein, computer-readable media does not include transitory computer-readable media (transitory media), such as modulated data signals and carrier waves.

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

[0078] Embodiments of this specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, and the like that perform specific tasks or implement specific abstract data types. Embodiments of this specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected via a communications network. In distributed computing environments, program modules may be located in local and remote computer storage media, including storage devices.

[0079] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between the various embodiments can be referenced across them. Each embodiment focuses on the differences from the other embodiments. In particular, since the system embodiments are generally similar to the method embodiments, their description is relatively simple. For relevant parts, reference can be made to the description of the method embodiments. Throughout this specification, reference to the terms "one embodiment," "some embodiments," "example," "specific example," or "some examples" means that the specific features, structures, materials, or characteristics described in conjunction with that embodiment or example are included in at least one embodiment or example of the embodiments in this specification. In this specification, the schematic representations of these terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, those skilled in the art may combine and integrate the different embodiments or examples, and features of different embodiments or examples, described in this specification, without conflict.

[0080] The foregoing is merely an embodiment of the present application and is not intended to limit the present application. For those skilled in the art, the present application may have various changes and variations. Any modifications, equivalent replacements, improvements, etc. made within the spirit and principles of the present application should all be included within the scope of the claims of the present application.

Claims

1. A method for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores, characterized in that: The two-phase fluid comprises an aqueous phase fluid and an oil phase fluid; and the method comprises: Collecting fluid parameters corresponding to at least one set of experimental samples; the experimental samples have different test factor attribute values ​​corresponding to the test factors; the fluid parameters include contact angle parameters between the water phase fluid and the oil phase fluid; screening a resistance control factor from test factors according to the fluid parameters; The two-phase fluid seepage resistance of the target formation is calculated based on the resistance control factor; the two-phase fluid seepage resistance includes the starting pressure and the flow resistance; the two-phase fluid seepage resistance of the target formation is calculated based on the resistance control factor, including: using the formula Calculate the starting pressure, where P f1 is the starting pressure, γ 12 is the interfacial tension between the wetting and non-wetting phases, r is the pore radius, θ R is the receding contact angle, θ A is the advancing contact angle; using the formula Formula for determining the difference between the cosine of the receding contact angle and the cosine of the advancing contact angle The formula Substitute into the formula Get the formula Where, P f2 is the flow resistance, k B is the Boltzmann constant, T is the temperature, r is the pore radius, λ is the average molecular transition distance, U is the moving speed of the three-phase contact line, K 0 is the average molecular transition frequency; using the formula Calculate flow resistance; The development plan is determined using the two-phase fluid seepage resistance.

2. The method according to claim 1, wherein The two-phase fluid seepage resistance includes a starting pressure; the starting pressure includes the difference in capillary forces corresponding to the water phase fluid and the oil phase fluid in a static state; when the two-phase fluid seepage resistance is the starting pressure, the fluid parameters include the surface tension parameters of the two-phase fluid; the surface tension parameters are obtained by the hanging drop method.

3. The method according to claim 2, wherein The experimental sample includes at least one of an octane-water system, a dodecane-water system, and a hexadecane-water system on a fluoropolymer coating on a dry mica sheet.

4. The method according to claim 1, wherein The two-phase fluid seepage resistance includes flow resistance; the flow resistance includes the difference in capillary forces corresponding to the water phase fluid and the oil phase fluid in a moving state; when the two-phase fluid seepage resistance is flow resistance, the fluid parameters include at least one of the molecular average transition frequency and the molecular average transition distance.

5. The method according to claim 4, wherein The experimental samples include at least one of a hexane-water system and a squalane-water system on a single crystal silicon wafer, and a C8-water system, a C12-water system and a C13-water system on apatite.

6. The method according to claim 1, wherein The step of screening the resistance control factor from the test factors according to the fluid parameters includes: Plot the contact angle parameters against the test factors; When the rate of change of the curve graph is greater than a change threshold, the test factor is determined to be a resistance control factor.

7. The method according to claim 1, wherein The contact angle parameters include an advancing contact angle and a receding contact angle; the advancing contact angle includes the contact angle formed when water displaces oil and water occupies part of the original space of oil, and the receding contact angle includes the contact angle formed when oil displaces water and oil occupies part of the original space of water; the resistance control factor includes at least one of a pore radius and a seepage velocity.

8. A device for determining a development plan based on the seepage resistance of two-phase fluid in shale nanopores, characterized in that: The two-phase fluid includes an aqueous phase fluid and an oil phase fluid; the device includes: A parameter acquisition module, configured to acquire fluid parameters corresponding to at least one set of experimental samples; the experimental samples corresponding to the test factors have different test factor attribute values; the fluid parameters include contact angle parameters between the water phase fluid and the oil phase fluid; a resistance control factor screening module, configured to screen resistance control factors from test factors according to the fluid parameters; A two-phase fluid seepage resistance calculation module is used to calculate the two-phase fluid seepage resistance of the target formation based on the resistance control factor; the two-phase fluid seepage resistance includes the starting pressure and the flow resistance; the calculation of the two-phase fluid seepage resistance of the target formation based on the resistance control factor includes: using the formula Calculate the starting pressure, where P f1 is the starting pressure, γ 12 is the interfacial tension between the wetting and non-wetting phases, r is the pore radius, θ R is the receding contact angle, θ A is the advancing contact angle; using the formula Formula for determining the difference between the cosine of the receding contact angle and the cosine of the advancing contact angle The formula Substitute into the formula Get the formula Where, P f2 is the flow resistance, k B is the Boltzmann constant, T is the temperature, r is the pore radius, λ is the average molecular transition distance, U is the moving speed of the three-phase contact line, K 0 is the average molecular transition frequency; Using the formula Calculate flow resistance; The development plan determination module is used to determine the development plan using the two-phase fluid seepage resistance.

Images