Three-phase numerical model of water migration in partially frozen geological media: model formulation, validation, and applications

2010 ◽  
Vol 15 (1) ◽  
pp. 69-85 ◽  
Author(s):  
Scott L. Painter
Author(s):  
J. Grasa ◽  
M. Sierra ◽  
N. Lauzeral ◽  
M.J. Muñoz ◽  
F.J. Miana-Mena ◽  
...  

SPE Journal ◽  
2020 ◽  
Vol 25 (06) ◽  
pp. 3265-3279
Author(s):  
Hamidreza Hamdi ◽  
Hamid Behmanesh ◽  
Christopher R. Clarkson

Summary Rate-transient analysis (RTA) is a useful reservoir/hydraulic fracture characterization method that can be applied to multifractured horizontal wells (MFHWs) producing from low-permeability (tight) and shale reservoirs. In this paper, we applied a recently developed three-phase RTA technique to the analysis of production data from an MFHW completed in a low-permeability volatile oil reservoir in the Western Canadian Sedimentary Basin. This RTA technique is used to analyze the transient linear flow regime for wells operated under constant flowing bottomhole pressure (BHP) conditions. With this method, the slope of the square-root-of-time plot applied to any of the producing phases can be used to directly calculate the linear flow parameter xfk without defining pseudovariables. The method requires a set of input pressure/volume/temperature (PVT) data and an estimate of two-phase relative permeability curves. For the field case studied herein, the PVT model is constructed by tuning an equation of state (EOS) from a set of PVT experiments, while the relative permeability curves are estimated from numerical model history-matchingresults. The subject well, an MFHW completed in 15 stages, produces oil, water, and gas at a nearly constant (measured downhole) flowing BHP. This well is completed in a low-permeability,near-critical volatile oil system. For this field case, application of the recently proposed RTA method leads to an estimate of xfk that is in close agreement (within 7%) with the results of a numerical model history match performed in parallel. The RTA method also provides pressure–saturation (P–S) relationships for all three phases that are within 2% of those derived from the numerical model. The derived P–S relationships are central to the use of other RTA methods that require calculation of multiphase pseudovariables. The three-phase RTA technique developed herein is a simple-yet-rigorous and accurate alternative to numerical model history matching for estimating xfk when fluid properties and relative permeability data are available.


2016 ◽  
Vol 711 ◽  
pp. 351-358 ◽  
Author(s):  
Nathalie Ramos ◽  
Guzel Shamsutdinova ◽  
Max A.N. Hendriks ◽  
Stefan Jacobsen

The last decades the concrete-ice abrasion process has been well known as a concrete surface degradation mechanism due to ice sliding. The topic is especially relevant for concrete gravity based structures in the Arctic offshore. The article presents a numerical model in which the onset of wear in the concrete-ice abrasion process is simulated. The simulations are performed on meso-scale, which means that concrete is modelled as a three-phase material in which paste, aggregates and the interface transition zone are distinguished. Lattice modelling is adopted for the numerical modeling. Hertzian contact theory which predicts excessive tensile stresses on the concrete surface due to sliding of ice asperities is used as an analytical basis for the numerical model. It was concluded that such model is able to capture both surface and subsurface cracking in the concrete.


Sensors ◽  
2021 ◽  
Vol 21 (8) ◽  
pp. 2643
Author(s):  
Qingqing Cao ◽  
Imad L. Al-Qadi

Ground-penetrating radar (GPR) has been used for asphalt concrete (AC) pavement density prediction for the past two decades. Recently, it has been considered as a method for pavement quality control and quality assurance. A numerical method to estimate asphalt pavement specific gravity from its dielectric properties was developed and validated. A three-phase numerical model considering aggregate, binder, and air void components was developed using an AC mixture generation algorithm. A take-and-add algorithm was used to generate the uneven air-void distribution in the three-phase model. The proposed three-phase model is capable of correlating pavement density and bulk and component dielectric properties. The model was validated using field data. Two methods were used to calculate the dielectric constant of the AC mixture, including reflection amplitude and two-way travel time methods. These were simulated and compared when vertical and longitudinal heterogeneity existed within the AC pavement layers. Results indicate that the reflection amplitude method is more sensitive to surface thin layers than the two-way travel time methods. Effect of air-void content, asphalt content, aggregate gradation, and aggregate dielectric constants on the GPR measurements were studied using the numerical model.


Sign in / Sign up

Export Citation Format

Share Document