Productivity Prediction Model of Perforated Horizontal Well Based on Permeability Calculation in Near-Well High Permeability Reservoir Area

To improve the productivity of oil wells,perforation technology is usually used to improve the productivity of horizontal wells in oilfield exploitation.After the perforation operation,the perforation channel around the wellbore will form a near-well high-permeability reservoir area with the penetration depth as the radius,that is,the formation has different permeability characteristics with the perforation depth as the dividing line.Generally,the permeability is measured by the permeability tester,but this approach has a high workload and limited application.In this paper,according to the reservoir characteristics of perforated horizontal wells,the reservoir is divided into two areas:the original reservoir area and the near-well high permeability reservoir area.Based on the theory of seepage mechanics and the formula of open hole productivity,the permeability calculation formula of near-well high permeability reservoir area with perforation parameters is deduced.According to the principle of seepage continuity,the seepage is regarded as the synthesis of two directions:the horizontal plane elliptic seepage field and the vertical plane radial seepage field,and the oil well productivity prediction model of the perforated horizontal well is established by partition.The model comparison demonstrates that the model is reasonable and feasible.To calculate and analyze the effect of oil well production and the law of influencing factors,actual production data of the oilfield are substituted into the oil well productivity formula.It can effectively guide the technical process design and effect prediction of perforated horizontal wells.

Two-dimensional horizontal visibility graph analysis of human brain aging on gray matter

Characterizing the trajectory of the healthy aging brain and exploring age-related structural changes in the brain can help deepen our understanding of the mechanism of brain aging.Currently,most structural magnetic resonance imaging literature explores brain aging merely from the perspective of morphological features,which cannot fully utilize the grayscale values containing important intrinsic information about brain structure.In this study,we propose the construction of two-dimensional horizontal visibility graphs based on the pixel intensity values of the gray matter slices directly.Normalized network structure entropy(NNSE)is then introduced to quantify the overall heterogeneities of these graphs.The results demonstrate a decrease in the NNSEs of gray matter with age.Compared with the middle-aged and the elderly,the larger values of the NNSE in the younger group may indicate more homogeneous network structures,smaller differences in importance between nodes and thus a more powerful ability to tolerate intrusion.In addition,the hub nodes of different adult age groups are primarily located in the precuneus,cingulate gyrus,superior temporal gyrus,inferior temporal gyrus,parahippocampal gyrus,insula,precentral gyrus and postcentral gyrus.Our study can provide a new perspective for understanding and exploring the structural mechanism of brain aging.

Risk Assessment of Deep-Water Horizontal X-Tree Installation

Due to the high potential risk and many influencing factors of subsea horizontal X-tree installation,to guarantee the successful completion of sea trials of domestic subsea horizontal X-trees,this paper established a modular risk evaluation model based on a fuzzy fault tree.First,through the analysis of the main process oftree down and combining the Offshore&Onshore Reliability Data(OREDA)failure statistics and the operation procedure and the data provided by the job,the fault tree model of risk analysis of the tree down installation was established.Then,by introducing the natural language of expert comprehensive evaluation and combining fuzzy principles,quantitative analysis was carried out,and the fuzzy number was used to calculate the failure probability of a basic event and the occurrence probability of a top event.Finally,through a sensitivity analysis of basic events,the basic events of top events significantly affected were determined,and risk control and prevention measures for the corresponding high-risk factors were proposed for subsea horizontal X-tree down installation.

Revolutionizing Tight Reservoir Production: A Novel Dual-Medium Unsteady Seepage Model for Optimizing Volumetrically Fractured Horizontal Wells

This study presents an avant-garde approach for predicting and optimizing production in tight reservoirs,employing a dual-medium unsteady seepage model specifically fashioned for volumetrically fractured horizontal wells.Traditional models often fail to fully capture the complex dynamics associated with these unconventional reservoirs.In a significant departure from these models,our approach incorporates an initiation pressure gradient and a discrete fracture seepage network,providing a more realistic representation of the seepage process.The model also integrates an enhanced fluid-solid interaction,which allows for a more comprehensive understanding of the fluid-structure interactions in the reservoir.This is achieved through the incorporation of improved permeability and stress coupling,leading to more precise predictions of reservoir behavior.The numerical solutions derived from the model are obtained through the sophisticated finite element method,ensuring high accuracy and computational efficiency.To ensure the model’s reliability and accuracy,the outcomes were tested against a real-world case,with results demonstrating strong alignment.A key revelation from the study is the significant difference between uncoupled and fully coupled volumetrically fractured horizontal wells,challenging conventional wisdom in the field.Additionally,the study delves into the effects of stress,fracture length,and fracture number on reservoir production,contributing valuable insights for the design and optimization of tight reservoirs.The findings from this study have the potential to revolutionize the field of tight reservoir prediction and management,offering significant advancements in petroleum engineering.The proposed approach brings forth a more nuanced understanding of tight reservoir systems and opens up new avenues for optimizing reservoir management and production.

