Calculation method and device for capillary force and relative permeability curve

[0085] Next, the technical scheme in the examples in this paper will be described in the following examples, and it is understood that the described embodiments are merely the embodiments of this Welcome, not all of the embodiments. Based on the embodiments herein, one of ordinary skill in the art does not have all other embodiments obtained under the pre-creative labor premise, it belongs to the scope of this paper.

[0086] The capillary force curve refers to a curve that changes the relationship between capillary force and the aqueous saturation of the core, relative to the permeability curve refers to the relative penetration rate of the oil phase and the aqueous saturation of the oil phase during the oil or water drift. Curve of change relationship, as well as a curve between the relative permeability of the aqueous phase and the aqueous saturation of the core.

[0087] Due to the multi-phase flow of the oil and gas reservoir core, the multiphase flow of water is an important issue in the field of oil and gas field development, which is mainly characterized by capillary and relative penetration curves in macro. In order to obtain capillary and relative permeability curve, the prior art mainly passed the accurate experiment of indoor core, and the experimental data was analyzed, and the capillary forces and relative permeability curves of the core were obtained.

[0088] However, due to the precision experiment of indoor core, the core of different penetration rates is considered for a few days or even months. This will greatly consume the energy of the researchers and lead to a decrease in research efficiency. Therefore, there is now an urgent need for a capillary force and a relative permeability curve to improve the efficiency of obtaining the core capillary forces and relative permeability curves.

[0089] In order to solve the above problems, the embodiments provide a method of calculating a capillary force and a relative permeability curve. figure 1 It is a schematic view showing a method of calculating a capillary force and a relative permeability curve provided herein, and the present specification provides a method of operation as described in Examples or Flowcharts, but based on conventional or uncustomary labor, it can include more Multiple or fewer operational steps. The sequence of steps listed in the examples is only one way in which a number of steps is performed, and does not mean unique execution order. When executed in the actual system or device product, it can be performed in parallel or parallel.

[0090] Refer figure 1 , A method of calculating capillary and relative penetration curve, including:

[0091] S101: Determining the cross section of the pore and the throat is a four-pointed core pore network model;

[0092] S102: Single-phase flow simulation calculation of the rock core pore network model, determine the absolute permeability of the core;

[0093]S103: The replacement simulation of the rock core pore network model is sequentially performed, respectively, determine the change in capillary force and core-saturation of the capillary force and the core-containing compaction of each of the alternative simulation calculations, and the effective penetration rate of the oil and the aqueous moisture Change relationship between saturation, changes in water phase effective penetration and moisture saturation of rock core;

[0094] S104: Calculation of the oil phase effective penetration and mandarin-containing relationship between the effects of the alternative simulation and the moisture saturation of the aqueous phase effective penetration and the moisture saturation of the core, respectively, the absolute permeability of the rock heart, respectively. In combination, it is determined that the relative permeability of the alternative simulation calculation of the oil phase and the variation of the degree of water saturation of the core, the relative penetration rate of the aqueous phase and the change of moisture saturation of the core.

[0095] Compared to the prior art, the indoor core precision experiment is performed, by establishing the core pore network model, and the model is used for multiple alternative simulation calculations, the method of obtaining the capillary force curve and the relative permeability curve is guaranteed to be accurate. The base consumption is short, and the research efficiency is high.

[0096] Capillary relative penetration curve refers to a curve that changes between capillary forces and aqueous saturation of rock core. The capillary force is directed in a direction towards the liquid concave surface, and its size is proportional to the surface tension of the liquid, which is inversely proportional to the radius of the tube, and is often shown in the formation capillary pore, which is often shown as two-phase non-miscible liquids (such as oil and water) bending interface. Pressure difference.

[0097] The relative permeability curve refers to a curve between the relative permeability of the oil phase and the aqueous saturation of the permeability of the oil phase, and the relative permeability of the aqueous phase and the aqueous saturation of the aqueous phase. Curve of the relationship. The relative penetration rate is the ratio of the effective penetration rate and absolute permeability of the absolute permeability when the multiphase fluid coexists.

[0098] Refer figure 2 For the core, there are several pores in the core, and the throat between the pore and the pores, in order to improve the accuracy of the simulation calculation, the cross section of the pore and the throat is four-star star, thereby determining the core Pore ​​network model. For the core pore network model, most of the positive six-faced body can be provided with an end face of the positive hexahedral to the end face, and the other end face opposite the entrance end surface is an exit end surface.

[0099] Further determine the number, central coordinates, volume, internal cut circular radius, cross-sectional area, length, and pores of the inlet pores or outlet pores (in the port on the end surface of the inlet (in the inlet end surface). The pores on the end face are the outlet pores). And determine the number, volume, inner cut circular radius, cross-sectional area, length, and the pore number communication with it. Further, using the chart discussion to determine the communication between the pores and the throat, remove the pores in the core pore network model, calculate the half angle of any four angle-shaped pores or throat:

[0100]

[0101] Among them, γ is a tetraight star pore or a throat half-angle, a pore or a cross-sectional area of ​​the throat, and R is a pore or an inner cutting radius of the throat.

[0102] The shape factor that can be calculated can be calculated:

[0103]

[0104] Where G is a shape factor of a pore or a throat cross section.

[0105] And the circumference of any hole or throat section:

[0106]

[0107] Where P is the perimeter of the pore or the throat section.

[0108] Refer image 3 In the embodiments of this paper, the rock core pore network model is single-phase flow simulation calculation, determine the absolute permeability of the core, including:

[0109] S301: Single-phase flow simulation calculation of the rock core pore network model, determine the conductivity coefficient of each pore in the rock core pore network model, and the conduction coefficient between the pores of each other;

[0110] S302: Determine the pressure of each pore based on the conductivity between the respective pores;

[0111] S303: The absolute permeability of the rock core is determined based on the pressure of each of the apertures and its corresponding conducting coefficients.

[0112] When the simulation of single-phase flow is calculated, the single-phase fluid is water or oil, generally in water to perform analog calculation of the mound heart permeability. First calculate the conductivity coefficient of any pore or throat:

[0113]

[0114] Where T is the pore or the conductivity of the throat, L is the length of the pore or the throat.

