计算最小走时和射线路径的界面网全局方法

用慢度分块均匀正方形模型将介质参数化，仅在正方形单元的边界上设置计算结点，这些结点构成界面网．根据Ｈｕｖｓｅｎｓ和Ｆｅｒｍａｔ原理，由不断扩张、收缩的波前点扫描代替波前面搜索，在波前点附近点的局部最小走时计算中对波前点之间的走时使用双曲线近似，通过比较确定最小走时和相应的次级源位置，记录在以界面网点位置为指针的３个一维数组中．借助这些数组通过向源搜索可计算任意点（包括界面网以外的点）上的全局最小走时和射线路径．这一方法不受介质慢度差异大小限制，占内存少，计算速度较快，适于走时反演和以Ｍａｓｌｏｖ射线理论为基础的波场计算．     

Meteors in the IAU Meteor Data Center on Hyperbolic Orbits

The hyperbolic meteor orbits among the 4,581 photographic and 62,906 radar meteors of the IAU MDC have been analysed using statistical methods. It was shown that the vast majority of hyperbolic orbits has been caused by the dispersion of determined velocities. The large proportion of hyperbolic orbits among the known meteor showers strongly suggests the hyperbolicity of the meteors is not real. The number of apparent hyperbolic orbits increases inversely proportional to the difference between the mean heliocentric velocity of meteor shower and the parabolic velocity limit. The number of hyperbolic meteors in the investigated catalogues does not, in any case, represent the number of interstellar meteors in observational data. The apparent hyperbolicity of these orbits is caused by a high spread in velocity determination, shifting a part of the data through the parabolic limit.     

Load testing of anchors for wire mesh and cable net rockfall slope protection systems

This study presents results of field tests conducted on anchors used to support wire mesh and cable net rockfall protection systems. The load transfer and failure characteristics of these anchors are different from those used in most civil applications in that loads are often applied transversely to the top of tendon rather than axially. The study included vertical as well as horizontal series of tests conducted on some anchors widely used in wire mesh and cable net rockfall protect systems. It was found that the deformation characteristics of these anchors under vertical loading are nonlinear. They are approximated by a hyperbolic formulation and used to calculate the ultimate capacity. Top-downward progressive cracking of the grout was observed during loading and influences the deformation characteristics of these anchors under horizontal loading. The anchors deflected excessively before they could attain their ultimate capacity in the horizontal direction. Based on the field tests, it appears that the deformation under horizontal loading in the systems can be limited by using an enlarged grout zone at the top.     

Strong, localised country-rock contamination and partial homogenisation of a mafic magma: An example from west central Sweden

A dolerite sill cutting slightly older basalt in west-central Sweden shows a strong chemical variation (54% < SiO₂ < 73%) within a restricted area (< 100 × 100 m²). The linear correlation among almost all elements is extremely high; in addition, _NdT is strongly correlated with the SiO₂ content. Least-square hyperbolic-ratio and three-element ratio modelling (common denominator) suggests that most of the chemical variation is explained by mixing and/or micro-mingling. In all, we test 407 hyperbolas, of which 402 are fitted. The five ratio pairs, which could not be fitted to a hyperbola using a least-square fitting procedure, have the ratio Th / Eu in common. Testing the goodness of fit is problematic for hyperbolic distributions; for comparing purposes we sum the distances to chords approximating the hyperbola. Mobile and immobile elements behave similarly, suggesting that no elements are lost or added from outside the system. The data suggests that already the most mafic of the analysed rocks is a mixture of the ‘normal’ dolerite and a siliceous crustal rock. A mafic magma intruded into the base of the crust, where it fractionated resulting in a decreased Mg number. The magma was then contaminated with country rocks in an intermediate magma chamber due to country rock melting; during mixing/mingling almost no fractionation took place. The contaminated rock suggests the presence of a fluid phase. This was probably a prerequisite for country-rock melting. Enrichment in some incompatible elements suggests that besides major mixing/mingling a thermochemical separation process has affected the most felsic rock enriching it in light rare earths and Zr.     

Sensitivity equations for hyperbolic conservation law-based flow models

The present paper focuses on the governing equations for the sensitivity of the variables to the parameters in flow models that can be described by one-dimensional scalar, hyperbolic conservation laws. The sensitivity is shown to obey a hyperbolic, scalar conservation law. The sensitivity is a conserved scalar except in the case of discontinuous flow solutions, where an extra, point source term must be added to the equations in order to enforce conservation. The propagation speed of the sensitivity waves being identical to that of the conserved variable in the original conservation law, the system of conservation laws formed by the original hyperbolic equation and the equation satisfied by the sensitivity is linearly degenerate. A consequence on the solution of the Riemann problem is that rarefaction waves for the variable of the original equation result in vacuum regions for the sensitivity. The numerical solution of the hyperbolic conservation law for the sensitivity by finite volume methods requires the implementation of a specific shock detection procedure. A set of necessary conditions is defined for the discretisation of the source term in the sensitivity equation. An application to the one-dimensional kinematic wave equation shows that the proposed numerical technique allows analytical solutions to be reproduced correctly. The computational examples show that first-order numerical schemes do not yield satisfactory numerical solutions in the neighbourhood of moving shocks and that higher-order schemes, such as the MUSCL scheme, should be used for sharp transients.     

