An Experimental and Numerical Study of In-Situ Steamdrive During Cyclic Steaming

1997 ◽  
Vol 12 (02) ◽  
pp. 144-150
Author(s):  
K. Dehghani ◽  
R.F. Meyer ◽  
H. Duran ◽  
M. Kumar ◽  
E.F. deZabala
Keyword(s):  
2021 ◽  
Author(s):  
Koki Oikawa ◽  
Hirotaka Saito ◽  
Seiichiro Kuroda ◽  
Kazunori Takahashi

<p>As an array antenna ground penetrating radar (GPR) system electronically switches any antenna combinations sequentially in milliseconds, multi-offset gather data, such as common mid-point (CMP) data, can be acquired almost seamlessly. However, due to the inflexibility of changing the antenna offset, only a limited number of scans can be obtained. The array GPR system has been used to collect time-lapse GPR data, including CMP data during the field infiltration experiment (Iwasaki et al., 2016). CMP data obtained by the array GPR are, however, too sparse to obtain reliable velocity using a standard velocity analysis, such as semblance analysis. We attempted to interpolate the sparse CMP data based on projection onto convex sets (POCS) algorithm (Yi et al., 2016) coupled with NMO correction to automatically determine optimum EM wave velocity. Our previous numerical study showed that the proposed method allows us to determine the EM wave velocity during the infiltration experiment.</p><p>The main objective of this study was to evaluate the performance of the proposed method to interpolate sparse array antenna GPR CMP data collected during the in-situ infiltration experiment at Tottori sand dunes. The interpolated CMP data were then used in the semblance analysis to determine the EM wave velocity, which was further used to compute the infiltration front depth. The estimated infiltration depths agreed well with independently obtained depths. This study demonstrated the possibility of developing an automatic velocity analysis based on POCS interpolation coupled with NMO correction for sparse CMP collected with array antenna GPR.</p>


2004 ◽  
Vol 22 (10) ◽  
pp. 3741-3749 ◽  
Author(s):  
A. Gonzalez-Esparza ◽  
A. Santillán ◽  
J. Ferrer

Abstract. We studied the heliospheric evolution in one and two dimensions of the interaction between two ejecta-like disturbances beyond the critical point: a faster ejecta 2 overtaking a previously launched slower ejecta 1. The study is based on a hydrodynamic model using the ZEUS-3-D code. This model can be applied to those cases where the interaction occurs far away from the Sun and there is no merging (magnetic reconnection) between the two ejecta. The simulation shows that when the faster ejecta 2 overtakes ejecta 1 there is an interchange of momentum between the two ejecta, where the leading ejecta 1 accelerates and the tracking ejecta 2 decelerates. Both ejecta tend to arrive at 1AU having similar speeds, but with the front of ejecta 1 propagating faster than the front of ejecta 2. The momentum is transferred from ejecta 2 to ejecta 1 when the shock initially driven by ejecta 2 passes through ejecta 1. Eventually the two shock waves driven by the two ejecta merge together into a single stronger shock. The 2-D simulation shows that the evolution of the interaction can be very complex and there are very different signatures of the same event at different viewing angles; however, the transferring of momentum between the two ejecta follows the same physical mechanism described above. These results are in qualitative agreement with in-situ plasma observations of "multiple magnetic clouds" detected at 1AU.


SPE Journal ◽  
2020 ◽  
Vol 25 (06) ◽  
pp. 3091-3110
Author(s):  
Ming Chen ◽  
Shicheng Zhang ◽  
Tong Zhou ◽  
Xinfang Ma ◽  
Yushi Zou

