Connectivity evaluation method and apparatus for rock pore channel, and server and medium

Abstract

A connectivity evaluation method and apparatus for a rock pore channel, and a server and a medium. The method comprises: acquiring mercury intrusion pressures and mercury intrusion volumes corresponding to the mercury intrusion pressures at at least two pressure change points for mercury entering a pore channel with a pore structure in rock to be evaluated; and on the basis of a first energy conservation formula, a second energy conservation formula, the mercury injection pressures and the mercury intrusion volumes, determining a connectivity coefficient of the pore channel, wherein the connectivity coefficient is used for evaluating the connectivity degree of the pore channel with the pore structure. Mercury intrusion pressures and mercury intrusion volumes are processed on the basis of the two energy conservation formulas, so as to enable the connectivity coefficient obtained on the basis of the processed mercury intrusion pressures and mercury intrusion volumes to be more accurate, thereby achieving more precise evaluation of the connectivity degree of pore channels in rock to be evaluated.

Description

Method and device for evaluating connectivity of rock pore channel, server and medium

[0001] The present application claims priority to the Chinese patent application No. 202410530514.1, filed on April 29, 2024, and entitled "Method and device for evaluating connectivity of rock pore channel, server and medium", the content of which is incorporated herein by reference in its entirety. TECHNICAL FIELD

[0002] The present application relates to the field of geological exploration and development, and in particular to a method and device for evaluating connectivity of rock pore channel, a server and a medium. BACKGROUND

[0003] The rock has a plurality of interconnected pore channels, and the interconnected pore channels form a pore structure.

[0004] Currently, the evaluation method for the connectivity of the pore structure in the rock uses traditional evaluation methods such as the semi-permeable barrier method, the mercury injection method and the centrifugal method. The traditional evaluation method has a large difference from the actual connectivity of the pore structure of the pore channel and has a low accuracy. SUMMARY

[0005] The embodiments of the present application provide a method and device for evaluating connectivity of rock pore channel, a server and a medium to improve the accuracy of evaluating the connectivity of the pore channel in the rock in the form of a pore structure.

[0006] In a first aspect, the embodiments of the present application provide a method for evaluating connectivity of rock pore channel, the method comprising: obtaining mercury entering a pore channel in a rock to be evaluated in the form of a pore structure, mercury injection pressure at at least two pressure change points and a mercury injection volume corresponding to the mercury injection pressure; determining a connectivity coefficient of the pore channel according to a first energy conservation formula, a second energy conservation formula, the mercury injection pressure and the mercury injection volume, the connectivity coefficient being used to evaluate a connectivity degree of the pore channel in the form of a pore structure; wherein the first energy conservation formula represents that work required for the mercury to enter the pore channel is equal to a sum of pressure work of the mercury injection pressure and volume work of the mercury injection volume, and the second energy conservation formula represents that the work required for the mercury to enter the pore channel is equal to a sum of surface energy generated by the mercury and the pore channel and interface energy generated by the mercury and air.

[0007] In some embodiments of the first aspect, determining the connectivity coefficient of the pore channel according to the first energy conservation formula, the second energy conservation formula, the mercury injection pressure and the mercury injection volume comprises: constructing a slope model of the pore channel in the form of a pore structure based on the first energy conservation formula and the second energy conservation formula; inputting the mercury injection pressure and the mercury injection volume into the slope model to output a slope discrete point set; and calculating the connectivity coefficient of the pore channel according to the slope discrete point set.

[0008] In some embodiments of the first aspect, constructing the slope model of the pore structure-shaped pore channel based on the first energy conservation formula and the second energy conservation formula comprises: constructing a hyperbolic model of the pore structure-shaped pore channel based on the first energy conservation formula and the second energy conservation formula, the hyperbolic model comprising a mercury intrusion pressure curve and a correlation curve, the mercury intrusion pressure curve representing a change relationship between the mercury intrusion pressure and the mercury intrusion volume, and the correlation curve representing a change relationship between the mercury intrusion pressure and the mercury intrusion volume obtained by converting the surface energy and the interfacial energy; and calculating the slope model based on the hyperbolic model.

[0009] In some embodiments of the first aspect, constructing the hyperbolic model of the pore structure-shaped pore channel based on the first energy conservation formula and the second energy conservation formula comprises: converting the equivalent relationship between the surface energy and the interfacial energy in the second energy conservation formula according to a capillary bundle equivalent assumption principle to obtain a third energy conservation formula representing a change relationship between the mercury intrusion pressure and the mercury intrusion volume; and constructing the hyperbolic model of the pore structure-shaped pore channel based on the third energy conservation formula and the first energy conservation formula.

[0010] In some embodiments of the first aspect, the slope model comprises a first tangent slope model of the mercury intrusion pressure curve and a second tangent slope model of the correlation curve; and inputting the mercury intrusion pressure and the mercury intrusion volume into the slope model to output a slope discrete point set, which comprises: inputting the mercury intrusion pressure and the mercury intrusion volume into the first tangent slope model to obtain a first discrete point set; inputting the mercury intrusion pressure and the mercury intrusion volume into the second tangent slope model to obtain a second discrete point set; and calculating the connectivity coefficient of the pore channel according to the slope discrete point set, which comprises: obtaining the connectivity coefficient of the pore channel according to the first discrete point set and the second discrete point set.

[0011] In some embodiments of the first aspect, obtaining the connectivity coefficient of the pore channel according to the first discrete point set and the second discrete point set comprises: determining an absolute value of a covariance between the first discrete point set and the second discrete point set as the connectivity coefficient of the pore channel.

[0012] In some embodiments of the first aspect, the connectivity coefficient of the pore channel is within a preset range, and the preset range is obtained by conversion according to a preset J function curve.

[0013] In a second aspect, a rock pore channel connectivity evaluation device is provided, which comprises,

[0014] The acquisition module is configured to acquire mercury entering a pore structure-shaped pore channel in a rock to be evaluated, mercury intrusion pressures at at least two pressure change points, and mercury intrusion volumes corresponding to the mercury intrusion pressures.

[0015] The calculation module is configured to determine a connectivity coefficient of the pore channel according to the first energy conservation formula, the second energy conservation formula, the mercury intrusion pressures, and the mercury intrusion volumes, and the connectivity coefficient is used to evaluate a connectivity degree of the pore structure-shaped pore channel.

[0016] wherein the first energy conservation formula represents that the work required for the mercury to enter the pore channel is equal to the sum of the pressure work of the mercury applied pressure and the volume work of the mercury entering volume, and the second energy conservation formula represents that the work required for the mercury to enter the pore channel is equal to the sum of the surface energy of the mercury and the pore channel and the interface energy of the mercury and air.

[0017] In some embodiments of the second aspect, the calculation module is specifically configured to:

[0018] construct a slope model of the pore channel in the pore structure based on the first energy conservation formula and the second energy conservation formula;

[0019] input the mercury injection pressure and the mercury entering volume into the slope model, and output a slope discrete point set;

[0020] calculate the connectivity coefficient of the pore channel according to the slope discrete point set.

[0021] In some embodiments of the second aspect, the calculation module is specifically configured to:

