Effect of stress-dependent microannulus aperture on well leakage

Debonding at the cement-casing interface is recognized as a principal failure mechanism leading to CO_(2) leakage in wells.This detachment gives rise to a microannulus,which notably possesses greater permeability than undamaged cement,undermining its sealing efficacy.Conventionally,the permeability of the microannulus is regarded as a uniform value throughout the well.However,fundamentally,a microannulus is one type of fracture,and its gap or aperture size is affected by the effective stress.In this work,we developed a unique experimental apparatus.This equipment facilitates the curing of cement inside a steel casing,the formation of a microannulus between the casing and the cement,and the investigation of the fluid flow dynamics along the microannulus under laboratory-replicated in situ conditions.The microannulus was formed by injecting fluid from one end of the setup,and receiving similar amount of fluid on the other end signified the development of the leakage channel.Additionally,strain gauges affixed to the casing’s external surface yielded key information on the microannulus’s opening and closure.We observed a noticeable decline in microannulus hydraulic aperture(or permeability)in relation to effective stress and an exponential equation fits their relationship.Our findings also indicate a distinct behavior when comparing liquid CO_(2) with water.Specifically,it is easier for liquid CO_(2) to create the microannulus.However,the hydraulic aperture range for this microannulus(0.7-6 mm)is considerably smaller than that created by water flow(2-17 mm).Finally,we integrated the stressdependent microannulus aperture size into the combined analysis of well mechanical integrity and well leakage.The outcomes consistently demonstrated that when factoring in the stress-dependent aperture sizes,the leakage rates are 3e5 times compared to a fixed aperture model.The traditional assumption of a constant aperture significantly underestimates fluid leakage risks.

A method to interpret fracture aperture of rock slope using adaptive shape and unmanned aerial vehicle multi-angle nap-of-the-object photogrammetry

The aperture of natural rock fractures significantly affects the deformation and strength properties of rock masses,as well as the hydrodynamic properties of fractured rock masses.The conventional measurement methods are inadequate for collecting data on high-steep rock slopes in complex mountainous regions.This study establishes a high-resolution three-dimensional model of a rock slope using unmanned aerial vehicle(UAV)multi-angle nap-of-the-object photogrammetry to obtain edge feature points of fractures.Fracture opening morphology is characterized using coordinate projection and transformation.Fracture central axis is determined using vertical measuring lines,allowing for the interpretation of aperture of adaptive fracture shape.The feasibility and reliability of the new method are verified at a construction site of a railway in southeast Tibet,China.The study shows that the fracture aperture has a significant interval effect and size effect.The optimal sampling length for fractures is approximately 0.5e1 m,and the optimal aperture interpretation results can be achieved when the measuring line spacing is 1%of the sampling length.Tensile fractures in the study area generally have larger apertures than shear fractures,and their tendency to increase with slope height is also greater than that of shear fractures.The aperture of tensile fractures is generally positively correlated with their trace length,while the correlation between the aperture of shear fractures and their trace length appears to be weak.Fractures of different orientations exhibit certain differences in their distribution of aperture,but generally follow the forms of normal,log-normal,and gamma distributions.This study provides essential data support for rock and slope stability evaluation,which is of significant practical importance.

Equivalent linear model for seismic damage evaluation of single-degree-of-freedom systems representing reinforced concrete structures considering cyclic degradation behavior

In this study,a novel equivalent damping ratio model that is suitable for reinforced concrete(RC)structures considering cyclic degradation behavior is developed,and a new equivalent linearization analysis method for implementing the proposed equivalent damping ratio model for use in seismic damage evaluation is presented.To this end,Ibarra’s peak-oriented model,which incorporates an energy-based degradation rule,is selected for representing hysteretic behavior of RC structure,and the optimized equivalent damping for predicting the maximum displacement response is presented by using the empirical method,in which the effect of cyclic degradation is considered.Moreover,the relationship between the hysteretic energy dissipation of the inelastic system and the elastic strain energy of the equivalent linear system is established so that the proposed equivalent linear system can be directly integrated with the Park-Ang seismic model to implement seismic damage evaluation.Due to the simplicity of the equivalent linearization method,the proposed method provides an efficient and reliable way of obtaining comprehensive insight into the seismic performance of RC structures.The verification demonstrates the validity of the proposed method.

