Acoustic field simulations of logging while drilling by cylindrical finite difference with variable grids

This study proposes an elastic finite difference(FD)time domain method with variable grids in three-dimensional cylindrical coordinates.The calculations will diverge and become less accurate by conventional cylindrical FD as the grid size gradually becomes more extensive with the increasing radius.To prevent grids from being too coarse in far fields,we compensate for the grid cell infl ation by refi ning the grid step in the azimuthal direction.The variable grid FD in the cylindrical coordinate systems has a higher effi ciency in solving acoustic logging while drilling(LWD)problems because the grid boundaries are consistent with those of the drill collar and the borehole.The proposed algorithm saves approximately 94%of the FD grids,80%of the computation time,and memory with a higher calculation accuracy than the FD on rectangular grids for the same models.We also calculate the acoustic LWD responses of the fl uid-fi lled borehole intersecting with fractures.Refl ections are generated at the fractures,which can be equivalent to an additional scattering source.The mode conversions between the collar and the Stoneley waves are revealed.The Stoneley spectra are more sensitive to the fracture.Finally,the logs in a heterogeneous formation with two refl ectors far from the borehole are modeled,and a means of estimating the azimuth of geological interfaces from refl ections is proposed.

Design of transit-time difference ultrasonic heat meter based on TDC-GP21

According to the transit-time difference measurement method,we proposed a design scheme of ultrasonic heat meter based on TDC-GP21.The measurement unit TDC-GP21 mainly completes transit-time measurement and inlet,outlet temperature measurement functions.Control unit and data processing unit based on MSP430F4152 of TI corporation complete functions including peripherals control,data analysis,temperature compensation algorithm,flow pattern compensation algorithm and low power consumption control.The design meets the Town Construction Professional Standard CJ 128-2007,and furthermore,some performances can be improved.

Non-contact Acoustic Contour Imaging Method Using Phase Difference in Extremely Shallow Underground

Buried object which exists at extremely shallow underground can be detected by using an acoustic vibration and a SLDV （scanning laser Doppler vibrometer）. Non-contact acoustic exploration can be realized by using air coupled sound. It was confirmed that discovery and the identification of a buried thing were possible using the OFR method （optimum frequency range method）. However, in this method, only the amplitude of the vibration speed spectrum is used and did not use the phase information. Therefore a new imaging method is proposed that used the phase information of the vibration speed spectrum. From the experimental results, the possibility of the outline extraction by using the phase difference between each scan point is confirmed. As a future task, the phase difference image will become clearer by changing the scan point density and long duration time of emitted wave.

3D numerical simulation in acoustic-elastic coupled media with staggered-grid finite-difference method

Acoustic-elastic coupled media is often encountered in most marine explorations, and accurate simulation of acoustic-elastic coupled media is of great significance. At present, the study of acoustic-elastic coupled media still assumes that the solid of the acoustic-elastic coupled media is isotropic, but this assumption is not in accordance with the actual situation. In this paper, we derive the solid media of acoustic-elastic coupled media from isotropic media to anisotropic media, and propose an acoustic-elastic coupled medium based ontransverse isotropic media with vertical symmetric axes(VTI) to improve the accuracy of forward modeling. Based on the relationship between the Thomsen parameter and the coefficient matrix of the anisotropic elastic wave equation, we transform the Thomson parameter into a velocity model with anisotropic properties. We use a staggered grid finite difference method to simulate the propagation of a wavefield in a three-dimensional acoustic-elastic coupled media. We obtain the snapshots of the wave field when the solid of the acoustic-elastic coupled media is an isotropic medium and a VTI media. When the solid of the acoustic-elastic coupled media is considered VTI media, we can observe the qP wave and qS wave that cannot be observed in the isotropic medium from the wave field snapshot. We can also find that the seismic records obtained by the method we use are more realistic. The algorithm proposed in this paper is of great significance for high-precision ocean numerical simulation.

Damage Processes of Polypropylene Fiber Reinforced Mortar in Different Fiber Content Revealed by Acoustic Emission Behavior

The performances of the cement-based materials can be improved by the incorporation of polypropylene fiber, but the damage processes become more complex with different fiber contents at the same time. The acoustic emission（AE） technology can achieve the global monitoring of internal damage in materials. The evolution process of failure mode and damage degree of polypropylene fiber reinforced mortar and concrete were analyzed by measuring the AE energy, RA value, AF value and b value. It was found that the cement matrix cracked on the initial stage, the cracks further developed on the medium stage and the fibers were pulled out on the last stage. The matrix cracked with minor injury cracks, but the fiber broke with serious damage cracks. The cumulative AE energy was proportional to the polypropylene fiber reinforced concrete and mortar＇s ductility. The damage mode and damage degree can be judged by identifying the damage stage obtained by the analysis of the AF value.