Development and application of a two-dimensional water quality model for the Daqinghe River Mouth of the Dianchi Lake

Water quality models are important tools to support the optimization of aquatic ecosystem rehabilitation programs and assess their efficiency. Basing on the flow conditions of the Daqinghe River Mouth of the Dianchi Lake, China, a two-dimensional water quality model was developed in the research. The hydrodynamics module was numerically solved by the alternating direction iteration （ADI） method. The parameters of the water quality module were obtained through the in situ experiments and the laboratory analyses that were conducted from 2006 to 2007. The model was calibrated and verified by the observation data in 2007. Among the four modelled key variables, i.e., water level, COD （in CODcr）, NH4＋-N and PO43-P the minimum value of the coefficient of determination （COD） was 0.69, indicating the model performed reasonably well. The developed model was then applied to simulate the water quality changes at a downstream cross-section assuming that the designed restoration programs were implemented. According to the simulated results, the restoration programs could cut down the loads of COD and PO43-P about 15%. Such a load reduction, unfortunately, would have very little effect on the NH4^＋-N removal. Moreover, the water quality at the outlet cross-section would be still in class V （3838-02）, indicating more measures should be taken to further reduce the loads. The study demonstrated the capability of water quality models to support aquatic ecosystem restorations.

Simulation of two-dimensional interior ballistics model of solid propellant electrothermal-chem ical launch with discharge rod plasma generator

Instead of the capillary plasma generator(CPG),a discharge rod plasma generator(DRPG)is used in the30 mm electrothermal-chemical(ETC)gun to improve the ignition uniformity of the solid propellant.An axisymmetric two-dimensional interior ballistics model of the solid propellant ETC gun(2D-IB-SPETCG)is presented to describe the process of the ETC launch.Both calculated pressure and projectile muzzle velocity accord well with the experimental results.The feasibility of the 2D-IB-SPETCG model is proved.Depending on the experimental data and initial parameters,detailed distribution of the ballistics parameters can be simulated.With the distribution of pressure and temperature of the gas phase and the propellant,the influence of plasma during the ignition process can be analyzed.Because of the radial flowing plasma,the propellant in the area of the DRPG is ignited within 0.01 ms,while all propellant in the chamber is ignited within 0.09 ms.The radial ignition delay time is much less than the axial delay time.During the ignition process,the radial pressure difference is less than 5 MPa at the place 0.025 m away from the breech.The radial ignition uniformity is proved.The temperature of the gas increases from several thousand K(conventional ignition)to several ten thousand K(plasma ignition).Compare the distribution of the density and temperature of the gas,we know that low density and high temperature gas appears near the exits of the DRPG,while high density and low temperature gas appears at the wall near the breech.The simulation of the 2D-IB-SPETCG model is an effective way to investigate the interior ballistics process of the ETC launch.The 2D-IB-SPETC model can be used for prediction and improvement of experiments.

Performance assessment of two-dimensional hydraulic models for generation of flood inundation maps in mountain river basins

Hydraulic models for the generation of flood inundation maps are not commonly applied in mountain river basins because of the difficulty in modeling the hydraulic behavior and the complex topography. This paper presents a comparative analysis of the performance of four twodimensional hydraulic models (HEC-RAS 2D, Iber 2D, Flood Modeller 2D, and PCSWMM 2D) with respect to the generation of flood inundation maps. The study area covers a 5-km reach of the Santa B-arbara River located in the Ecuadorian Andes, at 2330 masl, in Gualaceo. The model's performance was evaluated based on the water surface elevation and flood extent, in terms of the mean absolute difference and measure of fit. The analysis revealed that, for a given case, Iber 2D has the best performance in simulating the water level and inundation for flood events with 20- and 50-year return periods, respectively, followed by Flood Modeller 2D, HEC-RAS 2D, and PCSWMM 2D in terms of their performance. Grid resolution, the way in which hydraulic structures are mimicked, the model code, and the default value of the parameters are considered the main sources of prediction uncertainty.

Four-Beam Model for Vibration Analysis of a Cantilever Beam with an Embedded Horizontal Crack