[0115] And further determine the conduction coefficient between any two pores (pore I and the pore j):

[0116]

[0117] Where T ’ The conductivity coefficient between the pore I and the pore j, T i For the pore I conduction coefficient, T j For the conductivity of the pore J, T ij The conduction coefficient of the throat between the connecting pore I and the pore j.

[0118] The adjacent matrix can then be established by the method of the aperture and the throat, wherein the coefficient of the adjacent matrix is ​​the conductivity between any two pores. The adjacent matrix is ​​converted into a sparse matrix. Since the sparse matrix is ​​solved, the inlet end surface pressure value and the outlet end surface pressure value are set to any value, and the solution algorithm for the large-scale sparse matrix is ​​solved. After solving each The pressure of the pore.

[0119] That is to determine the traffic of the exit end:

[0120]

[0121] Where Q is the traffic of the exit end, P i For the pressure of pore I in communication with the exit end, P out Pressure for the export end surface, N C For the number of pores connected to the exit end surface.

[0122] Then, the absolute permeability is obtained by Darcy's law:

[0123]

[0124] Where K abs For the absolute permeability, μ is the viscosity of 1, the area of ​​the exit end surface of the s, L is the distance from the end face to the outlet end surface, P in The pressure of the entrance end surface.

[0125] Refer Figure 4 In this paper, the rock core pore network model is subjected to a variety of replacement simulation and simulation calculations, respectively, determine the change of capillary forces and core-containing water saturation of each alternative simulation calculation. The change relationship between equivalent permeability and moisture saturation of the core, the variation of the effective penetration rate of the aqueous phase and the moisture saturation of the core, including the following two steps:

[0126] 1 Step S401:

[0127] The simulation calculation of the oil-driven water, the simulation of the oil-driven water, and the change of the change of the oil-driven force and the moisture-saturation of the oil-driven, and the effect of effective penetration rate and the moisture saturation of the oil phase are related. Relationship between the effective permeability and the moisture saturation degree of the core;

[0128] 2 Step S402:

[0129] The simulation calculation of the rock core water-driven water drive oil is determined to determine the change of the change of the aqueous tube and the moisture saturation of the moisture and core water-saving, and the effective penetration rate of the oil is saturated with the effect of the oil phase. The change relationship of the degree, the change of the effective penetration rate of the aqueous phase and the moisture saturation of the core;

[0130] Alternatively, the simulation of the rock core water water drive oil is used to determine the simulation of the moisture core water-driven oil, the change of the change in capillary force and the moisture-saturation of the rocky, and the effective penetration rate and the core of the oil phase. Relationship between changes in water saturation, the relationship between the effective penetration rate of the aqueous phase and the aqueous saturation of the core;

[0131] Alternatively, the simulation calculation of the rigxin pore network model is mixed and wetted, and the simulation of mixed wetting rock core water drive oil is determined to calculate the change in capillary force and core-containing water saturation, and the oil phase effective penetration rate The relationship between the change in moisture saturation and the effect of water phase effective penetration and moisture saturation of the aqueous phase.

[0132] Refer Figure 5 with Image 6 , Figure 5 with Image 6 The deeper the color, the higher the saturation of the core, the higher the color, the higher the moisture-containing degree of the core. specific, Figure 5 The process shown in the oil-driven water is the completion process of the reservoir, and the reservoir is first filled with water, and the pore and the throat are filled in the initial state. First, the oil phase gradually replaces the aqueous phase, can gradually record the different capillary force corresponding to the core-containing saturation, the different oil phase effective penetration, the core-containing saturation, the different water phase effective penetration, the core of the core Water saturation. Then, if Image 6 The wettability according to the core is wet, the wet or mixed wet is determined to determine the simulation calculation of water-driven oil, oil and wet core water drive oil, or mixed wet rock, or mixed wetting rock, or mixed wet rock drive oil. Excellently replacing the oil phase, gradually record the different capillary force corresponding to the core-containing saturation, the different oil phase effective penetration, the core-containing saturation of the core, the effective penetration of different water phases, etc. .

[0133] Refer Figure 7 In this example, the simulation calculation of the oil-driven drive water is determined by the simulation of the oil-driven water, the simulation of the oil-driven, and the effect of the core-containing computation, the effective penetration rate of the oil phase. Change relationship between rock core with aquatic saturation, changes in the effective penetration rate of the aqueous phase and the moisture saturation of the core, including:

[0134] S701: Determines the throat of each throat, the first ciniry and saturation of the first mound heart under the initial state, the first capillary forces in the initial state, and the first capillary force and the first replacement step;

[0135] S702: Determination of the first phenantine in the initial state as the current first phenolic tube force;

[0136]S703: RC performs screening of all throats and pores, determines the throat of the throat of the throat of less than or equal to the current first capillary force, and all the pores communicated with the throat; the throat and the above The aqueous saturation degree of pores is utilized in the initial state, and the current first ciniry force corresponding to the current first capillary force is aqueous saturation; according to the current first core-containing water saturation, and each pore or The first oil phase conduction coefficient corresponding to the throat is respectively corresponding, and the first oil phase relative permeability and the first water phase relative permeability corresponding to the current first rock heart-containing water saturation; A capillary force is updated, and the current first capillary force and the first replacement step are determined to be updated, the current first pupping force, and determine the current first capillary force corresponding to the current first core Saturation, and the relative permeability of the first oil phase corresponding to the current first mound heart containing a relative penetration rate;

[0137] S704: The cycle ends until the current first capillary force is greater than the first capillary force in the termination state.

[0138] Refer Figure 8 Where the experimental value is actual experimental data for controlling the analog value, the simulation value is an experimental data in accordance with the calculation method of this article, which is compared, the simulation value is the same as the data trend of the experimental value, so the calculation method is used The resulting capillary forces and relative permeability curves are high. Specifically, since the pores and throat are filled with aqueous phase in the initial state, the first rock core containing aqueous saturation S w = 1, representing the aqueous saturation of the pore and the throat at this time. The oil phase enters the reservoir by the action of the hydrocarbon pressure, and the oil phase overcomes the capillary force to replace the water phase, and finally forms the reservoir, at which time the rock is water and wet. During this process, capillary is gradually increased, and the throat plays a major constraint. If the capillary force is greater than the throat of the throat, the oil phase will enter the corresponding throat in the respective throat. Running threaded threaded pressure value:

