Pressure Transient Analysis of Wells in the Fault-Karst Carbonate Reservoirs with Vertical Beads-on-String Structure: Case Studies in Shunbei Oilfield, Tarim Basin of Northwestern China

2021 ◽  
Author(s):  
Cao Wei ◽  
Shiqing Cheng ◽  
Yang Wang ◽  
Ruyuan Shang ◽  
Lang Zhu ◽  
...  

Abstract The newly discovered the fault-karst carbonate reservoirs in Tarim Basin are formed by large-scale tectonic fault activities and multiple-stage karstification. The ground outcrop, seismic reflection and well logging show that the large caves, vugs, tectonic fractures and matrix coexist in the reservoirs. The fractures interconnecting with caves in series to form vertical beads-on-string structure is the most common pattern. It is found that conventional models are difficult to match the recorded pressure data. To fill this gap, this work summarizes three kinds of patterns for this structure and presents three novel models to estimate formation properties. The physical models of multi-fracture-region multi-cave-region series connection are established by simplifying vertical beads-on-string structure. The corresponding mathematical models are developed, in which the flow in fracture regions obey Darcy's law, while the flow in cave regions obeys free flow. Importantly, the gravity is considered due to the flow along vertical direction. Then typical flow regimes are analyzed and sensitivity analysis is conducted. Our work shows that pressure-derivative curves show similar ‘concave’ characteristic for the cave storage regime and vug/matrix interporosity flow regime. The difference is that the pressure derivative for cave storage regime is unit slope, while this slope is not equal to unity for vug/matrix interporosity flow regime. Therefore, large cave and vug/matrix medium can be distinguished by slope of pressure derivative. More than that, the typical characteristic of vertical beads-on-string structure on type curve is that the cave storage regimes and linear flow regimes alternately appear. This characteristic helps the engineers identify vertical beads-on-string structure. A novel finding is that gravity effect could lead to unit-slope pressure and pressure derivative at late times, which is traditional recognized as the presence of a closed boundary when it may not be the case. Lastly, this methodology is applied to two cases from Shunbei Oilfield in which it is difficult to obtain good fitting quality and interpretation results using traditional methods. Besides conventional properties, the proposed methodology allows us to estimate other properties (e.g., cave height, cave radius), which are not readily obtained from conventional methods.

Author(s):  
Cao Wei ◽  
Shiqing Cheng ◽  
Gang Chen ◽  
Wenyang Shi ◽  
Jiaxin Wu ◽  
...  

Tarim Basin newly discovered the fault-karst carbonate reservoirs, which are formed by the large-scale tectonic fault activities and multiple-stage karstification. Four kinds of mediums coexist in the reservoirs, including the large cave, vug, tectonic fracture and matrix. The tectonic fractures interconnect with large caves in series to form the vertical beads-on-string structure, which is the most common connection pattern in reservoirs. To provide a well test method for evaluating this type of structure, this work firstly presents a multi-fracture-region multi-cave-region series connection physical model by simplifying vertical beads-on-string structure. We consider four kinds of mediums in the proposed physical model, including large caves, small vugs, high-angle tectonic fracture and rock matrix. The fracture regions mainly contain fracture, vug and matrix mediums. The cave regions contain cave medium. The corresponding mathematical model is also developed, in which the flow in fracture regions obeys the Darcy’s law, while the flow in cave regions is assumed to obey free flow. Furthermore, the gravity is taken into account because the flow is along the vertical direction. Then the typical flow regimes are analyzed and sensitivity analysis is conducted on crucial parameters. Results indicate that (a) the typical feature of vertical beads-on-string structure on type curves is that the cave storage regimes and linear flow regimes alternately appear; (b) the type curves will exhibit the cave storage regimes with unit-slope pressure derivative for the existence of large caves, which is different from the inter-porosity flow regimes for the existence of the vugs (slope ≠ 1); (c) the gravity effect could lead to unit-slope pressure and pressure derivative curves, which can be regarded as closed boundary in a peculiar sense; (d) gravity effect is difficult to be observed from well test curves with about 2-weeks test duration in real application. Finally, two cases from Shunbei Oilfield are interpreted to illustrate the practicability and feasibility of proposed method.


2016 ◽  
Vol 47 (S1) ◽  
pp. 149-160 ◽  
Author(s):  
Xijun Lai ◽  
Qiuhua Liang ◽  
Qun Huang ◽  
Jiahu Jiang ◽  
X. X. Lu

