An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

The First Global Map of Atmospheric Ammonia(NH_(3)) as Observed by the HIRAS/FY-3D Satellite

Atmospheric ammonia(NH_(3)) is a chemically active trace gas that plays an important role in the atmospheric environment and climate change. Satellite remote sensing is a powerful technique to monitor NH_(3) concentration based on the absorption lines of NH_(3) in the thermal infrared region. In this study, we establish a retrieval algorithm to derive the NH_(3)column from the Hyperspectral Infrared Atmospheric Sounder(HIRAS) onboard the Chinese Feng Yun(FY)-3D satellite and present the first atmospheric NH_(3) column global map observed by the HIRAS instrument. The HIRAS observations can well capture NH_(3) hotspots around the world, e.g., India, West Africa, and East China, where large NH_(3) emissions exist. The HIRAS NH_(3) columns are also compared to the space-based Infrared Atmospheric Sounding Interferometer(IASI)measurements, and we find that the two instruments observe a consistent NH_(3) global distribution, with correlation coefficient(R) values of 0.28–0.73. Finally, some remaining issues about the HIRAS NH_(3) retrieval are discussed.

Numerical Study of the Biomechanical Behavior of a 3D Printed Polymer Esophageal Stent in the Esophagus by BP Neural Network Algorithm

Esophageal disease is a common disorder of the digestive system that can severely affect the quality of life andprognosis of patients. Esophageal stenting is an effective treatment that has been widely used in clinical practice.However, esophageal stents of different types and parameters have varying adaptability and effectiveness forpatients, and they need to be individually selected according to the patient’s specific situation. The purposeof this study was to provide a reference for clinical doctors to choose suitable esophageal stents. We used 3Dprinting technology to fabricate esophageal stents with different ratios of thermoplastic polyurethane (TPU)/(Poly-ε-caprolactone) PCL polymer, and established an artificial neural network model that could predict the radial forceof esophageal stents based on the content of TPU, PCL and print parameter. We selected three optimal ratios formechanical performance tests and evaluated the biomechanical effects of different ratios of stents on esophagealimplantation, swallowing, and stent migration processes through finite element numerical simulation and in vitrosimulation tests. The results showed that different ratios of polymer stents had different mechanical properties,affecting the effectiveness of stent expansion treatment and the possibility of postoperative complications of stentimplantation.

3D forward modeling and response characteristics of low-resistivity overburden of the CFS-PML absorbing boundary for ground-well transient electromagnetic method

This study used the stable and convergent Dufort-Frankel method to differentially discretize the diffusion equation of the ground-well transient electromagnetic secondary field.The absorption boundary condition of complex frequency-shifted perfectly matched layer(CFS-PML)was used for truncation so that the low-frequency electromagnetic wave can be better absorbed at the model boundary.A typical three-dimensional(3D)homogeneous half-space model was established and a low-resistivity cube model was analyzed under the half-space condition.The response patterns and drivers of the low-resistivity cube model were discussed under the influence of a low-resistivity overburden.The absorption boundary conditions of CFS-PML significantly affected the low-frequency electromagnetic waves.For a low-resistivity cube around the borehole,its response curve exhibited a single-peak,and the extreme point of the curve corresponded to the center of the low-resistivity body.When the low-resistivity cube was directly below the borehole,the response curve showed three extreme values(two high and one low),with the low corresponding to the center of the low-resistivity body.The total field response of the low-resistivity overburden was stronger than that of the uniform half-space model due to the low-resistivity shielding effect of electromagnetic waves.When the receiving-transmitting distance gradually increased,the effect of the low-resistivity overburden was gradually weakened,and the response of the low-resistivity cube was strengthened.It was affected by the ratio of the overburden resistivity to the resistivity of the low-resistivity body.

A 3D shell-like approach using element-free Galerkin method for analysis of thin and thick plate structures

A new efficient meshless method based on the element-free Galerkin method is proposed to analyze the static deformation of thin and thick plate structures in this paper. Using the new 3D shell-like kinematics in analogy to the solid-shell concept of the finite element method, discretization is carried out by the nodes located on the upper and lower surfaces of the structures. The approximation of all unknown field variables is carried out by using the moving least squares （MLS） approximation scheme in the in-plane directions, while the linear interpolation is applied through the thickness direction. Thus, different boundary conditions are defined only using displacements and penalty method is used to enforce the essential boundary conditions. The constrained Galerkin weak form, which incorporates only dis- placement degrees of freedom （d.o.f.s）, is derived. A modified 3D constitutive relationship is adopted in order to avoid or eliminate some self-locking effects. The numeric efficiency of the proposed meshless formulation is illustrated by the numeric examples.