[0139]

[0140] in For the throat throat pressure value, σ is interface tension, θ r Contact angle in the process of oil-driven water, A α1 In order to wet the absence of the secondary angle area.

[0141] The calculation method of the unusual degree angle area of ​​the wet phase is as follows:

[0142]

[0143] The first capillary force in the initial state is the minimum throat pressure value in the throat pressure value obtained by the above formula (8), and the first capillary force in the termination state is through the above formula (8) The maximum threshold pressure value of all throat pressure values ​​obtained. Further, the first replacement step may be equal, or each step is not equal, and the specific form of the first replacement step is not limited. If each step is equal, the difference between the first capillary force and the first capillary force in the termination state can be taken, and the difference is performed, and the number of aliquots can be set according to the actual situation. For example, it can be equal to 10 times, 15 times, etc., the difference in the number of aliquots can be obtained is the first replacement step.

[0144] Determining the first wool in the initial state as the current first caps, then performs screening of all throat and pores, will be less than or equal to the current first capillary force, and with the throat All the pores of Tao Tong screens out, and the sipped throat and pore-containing water saturation is from the original S. w = 1 Update to:

[0145]

[0146] Where S w Water saturation for throat or pores, P c For the current first caps.

[0147] Of course, the aquatic saturation of the throat or pores that is not screened is S. w = 1, according to the filtered throat and the pore update, the water saturation of the throat and the pore-containing degree of the unsatisfactory, the current first ciniry force corresponding to the current first core containment is determined. :

[0148]

[0149] Where S w总 Current first core-containing water saturation for the current first capillary force, S w1 Water saturation after filtering the throat and pore update, V 1 The volume of the throat and the pore for screening, s w2 Water saturation of the throat and pores of unsatisfied, V 2 The volume of the throat and the pore, V is the total volume of all throat and pores.

[0150] Further, a first oil phase conduction coefficient corresponding to any pore or throat is calculated:

[0151] T o = T (1-S w ) 2 (12)

[0152] Where T o For the first oil phase conduction coefficient.

[0153] Further, the first aqueous phase conduction coefficient of any pores or throat is calculated:

[0154]

[0155] Where T w For the first aqueous phase conduction coefficient.

[0156] Finally, the relative penetration rate of the first oil corresponding to the current first mound heart is aqueous saturation:

[0157]

[0158] Among them, K ro For the first oil phase corresponding to the current first mound heart, K o The first oil phase corresponding to the current first mound is aqueous saturation.

[0159] It can be determined that the first aqueous phase relative permeability corresponding to the current first mound heart containing is:

[0160]

[0161] Among them, K rw For the first aqueous phase corresponding to the current first mound heart containing a relative penetration rate, k w The first water phase corresponding to the current first mound heart containing aqueous saturation.

[0162] Among them, the first oil phase effective penetration rate corresponding to the current first mound heart containing aqueous saturation. o Calculation method with the first water phase effective penetration rate corresponding to the current first mound heart-containing water saturation w All of the same, the first oil phase corresponding to the current first core-containing water saturation K o Calculation method as an example:

[0163] First, the first oil phase conduction coefficient T of the first oil phase corresponding to the throat is known. o , The conductivity coefficient between any two pores is obtained by the above formula (5), and T in the formula (5) i First oil phase conduction coefficient of pore I, T j First oil phase conduction coefficient of pore J, T ij The first oil phase conduction coefficient of the throat between the communication pore I and the pore j.

[0164] Then, the adjacent matrix can be established by the communication between the aperture and the throat, wherein the coefficient of the abutment matrix is ​​the first oil phase conduction coefficient between any two pores. The adjacent matrix is ​​converted into a sparse matrix. Because the sparse matrix is ​​solved, the inlet end surface pressure value and the outlet end surface pressure value are set to any value, and the solving algorithm for the large-scale sparse matrix is ​​solved, and after solving the current The pressure of the respective pores corresponding to the first mound heart containing water saturation.

[0165] Next, the first oil phase effective penetration of the current first mound heart containing water saturation is determined by the above formula (6) and (7). o. The corresponding change in formula (7) is:

[0166]

[0167] Where K o The first oil phase corresponding to the current first mound is aqueous saturation.

[0168] Similarly, the first water phase effective penetration ratio corresponding to the current first mound heart is aqueous w When the formula (5) - (7) is also modified, in which formula (7) is changed to:

[0169]

[0170] Where K w The first water phase corresponding to the current first mound heart containing aqueous saturation.

[0171] According to the same method, the current first capillary force and the first replacement step are determined to be the updated current first purse. For each updated current first caps, the current first capillary force corresponding to the current first core-containing water saturation is determined, and the first oil phase relative permeability corresponding to the current first mound heart containing water saturation and the first The aqueous phase is relatively permeable. The cycle is ended until the current first capillary force is greater than the first capillary tube in the termination state.

[0172] During this process, the current first rock is aqueous saturation as the abscissa, respectively, with the relative permeability of the current first pupping, the relative permeability of the first water phase, respectively, each horizontal vertical Several points under the coordinates, connect several points of each horizontal longitudinal coordinate, can constitute a first moratorial force curve, a first oil phase relative permeability curve, and a first aqueous phase relative permeability curve.

[0173] Refer Figure 9 In this paper, the simulation calculation of the rock core pore network model is used to determine the simulation of the moisture core water-driven oil, the change of the change of the catalysis and the variation of the core saturation of the rocky core, Change relationship between oil phase effective penetration and moisture saturation of the core, a variation of water-phase effective penetration and moisture saturation of the core, including:

[0174] S901: Determine the first card-breaking critical pressure value of each throat, the first collaborative filling critical pressure value of each hole, the second rock heart in the initial state, the second gruff in the initial state, termination The second hair tube force and the second replacement step size;

[0175] S902: Determination of the second hair stress in the initial state as the current second hair thinning tube force;