Influences of clean fracturing fluid viscosity and horizontal in-situ stress difference on hydraulic fracture propagation and morphology in coal seam

The viscosity of fracturing fluid and in-situ stress difference are the two important factors that affect the hydraulic fracturing pressure and propagation morphology. In this study, raw coal was used to prepare coal samples for experiments, and clean fracturing fluid samples were prepared using CTAB surfactant. A series of hydraulic fracturing tests were conducted with an in-house developed triaxial hydraulic fracturing simulator and the fracturing process was monitored with an acoustic emission instrument to analyze the influences of fracturing fluid viscosity and horizontal in-situ stress difference on coal fracture propagation. The results show that the number of branched fractures decreased, the fracture pattern became simpler, the fractures width increased obviously, and the distribution of AE event points was concentrated with the increase of the fracturing fluid viscosity or the horizontal in-situ stress difference. The acoustic emission energy decreases with the increase of fracturing fluid viscosity and increases with the increase of horizontal in situ stress difference. The low viscosity clean fracturing fluid has strong elasticity and is easy to be compressed into the tip of fractures, resulting in complex fractures. The high viscosity clean fracturing fluids are the opposite. Our experimental results provide a reference and scientific basis for the design and optimization of field hydraulic fracturing parameters.

Fracture identification based on remote detection acoustic reflection logging

Fracture identification is important for the evaluation of carbonate reservoirs. However, conventional logging equipment has small depth of investigation and cannot detect rock fractures more than three meters away from the borehole. Remote acoustic logging uses phase-controlled array-transmitting and long sound probes that increase the depth of investigation. The interpretation of logging data with respect to fractures is typically guided by practical experience rather than theory and is often ambiguous. We use remote acoustic reflection logging data and high-order finite-difference approximations in the forward modeling and prestack reverse-time migration to image fractures. First, we perform forward modeling of the fracture responses as a function of the fracture-borehole wall distance, aperture, and dip angle. Second, we extract the energy intensity within the imaging area to determine whether the fracture can be identified as the formation velocity is varied. Finally, we evaluate the effect of the fracture-borehole distance, fracture aperture, and dip angle on fracture identification.

Numerical simulation of an acoustic field generated by a phased arc array in a fluid-filled borehole

The acoustic tools widely used in borehole well logging and being developed in borehole acoustic reflection imaging do not have the function of azimuthal measurement due to a symmetric source, so they can not be used to evaluate the azimuthal character of borehole formation. In this paper, a 3D finite difference method was used to simulate the acoustic fields in a fluid-filled borehole generated by a traditional monopole source and a phased arc array. Acoustic waveforms were presented for both cases. The analysis of the simulated waveforms showed that different from the monopole source, the acoustic energy generated by the phased arc array transmitter mainly radiated to the borehole in a narrow azimuthal range, which was the key technique to implement azimuthal acoustic well logging. Similar to the monopole source, the waveforms generated by the phased arc array in the fluid-filled borehole also contain compressional （P） waves and shear （S） waves refracted along the borehole, which is the theoretical foundation of phased arc array acoustic well logging.

ACOUSTIC REFLECTION FROM A BAFFLED ELASTIC/PIEZOELECTRIC PLATE

A finite difference/boundary integral procedure to determine theacoustic reflected pressure from a fluid-loaded bi-laminated plate isdescribed. The bi-laminate is composed of a piezoelectric layer andan elastic layer in contact with the fluid, and is held by anacoustically hard baffle. In the numerical model, the fluid pressureat fluid/solid interface is replaced by a continuum of point sourcesweighted by the normal acceleration of the elastic plate, and thegoverning equation system is solved in the solid domain.

ACOUSTIC REFLECTION FROM A PERIODIC ELASTIC/PIEZOELECTRIC PLATE

The acoustic reflected pressure from a periodic elastic/piezoelectric laminated plate is studied for the purpose of acoustic reflection control.A finite difference/boundary integral procedure to determine the reflected pressure from the fluid-loaded plate is described.In the numerical model,a Green's function in the form of infinite sum is employed and a boundary integral is performed to replace the fluid pressure at fluid/solid interface by a continuum of point sources weighted by the normal acceleration of the elastic plate.The equation system is then solved only in the solid domain.It is demonstrated that an appropriate applied voltage potential across the piezoelectric layer has the effect of cancelling the fundamental propagating mode,and there,is no reflection for frequencies up to the cut-off frequency of the next propagating mode if the fundamental mode has been eliminated.