High-resolution photogrammetry to measure physical aperture of two separated rock fracture surfaces

Photogrammetry,reconstructing three-dimensional(3D)models from overlapping two-dimensional(2D)photos,finds application in rock mechanics and rock engineering to extract geometrical details of reconstructed objects,for example rock fractures.Fracture properties are important for determining the mechanical stability,permeability,strength,and shear behavior of the rock mass.Photogrammetry can be used to reconstruct detailed 3D models of two separated rock fracture surfaces to characterize fracture roughness and physical aperture,which controls the fluid flow,hydromechanical and shear behavior of the rock mass.This research aimed to determine the optimal number of scale bars required to produce high-precision 3D models of a fracture surface.A workflow has been developed to define the physical aperture of a fracture using photogrammetry.Three blocks of Kuru granite(25 cm×25 cm×10 cm)with an artificially induced fracture,were investigated.For scaling 3D models,321 markers were used as ground control points(GCPs)with predefined distances on each block.When the samples were wellmatched in their original positions,the entire block was photographed.Coordinate data of the GCPs were extracted from the 3D model of the blocks.Each half was surveyed separately and georeferenced by GCPs and merged into the same coordinate system.Two fracture surfaces were extracted from the 3D models and the vertical distance between the two surfaces was digitally calculated as physical aperture.Accuracy assessment of the photogrammetric reconstruction showed a 20-30 mm digital control distance accuracy when compared to known distances defined between markers.To attain this accuracy,the study found that at least 200 scale bars were required.Furthermore,photogrammetry was employed to measure changes in aperture under normal stresses.The results obtained from this approach were found to be in good agreement with those obtained using linear variable displacement transducers(LVDTs),with differences ranging from 1 mm to 8μm.

Synthetic Aperture Positioning: A Review

Obtaining precise position of interested emitters passively has wide applications in both civilian and military fields.Different from traditional parameter measurement and direct position determination(DPD)method,recently a new passive localization method based on synthetic aper-ture technique,named synthetic aperture positioning(SAP),has been proposed.The method com-pensates for the nonlinear phase produced by relative motion between the moving platform and the emitter,achieving coherent summation of intercepted signals.The SAP can obtain high-resolution and high-precision localization results at a low signal-to-noise ratio.This paper summarizes the research progress of SAP,including localization principles,spaceborne applications,and application scope analysis.Besides,the possible future outlook of SAP is considered.

High-Precision DOA Estimation Method Based on Synthetic Aperture Technique

The existing direction-of-arrival(DOA)estimation methods only utilize the current received signals,which are susceptible to noise.In this paper,a method for DOA estimation based on a motion platform is proposed to achieve high-precision DOA estimation by utilizing past and present signals.The concept of synthetic aperture is introduced to construct a linear DOA estima-tion model.A DOA fine-tuning method based on the linear model is proposed to eliminate the lin-ear DOA variation,achieving a non-coherent accumulation of DOA estimations.Moreover,the baseband modulation and the phase modulation caused by the range history are compensated to achieve the coherent accumulation of all the DOA estimations.Simulation results show that the proposed method can significantly improve the DOA estimated accuracy at low signal-to-noise ratios(SNR).

An efficient approach for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations

An efficient approach is proposed for the equivalent linearization of frame structures with plastic hinges under nonstationary seismic excitations.The concentrated plastic hinges,described by the Bouc-Wen model,are assumed to occur at the two ends of a linear-elastic beam element.The auxiliary differential equations governing the plastic rotational displacements and their corresponding hysteretic displacements are replaced with linearized differential equations.Then,the two sets of equations of motion for the original nonlinear system can be reduced to an expanded-order equivalent linearized equation of motion for equivalent linear systems.To solve the equation of motion for equivalent linear systems,the nonstationary random vibration analysis is carried out based on the explicit time-domain method with high efficiency.Finally,the proposed treatment method for initial values of equivalent parameters is investigated in conjunction with parallel computing technology,which provides a new way of obtaining the equivalent linear systems at different time instants.Based on the explicit time-domain method,the key responses of interest of the converged equivalent linear system can be calculated through dimension reduction analysis with high efficiency.Numerical examples indicate that the proposed approach has high computational efficiency,and shows good applicability to weak nonlinear and medium-intensity nonlinear systems.

A high-contrast imaging coronagraph for segmented-mirror large aperture telescopes using a spatial light modulator

The primary mirrors of current and future large telescopes always employ a segmented mirror configuration.The small but non-negligible gaps between neighboring segments cause additional diffraction,which restricts the performance of high-contrast coronagraph.To solve this problem,we propose a coronagraph system based on a single liquid crystal spatial light modulator(SLM).This spatial light modulator is used for amplitude apodization,and its feasibility and potential performance are demonstrated using a laboratory setup using the stochastic parallel gradient descent(SPGD)algorithm to control the spatial light modulator,which is based on point spread function(PSF)sensing and evaluation and optimized for maximum contrast in the discovery working area as a merit function.The system delivers a contrast in the order of 10−6,and shows excellent potential to be used in current and future large aperture telescopes,both on the ground and in space.

Modulated-ISRJ rejection using online dictionary learning for synthetic aperture radar imagery

In electromagnetic countermeasures circumstances,synthetic aperture radar(SAR)imagery usually suffers from severe quality degradation from modulated interrupt sampling repeater jamming(MISRJ),which usually owes considerable coherence with the SAR transmission waveform together with periodical modulation patterns.This paper develops an MISRJ suppression algorithm for SAR imagery with online dictionary learning.In the algorithm,the jamming modulation temporal properties are exploited with extracting and sorting MISRJ slices using fast-time autocorrelation.Online dictionary learning is followed to separate real signals from jamming slices.Under the learned representation,time-varying MISRJs are suppressed effectively.Both simulated and real-measured SAR data are also used to confirm advantages in suppressing time-varying MISRJs over traditional methods.

Analysis of secondary makeup characteristics of drill collar joint and prediction of downhole equivalent impact torque of Well SDTK1

Based on the three-dimensional elastic-plastic finite element analysis of the 8"(203.2 mm)drill collar joint,this paper studies the mechanical characteristics of the pin and box of NC56 drill collar joints under complex load conditions,as well as the downhole secondary makeup features,and calculates the downhole equivalent impact torque with the relative offset at the shoulder of internal and external threads.On the basis of verifying the correctness of the calculation results by using measured results in Well GT1,the prediction model of the downhole equivalent impact torque is formed and applied in the first extra-deep well with a depth over 10000 m in China(Well SDTK1).The results indicate that under complex loads,the stress distribution in drill collar joints is uneven,with relatively higher von Mises stress at the shoulder and the threads close to the shoulder.For 203.2 mm drill collar joints pre-tightened according to the make-up torque recommended by American Petroleum Institute standards,when the downhole equivalent impact torque exceeds 65 kN·m,the preload balance of the joint is disrupted,leading to secondary make-up of the joint.As the downhole equivalent impact torque increases,the relative offset at the shoulder of internal and external threads increases.The calculation results reveal that there exists significant downhole impact torque in Well SDTK1 with complex loading environment.It is necessary to use double shoulder collar joints to improve the impact torque resistance of the joint or optimize the operating parameters to reduce the downhole impact torque,and effectively prevent drilling tool failure.

An improved algorithm based on equivalent sound velocity profile method at large incident angle

