Method and apparatus for determining dynamic wetting angle

Abstract

The specification provides a method and device for determining a dynamic wetting angle. The method comprises: obtaining a capillary model radius, a displacement pressure difference, an initial wetting angle, and a first liquid film thickness; constructing a first correspondence relationship between a target included angle and the first liquid film thickness according to a fluid phase interface trajectory equation; constructing a second correspondence relationship between a separation pressure and the first liquid film thickness according to a separation pressure formula; determining a second liquid film thickness according to the capillary model radius, the displacement pressure difference, and the initial wetting angle; obtaining a target wetting angle according to the first correspondence relationship between the target included angle and the first liquid film thickness, the second correspondence relationship between the separation pressure and the first liquid film thickness, and the second liquid film thickness; detecting whether a difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold value; and taking the target wetting angle as the dynamic wetting angle in the case that the difference is less than the wetting angle threshold value. The method can accurately calculate the dynamic wetting angle in the two-phase seepage process of the capillary model.

Description

Technical Field

[0001] This specification pertains to the field of oil and gas reservoir seepage technology, and particularly relates to a method and apparatus for determining dynamic wetting angle. Background Technology

[0002] In oil and gas reservoir environments, the walls of reservoir pore throat structures exhibit different wettabilities (i.e., selective wetting) to different fluids due to the influence of mineral composition and liquid phase environment. This selective wetting of reservoir pore throat structures has a significant impact on the subsequent calculation of wetting angle.

[0003] In existing technologies, separation pressure theory is mostly used to characterize the fluid-structure interaction in microscale pore-throat structures. Based on this fluid-structure interaction, the wettability of the pore-throat structure wall is further obtained, and finally, the dynamic wetting angle is derived from the wettability. In this process, the physical model used is mineral-unstable fluid film-free fluid. The unstable fluid film in this model is not applicable to complex external seepage processes, thus leading to significant calculation errors.

[0004] There is currently no effective solution to the aforementioned technical problems. Summary of the Invention

[0005] This specification presents a method and apparatus for determining the dynamic wetting angle, which can accurately determine the dynamic wetting angle in the two-phase flow process of a capillary model.

[0006] The purpose of the embodiments in this specification is to provide a method for determining the dynamic wetting angle, including:

[0007] Obtain the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness;

[0008] Based on the fluid phase interface trajectory equation, a first correspondence between the target included angle and the first liquid film thickness is constructed.

[0009] Based on the separation pressure formula, a second correspondence between the separation pressure and the thickness of the first liquid film is established;

[0010] The thickness of the second liquid film is determined based on the capillary model radius, the displacement pressure difference, and the initial wetting angle.

[0011] The target wetting angle is obtained based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness.

[0012] Detect whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold;

[0013] If the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold, the target wetting angle is taken as the dynamic wetting angle.

[0014] Furthermore, in another embodiment of the method, after detecting whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold, the method further includes:

[0015] If the difference between the initial wetting angle and the target wetting angle is greater than or equal to the wetting angle threshold, the target wetting angle is taken as the initial wetting angle.

[0016] The thickness of the second liquid film is determined based on the capillary model radius, the displacement pressure difference, and the initial wetting angle.

[0017] The target wetting angle is obtained based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness.

[0018] Detect whether the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold.

[0019] Furthermore, in another embodiment of the method, the step of constructing a second correspondence between the separation pressure and the thickness of the first liquid film according to the separation pressure formula includes:

[0020] A second correspondence between the separation pressure and the thickness of the first liquid film is established according to the following separation pressure formula:

[0021]

[0022] Among them, Π h The separation pressure is represented by A, the Hamaker constant by h, the thickness of the first liquid film by λ, and the London wavelength by C. b N represents particle concentration. A y1 represents Avogadro's constant, k represents Boltzmann's constant, T represents the Kelvin temperature, y1 represents the surface potential of the first fluid surface, y2 represents the surface potential of the second fluid surface, κ represents the reciprocal of the Debye length, and A represents the surface potential of the second fluid surface. K h represents the structural force coefficient. s denoted by , cosh represents the hyperbolic cosine function, and sinh represents the hyperbolic sine function.

[0023] Furthermore, in another embodiment of the method, determining the second liquid film thickness based on the capillary model radius, the displacement pressure difference, and the initial wetting angle includes:

[0024] Calculate the equilibrium pressure based on the capillary model radius, the displacement pressure difference, and the initial wetting angle.

[0025] The thickness of the second liquid film is determined based on the equilibrium pressure.

[0026] Furthermore, in another embodiment of the method, calculating the equilibrium pressure based on the capillary model radius, the displacement pressure difference, and the initial wetting angle includes:

[0027] Calculate the equilibrium pressure using the following formula:

[0028]

[0029] in, ΔP represents the equilibrium pressure, σ represents the displacement pressure difference, θ represents the interfacial tension, θ represents the initial wetting angle, and r represents the capillary model radius.

[0030] Furthermore, in another embodiment of the method, obtaining the target wetting angle based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness includes:

[0031] Calculate the target wetting angle using the following formula:

[0032]

[0033] Where θ' represents the target wetting angle, α represents the target included angle, and cosα is determined based on the first correspondence between the target included angle and the thickness of the first liquid film. h Indicates separation pressure, Π h The thickness of the first liquid film is determined based on the second correspondence between the separation pressure and the thickness of the first liquid film, where h represents the thickness of the first liquid film and h0 represents the thickness of the second liquid film.

[0034] Furthermore, in another embodiment of the method, the wetting angle threshold is 0.05 degrees or 0.1 degrees.

[0035] Furthermore, in another embodiment of the method, obtaining the initial wetting angle includes: measuring the initial wetting angle based on a contact angle measuring instrument.

[0036] On the other hand, embodiments of this specification also provide a device for determining a dynamic wetting angle, including:

[0037] The acquisition module is used to acquire the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness;

[0038] The first construction module is used to construct a first correspondence between the target included angle and the first liquid film thickness based on the fluid phase interface trajectory equation;

[0039] The second construction module is used to construct a second correspondence between the separation pressure and the thickness of the first liquid film based on the separation pressure formula;

[0040] The first calculation module is used to determine the thickness of the second liquid film based on the capillary model radius, the displacement pressure difference, and the initial wetting angle.

[0041] The second calculation module is used to obtain the target wetting angle based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness;

[0042] The detection module is used to detect whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold.

[0043] The determining module is used to determine the target wetting angle as the dynamic wetting angle when the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold.

[0044] Furthermore, embodiments of this specification also provide a computer-readable storage medium storing computer instructions thereon, wherein the computer-readable storage medium executes the instructions to implement the above-described method for determining the dynamic wetting angle.

[0045] This specification provides a method for determining a dynamic wetting angle, which involves obtaining the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness; constructing a first correspondence between the target angle and the first liquid film thickness based on the fluid phase interface trajectory equation; constructing a second correspondence between the separation pressure and the first liquid film thickness based on the separation pressure formula; determining the second liquid film thickness based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; obtaining the target wetting angle based on the first correspondence between the target angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness; detecting whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold; and if the difference between the initial wetting angle and the target wetting angle is determined to be less than the wetting angle threshold, using the target wetting angle as the dynamic wetting angle.

[0046] Furthermore, when determining the thickness of the second liquid film based on the capillary model radius, the displacement pressure difference, and the initial wetting angle, the equilibrium pressure is calculated based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; and the thickness of the second liquid film is determined based on the equilibrium pressure. Attached Figure Description

[0047] To more clearly illustrate the embodiments of this specification, the accompanying drawings used in the embodiments will be briefly introduced below. The drawings described below are only some embodiments recorded in this specification. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0048] Figure 1 This is a flowchart illustrating an embodiment of a method for determining a dynamic wetting angle provided in this specification;

[0049] Figure 2 This is a schematic diagram of a capillary model provided in this manual;

[0050] Figure 3 This is a schematic diagram of a capillary model in equilibrium provided in this manual;

[0051] Figure 4 This is a diagram showing the predicted location of the seepage front in a capillary model in a specific scenario example of this manual;

[0052] Figure 5 This is a schematic diagram of the module structure of one embodiment of the dynamic wetting angle determination device provided in this specification;

[0053] Figure 6 This is a schematic diagram of the structural composition of a server provided in this manual. Detailed Implementation

[0054] To enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this specification, and not all embodiments. Based on the embodiments in this specification, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of this specification.

[0055] In oil and gas reservoir environments, the walls of reservoir pore throat structures exhibit different wettabilities (i.e., selective wetting) to different fluids due to the influence of mineral composition and liquid phase environment. This selective wetting of reservoir pore throat structures has a significant impact on the subsequent calculation of wetting angle.