[0176]S903: RC performs screening of all throats and pores, determining the throat of the first card and breaks the critical pressure value greater than the current second capillary force, and the first collaborative filling critical pressure value is greater than the pore of the current second capillary force; Using the throat of the thief and the pore-containing water saturation of the pore, the primary state of the second core containing water saturation is determined to determine the current second capillary force corresponding to the current second mound heart containing water saturation; according to the current second The core saturation of the core is aqueous saturation, and each pore or throat respectively corresponds to the second oil phase conduction coefficient and the second aqueous phase conduction coefficient, determines the second oil phase relative permeability and second water corresponding to the current second core-containing water saturation. Phase relative penetration rate; update the current second piratory force, determine the difference between the current second capillary force to the second replacement step, and determines the current first-grade-based tube force after the update. The current second mound heart-containing saturation corresponding to the dichlor fine tube, and the relative penetration rate of the second oil phase corresponding to the current second rock core containing the relative penetration rate;

[0177] S904: The cycle ends until the current second capillary force is smaller than the second capillary force in the termination state.

[0178] Refer Figure 10 Where the experimental value is actual experimental data for controlling the analog value, the simulation value is an experimental data in accordance with the calculation method of this article, which is compared, the simulation value is the same as the data trend of the experimental value, so the calculation method is used The resulting capillary forces and relative permeability curves are high. Specifically, since the simulation of the moisture core water drive oil process is over the end of the oil-driven water, the existing pore and throat are filled, and the wet lag effect is considered. Therefore, the second core saturation of the second rock in the initial state is a first mound heart containing aqueous saturation at the end of the oil-driven water.

[0179] The first card-breaking critical pressure value when the interruption of the stream is:

[0180]

[0181] in For the first card, the critical pressure value, σ is the interface tension, θ a Contact angles during the water drive.

[0182] The first collaborative filling critical pressure value of any hole is:

[0183]

[0184] in For the first collaborative filling critical pressure value, A 1 = 0, X is a random number of 0 to 1, and n is the number of oil-containing throatings connected to the corresponding pore.

[0185] When the interrupt of the oil phase occurs, the throat will be filled by the aqueous phase, and then further determines whether or not the pores communicated with the throat is filled with aqueous phase.

[0186] It can be set to the first state in the initial state to the maximum threshold pressure value in all throat of the above formula (8), and the second hair thin tube force in the termination state is set. Further, the second replacement step may be equal, or each step is not equal, and this paper is not limited to the specific form of the second replacement step. If each step is equal, the difference between the second capillary force in the termination state and the second capillary force in the initial state can be taken, and the difference value is equal to the difference, the number of aliquots can be set according to the actual situation. For example, it can be divided into 10 times, 15 times, etc., the difference between the number of aliquots is divided is the second substation step.

[0187] Determining the second packing tube force in the initial state is the current second caps, and then performs screening of all throat and pores, and the first card-down critical pressure value is greater than the throat of the current second capillary force, and the first A collaborative filling critical pressure value is larger than the current second capillary force, is screened. The aqueous saturation of the screen-out throat and the pore is updated to S. w = 1.

[0188] According to the above formula (11), aqueous saturation after the screening of the throat and the pore is updated and the laryngeal saturation of the undulated throat and the pore is determined to determine the current second capillary force corresponding to the current second capillary force. Spend. Equation (11) w总 The current second cinnitational aqueous saturation is the current second capillary force.

[0189] Further, the second oil phase conduction coefficient corresponding to any aperture or throat is calculated by the above formula (12), and T in the formula (12) o For the second oil phase conduction coefficient.

[0190] Further, the second aqueous phase conduction coefficient of any hole or throat is calculated by the formula (13), and T in the formula (13) w For the second aqueous phase conduction coefficient.

[0191] Finally, the second oil phase relative permeability corresponding to the current second rock heart containing aqueous saturation is determined by the formula (14). Equation (14), K ro Relative penetration rate corresponding to the current second rock heart containing water saturation, k o The second oil phase corresponding to the current second rock is aqueous saturation.

[0192] The second aqueous phase relative permeability corresponding to the current second mound heart-containing water saturation is determined using the formula (15). Formula (15) k rw The relative penetration rate corresponding to the current second core-containing water saturation degree, k w The second water phase corresponding to the current second rock heart is aqueous permeability.

[0193] Among them, the second oil phase effective penetration ratio corresponding to the current second rock is aqueous saturation. o Calculation method with the second aqueous phase effective penetration rate corresponding to the current second core-containing water saturation w All of the same, the second oil phase effective penetration rate corresponding to the current second rock heart water saturation o Calculation method as an example:

[0194] First, the second oil phase conduction coefficient T of any of the caution or throat is known. o , The conductivity coefficient between any two pores is obtained by the above formula (5), and T in the formula (5) i For the second oil phase conduction coefficient of pore I, T j For the second oil phase conduction coefficient of pore J, T ij The second oil phase conduction coefficient of the throat between the communication hole I and the pore j is connected.

[0195] Then, the adjacent matrix can be established by the communication between the chart theory, with the communication between the pores and the throat, wherein the coefficient of the abutment matrix is ​​the second oil phase conduction coefficient between any two pores. The adjacent matrix is ​​converted into a sparse matrix. Because the sparse matrix is ​​solved, the inlet end surface pressure value and the outlet end surface pressure value are set to any value, and the solving algorithm for the large-scale sparse matrix is ​​solved, and after solving the current The pressure of the respective pores corresponding to the second core containing water saturation.

[0196] Next, the second oil phase effective penetration of the current second rock heart containing water saturation is determined by the above formula (6) and (7). o. Where the formula (7) corresponding to the formula (16), in the formula (16) k o The second oil phase corresponding to the current second rock is aqueous saturation.

[0197] Similarly, the second aqueous phase effective penetration rate corresponding to the current second rock core containing water saturation is calculated. w When the formula (5) - (7) is also modified, in which the equation (7) corresponding to the formula (17), the formula (17) k w The second water phase corresponding to the current second rock heart is aqueous permeability.