With the development of ultra-wide coverage technology,multibeam echo-sounder(MBES)system has put forward higher requirements for localization accuracy and computational efficiency of ray tracing method.The classical equivalent sound speed profile(ESSP)method replaces the measured sound velocity profile(SVP)with a simple constant gradient SVP,reducing the computational workload of beam positioning.However,in deep-sea environment,the depth measurement error of this method rapidly increases from the central beam to the edge beam.By analyzing the positioning error of the ESSP method at edge beam,it is discovered that the positioning error increases monotonically with the incident angle,and the relationship between them could be expressed by polynomial function.Therefore,an error correction algorithm based on polynomial fitting is obtained.The simulation experiment conducted on an inclined seafloor shows that the proposed algorithm exhibits comparable efficiency to the original ESSP method,while significantly improving bathymetry accuracy by nearly eight times in the edge beam.

Real-Time Fracture Aperture Identification Using Mud Loss Data and Solution for LCM Combination

Managing server lost circulation is a major challenge of drilling operation in naturally fractured formations and it causes much nonproductive rig time. When encountered with loss, the fracture aperture intersecting the wellbore is not well-identified in time, which has a significant impact on the decision of drilling operation and the undesired result of loss curing. Therefore, the onset of fracture is identified in a timely manner and evaluated comprehensively to formulate an appropriate strategy over time. However, the mud loss date, which is the primary source of information retrieved from the drilling process, was not properly used in real-time prediction of fracture aperture. This article provides a detailed mathematical study to discuss the mechanism of mud invasion in the near-wellbore region and prediction of fracture aperture. The fracture aperture can be calculated from mud-loss data by solving a cubic equation with input parameters given by the well radius, the overpressure ratio, and the maximum mud-loss volume. It permits the proper selection of loss-circulation material (LCM) with respect to particle size distribution and fiber usage. The case study illustrates the applicability of this methodology with a discussion on LCM particle distribution in different scenarios and the result demonstrates the outcome of inappropriate LCM usage and the advantages of the novel fiber-based LCM treatment.

Application of improved equivalent edge currents in synthetic aperture radar imaging

Target modeling and scattering function calculating are important prerequisites and groundwork for the synthetic aperture radar（SAR） imaging simulation.According to the difficult problems that normal methods cannot calculate the scattering function of electrically large object under the condition to wideband,an effective method of improved equivalent edge currents is presented and applied to SAR imaging simulation for the first time.This method improves calculating velocity and has relatively high precision.The concrete steps of applying the method are given.By way of the simulation experiment,the effectiveness of the method is verified.

Development of short prototype of dual aperture quadrupole magnet for CEPC ring

Main quadrupole magnets are critical for the Circular Electron and Positron Collider(CEPC)and are specifically designed as dual aperture quadrupole(DAQ)magnets.However,the field crosstalk between the two apertures presents challenges.As the CEPC will work at four beam energies of Z,W,Higgs and ttbar mode,the DAQ magnets will operate at four field gradients spanning from 3.18 to 12.63 T/m.The first short quadrupole magnet prototype with the bore diameter of 76 mm and magnetic length of 1.0 m revealed the problems of large magnetic field harmonics and a magnetic center shift within the beam energy range.Accordingly,a compensation method was proposed in this work to solve the field crosstalk effect.By adjusting the gap height at the middle of the two apertures,the field harmonics and magnetic center shift are significantly reduced.After optimization,the short prototype was modified using a new scheme.The field simulations are validated from the magnetic measurement results.Further,the multipole field meets the requirements of the four beam energies.The detailed magnetic field optimization,field harmonics adjustment,and measurement results are presented herein.

An abnormal carbody swaying of intercity EMU train caused by low wheel-rail equivalent conicity and damping force unloading of yaw damper