[0022] construct a hyperbolic curve model of the pore channel in the pore structure based on the first energy conservation formula and the second energy conservation formula, the hyperbolic curve model including a mercury injection curve and a correlation curve, the mercury injection curve representing a change relationship between the mercury injection pressure and the mercury entering volume, and the correlation curve representing a change relationship between the mercury injection pressure and the mercury entering volume converted from the surface energy and the interface energy;

[0023] calculate the slope model based on the hyperbolic curve model.

[0024] In some embodiments of the second aspect, the calculation module is specifically configured to:

[0025] convert the equivalent relationship between the surface energy and the interface energy in the second energy conservation formula according to the capillary bundle equivalent assumption principle, obtain a change relationship between the mercury injection pressure and the mercury entering volume, and obtain a third energy conservation formula;

[0026] construct a hyperbolic curve model of the pore channel in the pore structure based on the third energy conservation formula and the first energy conservation formula.

[0027] In some embodiments of the second aspect, the slope model includes a first tangent slope model of the mercury injection curve, and a second tangent slope model of the correlation curve;

[0028] The calculation module inputs the mercury injection pressure and the mercury entering volume into the slope model, and outputs the slope discrete point set, and is specifically configured to:

[0029] input the mercury injection pressure and the mercury entering volume into the first tangent slope model to obtain a first discrete point set;

[0030] inputting the mercury injection pressure and the mercury injection volume into the second tangent slope model to obtain a second discrete point set;

[0031] The computing module is configured to calculate the connectivity coefficient of the pore channel according to the slope discrete point set, and specifically configured to:

[0032] The connectivity coefficient of the pore channel is obtained according to the first discrete point set and the second discrete point set.

[0033] In some embodiments of the second aspect, the computing module is specifically configured to:

[0034] The absolute value of the covariance between the first discrete point set and the second discrete point set is determined as the connectivity coefficient of the pore channel.

[0035] In some embodiments of the second aspect, the connectivity coefficient of the pore channel is within a preset range, and the preset range is obtained according to a preset J function curve.

[0036] The third aspect provides a server, comprising:

[0037] a processor, a memory, and a communication interface;

[0038] The memory is configured to store executable instructions of the processor;

[0039] The processor executes the computer executable instructions stored in the memory to perform the following steps:

[0040] obtaining the mercury injection pressure and the mercury injection volume at at least two pressure change points of the mercury entering the pore channel in the rock in a pore structure;

[0041] determining the connectivity coefficient of the pore channel according to the first energy conservation formula, the second energy conservation formula, the mercury injection pressure and the mercury injection volume, and the connectivity coefficient is used to evaluate the connectivity degree of the pore channel in the pore structure;

[0042] The first energy conservation formula represents that the work required for the mercury to enter the pore channel is equal to the sum of the pressure work of the mercury injection pressure applied to the mercury and the volume work of the mercury injection volume, and the second energy conservation formula represents that the work required for the mercury to enter the pore channel is equal to the sum of the surface energy generated by the mercury and the pore channel and the interface energy generated by the mercury and air.

[0043] In some embodiments of the third aspect, the processor executing the computer executable instructions stored in the memory is further configured to perform the following steps:

[0044] constructing a slope model of the pore channel in the pore structure based on the first energy conservation formula and the second energy conservation formula;

[0045] inputting the mercury injection pressure and the mercury injection volume into the slope model to output a slope discrete point set;

[0046] The connectivity coefficient of the pore channel is calculated according to the slope discrete point set.

[0047] In some embodiments of the third aspect, the computer-executable instructions stored in the memory that the processor executes further cause the processor to perform the following steps:

[0048] A hyperbolic model of the pore channel in the pore structure is constructed based on the first energy conservation formula and the second energy conservation formula, the hyperbolic model including a mercury injection curve and a correlation curve, the mercury injection curve representing a change relationship between a mercury injection pressure and a mercury injection volume, and the correlation curve representing a change relationship between the mercury injection pressure and the mercury injection volume obtained by converting the surface energy and the interfacial energy;

[0049] A slope model is calculated based on the hyperbolic model.

[0050] In some embodiments of the third aspect, the computer-executable instructions stored in the memory that the processor executes further cause the processor to perform the following steps:

[0051] An equivalent relationship between the surface energy and the interfacial energy in the second energy conservation formula is converted according to a capillary bundle equivalent assumption principle, to obtain a change relationship between the mercury injection pressure and the mercury injection volume, and to obtain a third energy conservation formula;

[0052] A hyperbolic model of the pore channel in the pore structure is constructed based on the third energy conservation formula and the first energy conservation formula.

[0053] In some embodiments of the third aspect, the slope model includes a first tangent slope model of the mercury injection curve, and a second tangent slope model of the correlation curve.

[0054] The computer-executable instructions stored in the memory that the processor executes further cause the processor to perform the following steps:

[0055] The mercury injection pressure and the mercury injection volume are input into the first tangent slope model to obtain a first discrete point set.

[0056] The mercury injection pressure and the mercury injection volume are input into the second tangent slope model to obtain a second discrete point set.

[0057] A connectivity coefficient of the pore channel is calculated according to the slope discrete point set, including:

[0058] The connectivity coefficient of the pore channel is obtained according to the first discrete point set and the second discrete point set.

[0059] In some embodiments of the third aspect, the computer-executable instructions stored in the memory that the processor executes further cause the processor to perform the following steps:

[0060] An absolute value of a covariance between the first discrete point set and the second discrete point set is determined as the connectivity coefficient of the pore channel.

[0061] In some embodiments of the third aspect, the connectivity coefficient of the pore channel is within a preset range, and the preset range is converted according to a preset J function curve.

[0062] The fourth aspect provides a readable storage medium, which stores a computer program, and the computer program is executed by a processor to implement the method of the first aspect.

[0063] The fifth aspect provides a computer program, which includes a computer program stored in a computer readable storage medium, at least one processor can read the computer program from the computer readable storage medium, and the at least one processor executes the computer program to implement the method of the first aspect.

[0064] The embodiments of the present application provide a rock pore channel connectivity evaluation method, device, server and medium. When mercury enters the pore channel in a gap structure of a rock to be evaluated, the mercury injection pressure of at least two pressure change points is obtained, and the mercury injection volume corresponding to the mercury injection pressure is obtained. The connectivity coefficient of the pore channel of the rock to be evaluated is obtained according to the mercury injection pressure, the mercury injection volume, a first energy conservation formula and a second energy conservation formula, so as to evaluate the connectivity degree of the pore channel of the rock to be evaluated. The mercury injection pressure and the mercury injection volume are processed in combination with the two energy conservation formulas, so that the connectivity coefficient obtained according to the processed mercury injection pressure and the mercury injection volume is more accurate, and the evaluation of the connectivity degree of the pore channel of the rock to be evaluated is more accurate. BRIEF DESCRIPTION OF DRAWINGS

[0065] In order to more clearly illustrate the technical solutions of the embodiments of the present application or the prior art, the drawings required to be used in the embodiments or the prior art description will be briefly introduced. Obviously, the drawings in the following description are some embodiments of the present application, and other drawings can also be obtained by those skilled in the art without creative labor.

[0066] Fig. 1 is a schematic diagram of an application scenario of the rock pore channel connectivity evaluation method provided by the embodiments of the present application;

[0067] Fig. 2 is a flowchart of the rock pore channel connectivity evaluation method provided by the embodiments of the present application;

[0068] Fig. 3 is another schematic diagram of an application scenario of the rock pore channel connectivity evaluation method provided by the embodiments of the present application;