[0198] According to the same method, the difference between the current second capillary force and the second replacement step is determined as the updated current second hairy tube force. For each updated, the current second wool tube force is determined to determine the current second capillary force corresponding to the current second core-containing water saturation, and the second oil phase relative penetration rate corresponding to the current second core-containing water saturation and second The aqueous phase is relatively permeable. The cycle ends until the current second capillary force is less than the second capillary force in the termination state.

[0199] In this process, in the current second rock-containing, the water saturation is the abscissa, respectively, the relative penetration rate of the current second hair, the relative penetration rate, the second water phase, respectively, each horizontal vertical Several points under the coordinates, connect several points of each horizontal longitudinal coordinate can constitute a second hair tube force curve, a second oil phase relative permeability curve, and a second aqueous phase relative permeability curve.

[0200] Refer Figure 11 In this paper, the simulation calculation of the rock core pore network model is used to determine the simulation of the moisturine core water drive oil, the change of the change of capillary force and the moisture-saturation of the core, and the relationship between the level of water saturation. Change relationship between oil phase effective penetration and moisture saturation of the core, a variation of water-phase effective penetration and moisture saturation of the core, including:

[0201] S1101: Determine the first piston-filled critical pressure value of each throat, the third mound heart in the initial state, the third capsule tube force in the initial state, the third wool tube force and the third substation in the termination state Step size;

[0202] S1102: Determination of the third caps in the initial state as the current third phenolic tube force;

[0203] S1103: RC performs screening of all throats and pores, determines the throat of the first piston filling critical pressure is less than or equal to the current third capillary force, and all of the pores communicating with the throat; The throat and the aqueous saturation of the throat and the aperture are updated in the initial state, the current third capillary force corresponding to the current third capillary force is aqueous saturation; according to the current third mound heart containing water saturation, And the third oil phase conductivity coefficient corresponding to each pore or throat, determines the relative permeability of the current third mound heart-containing water saturation and the relative permeability of the third water phase; The current third capillary force is updated, and the current third capillary force is determined as the difference between the third replacement step, and the current third capillary force corresponds to the current third capillary force. The current third mound heart contains water saturation, and the relative penetration rate of the third oil phase corresponding to the current third mound heart is ambient and the relative permeability of the third water phase;

[0204] S1104: The cycle ends until the current third capillary force is smaller than the third capillary force in the termination state.

[0205] Refer Figure 12 Where the experimental value is actual experimental data for controlling the analog value, the simulation value is an experimental data in accordance with the calculation method of this article, which is compared, the simulation value is the same as the data trend of the experimental value, so the calculation method is used The resulting capillary forces and relative permeability curves are high. Specifically, for the simulation of the hydrophilic driven oil driver oil process, the third mound heart in the initial state is aqueous saturation at the end of the oil-driven water.

[0206] when Hour,

[0207] The first piston filling critical pressure value of either the Throat is:

[0208]

[0209] in For the first piston filling critical pressure value, σ is interface tension, θ a Contact angles during the water drive.

[0210]

[0211] Among them α2 In order to wet the absence of the secondary angle area.

[0212] when When the pore or throat occurs, there is no oil film, the water phase is completely filled, and the corresponding first piston filling critical pressure value is:

[0213]

[0214] The third capillary force in the initial state is 0, the third capillary force in the termination state is the minimum value of the first piston filling critical pressure value of any of the formula (20) and the formula (22).

[0215] Further, the third replacement step may be equal, or each step is not equal, and this paper does not limit the specific form of the third replacement step. If each step is equal, the difference between the third capillary force in the termination state and the third capillary force in the initial state can be taken, and the difference value is equal to the difference, and the number of aliquots can be set according to the actual situation. For example, it can be equal to 10 times, 15 times, etc., the difference between the number of aliquots is the third replacement step.

[0216] Determining the third caps of the initial state as the current third wool tube force, then circulating to perform screening of all throat and pores, pulling the first piston filling critical pressure value of less than or equal to the current third capillary force And all the pores communicated with the throat are screened. The aqueous saturation of the screen-out throat and the pore is updated to S. w = 1.

[0217] According to the above formula (11), aqueous saturation of the throat and the pore update and aqueous saturation of the throat and the pore, the current third capillary force corresponding to the current third capillary force is determined to determine the current third capillary solution. Spend. Equation (11) w总 The current third capillary is aqueous saturation for the current third capillary force.

[0218] Further, a third oil phase conduction coefficient corresponding to any aperture or throat is calculated by the above formula (12), in the formula (12) o For the third oil phase conduction coefficient.

[0219] Further, the third aqueous phase conduction coefficient of any hole or throat is calculated using the formula (13), and T in the formula (13) w For the third water phase conduction coefficient.

[0220] Finally, the relative penetration of the current third mound heart-containing water saturation can be determined using the formula (14). Equation (14), K ro The relative penetration rate corresponding to the current third mound heart containing aqueous saturation, k o A third oil phase effective penetration for the current third mound heart-containing water saturation.

[0221] Using the formula (15), the relative permeability of the current third mound heart-containing water saturation is determined. Formula (15) k rw Relative penetration rate corresponding to the current third mound heart containing water saturation, k w A third water phase corresponding to the current third mound heart containing a third water phase.

[0222] Among them, the current third mound heart is aqueous saturation corresponding to the third oil phase effective penetration. o Calculation method with the third water phase effective penetration ratio corresponding to the current third mound heart-containing water saturation w All of the same, the three oil phase effective penetration rate corresponding to the current third mound heart containing water saturation o Calculation method as an example:

[0223] First, the third oil phase conduction coefficient T of a polar or throat is known. o , The conductivity coefficient between any two pores is obtained by the above formula (5), and T in the formula (5) i The third oil phase conductor coefficient of pore I, T j Third oil phase conduction coefficient of pore J, T ij The third oil phase conduction coefficient of the throat between the communication pores I and the pore j.

[0224] Then, the adjacent matrix can be established by the communication between the aperture and the throat, wherein the coefficient of the abutment matrix is ​​the third oil phase conduction coefficient between any two pores. The adjacent matrix is ​​converted into a sparse matrix. Because the sparse matrix is ​​solved, the inlet end surface pressure value and the outlet end surface pressure value are set to any value, and the solving algorithm for the large-scale sparse matrix is ​​solved, and after solving the current The pressure of the respective pores of the third mound heart is aqueous saturation.