[0056] Considering that existing technologies mostly employ separation pressure theory to characterize fluid-structure interaction in microscale pore-throat structures, and then further derive the wettability of the pore-throat structure wall based on fluid-structure interaction, and finally obtain the dynamic wetting angle based on the wettability. In the above process, the physical model used is mineral-unstable fluid film-free fluid. The unstable fluid film in this model is not applicable to complex external seepage processes and is easily unstable due to external factors, thus leading to large calculation errors.

[0057] In view of the above-mentioned problems in existing methods and the specific reasons for these problems, this application introduces a method for determining the dynamic wetting angle based on the separation pressure and displacement pressure difference, so as to accurately obtain the dynamic wetting angle.

[0058] Based on the above approach, this specification proposes a method for determining the dynamic wetting angle. First, the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness are obtained. Based on the fluid interface trajectory equation, a first correspondence is established between the target angle and the first liquid film thickness. Based on the separation pressure formula, a second correspondence is established between the separation pressure and the first liquid film thickness. The second liquid film thickness is determined based on the capillary model radius, the displacement pressure difference, and the initial wetting angle. Then, the target wetting angle is obtained based on the first correspondence between the target angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness. Finally, it is detected whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold. If the difference between the initial wetting angle and the target wetting angle is determined to be less than the wetting angle threshold, the target wetting angle is taken as the dynamic wetting angle.

[0059] See Figure 1 As shown in the embodiments of this specification, a method for determining the dynamic wetting angle is provided. In specific implementation, this method may include the following:

[0060] S101: Obtain the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness.

[0061] In some embodiments, the method for determining the dynamic wetting angle described above employs a "mineral-stable fluid film-free fluid" physical model and is implemented based on a capillary model. (See [reference needed]). Figure 2 As shown, Figure 2This is a schematic diagram of a capillary model. The wall of the capillary model is an ideal surface, and inside the model are a first fluid and a second fluid. The first fluid is the displacing phase fluid, which adsorbs onto the wall of the capillary model and forms a microscopic liquid film (liquid film) between the second fluid and the wall of the capillary model. The second fluid is the displaced phase fluid, and the interface between the first and second fluids is a circular arc. The macroscopic contact angle between the first fluid and the wall of the capillary model is the initial wetting angle, denoted as θ.

[0062] In some embodiments, during the process of the first fluid displacing the second fluid, the first fluid and the second fluid gradually reach a state of force equilibrium, such as... Figure 3 As shown, a force analysis is performed on the phase interface in equilibrium, where the phase interface is simultaneously subjected to the displacement pressure difference ΔP and the separation pressure Π. h The two forces act in opposite directions. The angle between the displacement pressure difference and the horizontal plane is called the target angle, denoted by α. The thickness of the first liquid film is the vertical distance between the second fluid and the wall of the capillary model in the initial state, denoted by h. The thickness of the second liquid film is the distance between the second fluid and the wall of the capillary model in the equilibrium state, denoted by h0. The interfacial tension is denoted by σ, and the interfacial tension can be obtained by the hanging ring method. Figure 3 Within the rectangular area, a cross-sectional view of the capillary model is shown in the vertical direction. The interior of the capillary model wall is covered with a thin film of liquid, and the radius of the capillary model is denoted as r.

[0063] In some embodiments, the initial wetting angle can be obtained by measuring the initial wetting angle using a contact angle measuring instrument.

[0064] S102: Based on the fluid phase interface trajectory equation, construct the first correspondence between the target included angle and the first liquid film thickness.

[0065] In some embodiments, when the initial wetting angle and the capillary model radius are fixed, the fluid phase interface trajectory equation is also fixed. Therefore, the first correspondence between the target angle and the first liquid film thickness can be constructed through the fluid phase interface trajectory equation.

[0066] In some embodiments, a first correspondence is established between the target included angle and the thickness of the first liquid film, which may specifically be: establishing a first correspondence between cosα and h.

[0067] S103: Based on the separation pressure formula, construct a second correspondence between the separation pressure and the thickness of the first liquid film.

[0068] In some embodiments, a second correspondence between the separation pressure and the thickness of the first liquid film is established according to the separation pressure formula. Specifically, this may include:

[0069] A second correspondence between the separation pressure and the thickness of the first liquid film is established according to the following separation pressure formula:

[0070]

[0071] Among them, Π h The separation pressure is represented by A, the Hamaker constant by h, the thickness of the first liquid film by λ, and the London wavelength by C. b N represents particle concentration. A y1 represents Avogadro's constant, k represents Boltzmann's constant, T represents the Kelvin temperature, y1 represents the surface potential of the first fluid surface, y2 represents the surface potential of the second fluid surface, κ represents the reciprocal of the Debye length, and A represents the surface potential of the second fluid surface. K h represents the structural force coefficient. s denoted by , cosh represents the hyperbolic cosine function, and sinh represents the hyperbolic sine function.

[0072] S104: Determine the thickness of the second liquid film based on the capillary model radius, the displacement pressure difference, and the initial wetting angle.

[0073] In some embodiments, the thickness of the second liquid film is determined based on the capillary model radius, the displacement pressure difference, and the initial wetting angle. Specifically, this may include:

[0074] S1: Calculate the equilibrium pressure based on the capillary model radius, the displacement pressure difference, and the initial wetting angle;

[0075] S2: Determine the thickness of the second liquid film based on the equilibrium pressure.

[0076] In some embodiments, equilibrium pressure refers to the separation pressure under equilibrium conditions.

[0077] In some embodiments, the equilibrium pressure is calculated based on the capillary model radius, the displacement pressure difference, and the initial wetting angle, including:

[0078] Calculate the equilibrium pressure using the following formula:

[0079]

[0080] Among them, Π h0 ΔP represents the equilibrium pressure, σ represents the displacement pressure difference, θ represents the interfacial tension, θ represents the initial wetting angle, and r represents the capillary model radius.

[0081] In some embodiments, after obtaining the equilibrium pressure, the thickness of the second liquid film can be calculated according to the separation pressure formula.

[0082] S105: The target wetting angle is obtained based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness.

[0083] In some embodiments, the target wetting angle is obtained based on a first correspondence between the target included angle and the first liquid film thickness, a second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness, including:

[0084] Calculate the target wetting angle using the following formula:

[0085]

[0086] Where θ' represents the target wetting angle, α represents the target included angle, and cosα is determined based on the first correspondence between the target included angle and the thickness of the first liquid film. h Indicates separation pressure, Π h The thickness of the first liquid film is determined based on the second correspondence between the separation pressure and the thickness of the first liquid film, where h represents the thickness of the first liquid film and h0 represents the thickness of the second liquid film.

[0087] S106: Detect whether the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold.

[0088] In some embodiments, the aforementioned wetting angle threshold can be 0.05 degrees or 0.1 degrees. It should be noted, however, that other wetting angle thresholds can be selected depending on the actual application scenario, and this specification does not limit such selection.

[0089] S107: If the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold, the target wetting angle is taken as the dynamic wetting angle.

[0090] In some embodiments, after obtaining the dynamic wetting angle, the position of the seepage front in the capillary model can also be predicted based on the dynamic wetting angle.

[0091] Specifically, the two-phase flow process in the capillary model conforms to the Wash-burn equation:

[0092]

[0093] Where t represents time, μ1 represents the viscosity of the first fluid, μ2 represents the viscosity of the second fluid, and L represents the position of the seepage front in the capillary model. c This represents the total length of the capillary model.

[0094] Introducing the separation pressure, the above formula can be transformed into the following form:

[0095]

[0096] Among them, Π eff Indicates the effective separation pressure.

[0097] The above formula can be used to obtain the position of the seepage front in the capillary model at different times.

[0098] In some embodiments, after detecting whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold, the method further includes:

[0099] S1: If the difference between the initial wetting angle and the target wetting angle is greater than or equal to the wetting angle threshold, the target wetting angle is taken as the initial wetting angle;

[0100] S2: Determine the thickness of the second liquid film based on the capillary model radius, the displacement pressure difference, and the initial wetting angle;

[0101] S3: The target wetting angle is obtained based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness;

[0102] S4: Detect whether the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold.

[0103] In some embodiments, if the difference between the initial wetting angle and the target wetting angle is determined to be greater than or equal to a wetting angle threshold, the target wetting angle is taken as the initial wetting angle; and the second liquid film thickness is determined again based on the initial wetting angle, the capillary model radius, and the displacement pressure difference; and the target wetting angle is calculated again based on the second liquid film thickness. It is then checked whether the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold. If the difference between the initial wetting angle and the target wetting angle is determined to be greater than or equal to the wetting angle threshold, the above operation is repeated until the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold. At this point, the target wetting angle is taken as the dynamic wetting angle.