[0069] Fig. 4 is another flowchart of the rock pore channel connectivity evaluation method provided by the embodiments of the present application;

[0070] Fig. 5 is another schematic diagram of an application scenario of the rock pore channel connectivity evaluation method provided by the embodiments of the present application;

[0071] Fig. 6 is a schematic diagram of a mercury injection curve of a rock pore channel obtained based on a traditional evaluation method;

[0072] Fig. 7 is a schematic diagram of a structure of an embodiment of a rock pore channel connectivity evaluation device provided by the present application;

[0073] Fig. 8 is a schematic diagram of a structure of a server provided by an embodiment of the present application. DETAILED DESCRIPTION

[0074] To make the objectives, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments made by those of ordinary skill in the art under the inspiration of the present application belong to the scope of protection of the present application.

[0075] The terms "first", "second", "third", "fourth" and the like (if any) in the specification and claims of the present application and the above drawings are used to distinguish similar objects, and do not necessarily indicate a specific order or a chronological sequence. It should be understood that the data thus used can be interchanged under appropriate circumstances, so that the embodiments of the present application described herein can be implemented in an order other than that illustrated or described herein. In addition, the terms "include" and "have" and any variations thereof are intended to cover non-exclusive inclusion, for example, a process, method, system, product or device that includes a series of steps or units does not necessarily have to be limited to those steps or units clearly listed, but can include other steps or units not clearly listed or inherent to the process, method, product or device.

[0076] The rock has a plurality of interconnected pore channels, and the interconnected pore channels form a pore structure.

[0077] At present, the evaluation method for the connectivity of the pore structure-shaped pore channels in the rock uses traditional evaluation methods such as the semi-permeable barrier method, the mercury injection method and the centrifugal method.

[0078] In the traditional evaluation method, the pore channels are equivalent to capillary tubes, and a capillary pressure curve of the capillary tubes (the curved surface additional pressure of the fluid level rising or falling in the pore channels is equivalent to the capillary pressure) is calculated, and the connectivity of the pore channels is evaluated according to the capillary pressure curve.

[0079] When the fluid includes water and oil, the oil and water phases flow in the pore channel. In a high-permeability oil reservoir, a traditional evaluation method only considers the capillary pressure of the water-oil interface reaching an equilibrium state, and ignores the dynamic effect of the capillary pressure of the water-oil interface not reaching the equilibrium state, so that the calculated capillary pressure curve of the high-permeability oil reservoir is not accurate enough, and the connectivity evaluation of the pore channel is not accurate enough.

[0080] In a pore channel of a low-permeability oil reservoir, due to the small pore throat of the pore channel, the permeability of the pore channel is low, and the water-oil seepage speed in the pore channel is low. Compared with a high-permeability oil reservoir, the dynamic effect of the capillary pressure is more obvious when the water-oil interface does not reach an equilibrium state. The calculated capillary pressure curve of the low-permeability oil reservoir is even more inaccurate, and the connectivity evaluation of the pore channel is even more inaccurate.

[0081] In addition, the traditional evaluation method is derived from clastic rocks such as conglomerate and sandstone, and the analysis object is relatively single, and the applicability is relatively narrow. For example, for a Middle East pore channel type carbonate rock, due to its unique sedimentary environment, diagenesis and pore throat composition of the pore channel, the traditional evaluation method is used to evaluate the Middle East pore channel type carbonate rock, and the calculated capillary pressure curve shows that the stable stage time of the curve is short, the curve shape is complex, and the characteristic section is not obvious. The calculated pore throat distribution of the Middle East pore channel type carbonate rock has a large dispersion degree, the permeability contribution value is dispersed, the connectivity evaluation is not accurate, and the difference with the actual situation of the Middle East pore channel type carbonate rock is large.

[0082] In view of this, the application provides a rock pore channel connectivity evaluation method, device, server and medium to improve the accuracy of evaluating the connectivity of the pore channel in the rock in a pore structure.

[0083] Please refer to FIG. 1, which is an application scenario of the rock pore channel connectivity evaluation method of the embodiment of the application. The application scenario includes a control device 10, a mercury injection device 20, and an evaluation object of the control device 10, which is a rock to be evaluated. The rock to be evaluated has a pore channel in a pore structure. The connectivity evaluation method is applied to the control device 10. Specifically, the control device 10 controls the mercury injection device 20 to press mercury into the pore channel. The control device 10 acquires the mercury injection pressure from the mercury injection device 20 and the mercury injection volume corresponding to the mercury injection pressure, and then determines the connectivity of the pore channel according to the mercury injection pressure, the mercury injection volume, a first energy conservation formula and a second energy conservation formula.

[0084] In some embodiments, the control device can be integrated into the mercury injection device.

[0085] In some embodiments, the control device comprises a server, a mobile phone, a tablet computer, a desktop computer, a laptop computer, a handheld computer, a notebook computer, an ultra-mobile personal computer (UMPC), and the like. In some embodiments, the mercury injection device comprises a mercury injection apparatus.

[0086] It should be noted that FIG. 1 is only a schematic diagram of an application scenario provided by the embodiments of the present application, and the actual forms of various devices included in FIG. 1 are not limited, and the interaction mode between the devices in FIG. 1 is not limited, and in the specific application of the scheme, the actual needs can be set.

[0087] In the following, the technical scheme of the present application is described in detail through specific embodiments. It should be noted that the following specific embodiments can be combined with each other, and the same or similar concepts or processes can not be described in detail in some embodiments.

[0088] FIG. 2 is a flowchart of an embodiment of a rock pore connectivity evaluation method provided by the present application. The execution subject of the embodiments of the present application can be a control device. The method in the present embodiment can be realized by software, hardware or a combination of software and hardware. As shown in FIG. 2, the connectivity evaluation method specifically includes the following steps:

[0089] Step S210: obtaining mercury entering a pore structure-shaped pore in the rock to be evaluated, mercury injection pressure at at least two pressure change points, and mercury injection volume corresponding to the mercury injection pressure.

[0090] As shown in FIG. 3, the rock to be evaluated can be a reservoir rock. The space for fluid flow in the reservoir rock is a complex small pore formed by some curved, different sizes and tortuously connected small pores.

[0091] The mercury injection pressure refers to the force of the mercury injection device to press the mercury into the pore. The mercury injection volume refers to the volume of the mercury entering the pore. The pressure change point refers to the point where the mercury injection pressure changes.

[0092] Taking FIG. 1 as an example, the mercury injection device 20 applies pressure to the mercury to generate the mercury injection pressure. With the application of the mercury injection pressure, the mercury enters the pore to generate the mercury injection volume. Specifically, when the mercury injection device 20 applies pressure to the mercury, the pressure instrument of the mercury injection device 20 can sense the mercury injection pressure, and then the mercury injection device 20 can convert the mercury injection volume according to the mercury injection pressure. That is, each mercury injection pressure has a corresponding mercury injection volume.

[0093] It can be understood that, based on the characteristics of mercury material, the greater the pressure of mercury injection, the greater the volume of mercury injection. Therefore, in the process of mercury entering the pore channel, the pressure of mercury injection will continue to increase, and the point of pressure change is called the pressure change point. From the mercury entering the pore channel to the complete saturation of the pore channel, there are multiple pressure change points. Therefore, the control device can obtain the mercury injection pressure of at least two pressure change points and the mercury injection volume corresponding to the mercury injection pressure.