As one of the main failure modes, embedded cracks occur in beam structures due to periodic loads. Hence it is useful to investigate the dynamic characteristics of a beam structure with an embedded crack for early crack detection and diagnosis. A new four-beam model with local flexibilities at crack tips is developed to investigate the transverse vibration of a cantilever beam with an embedded horizontal crack; two separate beam segments are used to model the crack region to allow opening of crack surfaces. Each beam segment is considered as an Euler-Bernoulli beam. The governing equations and the matching and boundary conditions of the four-beam model are derived using Hamilton＇s principle. The natural frequencies and mode shapes of the four-beam model are calculated using the transfer matrix method. The effects of the crack length, depth, and location on the first three natural frequencies and mode shapes of the cracked cantilever beam are investigated. A continuous wavelet transform method is used to analyze the mode shapes of the cracked cantilever beam. It is shown that sudden changes in spatial variations of the wavelet coefficients of the mode shapes can be used to identify the length and location of an embedded horizontal crack. The first three natural frequencies and mode shapes of a cantilever beam with an embedded crack from the finite element method and an experimental investigation are used to validate the proposed model. Local deformations in the vicinity of the crack tips can be described by the proposed four-beam model, which cannot be captured by previous methods.

Model tests and numerical analyses on horizontal impedance functions of inclined single piles embedded in cohesionless soil

Horizontal impedance functions of inclined single piles are measured experimentally for model soil-pile systems with both the effects of local soil nonlinearity and resonant characteristics.Two practical pile inclinations of 5掳 and 10掳 in addition to a vertical pile embedded in cohesionless soil and subjected to lateral harmonic pile head loadings for a wide range of frequencies are considered.Results obtained with low-to-high amplitude of lateral loadings on model soil-pile systems encased in a laminar shear box show that the local nonlinearities have a profound impact on the horizontal impedance functions of piles.Horizontal impedance functions of inclined piles are found to be smaller than the vertical pile and the values decrease as the angle of pile inclination increases.Distinct values of horizontal impedance functions are obtained for the ＇positive＇ and ＇negative＇ cycles of harmonic loadings,leading to asymmetric force-displacement relationships for the inclined piles.Validation of these experimental results is carried out through three-dimensional nonlinear finite element analyses,and the results from the numerical models are in good agreement with the experimental data.Sensitivity analyses conducted on the numerical models suggest that the consideration of local nonlinearity at the vicinity of the soil-pile interface influence the response of the soil-pile systems.

Modeling Arctic Intermediate Water: The effects of Neptune parameterization and horizontal resolution

Arctic Intermediate Water （AIW）, advected from the North Atlantic Ocean, has a potential influence on climate in the Arctic region, but is poorly simulated in coarse resolution models. In this study, a coupled ice-ocean model is used to investigate features of AIW by conducting two sensitivity experiments based on Neptune parameterization and horizontal resolution. The re- suits show that both experiments improve the modeling of temperature profiles in the western Eurasian Basin, mainly as a result of more realistic volume and heat transport through the Fram Strait. Topographical flows are well reproduced using Neptune parame- terization or a finer horizontal resolution. In the eddy-permitting model with relatively higher resolution, the velocity field is more realistic than in the Neptune parameterization model, and complex inflow and outflow belts of barotropic structure are well repro- duced. The findings of this study suggest that increased model resolution, as provided by an eddy-resolving model, is needed to reproduce realistic circulation and thermohaline structure in the Arctic, since the Rossby radius of deformation is only several kilometers in the Arctic Ocean. This paper focuses on the external heat input rather than internal mixing process, and obtains a conclusion that the heat input from the Fram Strait is a main factor to reproduce AIW in the Eurasian Basin successfully, at least for the western part.

Horizontal structure of convergent wind shear associated with sporadic E layers over East Asia

At present,the main detection instruments for observing sporadic E(Es)layers are ground-based radars,dense networks of ground-based global navigation satellite system(GNSS)receivers,and GNSS radio occultation,but they cannot capture the whole picture of the horizontal structure of Es layers.This study employs the Whole Atmosphere Community Climate Model with thermosphere and ionosphere eXtension model(WACCM-X 2.1)to derive the horizontal structure of the ion convergence region(HSICR)to explore the shapes of the large-scale Es layers over East Asia for the period from June 1 to August 31,2008.The simulation produced the various shapes of the HSICRs elongated in the northwest-southeast,northeast-southwest,or composed of individual small patches.The close connection between Es layer critical frequency(foEs)and vertical ion convergence indicates that the HSICR is a good candidate for revealing and explaining the horizontal structure of the large-scale Es layers.

Analysis of Wellbore Flow in Shale Gas Horizontal Wells