The full operation of the Three Gorges Dam (TGD) has altered the downstream natural flow regime. Flow regime changes have resulted in profound influences on the utility of water resources and hence a large area with a riparian ecosystem including China's two largest freshwater lakes in the Middle Yangtze. Because of complicated flow regimes in this large-scale river–lake system, the TGD's impacts on flow regimes are highly heterogeneous and require to be carefully addressed. To better understand them, we estimated water level and discharge changes solely induced by the TGD from 2006 to 2011 using a hydrodynamic model that facilitates the separation of the TGD's contribution to flow regimes. Results indicated that water regulation of the TGD caused profound impacts on the flow regimes of the Middle Yangtze. In the impoundment period from mid-September to October, rapid and significant decline of the water discharge downstream the TGD produced a prolonged dry season that occurred around 10 days earlier than before. Our analysis elucidated a pattern of recent changes in the flow regimes caused by the TGD. The findings are useful for addressing the TGD-induced environmental issues, optimizing the TGD's operation, and generating adaptive management strategy for the complex river–lake ecosystem.


2021 ◽  
Vol 7 (1) ◽  
Author(s):  
Tao Yue ◽  
Da Zhao ◽  
Duc T. T. Phan ◽  
Xiaolin Wang ◽  
Joshua Jonghyun Park ◽  
...  

AbstractThe vascular network of the circulatory system plays a vital role in maintaining homeostasis in the human body. In this paper, a novel modular microfluidic system with a vertical two-layered configuration is developed to generate large-scale perfused microvascular networks in vitro. The two-layer polydimethylsiloxane (PDMS) configuration allows the tissue chambers and medium channels not only to be designed and fabricated independently but also to be aligned and bonded accordingly. This method can produce a modular microfluidic system that has high flexibility and scalability to design an integrated platform with multiple perfused vascularized tissues with high densities. The medium channel was designed with a rhombic shape and fabricated to be semiclosed to form a capillary burst valve in the vertical direction, serving as the interface between the medium channels and tissue chambers. Angiogenesis and anastomosis at the vertical interface were successfully achieved by using different combinations of tissue chambers and medium channels. Various large-scale microvascular networks were generated and quantified in terms of vessel length and density. Minimal leakage of the perfused 70-kDa FITC-dextran confirmed the lumenization of the microvascular networks and the formation of tight vertical interconnections between the microvascular networks and medium channels in different structural layers. This platform enables the culturing of interconnected, large-scale perfused vascularized tissue networks with high density and scalability for a wide range of multiorgan-on-a-chip applications, including basic biological studies and drug screening.


Author(s):  
James J. Bell ◽  
David K.A. Barnes

Sponge communities were sampled at 3 m depth intervals at six sites experiencing different flow regimes at Lough Hyne, Ireland. Sponges were identified and classified within the following morphological groups: encrusting, massive, globular, pedunculate, tubular, flabellate, arborescent, repent and papillate morphological types on both vertical (≈90°) and inclined (≈45°) surfaces.Differences in the proportional abundance of the sponge body forms and density (sponge m−2) were observed between sites and depths. The density of sponges increased with depth at sites of slight to moderate current flow, but not at the site where current flow was turbulent. Morphological diversity of sponge communities decreased with increasing current flow due to the removal of delicate forms such as pedunculate and arborescent shaped sponges. Massive and encrusting morphologies dominated at the high-energy sites (fast and turbulent flow regimes) due to a high basal area to volume ratio, which prevents removal from cliff surfaces. However, pedunculate, papillate and arborescent types dominated at the low current sites as these shapes may help to prevent the settlement of sediment on sponge surfaces. Bray–Curtis Similarity analysis and Correspondence Analysis were used to distinguish between five different morphological communities.


2002 ◽  
Vol 464 ◽  
pp. 81-111 ◽  
Author(s):  
OLEG ZIKANOV ◽  
DONALD N. SLINN ◽  
MANHAR R. DHANAK

We present the results of large-eddy simulations (LES) of turbulent thermal convection generated by surface cooling in a finite-depth stably stratified horizontal layer with an isothermal bottom surface. The flow is a simplified model of turbulent convection occurring in the warm shallow ocean during adverse weather events. Simulations are performed in a 6 × 6 × 1 aspect ratio computational domain using the pseudo-spectral Fourier method in the horizontal plane and finite-difference discretization on a high-resolution clustered grid in the vertical direction. A moderate value of the Reynolds number and two different values of the Richardson number corresponding to a weak initial stratification are considered. A version of the dynamic model is applied as a subgrid-scale (SGS) closure. Its performance is evaluated based on comparison with the results of direct numerical simulations (DNS) and simulations using the Smagorinsky model. Comprehensive study of the spatial structure and statistical properties of the developed turbulent state shows some similarity to Rayleigh–Bénard convection and other types of turbulent thermal convection in horizontal layers, but also reveals distinctive features such as the dominance of a large-scale pattern of descending plumes and strong turbulent fluctuations near the surface.