[0094] Step S220: determining the connectivity coefficient of the pore channel according to the first energy conservation formula, the second energy conservation formula, the mercury injection pressure and the mercury injection volume.

[0095] The connectivity coefficient of the pore channel is used to evaluate the connectivity degree of the pore channel in the form of a void structure. As shown in FIG. 1, assuming that each pore channel can be wetted by a fluid, the connectivity coefficient of the pore channel of the rock to be measured is 1, indicating that the pore channel of the rock to be measured is completely connected.

[0096] The first energy conservation formula represents that the work required for mercury to enter the pore channel is equal to the sum of the pressure work of the applied mercury injection pressure and the volume work of the mercury injection volume.

[0097] Exemplarily, the first energy conservation formula (1) is as follows: dW = PdV + VdP (1)

[0098] Wherein, dW represents the work required for mercury to enter the pore channel, P represents the mercury injection pressure, V represents the mercury injection volume, PdV represents the volume work of the mercury injection volume, the volume work represents the work done when the mercury injection volume changes, VdP represents the pressure work of the mercury injection pressure, and the pressure work represents the work done when the mercury injection pressure changes.

[0099] The second energy conservation formula represents that the work required for mercury to enter the pore channel is equal to the sum of the surface energy generated by the mercury and the pore channel and the interfacial energy generated by the mercury and the air.

[0100] Exemplarily, the second energy conservation formula (2) is as follows: dW = σd(A sl -A sv )+ σdS (2)

[0101] The difference from formula (1) is that σ represents the mercury interfacial tension, A sl represents the contact area of the mercury and the wall surface of the pore channel. A sv represents the contact area of the liquid mercury and the mercury vapor, and S represents the interface area of the mercury vapor and the liquid mercury.

[0102] Since the relationship between mercury and rock is non-wetting, formula (2) can be converted into formula (3) according to Young's equation. That is, formula (3) is the second energy conservation formula.

[0103] dW = σdA + σdS (3)

[0104] wherein A and A sl The meanings are the same, and both represent the contact area of mercury and the wall surface of the pore. The difference from formula (1) is that σdA represents the surface energy increased when the mercury spreads on the wall surface of the pore. In the embodiments of the present application, σdA is simply referred to as the surface energy generated by the mercury and the pore. σdS represents the interface energy possessed by the curved interface at the liquid mercury-mercury vapor interface. It can be understood that the curved interface refers to the interface where the mercury contacts the air, and therefore, σdS is simply referred to as the interface energy generated by the mercury and the air in the embodiments of the present application.

[0105] It can be understood that based on formula (1) and formula (3), the terms are eliminated to obtain a data model for preprocessing the mercury intrusion pressure and the mercury intrusion volume (that is, the variables of the data model are the mercury intrusion pressure and the mercury intrusion volume). Then, the mercury intrusion pressure and the mercury intrusion volume are input into the data model, and the processed mercury intrusion pressure and the mercury intrusion volume are output. The connectivity coefficient of the pore can be calculated according to the processed mercury intrusion pressure and the mercury intrusion volume.

[0106] It can be understood that in the above technical solution, when the mercury enters the pore with a void structure in the rock to be evaluated, the mercury intrusion pressure of at least two pressure change points and the mercury intrusion volume corresponding to the mercury intrusion pressure are obtained, and the connectivity coefficient of the pore of the rock to be evaluated is obtained according to the mercury intrusion pressure, the mercury intrusion volume, the first energy conservation formula and the second energy conservation formula, so as to evaluate the connectivity degree of the pore of the rock to be evaluated. The mercury intrusion pressure and the mercury intrusion volume are processed in combination with the two energy conservation formulas, so that the connectivity coefficient obtained according to the processed mercury intrusion pressure and the mercury intrusion volume is more accurate, and the evaluation of the connectivity degree of the pore of the rock to be evaluated is more accurate.

[0107] Referring to FIG. 4, in some embodiments, step S220 includes: determining the connectivity coefficient of the pore according to the first energy conservation formula, the second energy conservation formula, the mercury intrusion pressure and the mercury intrusion volume, including: constructing a slope model of the pore with a pore structure based on the first energy conservation formula and the second energy conservation formula, then inputting the mercury intrusion pressure and the mercury intrusion volume into the slope model, outputting a set of slope discrete points, and finally calculating the connectivity coefficient of the pore according to the set of slope discrete points.

[0108] Specifically, the step S410 includes the following steps:

[0109] Step S410: According to the capillary bundle equivalent assumption principle, the equivalent relationship between the surface energy and the interface energy in the second energy conservation formula is converted to obtain the change relationship between the mercury intrusion pressure and the mercury intrusion volume, and a third energy conservation formula is obtained.

[0110] The capillary bundle equivalent assumption principle is that a pore channel in a void structure is equivalent to a multi-root capillary tube with a rough surface and a variable cross-section. Reservoir rock can be regarded as a multi-dimensional interconnected capillary network, and therefore, the capillary bundle model can be used to simplify the capillary network into a model composed of capillary tubes with different sizes in parallel.

[0111] According to the capillary bundle equivalent assumption principle, after the pore channel is equivalent to a multi-root capillary tube with an approximately equal diameter, the volume sum of all the pore channels in the rock to be measured is formula (4), and the area sum of the inner surfaces of all the pore channels is formula (5) and (6).

[0112] wherein V represents the mercury injection volume, from the perspective of filling the pore channel with mercury, V can also represent the volume sum of all the pore channels, π represents the circle constant, r represents the pore channel radius corresponding to the mercury injection pressure, and the mercury injection pressure sensed by the pressure instrument can be converted into the corresponding pore channel radius. N(r) represents the number of pore channels with a radius of r, L represents the length of the pore channel, and dr represents the differential radius r.

[0113] wherein the difference from formula (4) is that A represents the contact area between the mercury and the wall surface of the pore channel, and from the perspective of filling the pore channel with mercury, A can also represent the area sum of the inner surfaces of all the pore channels. S = 2πr 2 N(r) (6)

[0114] wherein the difference from formula (4) is that S represents the area sum of the inner surfaces of all the pore channels.

[0115] It can be understood that formula (5) and (6) calculate the area sum of the inner surfaces of all the pore channels in different ways.

[0116] At the same time, the relationship between the mercury injection pressure and the pore channel radius is formula (7).

[0117] By combining formula (4), (5), (6) and (7), according to the definitions of the variables in the formula, the second energy conservation formula (2) can be converted into the third energy conservation formula (8).

[0118] It can be understood that formula (8) converts the equivalent relationship between the surface energy σdA and the interfacial energy σdS in the second energy conservation formula, and obtains the change relationship between the mercury injection pressure and the mercury injection volume.

[0119] Step S420: based on the third energy conservation formula and the first energy conservation formula, a hyperbolic model of the pore channel is constructed.

[0120] Specifically, based on the energy conservation principle, since the required work dW done by the mercury in the first energy conservation formula (1) and the second energy conservation formula (3) is consistent, the required work dW and the mercury interfacial tension σ can be eliminated to obtain formula (9):

[0121] It can be seen from formula (9) that the change rate of the mercury vapor and liquid mercury interfacial area S with the contact area A of the mercury and the wall of the pore is equal to the tangent slope of the log-log coordinate pressure mercury curve.

