单程波方程偏移算子高阶辛格式

针对地震偏移算法中单程波算子中e指数项的近似，提出一种高阶辛格式算法. 与非辛格式的近似、一阶辛格式近似和两步法进行数值对比表明，文中提出的高阶辛格式具有更高的精度. 脉冲响应的计算也表明文中的近似方法具有降低频散效应和适应横向变速的能力. 用Marmousi模型验证了对复杂构造的成像的能力，实际资料计算表明本文方法适合于较大深度的稳定准确成像.     

Interaction of water waves with vertical cylinders using null-field integral equations

The scattering of water waves by bottom-mounted vertical circular cylinders is solved by using the null-field integral equations in conjunction with degenerate kernels and Fourier series to avoid calculating the Cauchy and Hadamard principal values. In the implementation, the null-field point can be exactly located on the real boundary owing to the introduction of degenerate kernels for fundamental solutions. An adaptive observer system of polar coordinates is considered to fully employ the properties of degenerate kernels. For the hypersingular equation, vector decomposition for the radial and tangential gradients is carefully considered. This method can be seen as a semi-analytical approach since errors attribute from the truncation of Fourier series. Neither hypersingularity in the Burton and Miller approach nor the CHIEF concept was required to deal with the problem of irregular frequencies. Five advantages of free of calculating principal value, well-posed algebraic system, convergence rate of exponential order, meshfree and elimination of boundary-layer effect, are achieved by using the present approach. Numerical results are given for the forces and free-surface elevation around the circular boundaries. Also, the near-trapped behavior arisen from the physical resonance is detected. A general-purpose program for water wave impinging several circular cylinders with arbitrary number, radii, and positions was developed. Several examples of water wave structure interaction by vertical circular cylinders were demonstrated to see the validity of the present formulation.     

Second-species solutions in the circular and elliptic restricted three-body problem

An analytical proof of the existence of some kinds of periodic orbits of second species of Poincaré, both in the Circular and Elliptic Restricted three-body problem, is given for small values of the mass parameter. The proof uses the asymptotic approximations for the solutions and the matching theory developed by Breakwell and Perko. In the paper their results are extended to the Elliptic problem and applied to prove the existence of second-species solutions generated by rectilinear ellipses in the Circular problem and nearly-rectilinear ones in the Elliptic case.     

Finite cube elements method for calculating gravity anomaly and structural index of solid and fractal bodies with defined boundaries

In this paper, we present the finite cube elements method (FCEM); a novel numerical tool for calculating the gravity anomaly g and structural index SI of solid models with defined boundaries and variable density distributions, tilted or in normal position (e.g. blocks, faulted blocks, cylinders, spheres, hemispheres, triaxial ellipsoids). Extending the calculation to fractal objects, such as Menger sponges of different orders and bodies defined by polyhedrons, demonstrates the robustness of FCEM. In addition, approximating the cube element by a sphere of equal volume makes the calculation of gravitation and related derivatives much simpler. In gravity modelling of a sphere, cubes with edges of 100 m and 200 m achieve a good compromise between running time and overall error.
Displaying the distribution of SI of the studied models on contour maps and profiles will have a strong impact on the forward and inverse modelling of potential field data, especially for Euler deconvolution.
For Menger sponges, plots of gravity elements g and its derivatives show similar patterns independent of fractal order. Moreover, both the pattern and magnitude of SI are independent of fractal order, allowing the use of SI as a new invariant measure for fractal objects. However, SI pattern and magnitude strongly depend on the depth to the buried bodies as do other elements
In this study, we also present a new type of plot; the structural index against distance variation diagrams from which we extract the three critical SI ( CSI ) values, one per axis. The inversion of gravity anomaly data at CSI values gives the optimal mean location of the buried body.     

Instantaneous deformation from continuous GPS: contributions from quasi-periodic loads

Continuous GPS (CGPS) coordinate time-series are known to experience repeating deformation signals with seasonal and other periods. It is unlikely that these signals represent perfect sinusoids with temporally constant amplitude. We develop an analysis method that accommodates temporal variations in the amplitudes of sinusoidal signals. We apply the method to simulated coordinate time-series to numerically explore the potential consequences of neglecting decadal variation in amplitude of annual motions on the residual-error spectra of CGPS measurements, as well as potential bias in estimates for secular site velocity. We find that secular velocity bias can be appreciable for shorter time-series, and that residual-error time-series of longer duration may contain significant power in a broad band centred on semi-annual frequency if temporal variation in the amplitude of annual motions is not accounted for in the model used to reduce the observations to residuals. It may be difficult to differentiate the bandpass filtered signature of mismodelled loading signals from power-law noise, using residual-error spectra for shorter time-series. We provide an example application to a ∼9-yr coordinate time-series for a CGPS station located in southern California at Carbon Creek Control Structure (CCCS), which is known to experience large amplitude seasonal motions associated with the Santa Ana aquifer system.     

一种机载雷达点云数据的快速分类方法

机载激光雷达扫描技术可以以点云形式快速获取地形表面高精度三维信息。基于激光雷达扫描数据及建筑物本身的拓扑信息就可以对建筑物进行精确的重建,而重建中最关键的技术是对点云数据进行分类,进而进行地物识别。对大规模三维点云数据进行快速分类,提出一种采用区域分割结合基于最小二乘平差的多项式拟合方法,将大量离散的三维点云分割后进行多项式拟合,并将二维数据分类转化为一维数据分类。在分类的基础上,将建筑物几何规则作为约束条件提取了房屋边缘。实验分析表明,该方法既能去除多余噪声,又能有效保留特征点,分类的总误差率低于3%。     

库车坳陷复杂高陡构造地震成像研究

复杂构造地震成像主要取决于叠前地震数据品质、偏移速度可靠性和偏移算子成像精度. 库车坳陷异常复杂的近地表条件导致极低信噪比的地震采集数据. 该区逆冲推覆高陡构造刺穿盐体大面积分布, 盐层厚度变化大、顶底面形态复杂, 盐下断裂带破碎、小断块发育, 形成异常复杂的地震成像问题. 本文重点研究三个关键环节:(1)精细的叠前地震预处理研究: 根据该区地震地质复杂性和地震资料特征, 采用一些新的方法技术和技术组合从振幅与时移的大、中、小尺度变化三个层次来解决资料信噪比问题, 重建深部反射信号; (2)三级偏移速度分析研究:利用库车坳陷盐刺穿逆冲推覆构造建模理论及变速成图配套技术解决叠前时间偏移速度场时深转换问题,利用井约束低频速度地震迭代反演技术解决连井层速度场与偏移速度场的融合问题,实现从DMO速度分析、叠前时间偏移速度分析到叠前深度偏移速度分析的有机衔接,建立拓扑结构相对保持的叠前深度偏移速度模型;(3)基于退化Fourier偏移算子的半解析波动方程叠前时间和深度偏移研究, 极大地改善了地震偏移过程中高波数波的成像问题. 通过对库车坳陷大北、博孜、却勒、西秋4和西秋10等复杂高陡构造的叠前时间和深度偏移地震成像处理,取得了较好的应用效果.     

Part II. Discontinuous Galerkin method applied to a single phase flow in porous media

Discontinuous Galerkin numerical simulations of single phase flow problem are described in this paper. The simulations show the advantages of using discontinuous approximation spaces. hp convergence results are obtained for smooth solutions. Unstructured meshes and unsmooth solutions are also considered.