Full-waveform Velocity Inversion Based on the Acoustic Wave Equation

Full-waveform velocity inversion based on the acoustic wave equation in the time domain is investigated in this paper. The inversion is the iterative minimization of the misfit between observed data and synthetic data obtained by a numerical solution of the wave equation. Two inversion algorithms in combination with the CG method and the BFGS method are described respectively. Numerical computations for two models including the benchmark Marmousi model with complex structure are implemented. The inversion results show that the BFGS-based algorithm behaves better in inversion than the CG-based algorithm does. Moreover, the good inversion result for Marmousi model with the BFGS-based algorithm suggests the quasi-Newton methods can provide an important tool for large-scale velocity inversion. More computations demonstrate the correctness and effectives of our inversion algorithms and code.

Effect of the fluctuant acoustic channel on the gain of a linear array in the ocean waveguide

The inhomogenous ocean waveguide,which leads the amplitude and phase of the signal arriving at a hydrophone array to fluctuate,is one of the causes that make the array gain deviate from its ideal value.The relationship between the array gain and the fluctuant acoustic channel is studied theoretically.The analytical expression of the array gain is derived via an acoustic channel transfer function on the assumption that the ambient noise field is isotropic.The expression is expanded via the Euler formula to give an insight into the effect of the fluctuant acoustic channel on the array gain.The result demonstrates that the amplitude fluctuation of the acoustic channel transfer functions has a slight effect on the array gain;however,the uniformity of the phase difference between the weighting coefficient and the channel transfer function on all the hydrophones in the array is a major factor that leads the array gain to further deviate from its ideal value.The numerical verification is conducted in the downslope waveguide,in which the gain of a horizontal uniform linear array(HLA)with a wide-aperture operating in the continental slope area is considered.Numerical result is consistent with the theoretical analysis.

Explosions and seismic phenomena based on exciting of acoustic-electromagnetic waves

During earthquakes and strong underground explosions it is possible to observe two different effects. The first one is connected with the destruction of media, and this causes acoustic and later hybrid acoustic-electromagnetic waves in an epicenter in the atmosphere and in the ionosphere. Another one is connected with cracks in crystals of rocks, which seems more interesting, because it is possible to recognize the natural earthquakes and industrial explosions. In the first effects due to nonlinear elastic modules the acoustic waves move through the lithosphere and transform their spectra from VLF (very low frequencies ~ 1 - 10 kHz) at the depth of about 30 km into the lower part of ELF (extremely low frequencies, ~ 3 Hz - 1 kHz) on the Earth’s surface, then they pass the atmosphere and penetrate into the ionosphere. During the nonlinear acoustic passage through the atmosphere and the ionosphere, the spectrum transforms from ELF into ULF (ultra low frequencies, < 1 Hz) range. In this review article the classification of spectrum and analysis of two cases of the destruction of rocks in the lithosphere is presented. The rocks possess piezoelectric and piezomagnetic properties. In this case the electromagnetic emission is excited by the fracturing in plates of crystals. The difference of emission from piezoelectric and magnetic plates in cases of industrial explosions and natural seismic events including volcanic phenomena gives a possibility to analyze the method of its identification. The consideration is based on the model of the plate of a finite size with an uniformly moving crack.

Efficient solution of large-scale matrix of acoustic wave equations in 3D frequency domain

In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.

Numerical Simulation of Acoustic-Gravity Waves Propagation in a Heterogeneous Earth-Atmosphere Model with Wind in the Atmosphere

A numerical-analytical solution for seismic and acoustic-gravity waves propagation is applied to a heterogeneous “Earth-Atmosphere” model. Seismic wave propagation in an elastic half-space is described by a system of first order dynamic equations of elasticity theory. Propagation of acoustic-gravity waves in the atmosphere is described by the linearized Navier-Stokes equations with the wind. The algorithm proposed is based on the integral Laguerre transform with respect to time, the finite integral Fourier transform along the spatial coordinate with the finite difference solution of the reduced problem.

轴承故障声发射信号的涡旋聚焦传输方法

针对滚动轴承内部供电及电信号传输难题,研究了基于声发射的轴承智能感知技术。首先,基于故障轴承转动振动特性构建声信号时序相差涡旋传输模型。然后,设计了环形中空式声超表面结构,研究了声涡旋聚焦机制及超表面透射性能,并仿真分析了故障类型与声场涡旋拓扑荷数的关联关系。最后,通过故障轴承的模拟与试验装置验证了轴承故障声发射信号的涡旋聚焦传输方法。

Finite difference time domain analysis of two-dimensional surface acoustic wave piezoelectric phononic crystals at radio frequency