Theflow behavior of shale gas horizontal wells is relatively complex,and this should be regarded as the main reason for which conventional pipeflow models are not suitable to describe the related dynamics.In this study,numerical simulations have been conducted to determine the gas-liquid distribution in these wells.In particular,using the measuredflow pressure data related to 97 groups of shale gas wells as a basis,9 distinct pipeflow models have been assessed,and the models displaying a high calculation accuracy for different water-gas ratio(WGR)ranges have been identified.The results show that:(1)The variation law of WGR in gas well satisfies a power function relation.(2)The well structure is the main factor affecting the gas-liquid distribution in the wellbore.(3)The Beggs&Brill,Hagedorn&Brown and Gray models exhibit a high calculation accuracy.

HIGH RESOLUTION PARAMETRIC MODELLING FOR TWO-DIMENSIONAL RADAR TARGET USING PRONY ALGORITHM

On the conditions of low-resolution radar, a parametric model for two-dimensional radar target is described here according to the theory of electromagnetic scattering and the geometrical theory of diffraction. A high resolution estimation algorithm to extract the model parameters is also developed by building the relation of the scattering model and Prony model. The analysis of Cramer-Rao bound and simulation show that the method here has better statistical performance. The simulated analysis also indicates that the accurate extraction of the diffraction coefficient of scattering center is restricted by signal to noise ratio, radar center frequency and radar bandwidth.

Rock physics and seismic modeling ofshale reservoirs with horizontal fractures

The presence of horizontal fractures enhances seismic anisotropy of shales. Calculation based on the effective medium theory indicates that horizontal fractures have little effects on velocities along the direction parallel to fractures,but can significantly reduce velocities along the direction normal to fractures. Seismic responses of shales with horizontal fractures are calculated based on the reflector model and the anisotropic propagator matrix method,in which the reflections are a combination of the contrast in impedance due to the variations in fracture density,anisotropic propagation of waves within the shales,and the tuning and interferences associated with layer thickness. Calculated results indicate that seismic reflections are sensitive to reservoir layer thickness and fracture density. Anisotropic propagation alters amplitudes and phases of reflections. It corresponds to higher reflection amplitudes for the case of surrounding sandstone with higher velocity because the increase in fracture density increases the contrast in impedance between the shale and sandstone. In contrast,the surrounding sandstone with lower velocity corresponds to lower reflection amplitudes for the increase in fracture density.

Aseismic negative dislocation model and deformation analysis of crustal horizontal movement during 1999--2001 in the northeast margin of Qinghai-Xizang block

Through numerical simulation for GPS data, aseismic negative dislocation model for crustal horizontal movement during 1999~2001 in the northeast margin of Qinghai-Xizang block is presented, combined with the spatial distribution of apparent strain field in this area, the characteristics of motion and deformation of active blocks and their boundary faults, together with the place and intensity of strain accumulation are analyzed. It is shown that: a) 9 active blocks appeared totally clockwise motion from eastward by north to eastward by south. Obvious sinistral strike-slip and NE-NEE relative compressive motion between the blocks separated by Qilianshan-Haiyuan fault zone was discovered; b) 20 fault segments (most of them showed compression) locked the relative motion between blocks to varying degrees, among the total, the mid-east segment of Qilianshan fault (containing the place where it meets Riyueshan-Lajishan fault) and the place where it meets Haiyuan fault and Zhuanglanghe fault, more favored accumulation of strain. Moreover, the region where Riyueshan-Lajishan fault meets north boundary of Qaidam block may have strain accumulation to some degree. c) Obtained magnitude of block velocities and locking of their boundaries were less than relevant results for observation in the period of 1993~1999.

A Numerical Algorithm Based on Quadratic Finite Element for Two-Dimensional Nonlinear Time Fractional Thermal Diffusion Model

In this article,a high-order scheme,which is formulated by combining the quadratic finite element method in space with a second-order time discrete scheme,is developed for looking for the numerical solution of a two-dimensional nonlinear time fractional thermal diffusion model.The time Caputo fractional derivative is approximated by using the L2-1formula,the first-order derivative and nonlinear term are discretized by some second-order approximation formulas,and the quadratic finite element is used to approximate the spatial direction.The error accuracy O(h3+
t2)is obtained,which is verified by the numerical results.

Global Two-Dimensional Chemistry Model and Simulation of Atmospheric Chemical Composition