杨梅素-3-O-β-D-乳糖的合成方法优化

杨梅素是一种主要用于抗炎的多酚类黄酮化合物,为改善其水溶性和稳定性,合成了杨梅素-3-O-β-D-乳糖,并开发了一条高效的合成路线。首先以杨梅苷为起始原料,对杨梅苷的5位、7位、3′位、4′位和5′位酚羟基进行保护,形成五苄基保护的杨梅苷。然后脱除3位鼠李糖,再与乙酰基保护的乳糖溴苷反应,最后通过依次脱去乙酰基和苄基得到目标化合物。经过路径优化,最终以更高效的方法合成了杨梅素-3-O-β-D-乳糖,减少了副产物的生成,提高了原合成路线总收率,此方法的总产率为44.5%。该路径得到了2个新型的杨梅素衍生物中间体。目标化合物的结构均经过1 H NMR,13 C NMR和HR-MS(ESI)确证。

Reconstruction of 3-D digital cores using a hybrid method

A 3-D 数字核心描述岩石的毛孔空间微观结构。X 光检查微 CT 扫描是最精确、直接却昂贵的方法获得 3-D 数字核心。在这研究，我们建议联合沉积模拟的一个混合方法并且退火模仿(SA ) 产生 3-D 的方法数字核心基于岩石的 2-D 图象。有沉积模拟的方法开始造 3-D 数字核心，它是为 SA 方法的起始的配置。我们用 SA 方法更新起始的数字核心匹配 2-D 岩石图象的自相关功能并且最后造最后的 3-D 数字核心。与典型 SA 方法相比，混合方法显著地减少了计算时间。本地孔理论被用于份量上比较重建的 3-D 有 X 光检查微 CT 3-D 图象的数字核心。结果显示 3-D 数字核心由混合方法重建了有同质和 X 光检查微 CT 图象的类似于那些的几何连接。形成因素和重建的 3-D 的渗透数字核心用有限元素方法被估计(女性) 并且格子 Boltzmann 方法(LBM ) 分别地。模仿的结果在对试验性的大小的好同意。模拟结果的比较建议数字核心由更多仔细独自比典型 SA 方法反映真运输性质的混合方法重建了。

High Efficient Methods of Content-based 3D Model Retrieval

Content-based 3D model retrieval is of great help to facilitate the reuse of existing designs and to inspire designers during conceptual design. However, there is still a gap to apply it in industry due to the low time efficiency. This paper presents two new methods with high efficiency to build a Content-based 3D model retrieval system. First, an improvement is made on the ＂Shape Distribution （D2）＂ algorithm, and a new algorithm named ＂Quick D2＂ is proposed. Four sample 3D mechanical models are used in an experiment to compare the time cost of the two algorithms. The result indicates that the time cost of Quick D2 is much lower than that of D2, while the descriptors extracted by the two algorithms are almost the same. Second, an expandable 3D model repository index method with high performance, namely, RBK index, is presented. On the basis of RBK index, the search space is pruned effectively during the search process, leading to a speed up of the whole system. The factors that influence the values of the key parameters of RBK index are discussed and an experimental method to find the optimal values of the key parameters is given. Finally, ＂3D Searcher＂, a content-based 3D model retrieval system is developed. By using the methods proposed, the time cost for the system to respond one query online is reduced by 75% on average. The system has been implemented in a manufacturing enterprise, and practical query examples during a case of the automobile rear axle design are also shown. The research method presented shows a new research perspective and can effectively improve the content-based 3D model retrieval efficiency.

Numerical simulation of standing wave with 3D predictor-corrector finite difference method for potential flow equations

A three-dimensional （3D） predictor-corrector finite difference method for standing wave is developed. It is applied to solve the 3D nonlinear potential flow equa- tions with a free surface. The 3D irregular tank is mapped onto a fixed cubic tank through the proper coordinate transform schemes. The cubic tank is distributed by the staggered meshgrid, and the staggered meshgrid is used to denote the variables of the flow field. The predictor-corrector finite difference method is given to develop the difference equa- tions of the dynamic boundary equation and kinematic boundary equation. Experimental results show that, using the finite difference method of the predictor-corrector scheme, the numerical solutions agree well with the published results. The wave profiles of the standing wave with different amplitudes and wave lengths are studied. The numerical solutions are also analyzed and presented graphically.

3-D Geological Modeling-Concept,Methods and Key Techniques