The finite-difference time-domain （FDTD） method is proposed for analyzing the surface acoustic wave （SAW） propagation in two-dimensional （2D） piezoelectric phononic crystals （PCs） at radio frequency （RF）, and also experiments are established to demonstrate its analysis result of the PCs＇ band gaps. The FDTD method takes the piezoelectric effect of PCs into account, in which periodic boundary conditions are used to decrease memory/time consumption and the perfectly matched layer boundary conditions are adopted as the SAW absorbers to attenuate artificial reflections. Two SAW delay lines are established with/without piezoelectric PCs located between interdigital transducers. By removing several echoes with window gating function in time domain, delay lines transmission function is achieved. The PCs＇ transmission functions and band gaps are obtained by comparing them in these two delay lines. When Aluminum/128°YX-LiNbO3 is adopted as scatter and substrate material, the PCs＇ band gap is calculated by this FDTD method and COMSOL respectively. Results show that computational results of FDTD agree well with experimental results and are better than that of COMSOL.

An optimized staggered variable-grid finite-difference scheme and its application in cross-well acoustic survey

In this paper, an optimized staggered variable-grid finite-difference (FD) method is developed in veloc- ity-stress elastic wave equations. On the basis of the dispersion-relation-preserving (DRP), a fourth-order finite-difference operator on non-uniform grids is constructed. The proposed algorithm is a continuous variable-grid method. It does not need interpolations for the field variables between re- gions with the fine spacing and the coarse one. The accuracy of the optimized scheme has been veri- fied with an analytical solution and a regular staggered-grid FD method for the eighth order accuracy in space. The comparisons of the proposed scheme with the variable-grid FD method based on Taylor series expansion are made. It is demonstrated that this optimized scheme has less dispersion errors than that with Taylor's series expansion. Thus, the proposed scheme uses coarser grids in numerical simulations than that constructed by the Taylor's series expansion. Finally, the capability of the opti- mized FD is demonstrated for a complex cross-well acoustic simulation. The numerical experiment shows that this method greatly saves storage requirements and computational time, and is stable.

VTI介质修正声学近似qP波波动方程与模拟

各向异性介质中纵横波通常耦合在一起传播,纵横波解耦是地震波传播理论研究的重要内容。经典的声学近似通过设置垂向qSV波速度V_(S0)为0来解耦qP波场,但是存在退化qSV波等问题。本研究将垂向qSV波速度V_(S0)考虑成波数k_(x)、k_(z)和各向异性参数ε、δ的函数,对声学近似进行修正,推导了VTI介质修正声学近似qP波频散关系和波动方程。VTI介质修正声学近似qP波波动方程包含椭圆项和非椭圆项两部分,因此采用混合有限差分/伪谱算法进行求解,即采用有限差分法求解椭圆项、伪谱法求解非椭圆项。频散关系分析和数值示例表明,基于修正声学近似的qP波波动方程不包含退化qSV波,是纯qP波方程,与弹性波方程模拟结果吻合较好,且具有较高精度;该方程在ε≥δ和ε<δ的VTI介质中均是稳定的,并且精度高于声学近似模拟结果。

基于深度神经网络的7065铝合金厚板应力检测模型

针对工业生产中传统超声应力检测法对铝合金厚板在不同拉伸率和不同温度条件下存在的测量误差的问题,以7065铝合金厚板为实验对象,提出一种在不同拉伸率和不同温度条件下的基于树突神经网络的应力预测模型与传统超声检测法融合的应力检测模型,然后使用改进的GSA-GRNN对该应力检测模型进行温度补偿。以南南铝公司生产的7065铝合金厚板为研究对象,使用恒温槽为超声检测提供恒温环境,分别对不同拉伸率、不同温度下的7065铝合金厚板进行超声检测,将声时差、拉伸率作为输入参数,应力作为输出参数,创建一个基于树突神经网络的应力检测模型,然后将应力检测模型的输出作为输入,使用改进的GSA-GRNN建立温度补偿模型对应力检测模型进行温度补偿。研究结果表明:融合了传统超声声时差的检测模型均方根误差为0.84636,相关系数为0.99743,和其他神经网络模型对比,该模型拥有更好的精度;在对该模型进行温度补偿后,模型的应力均方根误差和相关系数分别可以达到0.78848和0.99844,模型的精度得到了进一步的提升。证明基于数据驱动的神经网络融合传统超声检测可以有效降低检测误差,同时省去传统检测方法人工计算应力的时间,提高了检测效率。研究结果可以为基于数据驱动的应力检测模型提供进一步的优化参考。