Constitutive law for well limit condition for element WEPRO
Implemented by: F. Collin, 2001
This constitutive law is used for non-linear flow analysis of solids and allows to model fluids exchanges between a well and the reservoir.
The retention curve allows to compute the relative permeabilities and to see which fluids is mobile.
Prepro: LWPRO.F
Lagamine: WPRO.F
Plane stress state | NO |
Plane strain state | YES |
Axisymmetric state | YES |
3D state | NO |
Generalized plane state | NO |
Line 1 (2I5, 60A1) | |
---|---|
IL | Law number |
ITYPE | 197 |
COMMENT | Any comment (up to 60 characters) that will be reproduced on the output listing |
Line 1 (4I5) | |
---|---|
ISR | Formulation index for $S_w$ (see Appendix 8) |
IKW | Formulation index for $k_{r,w}$ (see Appendix 8) |
IKP | Formulation index for $k_{r,p}$ (see Appendix 8) |
IINJ | Index to specify the type of well |
= 0 : production well | |
= 1 : injection well |
Line 1 (1G10.0) | |
---|---|
TFACT | Soil geometric factor |
Line 2 (4G10.0) | |
VISCW | Water dynamic viscosity |
RHOW | Water density |
VISCP | Oil dynamic viscosity |
RHOP | Oil density |
Line 3 (7G10.0) | |
CSR1 | 1st coefficient of the function $S_w$ |
CSR2 | 2nd coefficient of the function $S_w$ |
CSR3 | 3rd coefficient of the function $S_w$ |
CSR4 | 4th coefficient of the function $S_w$ |
SRES | Residual saturation degree (=$S_{res}$) |
SRFIELD | Field saturation degree ($S_{r,field}$) |
AIREV | Air entry value [Pa] |
Line 4 (3G10.0) | |
CKW1 | 1st coefficient of the function $k_{r,w}$ |
CKW2 | 2nd coefficient of the function $k_{r,w}$ |
CKW3 | 3rd coefficient of the function $k_{r,w}$ |
Line 5 (2G10.0) | |
CKP1 | 1st coefficient of the function $k_{r,p}$ |
CKP2 | 2nd coefficient of the function $k_{r,p}$ |
Line 6 (1G10.0) | |
KRMIN | Minimum value of $k_r$ |
2
SIG(1) | Current value of water flow |
SIG(2) | Current value of oil flow |
5
Q(1) | none |
Q(2) | Water saturation degree |
Q(3) | Element jacobian |
Q(4) | Well pressure at the integration point |
Q(5) | Generalised reservoir pressure at the integration point |
For IINJ = 0 (production well), the fluid flow has the following expression : \[SIG(1) = f_w=T^{well}.\frac{k_{rw}}{\mu_w}(p_w-p_{well})\]\[SIG(2) = f_o=T^{well}.\frac{k_{ro}}{\mu_o}(p_o-p_{well})\]
For IINJ = 1 (injection well), the water injection flow has the following expression : \[SIG(1) = f_w = T^{well}.\left(\frac{k_{rw}}{\mu_w}+\frac{k_{ro}}{\mu_o}\right).(p_w-p_{well})\]\[SIG(2)=0\]