[0225] Next, the third oil phase effective penetration of the current third mound heart containing water saturation is determined by the above formula (6) and (7). o. Where the formula (7) corresponding to the formula (16), in the formula (16) k o A third oil phase effective penetration for the current third mound heart-containing water saturation.

[0226] Similarly, the third water phase effective penetration rate corresponding to the current third mound heart containing water saturation is K. w When the formula (5) - (7) is also modified, in which the equation (7) corresponding to the formula (17), the formula (17) k w A third water phase corresponding to the current third mound heart containing a third water phase.

[0227] According to the same method, the difference between the current third capillary force and the third replacement step is determined as the updated current third wool tube force. For each updated current third capsule, it is necessary to determine the current third capillary forces corresponding to the current third mound heart containing water-containing degree, and the relative penetration rate corresponding to the current third mound heart-containing water saturation. The aqueous phase is relatively permeable. The cycle is ended until the current third capillary force is less than the third capillary force in the termination state.

[0228] In this process, the current third mound heart is aqueous saturation as the abscissa, respectively, the relative permeability of the current third wool, the relative permeability of the third water, respectively, each horizontal vertical Some points under the coordinates, connect several points of each horizontal longitudinal coordinate, which can constitute a third capillary force curve, a third oil phase relative permeability curve, and a third aqueous phase relative permeability curve.

[0229] Refer Figure 13 In this example, the simulation calculation of the rigxin pore network model is mixed and wetted, and the simulation of mixed wetting rock is determined by the simulation of the mixed wetting rock drive oil. Relationship, oil phase effective penetration rate and variation of moisture saturation of the core, the change in water phase effective penetration and the moisture saturation of the core, including:

[0230] S1301: Determine the second card-breaking critical pressure value of each wet thoric, a second synergistic filling critical pressure value of each water and damp pore, a second piston-filled critical pressure value of each wet throat, initial In the state, the fourth wool tube forces in the initial state, the fourth purse in the abstain state, the fourth cap, the fourth hair tube force, and the fourth replacement step;

[0231] S1302: Determine the fourth massive tube for the initial state as the current fourth phenolic tube force;

[0232] S1303: RC performs screening of all throats and pores; it is judged whether or not the current first capillary force is greater than the fourth wool stiff force in the abort state; if the current first capillary force is greater than the inventive state The tetra-wool is tube, then the throat of the second card-off critical pressure value is greater than the current fourth capillary force, or the second synergistic filling critical pressure value is greater than the pore of the current fourth capillary force; if the current first hair tube The force is smaller than or equal to the fourth wool tubing force in the abstain state, and the throat in which the second piston filling critical pressure value is less than or equal to the current fourth capillary force, and all the pores communicating with the throat; Screening the throat and the aqueous saturation degree of the pore, the fourth rock heart containing the fourth mound heart, and determines the current fourth wool-containing power-saturation corresponding to the current fourth phenolic, according to the current fourth core. The aquatic saturation, and the fourth oil phase conduction coefficient and the fourth water phase conduction coefficient of each pore or throat, determine the fourth oil phase relative permeability and fourth water phase corresponding to the current fourth mound heart-containing water saturation. Relative penetration; update the current fourth hair tube force, the difference between the current fourth capillary force and the fourth replacement step is determined as the update, the current fourth purse, one by one, determines the current fourth Current fourth core-containing saturation corresponding to capillary force, and the relative penetration rate of the fourth oil phase corresponding to the current fourth core-containing water saturation;

[0233] S1304: The cycle ends until the current fourth phenolic tube force is smaller than the fourth piking force in the termination state.

[0234] Specifically, since the core is a mixed moisturizing core, the pore, which can set the percentage of the pore in the core as needed, and the throat that is connected is also a wet throat, and the remaining pores are a water and humid hole. The throat with the throat is a water and humid throat, with a percentage of 50%, 60%, and the like. The oil and humidifier and the wet throat can be carried out in accordance with step S501-S504 described above, in which the water and damp holes and water moisture throat can be carried out in accordance with step S401-S404 described above.

[0235] Refer Figure 14 Where the experimental value is actual experimental data for controlling the analog value, the simulation value is an experimental data in accordance with the calculation method of this article, which is compared, the simulation value is the same as the data trend of the experimental value, so the calculation method is used The resulting capillary forces and relative permeability curves are high. Specifically, since the simulation of the moisture core water drive oil process is over the end of the oil-driven water, the existing pore and throat are filled, and the wet lag effect is considered.

[0236] By determining the second card-off critical pressure value of any aqueous throat when any aqueous throat is determined by the above formula (18), in the formula (18) For the second card-off critical pressure value.

[0237] The second synergistic filling critical pressure value of any aqueous and damp pore is determined by the above formula (19), in the formula (19) For the second collaborative filling critical pressure value.

[0238] When the contact angle during water driving oil Hour,

[0239] The second piston filling critical pressure value of any wet throat is determined by the formula (20), in the formula For the second piston filling critical pressure value

[0240] When the contact angle during water driving oil When the pore or throat occurs, there is no oil film, and the aqueous phase is fully filled, and the second piston filling critical pressure value of any wet throat is determined by equation (22).

[0241] The fourth mound heart in the initial state is aqueous saturation of the first mound heart in the end of the oil-driven water. It can be set to the first state, the maximum throat pressure value in all throat pressure values ​​obtained by the above formula (8), and the fourth pneumatic tube force in the stop state is 0, and the setting is terminated. The fourth wooling tube force in the state is the minimum value in the second piston filling critical pressure value of any wet throat.

[0242]The fourth replacement step may be equal, or each step is not equal, and this paper is not limited to the specific form of the fourth replacement step. If each step is equal, it can take the difference between the fourth capillary force in the termination state and the fourth cap of the fourth cap, etc., the difference is equal to the difference, and the number of aliquots can be set according to the actual situation. For example, it can be divided into 10 times, 15 times, etc., the difference in the number of aliquots is divided into the fourth substation.