Summary Creating uniform multiple fractures is a challenging task due to reservoir heterogeneity and stress shadow. Limited-entry perforation and in-stage diversion are commonly used to improve multifracture treatments. Many studies have investigated the mechanism of limited-entry perforation for multifracture treatments, but relatively few have focused on the in-stage diversion process. The design of in-stage diversion is usually through trial and error because of the lack of a simulator. In this study, we present a fully coupled planar 2D multifracture model for simulating the in-stage diversion process. The objective is to evaluate flux redistribution after diversion and optimize the dosage of diverters and diversion timing under different in-stage in-situ stress difference. Our model considers ball sealer allocation and solves flux redistribution after diversion through a fully coupled multifracture model. A supertimestepping explicit algorithm is adopted to solve the solid/fluid coupling equations efficiently. Multifracture fronts are captured by using tip asymptotes and an adaptive time-marching approach. The modeling results are validated against analytical solutions for a plane-strain Khristianovic-Geertsma de Klerk (KGD) model. A series of numerical simulations are conducted to investigate the multifracture growth under different in-stage diversion operations. Parametric studies reveal that the in-stage in-situ stress difference is a critical parameter for diversion designs. When the in-situ stress difference is larger than 2 MPa, the fracture in the high-stress zone can hardly be initiated before diversion for a general fracturing design. More ball sealers are required for the formations with higher in-stage in-situ stress difference. The diverting time should be earlier for formations with high in-stage stress differences as well. Adding more perforation holes in the zone with higher in-situ stress is recommended to achieve even flux distribution. The results of this study can help understand the multifracture growth mechanism during in-stage diversion and optimize the diversion design timely.


2020 ◽  
Author(s):  
Nathalie Casas ◽  
Guilhem Mollon ◽  
Ali Daouadji

<p>How do earthquakes start? What are the parameters influencing fault evolutions? What are the local parameters controlling the seismic or aseismic character of slip?</p><p>To predict the dynamic behaviour of faults, it is important to understand slip mechanisms and their source. Lab or in-situ experiments can be very helpful, but tribological experience has shown that it is complicated to install local sensors inside a mechanical contact, and that they could disturb the behaviour of the sheared medium. Even with technical improvements on lab tools, some interesting data regarding gouge kinematics and rheology remains very difficult or impossible to obtain. Numerical modelling seems to be another way of understanding physics of earthquakes.</p><p>Fault zone usually present a granular gouge, coming from the wear material of previous slips. That is why, in this study, we present a numerical model to observe the evolution and behaviours of fault gouges. We chose to focus on physics of contacts inside a granular gouge at a millimetre-scale, studying contact interactions and friction coefficient between the different bodies. In order to get access to this kind of information, we implement a 2D granular fault gouge with Discrete Element Modelling in the software MELODY (Mollon, 2016). The gouge model involves two rough surfaces representing the rock walls separated by the granular gouge.</p><p>One of the interests of this code is its ability to represent realistic non-circular grain shapes with a Fourier-Voronoï method (Mollon et al., 2012). As most of the simulations reported in the literature use circular (2D)/spherical (3D) grains, we wanted to analyse numerically the contribution of angular grains. We confirm that they lead to higher friction coefficients and different global behaviours (Mair et al., 2002), (Guo et al., 2004).</p><p>In a first model, we investigate dry contacts to spotlight the influence of inter-particular cohesion and small particles on slip behaviour and static friction. A second model is carried out to observe aseismic and seismic slips occurring within the gouge. As stability depends on the interplay between the peak of static friction and the stiffness of the surrounding medium, the model includes the stiffness of the loading apparatus on the rock walls.</p><p>The work presented here focuses on millimetre-scale phenomena, but the employed model cannot be extended to the scale of the entire fault network, for computational cost reasons. It is expected, however, that it will lead to a better understanding of local behaviours that may be injected as simplified interface laws in larger-scale simulations.</p>


2021 ◽  
pp. 179011
Author(s):  
Linglong Wang ◽  
Xuecheng Wu ◽  
Xiang Gao ◽  
Yingchun Wu ◽  
Kefa Cen

2011 ◽  
Vol 117-119 ◽  
pp. 1723-1727 ◽  
Author(s):  
Jun Qi Wang

Deeply buried tunnels usually lie in high stress fields, whose horizontal stress which is not uniform is far larger than vertical stress, and their stability is dominated by the original in-situ stresses. With three-dimensional nonlinear finite element method, the axis orientation effects of tunnel on the displacement and stability of two types of surrounding rocks are studied systematically for one water diversion project. The tunnel lies in different original stress fields whose maximum horizontal principal stress is parallel with or perpendicular to the axis and lies in different kinds of rocks. The numerical analysis results show that the plastic zones develop in side wall of tunnel mostly when the horizontal maximum principal stress is parallel with the tunnel axis while the plastic zones distribute in the top and bottom of tunnel when the horizontal maximum principal stress is perpendicular to the tunnel axis. It is concluded that the principle of tunnel axis should be parallel with horizontal maximum principal stress regulated by the “specification for design of hydraulic tunnel” is not available for the stability of tunnel always.


Sign in / Sign up

Export Citation Format

Share Document