Joint 3D traveltime calculation based on fast marching method and wavefront construction

3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.

基于FAST和Sobol指数法的雷达系统效能敏感性分析

为了以较小的代价快速提升对空探测雷达系统的综合效能,首先搭建了对空探测雷达系统综合效能的评估指标体系,在此基础上完成了综合效能评估模型的建立。选择两种基于方差原理的全局敏感性分析方法,即傅里叶幅度敏感性测试(Fourier amplitudes sensitivity test,FAST)和Sobol指数法,分别从频域和时域对建立的评估模型进行敏感性分析。在基于FAST进行敏感性分析时,提出了21维输入指标所对应的特征频率集的选取方法;在基于Sobol指数法进行敏感性分析时,考虑了二阶及高阶敏感性指数的影响。对比FAST和Sobol指数法的分析结果,两种方法获得的敏感性分析结果具有较好的一致性。其中,功率和频率为对空探测雷达系统的主要影响指标,为雷达系统的研发和升级提供了可靠性依据。

基于Fast Marching Method的肝脏CT图像序列自动分割

快速行进法(fast marching method,FMM)已被证明在图像分割方面具有优势,在此基础上提出了一个混合分割的算法。这个方法加入了图像分割后处理步骤,成功解决了活体肝脏CT系列图像自动分割问题。首先是通过滤波去噪等处理得到速率系数图像,然后根据CT图像相邻层间的相似特点计算FMM所需的参数进行图像分割,最后使用开运算修正肝脏边缘。整个序列分割过程只需用户定义一个种子点,减少了人工干预,从而提高了效率和准确性。

Fast precise integration method for hyperbolic heat conduction problems

A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.

Fast Marching Method for Microseismic Source Location in Cavern-Containing Rockmass: Performance Analysis and Engineering Application

Microseismic(MS)event locations are vital aspect of MS monitoring technology used to delineate the damage zone inside the surrounding rock mass.However,complex geological conditions can impose significantly adverse effects on the final location results.To achieve a high-accuracy location in a complex cavern-containing structure,this study develops an MS location method using the fast marching method(FMM)with a second-order difference approach(FMM2).Based on the established velocity model with three-dimensional(3D)discrete grids,the realization of the MS location can be achieved by searching the minimum residual between the theoretical and actual first arrival times.Moreover,based on the calculation results of FMM2,the propagation paths from the MS sources to MS sensors can be obtained using the linear interpolation approach and the Runge–Kutta method.These methods were validated through a series of numerical experiments.In addition,our proposed method was applied to locate the recorded blasting and MS events that occurred during the excavation period of the underground caverns at the Houziyan hydropower station.The location results of the blasting activities show that our method can effectively reduce the location error compared with the results based on the uniform velocity model.Furthermore,the obtained MS location was verified through the occurrence of shotcrete fractures and spalling,and the monitoring results of the in-situ multipoint extensometer.Our proposed method can offer a more accurate rock fracture location and facilitate the delineation of damage zones inside the surrounding rock mass.

基于Kano-FAST的瓦楞纸激光印刷设备设计研究

目的推动印刷行业朝更高效的方向迈进,提高印刷设备的易操作性和视觉识别性。方法将Kano和FAST(Function Analysis System Technique)模型引入瓦楞纸激光打印设备设计的前期应用需求分析中,通过问卷的方法获取用户的基本要求,并划分为几个子类型,进而建立Kano的二维功能属性模型。采用FAST法建立功能树,辅助使用Kano模型,从而更精准地分析用户需求,并更好地根据其需求进行优化设计。结果综合运用设计原理,针对性地挖掘瓦楞纸激光印刷设备在造型识别性、操作易用性、生产安全性上存在的问题,进而输出更优解。结论该设计方法的引入有助于为同类型的印刷设备设计提供参考,并引起更多相关厂家的重视,推动印刷行业向更积极的方向发展。

Calculation of all-time apparent resistivity of large loop transient electromagnetic method with very fast simulated annealing

In large loop transient electromagnetic method(TEM),the late time apparent resistivity formula cannot truly reflect the geoelectric model,thus it needs to define the all-time apparent resistivity with the position information of measuring point.Utilizing very fast simulated annealing(VFSA) to fit the theoretical electromagnetic force(EMF) and measured EMF could obtain the all-time apparent resistivity of the measuring points in rectangular transmitting loop.The selective cope of initial model of VFSA could be confirmed by taking the late time apparent resistivity of transient electromagnetic method as the prior information.For verifying the correctness,the all-time apparent resistivities of the geoelectric models were calculated by VFSA and dichotomy,respectively.The results indicate that the relative differences of apparent resistivities calculated by these two methods are within 3%.The change of measuring point position has little influence on the tracing pattern of all-time apparent resistivity.The first branch of the curve of all-time apparent resistivity is close to the resistivity of the first layer medium and the last branch is close to the resistivity of the last layer medium,which proves the correctness of the arithmetics proposed.

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.