[0243] The fourth cap of the initial state is determined as the current fourth phenoliostatic force. Since the core is a mixed wet rock heart, the current first porrontal tube force is subjected to the fourth-raised step size in the initial state, and the water and wet pore and water wet throat are replaced by the fourth drive. The sip is the throat of the second card-off critical pressure value greater than the current fourth capillary force, or the second synergistic filling critical pressure value is greater than the pores of the current fourth capillary force. Until the current fourth hair tube force is reduced to the fourth porrontal fine tube force in the abort state, the water and damp pores and the water and wet throat are replaced, and the oil and wet pores and wet throat are started: screen out of the second piston The filling critical pressure value is less than or equal to the current fourth capillary force, and all the pores communicated with the throat.

[0244] As the replacement, during each step of the replacement process, the water saturation of the throat and the pores of the screen is updated to S. w = 1.

[0245] According to the above formula (11), the screen-selected throat and the pore-containing moisture saturation and the unpackable throat and the pore-containing saturation are determined to determine the current fourth federal heart-containing water saturation. Spend. Equation (11) w总 Current fourth mound heart with the current fourth phenoliostatic force containing water saturation.

[0246] Further, the fourth oil phase conductivity coefficient corresponding to any aperture or throat is calculated by the above formula (12), and T in the formula (12) o For the fourth oil phase conduction coefficient.

[0247] Further, the fourth water phase conduction coefficient of any pore or throat is calculated by the formula (13), and T in the formula (13) w For the fourth water phase conduction coefficient.

[0248] Finally, the secondary oil phase relative penetration rate corresponding to the current fourth rock heart is used in the formula (14). Equation (14), K ro Relative penetration rate corresponding to the current fourth core-containing water saturation, k o An effective permeability of the fourth oil corresponding to the current fourth mound heart containing water saturation.

[0249] The fourth water phase relative permeability corresponding to the current fourth core-containing water saturation is determined using the formula (15). Formula (15) k rw For the fourth water phase relative penetration rate corresponding to the current fourth core-containing water saturation, k w A fourth water phase effective penetration rate corresponding to the current fourth mound heart water saturation.

[0250] Among them, the current fourth rock is aqueous saturation of the fourth oil phase effective penetration rate K o Calculation method with the fourth water phase effective penetration rate corresponding to the current fourth rock heart-containing water saturation w All of the same, the fourth oil phase effective penetration ratio corresponding to the current fourth core-containing water saturation o Calculation method as an example:

[0251] First, the fourth oil phase conduction coefficient T of a polar or throat is known. o , The conductivity coefficient between any two pores is obtained by the above formula (5), and T in the formula (5) i The fourth oil phase conduction coefficient of pore I, T j The fourth oil phase conduction coefficient of pore J, T ij The fourth oil phase conduction coefficient of the throat between the communication pore I and the pore j.

[0252] Then, the adjacent matrix can be established by the communication between the chart theory, with the communication between the pores and the throat, wherein the coefficient of the abutment matrix is ​​the fourth oil phase conductivity coefficient between any two pores. The adjacent matrix is ​​converted into a sparse matrix. Because the sparse matrix is ​​solved, the inlet end surface pressure value and the outlet end surface pressure value are set to any value, and the solving algorithm for the large-scale sparse matrix is ​​solved, and after solving the current The pressure of the respective pores corresponding to the fourth mound heart containing water saturation.

[0253] Next, the fourth oil phase effective penetration of the current fourth core-containing water saturation is determined by the above formula (6) and (7). o. Where the formula (7) corresponding to the formula (16), in the formula (16) k o An effective permeability of the fourth oil corresponding to the current fourth mound heart containing water saturation.

[0254] Similarly, calculating the fourth water phase effective penetration rate corresponding to the current fourth core-containing water saturation w When the formula (5) - (7) is also modified, in which the equation (7) corresponding to the formula (17), the formula (17) k w A fourth water phase effective penetration rate corresponding to the current fourth mound heart water saturation.

[0255] According to the same method, the difference between the current fourth capillary force and the fourth replacement step is determined to be updated, the current fourth wool tube force. For each updated current fourth phenoli, it is necessary to determine the current fourth milky-containing water saturation corresponding to the current fourth phenoliosis, and the fourth oil phase relative penetration rate corresponding to the current fourth mound heart containing water saturation and fourth The aqueous phase is relatively permeable. The circulation ends the cycle until the current fourth wool tube force is less than the fourth capsule tube force in the termination state.

[0256] In this process, the current fourth core-containing water saturation is the abscissa, respectively, the relative permeability of the current fourth wool, the fourth oil phase, respectively, and the fourth water phase relative penetration is ordinary counted, each horizontal vertical Several points under the coordinates, connect several points of each horizontal longitudinal coordinate, which can constitute the fourth pelletizer curve, the fourth oil phase relative permeability curve, and the fourth water phase relative permeability curve.

[0257] Based on the above described capillary forces and relative permeability curves, the embodiments provide a capillary force and a computing device of the relative permeability curve. The apparatus described may include a system (including distributed systems), software (applications), modules, components, servers, clients, and the like and combine the necessary implementation hardware. Based on the same innovation concept, the apparatus in one or more embodiments provided herein is as described below. Since the implementation scheme of the problem solves the problem is similar, the implementation of the embodiments of the embodiments will be described in the implementation and repetition of the foregoing method. The term "unit" or "module" as used below can achieve a combination of software and / or hardware for predetermined functions. Although the apparatus described in the following examples is preferably implemented in software, the implementation of hardware, or combinations of software and hardware may also be conceived.

[0258] specifically, Figure 15 It is a schematic diagram of a module structure of a computing device of a capillary force and a relative permeability curve provided herein, referring to Figure 15 As shown, a computing device of a capillary force and a relative permeability curve provided herein includes: a model determination module 100, an absolute permeability determining module 200, an alternative to analog calculation module 300, a relationship determining module 400.

[0259] Model determination module 100: Determine the cross section of the pore and the throat of the four-stranded rock pore network model;

[0260] Absolute Permeability Determination Module 200: Single-phase flow simulation calculation of the rock core pore network model, determine the absolute permeability of the core;