Advected fields in maps: II. Dynamo action in the stretch–fold–shear map

The growth of magnetic field is considered in the stretch–fold–shear map in the limit of weak diffusion. Numerical results are given for insulating, perfectly conducting and periodic boundary conditions. The resulting eigenvalue branches and magnetic fields are related to eigenvalue branches for perfect dynamo action, obtained for zero diffusion using a complex variable formulation.

The effect of diffusion on these perfect dynamo modes depends on their structure, growth rate and the diffusive boundary conditions employed. In some cases, the effect of diffusion is a small perturbation, giving a correction going to zero in the limit of weak diffusion, with a scaling exponent given analytically. In other cases weak diffusion can entirely destroy a perfect dynamo branch. Diffusive boundary layers can also generate entirely new branches.

These different cases are elucidated, and within the framework of the asymptotic approximations used (which do not constitute a rigorous proof), it is seen that for all three boundary conditions employed, the stretch–fold–shear map is a fast dynamo.     

Voronoi图和双曲线联合方法在地震快速定位中的应用

基于Voronoi图和地震震中近似满足双曲线分布,提出了一种联合应用V图和双曲线的地震预警快速定位方法.并对青海地震台网记录的地震进行了重新定位处理.研究结果发现,该地震定位方法能满足地震速报的精度要求,在时效性上也可满足地震预警的需要.该方法在地震预警方面有应用前景.     

An approximate Riemann solver for sensitivity equations with discontinuous solutions

The sensitivity of a model output (called a variable) to a parameter can be defined as the partial derivative of the variable with respect to the parameter. When the governing equations are not differentiable with respect to this parameter, problems arise in the numerical solution of the sensitivity equations, such as locally infinite values or instability. An approximate Riemann solver is thus proposed for direct sensitivity calculation for hyperbolic systems of conservation laws in the presence of discontinuous solutions. The proposed approach uses an extra source term in the form of a Dirac function to restore sensitivity balance across the shocks. It is valid for systems such as the Euler equations for gas dynamics or the shallow water equations for free surface flow. The method is first detailed and its application to the shallow water equations is proposed, with some test cases such as dike- or dam-break problems with or without source terms. An application to a two-dimensional flow problem illustrates the superiority of direct sensitivity calculation over the classical empirical approach.     

地震勘探时空域双曲滤波因子的特点与提高信噪比的能力

为改善强随机噪声背景中地震勘探资料信噪比（S/N）,提出了混沌振子算法;该算法可实现对未知同相轴的确定. 本文的时空域双曲滤波（hyperbolic time\|distance relation filter,HTDF）是上述算法的后续处理,即处理得到存在的弱同相轴.HTDF的理论基础是最小平方滤波和混沌振子算法. 为说明滤波过程与效果,对于包括一个反射同相轴的地震记录进行了多道滤波因子计算. 结果表明,与一个输出道相应的双曲滤波因子组具有复杂的时频域结构,突出的特点是振幅谱“波动式陷频”表现;弥补或消除这个不足的方法是与输出道相应的多个滤波因子之间存在的振幅谱“互补”性.利用相同随机噪声背景下的共炮点理论记录,经与域双曲速度滤波、褶积滤波处理方法比较,在振幅谱、子波初至、零相位性、S/N改善等方面,HTDF的滤波效果优于另外两种.最后探讨了该滤波技术的机理,包括同相轴内各子波间的相似性,以及对组内因子实施的多次有效相加处理.     

潮汕坳陷多次波组合压制技术研究

潮汕坳陷中生界发育,是我国海上重要的油气勘探后备区。但该区多次波非常发育,多次波识别和压制一直是该区地震资料处理的关键和难点。合理地选取有效的多次波压制方法是获得高品质地震剖面的关键所在。针对原始资料中多次波异常发育的情况,文章分析了该区多种类型多次波的特点,尝试了零偏移距截距时间-射线参数(τ-p)域预测反褶积、自由界面多次波衰减法(Surface-Related Multiple Elimination,SRME)和双曲线Radon变换,对多次波进行组合压制,克服了仅依靠一种手段难以压制多种类型多次波的局限。从资料处理效果看,τ-p域预测反褶积方法对短周期的鸣震和微曲多次波的压制效果较好,SRME有效地压制自由表面多次波,且在近炮检距范围内压制多次波效果很好,而双曲线Radon变换是衰减长周期的层间多次波和未被充分压制的自由表面多次波很好的选择。组合压制多次波后的剖面质量明显得到改善。