[0122] Based on formula (9) and the third energy conservation formula (8), a hyperbolic model (10) with pressure mercury pressure and mercury injection volume as variables can be obtained.

[0123] It can be understood that the hyperbolic model (10) includes a formula that can represent the pressure mercury curve and a formula that can represent the correlation curve. The pressure mercury curve represents the change relationship between the pressure mercury pressure P and the mercury injection volume V. The correlation curve represents the change relationship between the pressure mercury pressure P and the mercury injection volume V converted from the surface energy σdA and the interfacial energy σdS.

[0124] Step S430: obtaining a slope model based on the logarithm of the hyperbolic model.

[0125] It can be understood that the hyperbolic model includes the pressure mercury curve and the correlation curve. Since the pressure mercury curve can represent the change relationship between the pressure mercury pressure P and the mercury injection volume V, and the correlation curve represents the change relationship between the pressure mercury pressure P and the mercury injection volume V converted from the surface energy σdA and the interfacial energy σdS. Therefore, the change trend of the pressure mercury curve and the correlation curve should be approximately the same. In order to facilitate the analysis of the change trend, the hyperbolic model (10) is taken logarithm to obtain formula (11).

[0126] The formula (11) can be simplified to obtain formula (12). log k2=-log k1+C(12)

[0127] Wherein, k2 represents the first tangent slope of the point on the log-log pressure mercury curve, k1 represents the second tangent slope of the point on the semi-logarithmic curve of the inverse of the pressure mercury pressure and the mercury injection volume, and it can be understood that the semi-logarithmic curve of the inverse of the pressure mercury pressure and the mercury injection volume also refers to the correlation curve of the present application. C represents a constant.

[0128] Based on formula (11) and (12), a slope model can be constructed, which includes a first slope model (13) and a second slope model (14).

[0129] Wherein, k2i represents the slope of the tangent line corresponding to the i-th pressure change point on the double logarithmic intrusion mercury curve, P i+1 represents the intrusion mercury pressure corresponding to the i+1-th pressure change point, P i-1 represents the intrusion mercury pressure corresponding to the i-1-th pressure change point, V i+1 represents the intrusion mercury volume corresponding to the i+1-th pressure change point, P i-1 represents the intrusion mercury volume corresponding to the i-1-th pressure change point.

[0130] wherein the difference from formula (13) is that k 1i represents the slope of the tangent line corresponding to the i-th pressure change point on the semi-logarithmic curve of the reciprocal of the intrusion mercury pressure versus the intrusion mercury volume (i.e., the correlation curve).

[0131] It can be understood that, in an ideal state, when the connectivity of the pore structure-shaped pore channels is 100%, the equation relationship of formula (12) is absolutely established, that is, the slopes of the intrusion mercury curve and the correlation curve are in a linear correlation relationship. In the actual pore structure of the rock, the connectivity is not 100%, and there is a certain deviation from the ideal state. The correlation of the two curves will also change accordingly. Therefore, the actually obtained intrusion mercury pressure and intrusion mercury volume are input into the slope model (including the first slope model (13) and the second slope model (14)), preprocessed, and the preprocessed intrusion mercury pressure and intrusion mercury volume are output. The connectivity coefficient obtained based on the preprocessed intrusion mercury pressure and intrusion mercury volume can quantitatively represent the connectivity (connectivity degree) of the pore structure-shaped pore channels.

[0132] Step S440: inputting the intrusion mercury pressure and the intrusion mercury volume of at least two pressure change points into the slope model, and outputting a slope discrete point set.

[0133] The intrusion mercury pressure of at least two pressure change points, and the intrusion mercury volume corresponding to the intrusion mercury pressure are input into the first tangent slope model (13), and a first discrete point set {(V i ,p i )} is output. That is, k 2i ={(V i ,p i )}.

[0134] The intrusion mercury pressure of at least two pressure change points, and the intrusion mercury volume corresponding to the intrusion mercury pressure are input into the second tangent slope model (14), and a second discrete point set {(V i ,p i )} is output. That is, k 1i ={(V i ,p i )}.