Low-frequency carbody swaying phenomenon often occurs to railway vehicles due to hunting instability,which seriously deteriorates the ride comfort of passengers.This paper investigates low-frequency carbody swaying through experimental analysis and numerical simulation.In the tests,the carbody acceleration,the wheel-rail profiles,and the dynamic characteristics of dampers were measured to understand the characteristics of the abnormal carbody vibration and to find out its primary contributor.Linear and nonlinear numerical simulations on the mechanism and optimization measures were carried out to solve this carbody swaying issue.The results showed that the carbody swaying is the manifest of carbody hunting instability.The low equivalent conicity and the decrease of dynamic damping of the yaw damper are probably the cause of this phenomenon.The optimization measures to increase the equivalent conicity and dynamic damping of the yaw damper were put forward and verified by on-track tests.The results of this study could enrich the knowledge of carbody hunting and provide a reference for solving abnormal carbody vibrations.

Force Analysis of the Overconstrained Mechanisms Based on Equivalent Stiffness Considering Limb Axial Deformation

Overconstrained mechanism has the advantages of large bearing capacity and high motion reliability,but its force analysis is complex and difficult because the mechanism system contains overconstraints.Considering the limb axial deformation,taking typical 2SS+P and 7-SS passive overconstrained mechanisms,2SPS+P and 7-SPS active overconstrained mechanisms,and 2SPS+P and 7-SPS passive-input overconstrained mechanisms as examples,a new force analysis method based on the idea of equivalent stiffness is proposed.The equivalent stiffness matrix of passive overconstrained mechanism is derived by combining the force balance and deformation compatibility equations with consideration of axial elastic limb deformations.The relationship between the constraint wrench magnitudes and the external force,limb stiffness is established.The equivalent stiffness matrix of active overconstrained mechanism is derived by combining the force balance and displacement compatibility equations.Here,the relationship between the magnitudes of the actuated wrenches and the external force,limb stiffness is investigated.Combining with the equivalent stiffness of the passive overconstrained mechanism,an analytical relationship between the actuated forces of passive-input overconstrained mechanism and the output displacement,limb stiffness is explored.Finally,adaptability of the equivalent stiffness to overconstrained mechanisms is discussed,and the effect of the limb stiffness on overconstrained mechanisms force distribution is revealed.The research results provide a theoretical reference for the design,research and practical application of overconstrained mechanism.

Generalized-Extended-State-Observer and Equivalent-Input-Disturbance Methods for Active Disturbance Rejection: Deep Observation and Comparison

Active disturbance-rejection methods are effective in estimating and rejecting disturbances in both transient and steady-state responses.This paper presents a deep observation on and a comparison between two of those methods:the generalized extended-state observer(GESO)and the equivalent input disturbance(EID)from assumptions,system configurations,stability conditions,system design,disturbance-rejection performance,and extensibility.A time-domain index is introduced to assess the disturbance-rejection performance.A detailed observation of disturbance-suppression mechanisms reveals the superiority of the EID approach over the GESO method.A comparison between these two methods shows that assumptions on disturbances are more practical and the adjustment of disturbance-rejection performance is easier for the EID approach than for the GESO method.

Retrieving chlorophyll content and equivalent water thickness of Moso bamboo(Phyllostachys pubescens) forests under Pantana phyllostachysae Chao-induced stress from Sentinel-2A/B images in a multiple LUTs-based PROSAIL framework