[0104] By comprehensively considering the effects of fluid-solid interaction and displacement pressure difference in microscale seepage channels, the dynamic wetting angle can be accurately calculated in the two-phase seepage process of the capillary model.

[0105] In a specific scenario example, the method proposed in this application can be used to determine the dynamic wetting angle. The capillary model has a radius of 50 nm and is made of quartz. The first fluid is an aqueous surfactant solution, and the second fluid is a decane liquid. According to the separation pressure theory, the liquid film (i.e., the stable fluid film) is composed of the surfactant aqueous solution. The influence of the displacement pressure difference is temporarily disregarded, so the displacement pressure difference is set to 0 Pa. Based on the properties of the capillary model, the first fluid, and the second fluid, the Hamaker constant is taken as -0.045 × 10⁻²⁰ J. The calculated dynamic wetting angle in this scenario example is 14.95°. The dynamic wetting angle calculated using traditional methods is 19.70°. Since it is difficult to obtain the dynamic wetting angle using actual measurement methods, it is necessary to predict the position of the seepage front in the capillary model based on the dynamic wetting angle, in order to verify the reliability of the method proposed in this application based on the position of the seepage front in the capillary model. Furthermore, based on the dynamic wetting angle obtained by the method proposed in this application and the dynamic wetting angle obtained by the traditional method, the position of the seepage front in the capillary model is predicted, and the prediction results are as follows: Figure 4 As shown, Figure 4 The actual percolation curve in the model represents the position of the flow front in the capillary model obtained through actual measurement methods. Figure 4 It can be seen that by using the method proposed in this application to obtain the dynamic wetting angle and predict the position of the seepage front in the capillary model, the prediction result is closer to the actual seepage value. Therefore, the method proposed in this application has higher accuracy than the traditional method.

[0106] Based on the above method for determining the dynamic wetting angle, this specification also provides an embodiment of a device for determining the dynamic wetting angle, see reference. Figure 5 As shown, the device for determining the dynamic wetting angle specifically includes the following modules:

[0107] The acquisition module 501 is used to acquire the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness.

[0108] The first construction module 502 is used to construct a first correspondence between the target included angle and the first liquid film thickness based on the fluid phase interface trajectory equation;

[0109] The second construction module 503 is used to construct a second correspondence between the separation pressure and the thickness of the first liquid film according to the separation pressure formula;

[0110] The first calculation module 504 is used to determine the thickness of the second liquid film based on the capillary model radius, the displacement pressure difference, and the initial wetting angle.

[0111] The second calculation module 505 is used to obtain the target wetting angle based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness;

[0112] Detection module 506 is used to detect whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold.

[0113] The determining module 507 is used to determine the target wetting angle as the dynamic wetting angle when the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold.

[0114] In some embodiments, the acquisition module 501 described above can be used to acquire the initial wetting angle in the following manner: the initial wetting angle is obtained based on the measurement of a contact angle measuring instrument.

[0115] In some embodiments, the second construction module 503 described above can be specifically used to construct a second correspondence between the separation pressure and the thickness of the first liquid film according to the following separation pressure formula:

[0116]

[0117] Among them, Π h The separation pressure is represented by A, the Hamaker constant by h, the thickness of the first liquid film by λ, and the London wavelength by C. b N represents particle concentration. A y1 represents Avogadro's constant, k represents Boltzmann's constant, T represents the Kelvin temperature, y1 represents the surface potential of the first fluid surface, y2 represents the surface potential of the second fluid surface, κ represents the reciprocal of the Debye length, and A represents the surface potential of the second fluid surface. K h represents the structural force coefficient. s denoted by , cosh represents the hyperbolic cosine function, and sinh represents the hyperbolic sine function.

[0118] In some embodiments, the first calculation module 504 described above can be specifically used to calculate the equilibrium pressure based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; and to determine the thickness of the second liquid film based on the equilibrium pressure.

[0119] In some embodiments, the second calculation module 505 described above can be specifically used to calculate the target wetting angle according to the following formula:

[0120]