[0261] The analog simulation calculation module 300: The alternative to the rock core pore network model sequentially performs a variety of replacement types, which determine the change in capillary forces and core-containing water saturation of each replacement simulation calculation. The relationship between penetration and moisture-saturation degree of osmosis and moisture, water phase effective penetration and variation of moisture saturation of rock core;

[0262] Relationship determination module 400: The effect of the effective penetration calculation corresponding to each of the alternative simulation and the variation of the moisture saturation of the core, the water phase effective penetration and the variation of the moisture saturation of the core, respectively, respectively, respectively. The absolute permeability is combined to determine the relationship between the relative permeability of the oil phase corresponding to the alternative and simulation calculation of the alternative and madness and the change in the relationship between the relative permeability of the aqueous phase and the variation of the moisture saturation.

[0263] Refer Figure 16 As shown, based on the above-described capillary force and a relative permeability curve, a computer device 1602 is also provided herein, wherein the above method operates on computer device 1602. Computer device 1602 can include one or more processor 1604, such as one or more central processing unit (CPU), or graphics processor (GPU), each processing unit, can implement one or more hardware threads. Computer device 1602 can also include any memory 1606 for storing any kind of information such as code, setting, data, and the like, in a specific embodiment, computer program, memory 1606, can run on processor 1604, When the computer program is operated by the processor 1604, an instruction according to the method can be performed. Non-limiting, for example, memory 1606 can include any one or more combinations: any type of RAM, any type of ROM, flash device, hard disk, disc, etc. More generally, any memory can be used to store information. Further, any memory can provide volatility or non-volatile retaining. Further, any memory can represent the fixed or removable member of the computer device 1602. In one case, when the processor 1604 performs an associated instruction stored in any of the combinations of memory or memory, the computer device 1602 can perform any of the related instructions. Computer device 1602 also includes one or more drive mechanisms 1608 for interacting with any memory, such as hard disk drive mechanisms, optical disc drive mechanisms, and the like.

[0264] Computer device 1602 can also include an input / output module 1610 (I / O) for receiving various inputs (via input device 1612) and for providing various outputs (via output device 1614). A specific output mechanism can include presentation device 1616 and associated graphical user interface 1618 (GUI). In other embodiments, it is also possible to not include input / output module 1610 (I / O), input device 1612, and output device 1614, only one computer device in the network. Computer device 1602 can also include one or more network interfaces 1620 for exchanging data via one or more communication links 1622. One or more communication bus 1624 couple the components described above together.

[0265] Communication link 1622 can be implemented in any way, for example, by a local area network, a wide area network (e.g., a point-to-point connection, etc., or any combination thereof. Communication link 1622 can include any combination of hard-wired links, wireless links, routers, gateway functions, name servers, etc. disposed of any protocol or protocol.

[0266] Correspond to figure 1 , 3 The method in 4, 7, 9, 11, 13, which provides a computer readable storage medium that stores a computer program on the computer readable storage medium, which is executed when the processor is running The step of the method.

[0267]The embodiment of this paper also provides a computer readable instruction, wherein when the processor performs the instruction, the program makes the processor executes such as figure 1 , 3 The method shown in 4, 7, 9, 11, 13.

[0268] It should be understood that in various embodiments herein, the size of the serial numbers of the above processes does not mean that the execution order of each process should be determined in its function and inherent logic, and should not be applied to the implementation of this example. The process constitutes any defined.

[0269] It is also to be understood that in the examples, the terms "and / or" are merely a relationship of the associated object, indicating that there may be three relationships. For example, A and / or B, which may be represented: Alone A, while there are three situations of A and B, alone in b. In addition, the characters "/" herein generally indicate a "or" relationship of the associated object before and after.

[0270] One of ordinary skill in the art will appreciate that the unit and algorithm step described herein can be implemented in electronic hardware, computer software or binding of both electronic hardware, computer software or both, for clear explanation of hardware and software. Interchangeability, the composition and steps of each example have been generally described in terms of functions. These functions are executed in hardware or software, depending on the specific application and design constraints of the technical solution. Professional technicians can use different methods to implement the functions described in each particular application, but this implementation should not be considered exceeding the scope of this article.

[0271] Those skilled in the art will clearly understand that in order to describe convenient and concise, the specific operation of the above-described system, apparatus, and unit, can refer to the corresponding process in the foregoing method embodiment, and will not be described herein again.

[0272] In several embodiments provided herein, it should be understood that the disclosed systems, apparatus, and methods can be implemented in other ways. For example, the device embodiment described above is merely schematic, for example, the division of the unit, only one logic function division, and may have additional division, such as a plurality of units or components, may be combined or It can be integrated into another, or some features can be ignored, or not executed. Additionally, the coupling or direct coupling or communication connection between the displayed or discussed may be indirectly coupled or communication between some interfaces, devices or units, or is also electrically mechanically, mechanical or other forms.

[0273] The unit as the separation member may be or may not be physically separated, and the components displayed as the unit may be or may not be a physical unit, i.e., in one place, or can also be distributed to a plurality of network elements. The object of the present embodiment can be implemented in accordance with the actual needs to select some or all units.

[0274] Further, each functional unit in various embodiments will be integrated into one processing unit, or may be physically generated by each unit, or two or more units integrated in one unit. The above-described integrated units can be implemented in the form of hardware, or may be implemented in the form of a software functional unit.

[0275] The integrated unit can be stored in a computer readable storage medium if implemented in the form of a software functional unit and is used as a stand-alone product. Based on this understanding, the technical solutions of this paper essentially or to contribute to prior art, or all or part of the technical solution can be embodied in the form of software products, the computer software product is stored in a storage medium. These include several instructions to enable a computer device (which can be a personal computer, server, or network device, etc.) to perform all or part of the method described in various embodiments described herein. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, RAD-ONLY MEMORY), RAM, RAM, RANDOM Access Memory, disk, or disc or optical disk, etc. can store program code .

[0276] The specific embodiments are described in this paper to illustrate the principles and embodiments of this paper, and the description of the above embodiments is intended to help understand the methods of this paper and their core ideas; at the same time, for the general technicians in the art, according to the ideology of this article In terms of specific embodiments and applications, there will be changes, in summary, the contents of this specification should not be construed as limiting this paper.