Biochemical components of Moso bamboo(Phyllostachys pubescens)are critical to physiological and ecological processes and play an important role in the material and energy cycles of the ecosystem.The coupled PROSPECT with SAIL(PROSAIL)radiative transfer model is widely used for vegetation biochemical component content inversion.However,the presence of leaf-eating pests,such as Pantana phyllostachysae Chao(PPC),weakens the performance of the model for estimating biochemical components of Moso bamboo and thus must be considered.Therefore,this study considered pest-induced stress signals associated with Sentinel-2A/B images and field data and established multiple sets of biochemical canopy reflectance look-up tables(LUTs)based on the PROSAIL framework by setting different parameter ranges according to infestation levels.Quantitative inversions of leaf area index(LAI),leaf chlorophyll content(LCC),and leaf equivalent water thickness(LEWT)were derived.The scale conversions from LCC to canopy chlorophyll content(CCC)and LEWT to canopy equivalent water thickness(CEWT)were calculated.The results showed that LAI,CCC,and CEWT were inversely related with PPC-induced stress.When applying multiple LUTs,the p-values were<0.01;the R2 values for LAI,CCC,and CEWT were 0.71,0.68,and 0.65 with root mean square error(RMSE)(normalized RMSE,NRMSE)values of 0.38(0.16),17.56μg cm-2(0.20),and 0.02 cm(0.51),respectively.Compared to the values obtained for the traditional PROSAIL model,for October,R2 values increased by 0.05 and 0.10 and NRMSE decreased by 0.09 and 0.02 for CCC and CEWT,respectively and RMSE decreased by 0.35μg cm-2 for CCC.The feasibility of the inverse strategy for integrating pest-induced stress factors into the PROSAIL model,while establishing multiple LUTs under different pest-induced damage levels,was successfully demonstrated and can potentially enhance future vegetation parameter inversion and monitoring of bamboo forest health and ecosystems.

Aperture shape optimization in intensity-modulated radiation therapy planning

The gradient element of the aperture gradient map is utilized directly to generate the aperture shape without modulation.This process can be likened to choosing the direction of negative gradient descent for the generic aperture shape optimiza-tion.The negative gradient descent direction is more suitable under local conditions and has a slow convergence rate.To overcome these limitations,this study introduced conjugate gradients into aperture shape optimization based on gradient modulation.First,the aperture gradient map of the current beam was obtained for the proposed aperture shape optimiza-tion method,and the gradients of the aperture gradient map were modulated using conjugate gradients to form a modulated gradient map.The aperture shape was generated based on the modulated gradient map.The proposed optimization method does not change the optimal solution of the original optimization problem,but changes the iterative search direction when generating the aperture shape.The performance of the proposed method was verified using cases of head and neck cancer,and prostate cancer.The optimization results indicate that the proposed optimization method better protects the organs at risk and rapidly reduces the objective function value by ensuring a similar dose distribution to the planning target volume.Compared to the contrasting methods,the normal tissue complication probability obtained by the proposed optimization method decreased by up to 4.61%,and the optimization time of the proposed method decreased by 5.26%on average for ten cancer cases.The effectiveness and acceleration of the proposed method were verified through comparative experiments.According to the comparative experiments,the results indicate that the proposed optimization method is more suitable for clinical applications.It is feasible for the aperture shape optimization involving the proposed method.

Wide-spectrum optical synthetic aperture imaging via spatial intensity interferometry

High resolution imaging is achieved using increasingly larger apertures and successively shorter wavelengths.Optical aperture synthesis is an important high-resolution imaging technology used in astronomy.Conventional long baseline amplitude interferometry is susceptible to uncontrollable phase fluctuations,and the technical difficulty increases rapidly as the wavelength decreases.The intensity interferometry inspired by HBT experiment is essentially insensitive to phase fluctuations,but suffers from a narrow spectral bandwidth which results in a lack of effective photons.In this study,we propose optical synthetic aperture imaging based on spatial intensity interferometry.This not only realizes diffraction-limited optical aperture synthesis in a single shot,but also enables imaging with a wide spectral bandwidth,which greatly improves the optical energy efficiency of intensity interferometry.And this method is insensitive to the optical path difference between the sub-apertures.Simulations and experiments present optical aperture synthesis diffraction-limited imaging through spatial intensity interferometry in a 100 nm spectral width of visible light,whose maximum optical path difference between the sub-apertures reaches 69λ.This technique is expected to provide a solution for optical aperture synthesis over kilometer-long baselines at optical wavelengths.