[0121] Where θ' represents the target wetting angle, α represents the target included angle, and cosα is determined based on the first correspondence between the target included angle and the thickness of the first liquid film. hIndicates separation pressure, Π h The thickness of the first liquid film is determined based on the second correspondence between the separation pressure and the thickness of the first liquid film, where h represents the thickness of the first liquid film and h0 represents the thickness of the second liquid film.

[0122] In some embodiments, the above-described dynamic wetting angle determining device can further be used to: take the target wetting angle as the initial wetting angle when the difference between the initial wetting angle and the target wetting angle is greater than or equal to a wetting angle threshold; determine the second liquid film thickness based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; obtain the target wetting angle based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness; and detect whether the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold.

[0123] It should be noted that the units, devices, or modules described in the above embodiments can be implemented by computer chips or physical entities, or by products with certain functions. For ease of description, the above devices are described by dividing them into various modules according to their functions. Of course, in implementing this specification, the functions of each module can be implemented in one or more software and / or hardware, or the module that implements the same function can be implemented by a combination of multiple sub-modules or sub-units, etc. The device embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and there may be other division methods in actual implementation. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Furthermore, the coupling or direct coupling or communication connection between the devices or units shown or discussed can be through some interfaces, and the indirect coupling or communication connection between devices or units can be electrical, mechanical, or other forms.

[0124] This specification also provides a computer storage medium for a method of determining a dynamic wetting angle. The computer storage medium stores computer program instructions that, when executed, implement the following: obtaining the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness; constructing a first correspondence between the target angle and the first liquid film thickness based on the fluid interface trajectory equation; constructing a second correspondence between the separation pressure and the first liquid film thickness based on the separation pressure formula; determining a second liquid film thickness based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; obtaining a target wetting angle based on the first correspondence between the target angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness; detecting whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold; and if the difference between the initial wetting angle and the target wetting angle is determined to be less than the wetting angle threshold, using the target wetting angle as the dynamic wetting angle.

[0125] In this embodiment, the storage medium includes, but is not limited to, Random Access Memory (RAM), Read-Only Memory (ROM), cache, hard disk drive (HDD), or memory card. The memory can be used to store computer program instructions. The network communication unit can be an interface configured according to standards specified in the communication protocol for network connection communication.

[0126] In this embodiment, the specific functions and effects implemented by the program instructions stored in the computer storage medium can be explained in comparison with other implementation methods, and will not be repeated here.

[0127] This specification also provides a server, including a processor and a memory for storing processor-executable instructions. In specific implementations, the processor can perform the following steps according to the instructions: obtaining the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness; constructing a first correspondence between the target angle and the first liquid film thickness based on the fluid interface trajectory equation; constructing a second correspondence between the separation pressure and the first liquid film thickness based on the separation pressure formula; determining a second liquid film thickness based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; obtaining a target wetting angle based on the first correspondence between the target angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness; detecting whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold; and if the difference between the initial wetting angle and the target wetting angle is determined to be less than the wetting angle threshold, using the target wetting angle as the dynamic wetting angle.

[0128] To execute the above instructions more accurately, please refer to... Figure 6 As shown in the embodiments of this specification, another specific server is also provided, wherein the server includes a network communication port 601, a processor 602 and a memory 603, and the above structures are connected by internal cables so that the various structures can perform specific data interaction.

[0129] Specifically, the network communication port 601 can be used to obtain the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness.

[0130] The processor 602 can be specifically configured to: construct a second correspondence between the separation pressure and the thickness of the first liquid film according to the separation pressure formula; determine the thickness of the second liquid film according to the capillary model radius, the displacement pressure difference, and the initial wetting angle; obtain the target wetting angle according to the first correspondence between the target angle and the thickness of the first liquid film, the second correspondence between the separation pressure and the thickness of the first liquid film, and the thickness of the second liquid film; detect whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold; and if it is determined that the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold, use the target wetting angle as the dynamic wetting angle.

[0131] The memory 603 can be used to store the corresponding instruction program.

[0132] In this embodiment, the network communication port 601 can be a virtual port bound to different communication protocols, thereby enabling the sending or receiving of different data. For example, the network communication port can be a port responsible for web data communication, a port responsible for FTP data communication, or a port responsible for email data communication. Furthermore, the network communication port can also be a physical communication interface or communication chip. For example, it can be a wireless mobile network communication chip, such as GSM or CDMA; it can also be a Wi-Fi chip; or it can be a Bluetooth chip.