Trajectory tracking and obstacle avoidance method for robots based on fast terminal sliding mode

For accurate trajectory tracking and obstacle avoidance in finite time of a nonholonomic mobile robot,a trajectory tracking controller based on global fast terminal sliding mode method is proposed,which has the advantages of chattering-free and adjustable convergence time.First of all,the kinematics model of the robot is established in mobile carrier coordinates.Secondly,the global structure including terminal attractor and exponential convergence of the fast terminal sliding mode trajectory tracking controller is proved by Lyapunov stability theory,ensuring that the trajectory and heading angle tracking error converges to a smaller zero range in finite time.Finally,the artificial potential field obstacle avoidance method is introduced to make the robot not only track the reference trajectory strictly,but also avoid the obstacles.The simulation results show that the proposed method can achieve a stable tracking control in finite time for a given reference trajectory.

A fast explicit finite difference method for determination of wellhead injection pressure

A fast explicit finite difference method (FEFDM),derived from the differential equations of one-dimensional steady pipe flow,was presented for calculation of wellhead injection pressure.Recalculation with a traditional numerical method of the same equations corroborates well the reliability and rate of FEFDM.Moreover,a flow rate estimate method was developed for the project whose injection rate has not been clearly determined.A wellhead pressure regime determined by this method was successfully applied to the trial injection operations in Shihezi formation of Shenhua CCS Project,which is a good practice verification of FEFDM.At last,this method was used to evaluate the effect of friction and acceleration terms on the flow equation on the wellhead pressure.The result shows that for deep wellbore,the friction term can be omitted when flow rate is low and in a wide range of velocity the acceleration term can always be deleted.It is also shown that with flow rate increasing,the friction term can no longer be neglected.

Error Control Strategies for Numerical Integrations in Fast Collocation Methods

We propose two error control techniques for numerical integrations in fast multiscale collocation methods for solving Fredholm integral equations of the second kind with weakly singular kernels. Both techniques utilize quadratures for singular integrals using graded points. One has a polynomial order of accuracy if the integrand has a polynomial order of smoothness except at the singular point and the other has exponential order of accuracy if the integrand has an infinite order of smoothness except at the singular point. We estimate the order of convergence and computational complexity of the corresponding approximate solutions of the equation. We prove that the second technique preserves the order of convergence and computational complexity of the original collocation method. Numerical experiments are presented to illustrate the theoretical estimates.

A fast computational method for the landing footprints of space-to-ground vehicles

Fast computation of the landing footprint of a space-to-ground vehicle is a basic requirement for the deployment of parking orbits, as well as for enabling decision makers to develop real-time programs of transfer trajectories. In order to address the usually slow computational time for the determination of the landing footprint of a space-to-ground vehicle under finite thrust, this work proposes a method that uses polynomial equations to describe the boundaries of the landing footprint and uses back propagation(BP) neural networks to quickly determine the landing footprint of the space-to-ground vehicle. First, given orbital parameters and a manoeuvre moment, the solution model of the landing footprint of a space-to-ground vehicle under finite thrust is established. Second, given arbitrary orbital parameters and an arbitrary manoeuvre moment, a fast computational model for the landing footprint of a space-to-ground vehicle based on BP neural networks is provided.Finally, the simulation results demonstrate that under the premise of ensuring accuracy, the proposed method can quickly determine the landing footprint of a space-to-ground vehicle with arbitrary orbital parameters and arbitrary manoeuvre moments. The proposed fast computational method for determining a landing footprint lays a foundation for the parking-orbit configuration and supports the design of real-time transfer trajectories.

A wideband fast multipole boundary element method for half-space/plane-symmetric acoustic wave problems

This paper presents a novel wideband fast multipole boundary element approach to 3D half-space/planesymmetric acoustic wave problems.The half-space fundamental solution is employed in the boundary integral equations so that the tree structure required in the fast multipole algorithm is constructed for the boundary elements in the real domain only.Moreover,a set of symmetric relations between the multipole expansion coefficients of the real and image domains are derived,and the half-space fundamental solution is modified for the purpose of applying such relations to avoid calculating,translating and saving the multipole/local expansion coefficients of the image domain.The wideband adaptive multilevel fast multipole algorithm associated with the iterative solver GMRES is employed so that the present method is accurate and efficient for both lowand high-frequency acoustic wave problems.As for exterior acoustic problems,the Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem involved in the conventional boundary integral equation method.Details on the implementation of the present method are described,and numerical examples are given to demonstrate its accuracy and efficiency.