Author(s):  
Mahmud R. Amin ◽  
Nallamuthu Rajaratnam ◽  
David Z. Zhu

Abstract This work presents an analytical study of the flow and energy loss immediately downstream of rectangular sharp-crested weirs for free and submerged flows, using the theory of plane turbulent jets and the analysis of some relevant studies. The flow regimes downstream of the sharp-crested weir is characterized as the impinging jet and surface flow regimes. Based on the flow characteristics and the downstream tailwater depths, each flow regime is further classified, and the relative energy loss equation is developed. It is found that significant energy loss occurs for the regime of supercritical flow and the upper stage of impinging jet flow. The energy loss for the submerged flow regime is minimal.


1969 ◽  
Vol 9 (03) ◽  
pp. 293-300 ◽  
Author(s):  
J.E. Varnon ◽  
R.A. Greenkorn

Abstract This paper reports an investigation of unstable fingering in two-fluid flow in a porous medium to determine if lambda the dimensionless finger width, is unique For a viscous finger A is the ratio of finger width to the distance between the tips of the two trailing fingers adjacent to the leading finger. For a gravity finger lambda is defined as the ratio of finger width, to "height" of the medium perpendicular to hulk flow. This work confirms previous experiments and existing theory that for viscous fingering lambda approaches a value of 0.5 with increasing ratio of viscous to interfacial force. However, for a given fluid pair and given, medium, this ratio can he increased only by increasing the, velocity. Experiments on gas liquid systems show that the asymptotic value of lambda with velocity is not always 0.5. Apparently, for gas-liquid systems, the influence of the interfacial force cannot always he eliminated by increasing the velocity. For such systems lambda is a function of fluid pair and media permeability. If the gravity force normal to the hulk permeability. If the gravity force normal to the hulk flow is active, it damps out the viscous fingers except for an underlying or overlying finger. The dimensionless width of this gravity finger strongly depends on velocity and height of the medium, as well as the fluid and media properties. The existing experiments and theories are reviewed and the gravity, stable, and viscous flow regimes are described in view of these experiments and theories. The existence of a gravity-dominated unstable regime, a gravity-viscous balanced stable regime, and a viscous-anminated regime was demonstrated experimentally by increasing flow velocity bin a rectangular glass head model. Asymptotic values of the dimensionless finger width were determined in various-sized Hele-Shaw models with gravity perpendicular and parallel to flow. The dimensionless perpendicular and parallel to flow. The dimensionless finger width lambda was determined as a function of applied force, flow resistance, and fluid properties. The results are interpreted dimensionally. Some comments are made concerning possible scaling and meaningful extensions of theory to describe these regimes in three-dimensional flow. Previous description of unstable two-fluid flow in porous media is mainly restricted to studies of viscous-dominated instability. The direction of this study is to provide data and understanding to consider the more realistic problem of predicting flow in three dimensions that may result in instabilities that are combinations of all, four flow regimes. Introduction The unstable flow of two fluids is characterized by interface changes between the fluids as a result of changes in relative forces. In a given porous medium and for a given fluid pair the gravity force dominates flow at low displacement velocities. As the velocity increases the viscous forces begin to affect flow significantly, and eventually there is a balance between effects of the gravity and viscous forces. As velocity increases further, the viscous force dominates flow. In the plane parallel to gravity, four flow regimes result as the velocity is increased: a gravity-induced stable flow regime; a gravity-dominated unstable flow regime; a stable regime resulting from a balance between gravity and viscous forces; and a viscous-induced unstable flow regime. The gravity-induced stable regime is represented schematically in Fig. 1a. This general flow pattern persists with the displacing fluid contacting all of persists with the displacing fluid contacting all of the in-place fluid until the interface becomes parallel to the bulk flow. At this velocity a gravity finger forms, and the interface, is unstable in that the length of the gravity finger grows and the fluid behind the nose of the finger is practically nonmobile because of the small pressure gradient along the finger. The gravity-dominated unstable flow is shown schematically in Fig. 1b. As the injection rate is increased, the gravity finger thickens, perhaps until it spans the medium creating a stable interface where all of the in-place, fluid is again mobile. This regime would, not occur in the absence of gravity. It occurs due to the counter effects of the gravity and viscous forces (Fig. 1c). As the velocity of the displacing fluid increases, the viscous forces dominate, and, the interface breaks into viscous fingers (Fig. 1d). SPEJ p. 293


2016 ◽  
Author(s):  
Dominik Paprotny ◽  
Oswaldo Morales Nápoles