[0133] In this embodiment, the processor 602 can be implemented in any suitable manner. For example, the processor can take the form of a microprocessor or processor and a computer-readable medium storing computer-readable program code (e.g., software or firmware) executable by the (micro)processor, logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers, and embedded microcontrollers, etc. This specification is not limiting.

[0134] In this embodiment, the memory 603 may include multiple layers. In a digital system, anything that can store binary data can be a memory. In an integrated circuit, a circuit with storage function but no physical form is also called a memory, such as RAM, FIFO, etc. In a system, a storage device with a physical form is also called a memory, such as a memory stick, TF card, etc.

[0135] While this specification provides the steps of operation for the methods described in the embodiments or flowcharts, more or fewer steps may be included based on conventional or non-inventive means. The order of steps listed in the embodiments is merely one possible order of execution among many steps and does not represent the only possible order. In actual device or client product execution, the methods shown in the embodiments or drawings may be executed sequentially or in parallel (e.g., in a parallel processor or multi-threaded processing environment, or even a distributed data processing environment). The terms "comprising," "including," or any other variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, product, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, product, or apparatus. Without further limitations, the presence of other identical or equivalent elements in a process, method, product, or apparatus that includes said elements is not excluded. The terms "first," "second," etc., are used to denote names and do not indicate any particular order.

[0136] Those skilled in the art will also know that, besides implementing the controller using purely computer-readable program code, the same functions can be achieved by logically programming the method steps, making the controller function as logic gates, switches, application-specific integrated circuits (ASICs), programmable logic controllers (PLCs), and embedded microcontrollers. Therefore, such a controller can be considered a hardware component, and the devices within it used to implement various functions can also be considered structures within that hardware component. Alternatively, the devices used to implement various functions can be considered as both software modules implementing the method and structures within a hardware component.

[0137] This specification can be described in the general context of computer-executable instructions that are executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, classes, etc., that perform a specific task or implement a specific abstract data type. This specification can also be practiced in distributed computing environments, where tasks are performed by remote processing devices connected via a communication network. In distributed computing environments, program modules can reside in local and remote computer storage media, including storage devices.

[0138] As can be seen from the above description of the embodiments, those skilled in the art can clearly understand that this specification can be implemented by means of software plus necessary general-purpose hardware platforms. Based on this understanding, the technical solutions of this specification can essentially be embodied in the form of a software product. This computer software product can be stored in a storage medium, such as ROM / RAM, magnetic disk, optical disk, etc., and includes several instructions to cause a computer device (which may be a personal computer, mobile terminal, server, or network device, etc.) to execute the methods described in the various embodiments or some parts of the embodiments of this specification.

[0139] The various embodiments in this specification are described in a progressive manner. Similar or identical parts between embodiments can be referred to interchangeably. Each embodiment focuses on its differences from other embodiments. This specification can be used in numerous general-purpose or special-purpose computer system environments or configurations. Examples include: personal computers, server computers, handheld or portable devices, tablet devices, multiprocessor systems, microprocessor-based systems, set-top boxes, programmable electronic devices, network PCs, minicomputers, mainframe computers, and distributed computing environments including any of the above systems or devices, etc.

[0140] Although this specification has been described by way of examples, those skilled in the art will recognize that many variations and modifications are possible without departing from the spirit of this specification, and it is intended that the appended claims cover such variations and modifications without departing from the spirit of this specification.

Claims

1. A method for determining a dynamic wetting angle, characterized in that, include: Obtain the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness; Based on the fluid phase interface trajectory equation, a first correspondence between the target included angle and the first liquid film thickness is constructed. Based on the separation pressure formula, a second correspondence between the separation pressure and the thickness of the first liquid film is established; The second liquid film thickness is determined based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; this includes: calculating the equilibrium pressure based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; and determining the second liquid film thickness based on the equilibrium pressure. The target wetting angle is obtained based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness. Detect whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold; If the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold, the target wetting angle is taken as the dynamic wetting angle. Specifically, establishing a second correspondence between the separation pressure and the thickness of the first liquid film according to the separation pressure formula includes: establishing a second correspondence between the separation pressure and the thickness of the first liquid film according to the following separation pressure formula: in, Indicates separation pressure, A Represents the Hamaker constant. h Indicates the thickness of the first liquid film. Indicates London wavelength. Indicates particle concentration. Represents Avogadro's constant. k Represents the Boltzmann constant. Indicates Kelvin temperature. This represents the surface potential of the first fluid surface. This represents the surface potential of the second fluid surface. This represents the reciprocal of the length of Debye. Indicates the structural force coefficient. The characteristic decay length is represented by cosh, which represents the hyperbolic cosine function, and sinh represents the hyperbolic sine function. Based on the first correspondence between the target wetting angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness, the target wetting angle is obtained, including: Calculate the target wetting angle using the following formula: in, Indicates the target wetting angle. Indicates the included angle of the target. Determined based on the first correspondence between the target angle and the thickness of the first liquid film. Indicates separation pressure, Determined based on the second correspondence between separation pressure and the thickness of the first liquid film. h Indicates the thickness of the first liquid film. This indicates the thickness of the second liquid film.