Application of the extended Fourier amplitude sensitivity testing(FAST)method to inflated,axial stretched,and residually stressed cylinders

This paper is dedicated to applying the Fourier amplitude sensitivity test(FAST)method to the problem of mixed extension and inflation of a circular cylindrical tube in the presence of residual stresses.The metafunctions and the Ishigami function are considered in the sensitivity analysis(SA).The effects of the input variables on the output variables are investigated,and the most important parameters of the system under the applied pressure and axial force such as the axial stretch and the azimuthal stretch are determined.

Fast Measurement of Dielectric Conductivity for Space Application by Surface Potential Decay Method

Surface potential decay of polymers for electrical insulation can help to determine the dark conductivity for spacecraft charging analysis. Due to the existence of radiation-induced conductivity, it decays fast in the first few hours after irradiation and exponentially slowly for the remaining time. The measurement of dark conductivity with this method usually takes the slow part and needs a couple of days. Integrating the Fowler formula into the deep dielectric charging equations, we obtain a new expression for the fast decay part. The experimental data of different materials, dose rates and temperatures are fitted by the new expression. Both the dark conductivity and the radiation-induced conductivity are derived and compared with other methods. The result shows a good estimation of dark conductivity and radiation-induced conductivity in high-resistivity polymers, which enables a fast measurement of dielectric conductivity within about 600 rain after irradiation.

Power System State Estimation Solution With Zero Injection Constraints Using Modified Newton Method and Fast Decoupled Method in Polar Coordinate

如何保证零注入节点的注入功率在状态估计结果中严格为0是电力系统状态估计研究中的重要问题。在直角坐标下,由于零注入约束为线性约束,可使用修正牛顿法来有效地解决这一问题。因此,借鉴直角坐标下修正牛顿法的思路,提出了极坐标下的修正牛顿法和修正快速解耦估计。这些方法的计算流程与传统的极坐标下的牛顿法和快速解耦估计非常相似,计算速度与大权重法相当,同时能够保证零注入约束严格满足。仿真结果验证了所得结论。

AN AUTOMATIC SEGMENTATION METHOD FOR FAST IMAGING IN PET

A new segmentation method has been developed for PET fast imaging. The technique automatically segments the transmission images into different anatomical regions, it efficiently reduced the PET transmission scan time. The result shows that this method gives only 3 min-scan time which is perfect for attenuation correction of the PET images instead of the original 15-30 min-scan time. This approach has been successfully tested both on phantom and clinical data.

Scattering Analysis of Bi-anisotropic Materials Using Fast Dipole Method

Electromagnetic scattering from inhomogeneous three-dimensional( 3D) bi-anisotropic scatterers is formulated in terms of the volume integral equation( VIE) method. Based on the volume equivalence principle,the VIE is represented in terms of a pair of coupled bi-anisotropic polarized volume electric and magnetic flux densities. The VIE is solved using the method of moments( MoM) combined with tetrahedral mesh. Then the fast dipole method( FDM) based on the equivalent dipole method( EDM) is extended to analyze the scattering of bi-anisotropic media by solving the VIE. Finally,some numerical results are given to demonstrate the accuracy of the developed method for the scattering analysis of the bi-anisotropic media.

Modified slow-fast analysis method for slow-fast dynamical systems with two scales in frequency domain

A modified slow-fast analysis method is presented for the periodically excited non-autonomous dynamical system with an order gap between the exciting frequency and the natural frequency.By regarding the exciting term as a slow-varying parameter,a generalized autonomous fast subsystem can be defined,the equilibrium branches as well as the bifurcations of which can be employed to account for the mechanism of the bursting oscillations by combining the transformed phase portrait introduced.As an example,a typical periodically excited Hartley model is used to demonstrate the validness of the method,in which the exciting frequency is far less than the natural frequency.The equilibrium branches and their bifurcations of the fast subsystem with the variation of the slow-varying parameter are presented.Bursting oscillations for two typical cases are considered,which reveals that,fold bifurcation may cause the the trajectory to jump between different equilibrium branches,while Hopf bifurcation may cause the trajectory to oscillate around the stable limit cycle.

A fast computing method to distinguish the hyperbolic trajectory of an non-autonomous system

Attempting to find a fast computing method to DHT （distinguished hyperbolic trajectory）, this study first proves that the errors of the stable DHT can be ignored in normal direction when they are computed as the trajectories extend. This conclusion means that the stable flow with perturbation will approach to the real trajectory as it extends over time. Based on this theory and combined with the improved DHT computing method, this paper reports a new fast computing method to DHT, which magnifies the DHT computing speed without decreasing its accuracy.