[0135] It can be understood that the first discrete point set {(Vi i} is a set of mercury injection pressures preprocessed by the first tangent slope model and a set of mercury injection volumes corresponding to the mercury injection pressures, and reflects the relationship between the mercury injection pressures and the mercury injection volumes in the first tangent slope model. The second discrete point set {(V i i} is a set of mercury injection pressures preprocessed by the second tangent slope model and a set of mercury injection volumes corresponding to the mercury injection pressures, and reflects the relationship between the mercury injection pressures and the mercury injection volumes in the second tangent slope model.

[0136] Step S450: calculating the connectivity coefficient of the pore channel according to the slope discrete point set.

[0137] Specifically, the connectivity coefficient of the pore channel is obtained according to the first discrete point set and the second discrete point set. Exemplarily, the determination formula (15) of the connectivity coefficient is as follows:

[0138] wherein d represents the connectivity coefficient, k 2i represents the second discrete point set, k 2i represents the first discrete point set. represents the average value of the tangent slopes of all points on the double logarithmic mercury injection curve, that is, represents the average value of the first discrete point set. represents the average value of the tangent slopes of all points on the semi-logarithmic curve of the mercury injection pressure reciprocal and the mercury injection volume, that is, represents the average value of the second discrete point set. It can be understood that the formula (15) represents that the absolute value of the covariance between the first discrete point set and the second discrete point set is the connectivity coefficient of the pore channel with a pore structure.

[0139] In some embodiments, the connectivity coefficient of the pore channel is within a preset range, and the preset range is converted according to a preset J function curve.

[0140] After obtaining the connectivity coefficient, the preset range can be converted by the preset J function curve, and the slope model is optimized according to the preset range, so that the connectivity coefficient output by the optimized slope model is within the preset range. Exemplarily, according to the preset J function curve, a series of standard slopes are calculated, and a standard connectivity coefficient is calculated according to the standard slope. Taking the standard connectivity coefficient as a reference, a range is selected, and the range is the preset range, for example, the preset range is 0.8 to 1.

[0141] The connectivity coefficient obtained by the connectivity evaluation method of the embodiments of the present application is verified through the following experimental embodiments.

[0142] ​​Preparation of the experiment: 14 rocks to be evaluated are prepared. The 14 rocks to be evaluated are from the same formation of a non-homogeneous carbonate reservoir in an oilfield. The oil reservoir has complex geological conditions, various sedimentary and diagenetic processes, and strong heterogeneity.

[0143] The well names of the 14 rocks to be evaluated are W34, W110, W67, W60, W40, W104, W13, W31, W96, W18, W93, W82 and W50 respectively. The first value of each rock to be evaluated obtained by converting the porosity is approximately equal to the second value obtained by converting the permeability (the difference between the first value and the second value can be pre-set).

[0144] Experimental method: the control device executes the following steps:

[0145] a1: control the mercury injection device to press mercury into each rock to be evaluated.

[0146] a2: obtain the mercury injection pressure and the mercury injection volume of each rock to be evaluated at at least two pressure change points.

[0147] a3: input the mercury injection pressure and the mercury injection volume of step a2 into the slope model, and output the slope discrete point set.

[0148] a4: according to the slope discrete point set of step a3, the connectivity coefficient of each rock to be evaluated is calculated, and the connectivity coefficients of the 14 rocks to be evaluated are obtained.

[0149] The connectivity coefficients of the 14 rocks to be evaluated are shown in Table 1 below.

[0150] Table 1. Connectivity coefficient result table

[0151] It can be understood that in the case of approximate porosity and permeability, the closer the connectivity coefficient is to 1, the more accurate the evaluation of the connectivity degree. The connectivity coefficients of the 14 rocks to be evaluated in Table 1 are close to 1, which indicates that the connectivity evaluation method of the present application is accurate in evaluating the connectivity degree of the pore structure-like channels in the rock to be evaluated.

[0152] Then, the control device further executes the following steps:

[0153] b1: obtain the porosity and permeability of each rock to be evaluated.

[0154] b2: obtain the mercury injection curve according to the porosity and permeability of step b1.

[0155] b3: perform dimensionless processing on the mercury injection curve obtained in step b2 to obtain a processed curve, and obtain 14 processed curves.

[0156] Fig. 5 shows the 14 processed curves. It can be understood that the 14 processed curves in Fig. 5 are close to the preset J function curve. It is shown that the connectivity evaluation method of the embodiment of the present application is accurate in evaluating the connectivity degree of the pore structure-shaped pore channel in the rock to be evaluated.

[0157] Comparative experiment: According to the traditional evaluation method, the mercury injection pressure and the water saturation of the 14 rocks to be evaluated are obtained, and the mercury injection pressure and the water saturation are shown in Table 2. The mercury injection curve is drawn according to Table 2, and the drawing result is shown in Fig. 6.

[0158] Table 2. Mercury injection pressure and water saturation result table

[0159] It can be understood that the mercury injection curves obtained by the traditional evaluation method are quite different and have various curve shapes for the 14 rocks to be evaluated, and it is difficult to determine the connectivity of the pore structure-shaped pore channel by the conventional curve determination method.

[0160] It can be understood that the connectivity evaluation method of the embodiment of the present application processes the obtained mercury injection pressure and mercury injection volume by the slope model, and calculates the connectivity coefficient of the pore structure-shaped pore channel of the rock according to the processed mercury injection pressure and mercury injection volume. The type of obtained parameters is less, the calculation of the connectivity coefficient is convenient, and the connectivity coefficient has high accuracy.

[0161] The following is an apparatus embodiment of the present application, which can be used to execute the method embodiments of the present application. For details not disclosed in the apparatus embodiments of the present application, please refer to the method embodiments of the present application.

[0162] Fig. 7 is a structural schematic diagram of a rock pore channel connectivity evaluation apparatus embodiment one provided by the embodiment of the present application; the apparatus 30 can be integrated in the control apparatus 10 in the above-mentioned method embodiments, or can be realized by the control apparatus 10 in the above-mentioned method embodiments. As shown in Fig. 7, the connectivity evaluation apparatus 30 comprises an acquisition module 31 and a calculation module 32.

[0163] The acquisition module 31 is used to acquire the mercury entering the pore structure-shaped pore channel in the rock to be evaluated, the mercury injection pressure at at least two pressure change points, and the mercury injection volume corresponding to the mercury injection pressure.

[0164] The calculation module 32 is used to determine the connectivity coefficient of the pore channel according to the first energy conservation formula, the second energy conservation formula, the mercury injection pressure and the mercury injection volume, and the connectivity coefficient is used to evaluate the connectivity degree of the pore structure-shaped pore channel.

[0165] In some embodiments, the computing module 32 is specifically configured to determine the connectivity coefficient of the pore channel according to the first energy conservation formula, the second energy conservation formula, the mercury injection pressure and the mercury injection volume, comprising:

[0166] constructing a slope model of the pore channel in the form of a pore structure based on the first energy conservation formula and the second energy conservation formula; inputting the mercury injection pressure and the mercury injection volume into the slope model to output a set of slope discrete points; and calculating the connectivity coefficient of the pore channel according to the set of slope discrete points.

[0167] In some embodiments, the computing module 32 is specifically configured to construct a hyperbolic model of the pore channel in the form of a pore structure based on the first energy conservation formula and the second energy conservation formula, and calculate the slope model based on the logarithm of the hyperbolic model.

[0168] In some embodiments, the computing module 32 is specifically configured to convert the equivalent relationship between the surface energy and the interfacial energy in the second energy conservation formula according to the capillary bundle equivalent assumption principle to obtain a change relationship between the mercury injection pressure and the mercury injection volume, and obtain a third energy conservation formula; and construct a hyperbolic model of the pore channel in the form of a pore structure based on the third energy conservation formula and the first energy conservation formula.

[0169] In some embodiments, the computing module 32 is specifically configured to input the mercury injection pressure and the mercury injection volume into the first tangent slope model to obtain a first set of discrete points; input the mercury injection pressure and the mercury injection volume into the second tangent slope model to obtain a second set of discrete points; and obtain the connectivity coefficient of the pore channel according to the first set of discrete points and the second set of discrete points.

[0170] In some embodiments, the computing module 32 is specifically configured to determine the absolute value of the covariance between the first set of discrete points and the second set of discrete points as the connectivity coefficient of the pore channel.

[0171] The connectivity evaluation device 30 provided in the embodiment is used to execute the technical solutions in any of the method embodiments, and has similar implementation principles and technical effects, which will not be described here.

[0172] FIG. 8 is a structural schematic diagram of a server provided in an embodiment of the present application. As shown in FIG. 8, the server 40 comprises:

[0173] a processor 41, a memory 42, and a communication interface 43;

[0174] The memory 42 is used to store executable instructions of the processor 41.

[0175] The processor 41 executes the computer executable instructions stored in the memory 42 to perform the following steps:

[0176] The mercury entering the pore structure of the rock to be evaluated is obtained, and the mercury injection pressure and the mercury injection volume corresponding to the mercury injection pressure at at least two pressure change points are obtained;

[0177] According to the first energy conservation formula, the second energy conservation formula, the mercury injection pressure and the mercury injection volume, the connectivity coefficient of the pore channel is determined, and the connectivity coefficient is used to evaluate the connectivity degree of the pore channel in the pore structure;

[0178] The first energy conservation formula represents that the work required for the mercury to enter the pore channel is equal to the sum of the pressure work of the mercury injection pressure applied to the mercury and the volume work of the mercury injection volume, and the second energy conservation formula represents that the work required for the mercury to enter the pore channel is equal to the sum of the surface energy generated by the mercury and the pore channel and the interface energy generated by the mercury and the air.

[0179] In some embodiments, the processor executing the computer-executable instructions stored in the memory is further configured to perform the following steps:

[0180] Based on the first energy conservation formula and the second energy conservation formula, a slope model of the pore channel in the pore structure is constructed;

[0181] The mercury injection pressure and the mercury injection volume are input into the slope model, and a set of slope discrete points is output;

[0182] According to the set of slope discrete points, the connectivity coefficient of the pore channel is calculated.

[0183] In some embodiments, the processor executing the computer-executable instructions stored in the memory is further configured to perform the following steps:

[0184] Based on the first energy conservation formula and the second energy conservation formula, a hyperbolic model of the pore channel in the pore structure is constructed, the hyperbolic model includes a mercury injection pressure curve and an associated curve, the mercury injection pressure curve represents the change relationship between the mercury injection pressure and the mercury injection volume, and the associated curve represents the change relationship between the mercury injection pressure and the mercury injection volume converted from the surface energy and the interface energy;

[0185] Based on the hyperbolic model, a slope model is calculated.

[0186] In some embodiments, the processor executing the computer-executable instructions stored in the memory is further configured to perform the following steps:

[0187] According to the capillary bundle equivalent assumption principle, the equivalent relationship between the surface energy and the interface energy in the second energy conservation formula is converted to obtain the change relationship between the mercury injection pressure and the mercury injection volume, and a third energy conservation formula is obtained;

[0188] Based on the third energy conservation formula and the first energy conservation formula, a hyperbolic model of the pore channel in the pore structure is constructed.

[0189] In some embodiments, the slope model comprises a first tangent slope model of the mercury intrusion curve, and a second tangent slope model of the associated curve;

[0190] The processor executing the computer-executable instructions stored in the memory is further configured to perform the following steps:

[0191] The mercury intrusion pressure and the intrusion volume are input into the first tangent slope model to obtain a first discrete point set;

[0192] The mercury intrusion pressure and the intrusion volume are input into the second tangent slope model to obtain a second discrete point set;

[0193] The connectivity coefficient of the pore channel is calculated according to the slope discrete point set, comprising:

[0194] The connectivity coefficient of the pore channel is obtained according to the first discrete point set and the second discrete point set.

[0195] In some embodiments, the processor executing the computer-executable instructions stored in the memory is further configured to perform the following steps:

[0196] The absolute value of the covariance between the first discrete point set and the second discrete point set is determined as the connectivity coefficient of the pore channel.

[0197] In some embodiments, the connectivity coefficient of the pore channel is within a preset range, and the preset range is converted according to a preset J function curve.

[0198] Optionally, the memory 42 can be independent or integrated with the processor 41.

[0199] Optionally, when the memory 42 is independent of the processor 41, the server 40 can further comprise:

[0200] The bus 44, the memory 42 and the communication interface 43 are connected with the processor 41 through the bus 44 and complete the communication between each other, and the communication interface 43 is used for communicating with other devices.

[0201] Optionally, the communication interface 43 can be realized by a transceiver. The communication interface is used for realizing the communication between the database access device and other devices (such as a client, a read-write library and a read-only library). The memory can contain a random access memory (RAM) and can also include a non-volatile memory, such as at least one disk memory.

[0202] The bus 44 can be a peripheral component interconnect (PCI) bus or an extended industry standard architecture (EISA) bus, etc. The bus can be divided into an address bus, a data bus, a control bus, etc. For ease of representation, only one thick line is shown in the figure, but it does not mean that there is only one bus or only one type of bus.

[0203] The processor described above can be a general processor, including a central processing unit (CPU), a network processor (NP), etc.; can also be a digital signal processor (DSP), an application specific integrated circuit (ASIC), a field programmable gate array (FPGA) or other programmable logic device, a discrete gate or transistor logic device, a discrete hardware component.

[0204] The server is configured to implement the technical solutions of the control device in any of the preceding method embodiments, and has similar implementation principles and technical effects, which will not be described here.

[0205] The embodiments of the present application also provide a readable storage medium having a computer program stored thereon, and the computer program is configured to implement the technical solutions provided by any of the preceding embodiments when executed by a processor.

[0206] The embodiments of the present application also provide a computer program product comprising a computer program, and the computer program is configured to implement the technical solutions provided by any of the preceding method embodiments when executed by a processor.

[0207] Those skilled in the art can understand that all or part of the steps of the above-mentioned method embodiments can be completed by program instruction related hardware. The above-mentioned program can be stored in a computer readable storage medium. When the program is executed, the steps of the above-mentioned method embodiments are executed; and the above-mentioned storage medium includes ROM, RAM, magnetic disc or optical disc and various storage medium that can store program codes.

[0208] Finally, it should be pointed out that: the above embodiments are only used to illustrate the technical solutions of the present application, and not to limit them; although the present application has been described in detail with reference to the above embodiments, those skilled in the art should understand that: it can still modify the technical solutions recorded in the above embodiments, or make equivalent replacement for part or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A method of evaluating the connectivity of a rock formation, characterized by, The method comprises: acquiring mercury entering a pore structure-shaped pore channel in a rock to be evaluated, a mercury injection pressure at at least two pressure change points, and a mercury injection volume corresponding to the mercury injection pressure; determining a connectivity coefficient of the pore channel according to a first energy conservation formula, a second energy conservation formula, the mercury injection pressure, and the mercury injection volume, the connectivity coefficient being used to evaluate a connectivity degree of the pore structure-shaped pore channel; wherein the first energy conservation formula indicates that work required for mercury to enter a pore channel is equal to a sum of pressure work of a mercury injection pressure applied to the mercury and volume work of a mercury injection volume of the mercury, and the second energy conservation formula indicates that the work required for the mercury to enter the pore channel is equal to a sum of surface energy generated by the mercury and the pore channel and interface energy generated by the mercury and air.

2. The method of claim 1, wherein, The determining of the connectivity coefficient of the pore channel according to the first energy conservation formula, the second energy conservation formula, the mercury injection pressure, and the mercury injection volume comprises: constructing a slope model of the pore structure-shaped pore channel based on the first energy conservation formula and the second energy conservation formula; inputting the mercury injection pressure and the mercury injection volume into the slope model to output a slope discrete point set; calculating the connectivity coefficient of the pore channel according to the slope discrete point set.

3. The method of claim 2, wherein, The constructing of the slope model of the pore structure-shaped pore channel based on the first energy conservation formula and the second energy conservation formula comprises: constructing a hyperbolic curve model of the pore structure-shaped pore channel based on the first energy conservation formula and the second energy conservation formula, the hyperbolic curve model comprising a mercury injection pressure curve and a correlation curve, the mercury injection pressure curve representing a change relationship between the mercury injection pressure and the mercury injection volume, and the correlation curve representing a change relationship between the mercury injection pressure and the mercury injection volume converted according to the surface energy and the interface energy; calculating the slope model based on the hyperbolic curve model.

4. The method of claim 3, wherein, The constructing of the hyperbolic curve model of the pore structure-shaped pore channel based on the first energy conservation formula and the second energy conservation formula comprises: converting an equivalent relationship between the surface energy and the interface energy in the second energy conservation formula according to a capillary bundle equivalent assumption principle to obtain a third energy conservation formula representing a change relationship between the mercury injection pressure and the mercury injection volume; constructing the hyperbolic curve model of the pore structure-shaped pore channel based on the third energy conservation formula and the first energy conservation formula.

5. The method of claim 3, wherein, The slope model comprises a first tangent slope model of the mercury injection pressure curve and a second tangent slope model of the correlation curve. The inputting of the mercury injection pressure and the mercury injection volume into the slope model to output the slope discrete point set comprises: inputting the mercury injection pressure and the mercury injection volume into the first tangent slope model to obtain a first discrete point set; inputting the mercury injection pressure and the mercury injection volume into the second tangent slope model to obtain a second discrete point set; The calculating of the connectivity coefficient of the pore channel according to the slope discrete point set comprises: obtaining the connectivity coefficient of the pore channel according to the first discrete point set and the second discrete point set.

6. The method of claim 5, wherein, The obtaining of the connectivity coefficient of the pore channel according to the first discrete point set and the second discrete point set comprises: An absolute value of a covariance between the first discrete point set and the second discrete point set is determined as the connectivity coefficient of the pore channel.

7. The method according to any one of claims 1 to 6, characterized in that, The connectivity coefficient of the pore channel is within a preset range, and the preset range is converted according to a preset J function curve.

8. An apparatus for evaluating the connectivity of a rock formation, characterized by The connectivity evaluation device comprises, The acquisition module is configured to acquire mercury entering a pore channel in the form of a pore structure in the rock to be evaluated, pressure mercury pressures at at least two pressure change points, and mercury entering volumes corresponding to the pressure mercury pressures; The calculation module is configured to determine, according to a first energy conservation formula, a second energy conservation formula, the pressure mercury pressures, and the mercury entering volumes, a connectivity coefficient of the pore channel, the connectivity coefficient being used to evaluate a connectivity degree of the pore channel in the form of the pore structure. The first energy conservation formula indicates that work required for the mercury to enter the pore channel is equal to a sum of pressure work of the pressure mercury pressures applied to the mercury and volume work of the mercury entering volumes, and the second energy conservation formula indicates that the work required for the mercury to enter the pore channel is equal to a sum of surface energy generated by the mercury and the pore channel and interface energy generated by the mercury and air.

9. The apparatus of claim 8, wherein, The calculation module is specifically configured to: construct a slope model of the pore channel in the form of the pore structure based on the first energy conservation formula and the second energy conservation formula; input the pressure mercury pressures and the mercury entering volumes into the slope model, and output a slope discrete point set; determine, according to the slope discrete point set, the connectivity coefficient of the pore channel.

10. The apparatus of claim 9, wherein, The calculation module is specifically configured to: construct a hyperbolic curve model of the pore channel in the form of the pore structure based on the first energy conservation formula and the second energy conservation formula, the hyperbolic curve model including a mercury intrusion curve and a correlation curve, the mercury intrusion curve representing a change relationship between the pressure mercury pressures and the mercury entering volumes, and the correlation curve representing a change relationship between the pressure mercury pressures and the mercury entering volumes converted according to the surface energy and the interface energy; determine, based on the hyperbolic curve model, the slope model.

11. The apparatus of claim 10, wherein, The calculation module is specifically configured to: convert, according to a capillary bundle equivalent assumption principle, an equivalent relationship between the surface energy and the interface energy in the second energy conservation formula, to obtain a third energy conservation formula representing a change relationship between the pressure mercury pressures and the mercury entering volumes; construct the hyperbolic curve model of the pore channel in the form of the pore structure based on the third energy conservation formula and the first energy conservation formula.

12. The apparatus of claim 10, wherein, The slope model includes a first tangent slope model of the mercury intrusion curve and a second tangent slope model of the correlation curve; and the calculation module is specifically configured to: input the pressure mercury pressures and the mercury entering volumes into the first tangent slope model to obtain a first discrete point set; input the pressure mercury pressures and the mercury entering volumes into the second tangent slope model to obtain a second discrete point set; The calculation module determines, according to the slope discrete point set, the connectivity coefficient of the pore channel, and is specifically configured to: determine, according to the first discrete point set and the second discrete point set, the connectivity coefficient of the pore channel.

13. The apparatus of claim 12, wherein, The calculation module is specifically configured to: determine, according to the first discrete point set and the second discrete point set, an absolute value of a covariance between the first discrete point set and the second discrete point set as the connectivity coefficient of the pore channel.

14. The apparatus of claims 8-13, wherein, The communication coefficient of the pore channel is in a preset range, and the preset range is converted according to a preset J function curve.

15. A server, characterized by Comprise: Processor, memory, communication interface; The memory is used for storing executable instructions of the processor; Wherein, the processor executes the computer execution instructions stored in the memory to execute the following steps: Obtain the mercury entering the pore structure of the rock to be evaluated, the pressure mercury pressure at at least two pressure change points, and the mercury volume corresponding to the pressure mercury pressure; According to the first energy conservation formula, the second energy conservation formula, the pressure mercury pressure and the mercury volume, the communication coefficient of the pore channel is determined, and the communication coefficient is used for evaluating the communication degree of the pore structure of the pore channel; Wherein, the first energy conservation formula represents that the work required for mercury to enter the pore channel is equal to the sum of the pressure work of the pressure mercury pressure applied to the mercury and the volume work of the mercury volume, and the second energy conservation formula represents that the work required for mercury to enter the pore channel is equal to the sum of the surface energy generated by the mercury and the pore channel and the interface energy generated by the mercury and air.

16. The server of claim 15, wherein, The processor executing the computer execution instructions stored in the memory is also used to execute the following steps: Based on the first energy conservation formula and the second energy conservation formula, a slope model of the pore structure of the pore channel is constructed; The pressure mercury pressure and the mercury volume are input into the slope model, and a slope discrete point set is output; According to the slope discrete point set, the communication coefficient of the pore channel is calculated.

17. The server of claim 16, wherein, The processor executing the computer execution instructions stored in the memory is also used to execute the following steps: Based on the first energy conservation formula and the second energy conservation formula, a hyperbolic curve model of the pore structure of the pore channel is constructed, the hyperbolic curve model includes a pressure mercury curve and an associated curve, the pressure mercury curve represents the change relationship between the pressure mercury pressure and the mercury volume, and the associated curve represents the change relationship between the pressure mercury pressure and the mercury volume converted according to the surface energy and the interface energy; Based on the hyperbolic curve model, the slope model is calculated.

18. The server of claim 17, wherein, The processor executing the computer execution instructions stored in the memory is also used to execute the following steps: According to the capillary bundle equivalent assumption principle, the equivalent relationship between the surface energy and the interface energy in the second energy conservation formula is converted to obtain the change relationship between the pressure mercury pressure and the mercury volume, and a third energy conservation formula is obtained; Based on the third energy conservation formula and the first energy conservation formula, the hyperbolic curve model of the pore structure of the pore channel is constructed.

19. The server of claim 17, wherein, The slope model includes a first tangent slope model of the pressure mercury curve and a second tangent slope model of the associated curve; The processor executing the computer execution instructions stored in the memory is also used to execute the following steps: The pressure mercury pressure and the mercury volume are input into the first tangent slope model to obtain a first discrete point set; The pressure mercury pressure and the mercury volume are input into the second tangent slope model to obtain a second discrete point set; The communication coefficient of the pore channel calculated according to the slope discrete point set comprises: The connectivity coefficient of the pore channel is obtained according to the first discrete point set and the second discrete point set.

20. The server of claim 19, wherein, The processor further executes the computer-executed instructions stored in the memory to perform the following steps: The absolute value of the covariance between the first discrete point set and the second discrete point set is determined as the connectivity coefficient of the pore channel.

21. The server of claims 15-20, wherein, The connectivity coefficient of the pore channel is within a preset range, and the preset range is converted according to a preset J function curve.

22. A readable storage medium, having stored thereon a computer program, characterized in that, The computer program is executed by the processor to implement the method in any one of claims 1 to 7.

23. A computer program, characterized in that, The computer program comprises a computer program stored in a computer readable storage medium, at least one processor can read the computer program from the computer readable storage medium, and the at least one processor executes the computer program to implement the method in any one of claims 1 to 7.

Images