Abstract. Large-scale hydrological modelling of flood hazard requires adequate extreme discharge data. Models based on physics are applied alongside those utilizing only statistical analysis. The former requires enormous computation power, while the latter are most limited in accuracy and spatial coverage. In this paper we introduce an alternate, statistical approach based on Bayesian Networks (BN), a graphical model for dependent random variables. We use a non-parametric BN to describe the joint distribution of extreme discharges in European rivers and variables describing the geographical characteristics of their catchments. Data on annual maxima of daily discharges from more than 1800 river gauge stations were collected, together with information on terrain, land use and climate of catchments that drain to those locations. The (conditional) correlations between the variables are modelled through copulas, with the dependency structure defined in the network. The results show that using this method, mean annual maxima and return periods of discharges could be estimated with an accuracy similar to existing studies using physical models for Europe, and better than a comparable global statistical method. Performance of the model varies slightly between regions of Europe, but is consistent between different time periods, and is not affected by a split-sample validation. The BN was applied to a large domain covering all sizes of rivers in the continent, both for present and future climate, showing large variation in influence of climate change on river discharges, as well as large differences between emission scenarios. The method could be used to provide quick estimates of extreme discharges at any location for the purpose of obtaining input information for hydraulic modelling.


2017 ◽  
Vol 74 (11) ◽  
pp. 3515-3532 ◽  
Author(s):  
Shuguang Wang ◽  
Adam H. Sobel

Abstract A set of idealized cloud-permitting simulations is performed to explore the influence of small islands on precipitating convection as a function of large-scale wind speed. The islands are situated in a long narrow ocean domain that is in radiative–convective equilibrium (RCE) as a whole, constraining the domain-average precipitation. The island occupies a small part of the domain, so that significant precipitation variations over the island can occur, compensated by smaller variations over the larger surrounding oceanic area. While the prevailing wind speeds vary over flat islands, three distinct flow regimes occur. Rainfall is greatly enhanced, and a local symmetric circulation is formed in the time mean around the island, when the prevailing large-scale wind speed is small. The rainfall enhancement over the island is much reduced when the wind speed is increased to a moderate value. This difference is characterized by a change in the mechanisms by which convection is forced. A thermally forced sea breeze due to surface heating dominates when the large-scale wind is weak. Mechanically forced convection, on the other hand, is favored when the large-scale wind is moderately strong, and horizontal advection of temperature reduces the land–sea thermal contrast that drives the sea breeze. Further increases of the prevailing wind speed lead to strong asymmetry between the windward and leeward sides of the island, owing to gravity waves that result from the land–sea contrast in surface roughness as well as upward deflection of the horizontal flow by elevated diurnal heating. Small-amplitude topography (up to 800-m elevation is considered) has a quantitative impact but does not qualitatively alter the flow regimes or their dependence on wind speed.


2021 ◽  
Author(s):  
Aurore Lafond ◽  
Maurice Ringer ◽  
Florian Le Blay ◽  
Jiaxu Liu ◽  
Ekaterina Millan ◽  
...  

Abstract Abnormal surface pressure is typically the first indicator of a number of problematic events, including kicks, losses, washouts and stuck pipe. These events account for 60–70% of all drilling-related nonproductive time, so their early and accurate detection has the potential to save the industry billions of dollars. Detecting these events today requires an expert user watching multiple curves, which can be costly, and subject to human errors. The solution presented in this paper is aiming at augmenting traditional models with new machine learning techniques, which enable to detect these events automatically and help the monitoring of the drilling well. Today’s real-time monitoring systems employ complex physical models to estimate surface standpipe pressure while drilling. These require many inputs and are difficult to calibrate. Machine learning is an alternative method to predict pump pressure, but this alone needs significant labelled training data, which is often lacking in the drilling world. The new system combines these approaches: a machine learning framework is used to enable automated learning while the physical models work to compensate any gaps in the training data. The system uses only standard surface measurements, is fully automated, and is continuously retrained while drilling to ensure the most accurate pressure prediction. In addition, a stochastic (Bayesian) machine learning technique is used, which enables not only a prediction of the pressure, but also the uncertainty and confidence of this prediction. Last, the new system includes a data quality control workflow. It discards periods of low data quality for the pressure anomaly detection and enables to have a smarter real-time events analysis. The new system has been tested on historical wells using a new test and validation framework. The framework runs the system automatically on large volumes of both historical and simulated data, to enable cross-referencing the results with observations. In this paper, we show the results of the automated test framework as well as the capabilities of the new system in two specific case studies, one on land and another offshore. Moreover, large scale statistics enlighten the reliability and the efficiency of this new detection workflow. The new system builds on the trend in our industry to better capture and utilize digital data for optimizing drilling.


Sign in / Sign up

Export Citation Format

Share Document