3-D geological modeling plays an increasingly important role in Petroleum Geology, Mining Geology and Engineering Geology. The complexity of geological conditions requires different modeling methods in different situations. This paper summarizes the general concept of geological modeling; compares the characteristics of borehole-based modeling, cross-section based modeling and multi- source interactive modeling; analyses key techniques in 3-D geological modeling; and highlights the main difficulties and directions of future studies.

3D elastic waveform modeling with an optimized equivalent staggered-grid finite-difference method

Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.

Effect of heterogeneity on mechanical and micro-seismic behaviors of sandstone subjected to multi-level cyclic loading: A discrete element method investigation

In numerical simulation of the mechanical responses and acoustic emission(AE)characteristics of rocks under cyclic loading,the impacts of compositional heterogeneities of mineral grains have barely been considered.This will lead to a poor reproduction of rock’s behaviors in terms of stress-strain relationship and micro-seismic characteristics in numerical simulation.This work aims to analyze and reveal the impact of parameter heterogeneity on the rock’s fatigue and micro-seismic properties based on PFC3D.Two distribution patterns(uniform and Weibull distributions),are implemented to assign four critical parameters(i.e.tensile strength,cohesion,parallel bond stiffness and linear stiffness)for 32 sets of numerical schemes.The results show that the models with high heterogeneity of tensile strength and cohesion can better reproduce the stress-strain relationship as well as the patterns of cumulative AE counts and energy magnitude.The evolution of the proportion of three-level AE events in the laboratory test is consistent with the numerical results when the highly heterogeneous tensile strength and cohesion are distributed.The numerical results can provide practical guidance to the PFC-based modeling of rock heterogeneity when exposed to multi-level cyclic loading and AE monitoring.

3D Free-Standing Carbon Nanofibers Modified by Lithiophilic Metals Enabling Dendrite-Free Anodes for Li Metal Batteries

Li metal with high-energy density is considered as the most promising anode for the next-generation rechargeable Li metal batteries;however,the growth of Li dendrites seriously hinders its practical application.Herein,3D free-standing carbon nanofibers modified by lithiophilic metal particles(CNF/Me,Me=Sn,Fe,Co)are obtained in situ by the electrospinning method.Benefiting from the lithophilicity,the CNF/Me composite may effectively prevent the formation of Li dendrites in the Li metal batteries.The optimized CNF/Sn–Li composite electrode exhibits a stable cycle life of over 2350 h during Li plating/stripping.When matched with typical commercial LiFePO_(4)(LFP)cathode,the LFP//CNF/Sn–Li full cell presents a high initial discharge specific capacity of 139 mAh g^(−1)at 1 C,which remains at 146 mAh g^(−1)after 400 cycles.When another state-of-the-art commercial LiNi_(0.8)Co_(0.1)Mn_(0.1)O_(2)(NCM(811))cathode is used,the assembled NCM//CNF/Sn–Li full cell shows a large initial specific discharge capacity of 206 mAh g^(−1)at substantially enhanced 10 C,which keeps at the good capacity of 99 mAh g^(−1)after 300 cycles.These results are greatly superior to the counterparts with Li as the anodes,indicating the great potential for practical utilization of the advanced CNF/Sn–Li electrode.

3D Model Retrieval Method Based on Affinity Propagation Clustering

In order to improve the accuracy and efficiency of 3D model retrieval,the method based on affinity propagation clustering algorithm is proposed. Firstly,projection ray-based method is proposed to improve the feature extraction efficiency of 3D models. Based on the relationship between model and its projection,the intersection in 3D space is transformed into intersection in 2D space,which reduces the number of intersection and improves the efficiency of the extraction algorithm. In feature extraction,multi-layer spheres method is analyzed. The two-layer spheres method makes the feature vector more accurate and improves retrieval precision. Secondly,Semi-supervised Affinity Propagation ( S-AP) clustering is utilized because it can be applied to different cluster structures. The S-AP algorithm is adopted to find the center models and then the center model collection is built. During retrieval process,the collection is utilized to classify the query model into corresponding model base and then the most similar model is retrieved in the model base. Finally,75 sample models from Princeton library are selected to do the experiment and then 36 models are used for retrieval test. The results validate that the proposed method outperforms the original method and the retrieval precision and recall ratios are improved effectively.

Inner damage identification and residual strength assessment of a 3D printed tunnel with marble-like cementitious materials using piezoelectric transducers