2. The method according to claim 1, characterized in that, After detecting whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold, the method further includes: If the difference between the initial wetting angle and the target wetting angle is greater than or equal to the wetting angle threshold, the target wetting angle is taken as the initial wetting angle. The thickness of the second liquid film is determined based on the capillary model radius, the displacement pressure difference, and the initial wetting angle. The target wetting angle is obtained based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness. Detect whether the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold.

3. The method according to claim 1, characterized in that, The equilibrium pressure is calculated based on the capillary model radius, the displacement pressure difference, and the initial wetting angle, including: Calculate the equilibrium pressure using the following formula: in, Indicates balancing pressure, Indicates the displacement pressure difference. Indicates interfacial tension. Indicates the initial wetting angle. This represents the radius of the capillary model.

4. The method according to claim 1, characterized in that, The wetting angle threshold is 0.05 degrees or 0.1 degrees.

5. The method according to claim 1, characterized in that, Obtaining the initial wetting angle includes: measuring the initial wetting angle using a contact angle measuring instrument.

6. A device for determining a dynamic wetting angle, characterized in that, include: The acquisition module is used to acquire the capillary model radius, displacement pressure difference, initial wetting angle, and first liquid film thickness; The first construction module is used to construct a first correspondence between the target included angle and the first liquid film thickness based on the fluid phase interface trajectory equation; The second construction module is used to construct a second correspondence between the separation pressure and the thickness of the first liquid film based on the separation pressure formula; The first calculation module is used to determine the thickness of the second liquid film based on the capillary model radius, the displacement pressure difference, and the initial wetting angle. The second calculation module is used to obtain the target wetting angle based on the first correspondence between the target included angle and the first liquid film thickness, the second correspondence between the separation pressure and the first liquid film thickness, and the second liquid film thickness; The detection module is used to detect whether the difference between the initial wetting angle and the target wetting angle is less than a wetting angle threshold. The determining module is used to determine the target wetting angle as the dynamic wetting angle when the difference between the initial wetting angle and the target wetting angle is less than the wetting angle threshold. Specifically, the second construction module is used to construct a second correspondence between the separation pressure and the thickness of the first liquid film according to the following separation pressure formula: in, Indicates separation pressure, A Represents the Hamaker constant. h Indicates the thickness of the first liquid film. Indicates London wavelength. Indicates particle concentration. Represents Avogadro's constant. k Represents the Boltzmann constant. Indicates Kelvin temperature. This represents the surface potential of the first fluid surface. This represents the surface potential of the second fluid surface. This represents the reciprocal of the length of Debye. Indicates the structural force coefficient. The characteristic decay length is represented by cosh, which represents the hyperbolic cosine function, and sinh represents the hyperbolic sine function. The first calculation module is specifically used to: calculate the equilibrium pressure based on the capillary model radius, the displacement pressure difference, and the initial wetting angle; and determine the thickness of the second liquid film based on the equilibrium pressure. The second calculation module is specifically used to calculate the target wetting angle according to the following formula: in, Indicates the target wetting angle. Indicates the included angle of the target. Determined based on the first correspondence between the target angle and the thickness of the first liquid film. Indicates separation pressure, Determined based on the second correspondence between separation pressure and the thickness of the first liquid film. h Indicates the thickness of the first liquid film. This indicates the thickness of the second liquid film.

7. A computer-readable storage medium, characterized in that, It stores computer instructions that, when executed by a processor, implement the steps of the method according to any one of claims 1 to 5.

Images