A global two-dimensional zonally averaged chemistry model is developed to study the chemi-cal composition of atmosphere. The region of the model is from 90°S to 90°N and from the ground to the altitude of 20 km with a resolution of 5° x 1 km. The wind field is residual circulation calcu-lated from diabatic rate. 34 species and 104 chemical and photochemical reactions are considered in the model. The sources of CH4, CO and NOx, which are divided into seasonal sources and non-seasonal sources, are parameterized as a function of latitude and time. The chemical composi-tion of atmosphere was simulated with emission level of CH4, CO and NOx in 1990. The results are compared with observations and other model results, showing that the model is successful to simu-late the atmospheric chemical composition and distribution of CH4. Key words Global two-dimensional chemistry model - Atmospheric composition - Emission This work was supported by the State Key Program for basic research “ Climate Dynamics and Cli-mate Prediction Theory” (Pandeng-yu-21).The authors would like to express their thanks to the National Oceanic and Atmospheric Administration (NOAA), Climate Monitoring and Diagnostics Laboratory (CMDL), Carbon Cycle Group for providing the observational data of CO and CH4.

Simulation logging experiment and interpretation model of array production logging measurements in a horizontal well

The distributions of local velocity and local phase holdup along the radial direction of pipes are complicated because of gravity differentiation,and the distribution of fluid velocity fi eld changes along the gravity direction in horizontal wells.Therefore,measuring the mixture flow and water holdup is difficult,resulting in poor interpretation accuracy of the production logging output profile.In this paper,oil–water two-phase flow dynamic simulation logging experiments in horizontal oil–water two-phase fl ow simulation wells were conducted using the Multiple Array Production Suite,which comprises a capacitance array tool(CAT)and a spinner array tool(SAT),and then the response characteristics of SAT and CAT in diff erent fl ow rates and water cut production conditions were studied.According to the response characteristics of CAT in diff erent water holdup ranges,interpolation imaging along the wellbore section determines the water holdup distribution,and then,the oil–water two-phase velocity fi eld in the fl ow section was reconstructed on the basis of the fl ow section water holdup distribution and the logging value of SAT and combined with the rheological equation of viscous fl uid,and the calculation method of the oil–water partial phase fl ow rate in the fl ow section was proposed.This new approach was applied in the experiment data calculations,and the results are basically consistent with the experimental data.The total fl ow rate and water holdup from the calculation are in agreement with the set values in the experiment,suggesting that the method has high accuracy.

Comprehensive two-dimensional river ice model based on boundary-fitted coordinate transformation method

River ice is a natural phenomenon in cold regions, influenced by meteorology, geomorphology, and hydraulic conditions. River ice processes involve complex interactions between hydrodynamic, mechanical, and thermal processes, and they are also influenced by weather and hydrologic conditions. Because natural rivers are serpentine, with bends, narrows, and straight reaches, the commonly-used one-dimensional river ice models and two-dimensional models based on the rectangular Cartesian coordinates are incapable of simulating the physical phenomena accurately. In order to accurately simulate the complicated river geometry and overcome the difficulties of numerical simulation resulting from both complex boundaries and differences between length and width scales, a two-dimensional river ice numerical model based on a boundary-fitted coordinate transformation method was developed. The presented model considers the influence of the frazil ice accumulation under ice cover and the shape of the leading edge of ice cover during the freezing process. The model is capable of determining the velocity field, the distribution of water temperature, the concentration distribution of frazil ice, the transport of floating ice, the progression, stability, and thawing of ice cover, and the transport, accumulation, and erosion of ice under ice cover. A MacCormack scheme was used to solve the equations numerically. The model was validated with field observations from the Hequ Reach of the Yellow River. Comparison of simulation results with field data indicates that the model is capable of simulating the river ice process with high accuracy.

Real-time forward modeling and inversion of logging-while-drilling electromagnetic measurements in horizontal wells

Based on the pseudo-analytical equation of electromagnetic log for layered formation,an optimal boundary match method is proposed to adaptively truncate the encountered formation structures.An efficient integral method is put forward to significantly accelerate the convergence of Sommerfeld integral.By asymptotically approximating and subtracting the first reflection/transmission waves from the scattered field,the new Sommerfeld integral method has addressed difficulties encountered by the traditional digital filtering method,such as low computational precision and limited operating range,and realized the acceleration of the computation speed of logging-while-drilling electromagnetic measurements(LWD EM).By making use of the priori information from the offset/pilot wells and interactively adjusting the formation model,the optimum initial guesses of the inversion model is determined in order to predict the nearby formation boundaries.The gradient optimization algorithm is developed and an interactive inversion system for the LWD EM data from the horizontal wells is established.The inverted results of field data demonstrated that the real-time interactive inversion method is capable of providing the accurate boundaries of layers around the wellbore from the LWD EM,and it will benefit the wellbore trajectory optimization and reservoir interpretation.