Quantitative damage identification of surrounding rock is important to assess the current condition and residual strength of underground tunnels.In this work,an underground tunnel model with marble-like cementitious materials was first fabricated using the three-dimensional(3D)printing technique and then loaded to simulate its failure mode in the laboratory.Lead zirconate titanate piezoelectric(PZT)transducers were embedded in the surrounding rock around the tunnel in the process of 3D printing.A 3D monitoring network was formed to locate damage areas and evaluate damage extent during loading.Results show that as the load increased,main cracks firstly appeared above the tunnel roof and below the floor,and then they coalesced into the tunnel boundary.Finally,the tunnel model was broken into several parts.The resonant frequency and the peak of the conductance signature firstly shifted rightwards with loading due to the sealing of microcracks,and then shifted backwards after new cracks appeared.An overall increase in the root-mean-square deviation(RMSD)calculated from conductance signatures of all the PZT transducers was observed as the load(damage)increased.Damage-dependent equivalent stiffness parameters(ESPs)were calculated from the real and imaginary signatures of each PZT at different damage states.Satisfactory agreement between equivalent and experimental ESP values was achieved.Also,the relationship between the change of the ESP and the residual strength was obtained.The method paves the way for damage identification and residual strength estimation of other 3D printed structures in civil engineering.

3D Evaluation Method of Cutting Surface Topography of Carbon / Phenolic (C/Ph) Composite

3D evaluation method of cutting surface topography for C/Ph composites was established.The cutting surface was measured by Talyscan 150,using 3D non-contact measurement.Through the results of 2D and 3D roughness evaluating for C/Ph composite and Duralumin,the 2D evaluation method of the cutting surface topography of C/Ph composite loses a lot of information,the characteristics of the surface topography of C/Ph composite can be comprehensively and authentically evaluated only by 3D evaluation method.Furthermore,3D amplitude and spatial parameters were adopted to evaluate the surface.The results show that： the topography of the C/Ph composite is anisotropic,there are more valleys in the machined surface of C/Ph than that of duralumin,and there are not obvious feeding textures for C/Ph,which indicates the machining mechanism is different from the metal.In conclusion,the topography of the C/Ph composite cutting surface is anisotropic;the cutting surface of C/Ph composite needs 3D evaluation method.

An explicit method for numerical simulation of wave equations: 3D wave motion

In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order （M=1） recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.

The Application of Geostatistical Methods to Prepare the 3D Petrophysical Model of Oil Reservoir

In this study, the petrophysical parameters such as density, sonic, neutron, and porosity were investigated and presented in the 3D models. The 3D models were built using geostatistical method that is used to estimate studied parameters in the entire reservoir. For this purpose, the variogram of each parameter was determined to specify spatial correlation of data. Resulted variograms were non-monotonic. That shows anisotropy of structure. The lithology and porosity parameters are the main causes of this anisotropy. The 3D models also show that petrophysical data has higher variation in north part of reservoir than south part. In addition to, the west limb of reservoir shows higher porosity than east limb. The variation of sonic and neutron data are similar whereas the density data has opposed variation.

ANALYSIS OF SLAB EDGING BY A 3-D RIGID VISCO-PLASTIC FINITE ELEMENT METHOD

3-D rigid visco-plastic finite element method (FEM) is used in the analysisof metal forming processes, including strip and plate rolling, shape rolling, slab edging, specialstrip rolling. The shifted incomplete Cholesky decomposition of the stiffness matrix with thesolution of the equations for velocity increment by the conjugate gradient method is combined. Thistechnique, termed the shifted ICCG method, is then employed to solve the slab edging problem. Theperformance of this algorithm in terms of the number of iterations, friction variation, shiftedparameter psi and the results of simulation for processing parameters are analysed. Numerical testsand application of this technique verify the efficiency and stability of the shifted ICCG method inthe analysis of slab edging.

New Method for 3D Transient Eddy Current Field Calculation and Its Application in Magneto-Acoustic Tomography

A new method of 3D transient eddy current field calculation is proposed. The Maxwell equations with time component elimination （METCE） are derived under the assumption of magnetic quasi static approximation, especially for the sample of low conductivity. Based on METCE, we deduce a more efficient reconstruction algorithm of a 3D transient eddy current field. The computational burden is greatly reduced through the new algorithm, and the computational efficiency is improved. This new algorithm decompounds the space-time variables into two individual variables. The idea is to solve the spatial vector component firstly, and then multiply it by the corresponded time component. The iterative methods based on METCE are introduced to recover the distribution of conductivity in magneto-acoustic tomography. The reconstructed images of conductivity are consistent with the original distribution, which validate the new method.