共查询到20条相似文献，搜索用时 78 毫秒
1.
2.
Improving the accuracy of NMO corrections and of the corresponding interval velocities entails implementing a better approximation than the formula used since the beginning of seismic processing. The exact equations are not practical as they include many unknowns. The approximate expression has only two unknowns, the reflection time and the rms velocity, but becomes inaccurate for large apertures of the recording system and heterogeneous vertical velocities. Several methods of improving the accuracy have been considered, but the gains do not compensate for the dramatic increase in computing time. Two alternative equations are proposed: the first containing two parameters, the reflection time and the focusing time, is not valid for apertures much greater than is the standard formula, but has a much faster computing time and does not stretch the far traces; the other, containing three parameters, the reflection time, like focusing time and the tuning velocity, retains high frequencies for apertures about twice those allowed by the standard equation. Its computing time can be kept within the same limits. NMO equations, old and new, are designed strictly for horizontal layering, but remain reliable as long as the rays travel through the same layers in both the down and up directions. An equation, similar to Dix's formula, is given to compute the interval velocities. The entire scheme can be automated to produce interval-velocity sections without manual picking. 相似文献
3.
Thirteen equations based on the mass-transfer method for determining free water evaporation were expressed in seven generalized equations. These seven equations were then compared with pan evaporation at four climatological stations in north-western Ontario, Canada. The comparisons were based on monthly evaporation. Equations were compared by calibrating them on the entire data sets as well as by calibrating on part of the data and then verifying them on the remainder of the data. The results of comparison showed that all equations were in reasonable agreement with observed evaporation, and that the effect of wind velocity on monthly evaporation was marginal. However, when an equation with parameters obtained at one site was applied to compute evaporation at another site, the computed evaporation was not in good agreement with observed values. © 1997 by John Wiley & Sons, Ltd. 相似文献
4.
Liangguo Dong Zhongyi Fan Hongzhi Wang Benxin Chi Yuzhu Liu 《Geophysical Prospecting》2018,66(8):1503-1520
Reflection full waveform inversion can update subsurface velocity structure of the deeper part, but tends to get stuck in the local minima associated with the waveform misfit function. These local minima cause cycle skipping if the initial background velocity model is far from the true model. Since conventional reflection full waveform inversion using two‐way wave equation in time domain is computationally expensive and consumes a large amount of memory, we implement a correlation‐based reflection waveform inversion using one‐way wave equations to retrieve the background velocity. In this method, one‐way wave equations are used for the seismic wave forward modelling, migration/de‐migration and the gradient computation of objective function in frequency domain. Compared with the method using two‐way wave equation, the proposed method benefits from the lower computational cost of one‐way wave equations without significant accuracy reduction in the cases without steep dips. It also largely reduces the memory requirement by an order of magnitude than implementation using two‐way wave equation both for two‐ and three‐dimensional situations. Through numerical analysis, we also find that one‐way wave equations can better construct the low wavenumber reflection wavepath without producing high‐amplitude short‐wavelength components near the image points in the reflection full waveform inversion gradient. Synthetic test and real data application show that the proposed method efficiently updates the background velocity model. 相似文献
5.
SPH simulation of free surface flow over a sharp-crested weir 总被引：1，自引：0，他引：1
In this paper the numerical simulation of a free surface flow over a sharp-crested weir is presented. Since in this case the usual shallow water assumptions are not satisfied, we propose to solve the problem using the full weakly compressible Navier–Stokes equations with the Tait equation of state for water. The numerical method used consists of the new meshless Smooth Particle Hydrodynamics (SPH) formulation proposed by Ferrari et al. (2009) [8], that accurately tracks the free surface profile and provides monotone pressure fields. Thus, the unsteady evolution of the complex moving material interface (free surface) can been properly solved. The simulations involving about half a million of fluid particles have been run in parallel on two of the most powerful High Performance Computing (HPC) facilities in Europe. The validation of the results has been carried out analysing the pressure field and comparing the free surface profiles obtained with the SPH scheme with experimental measurements available in literature [18]. A very good quantitative agreement has been obtained. 相似文献
6.
为了提高AVO(amplitude versus offset)反演结果的精度和横向连续性,本文提出了一种新的AVO反演约束方法,该方法结合贝叶斯原理和卡尔曼滤波算法实现了对反演参数纵向和横向的同时约束.文章首先结合反演参数的纵向贝叶斯先验概率约束和反演参数的横向连续性假设建立了与卡尔曼滤波算法对应的AVO反演系统的数学模型,然后将该数学模型代入卡尔曼滤波算法框架,利用卡尔曼滤波算法实现了双向约束AVO反演.二维模型测试和实际数据测试结果表明,相对于单纯的纵向贝叶斯先验概率约束,双向约束能更准确地刻画参数的横向变化,得到更准确、横向连续性更好的反演结果. 相似文献
7.
It is increasingly recognized that effective river management requires a catchment scale approach. Sediment transport processes are relevant to a number of river functions but quantifying sediment fluxes at network scales is hampered by the difficulty of measuring the variables required for most sediment transport equations (e.g. shear stress, velocity, and flow depth). We develop new bedload and total load sediment transport equations based on specific stream power. These equations use data that are relatively easy to collect or estimate throughout stream networks using remote sensing and other available data: slope, discharge, channel width, and grain size. The new equations are parsimonious yet have similar accuracy to other, more established, alternatives. We further confirm previous findings that the dimensionless critical specific stream power for incipient particle motion is generally consistent across datasets, and that the uncertainty in this parameter has only a minor impact on calculated sediment transport rates. Finally, we test the new bedload transport equation by applying it in a simple channel incision model. Our model results are in close agreement to flume observations and can predict incision rates more accurately than a more complicated morphodynamic model. These new sediment transport equations are well suited for use at stream network scales, allowing quantification of this important process for river management applications. Copyright © 2017 John Wiley & Sons, Ltd. 相似文献
8.
通过讨论地球空间已有稳定核素内质子数与中子数的分布趋势,介绍了稳定核素的趋势分析方法及其有关周期性分布方程形式,给出了理论方程曲线与地球空间稳定核素实验数据分布点的对比结果,进而给出了稳定核素极限值和元素周期表中化学元素极限,以及其与正负粒子对的可能对应关系方程,包括位于电子中微子层面附近的粒子质量量级初步估计.随后通过建立真空物质能量状态的二个假设,及基于等效Binet方程,给出了与Einstein狭义相对论有关结论相融合的物质粒子以光速及超光速运动的质量及能量方程;作为推论,对这些方程与暗物质及暗能量的可能对应关系予以了初步探讨. 相似文献
9.
Representation and Analysis of the Earthquake Size Distribution: A Historical Review and Some New Approaches 总被引：7，自引：0，他引：7
T. Utsu 《Pure and Applied Geophysics》1999,155(2-4):509-535
—The size distribution of earthquakes has been investigated since the early 20th century. In 1932 Wadati assumed a power-law distribution n(E) = kE ?w for earthquake energy E and estimated the w value to be 1.7 ～ 2.1. Since the introduction of the magnitude-frequency relation by Gutenberg and Richter in 1944 in the form of log n(M) = a?bM, the spatial or temporal variation (or stability) of b value has been a frequently discussed subject in seismicity studies. The log n(M) versus M plots for some data sets exhibit considerable deviation from a straight line. Many modifications of the G-R relation have been proposed to represent such character. The modified equations include the truncated G-R equation, two-range G-R equation, equations with various additional terms to the original G-R equation. The gamma distribution of seismic moments is equivalent to one of these equations.¶In this paper we examine which equation is the most suitable to magnitude data from Japan and the world using AIC. In some cases, the original G-R equation is the most suitable, however in some cases other equations fit far better. The AIC is also a powerful tool to test the significance of the difference in parameter values between two sets of magnitude data under the assumption that the magnitudes are distributed according to a specified equation. Even if there is no significant difference in b value between two data sets (the G-R relation is assumed), we may find a significant difference between the same data sets under the assumption of another relation. To represent a character of the size distribution, there are indexes other than parameters in the magnitude-frequency distribution. The η value is one of such numbers. Although it is certain that these indexes vary among different data sets and are usable to represent a certain feature of seismicity, the usefulness of these indexes in some practical problems such as foreshock discrimination has not yet been established. 相似文献
10.
各向异性介质纯P波方程完全不受横波的干扰,在一定程度上可以减缓由于介质各向异性引起的数值不稳定,本文推导了具有垂直对称轴的横向各向同性(VTI)介质纯P波一阶速度-应力方程.由于纯P波方程存在一个分数形式的伪微分算子,无法直接采用有限差分法求解.针对该问题,本文采用伪谱法和高阶有限差分法联合求解波动方程,重点分析了混合法求解纯P波一阶速度-应力方程的稳定性问题,并给出了混合法求解纯P波方程的稳定性条件.数值模拟结果表明纯P波方程伪谱法和高阶有限差分混合法能够进行复杂介质的正演模拟,在强变速度、变密度的地球介质中仍然具有较好的稳定性. 相似文献
11.
The specific energy and specific force equations have many applications in open-channel flow problems. At present, these equations have analytical solution only for rectangular channels. Trial and error procedure also graphical solutions are the existing methods of solving these equations. No analytical solutions are available in the technical literature for these equations in trapezoidal and triangular channels because it is presumed that these equations are quintic equations. The inversion of such equations consists of finding the roots of quintic equations. In the current study for a given channel geometry and discharge, the subcritical (supercritical) depth is analytically found in terms of the other supercritical (subcritical) depth. For this purpose, by considering physically meaningful domains, a quintic equation has been reduced to a quartic equation. In the next step, this quartic equation has been converted to a resolvent cubic equation and two quadratic equations. This research shows these steps clearly to reach an acceptable physical analytic solution for water depth in trapezoidal and triangular channels. 相似文献
12.
采用τ -p变换研究各向同性纵波反射系数的振幅特性,通过分析纵波在τ -p域中的时差特性曲线信息,提出了一个和传统Shuey算法相似的新的二元数学模型,只需对目标反射层以及上中下三层进行τ -p域里的时差分析,而无须其他附加信息.结合模拟数据以及实际的测井数据,分别给出了详细的实验步骤及分析.实验结果证明此算法的反演结果和实际值吻合良好,可以应用于实际的工程分析.分析方法可以作为Shuey AVO的替代,同时具有相似的准确度,结果令人满意. 相似文献
13.
Ito equations are derived for simple stochastic cellular automaton with parameters describing efficiencies for avalanche triggering
and cell occupation. Analytical results are compared with the numerical one obtained from the histogram method. Good agreement
for various parameters supports the wide applicability of the Ito equation as a macroscopic model of some cellular automata
and complex natural phenomena which manifest random energy release. Also, the paper is an example of effectiveness of histogram
procedure as an adequate method of nonlinear modeling of time series. 相似文献
14.
Pure acoustic wave propagation in transversely isotropic media by the pseudospectral method 总被引：4，自引：0，他引：4
Seismic wave propagation in transversely isotropic (TI) media is commonly described by a set of coupled partial differential equations, derived from the acoustic approximation. These equations produce pure P‐wave responses in elliptically anisotropic media but generate undesired shear‐wave components for more general TI anisotropy. Furthermore, these equations suffer from instabilities when the anisotropy parameter ε is less than δ. One solution to both problems is to use pure acoustic anisotropic wave equations, which can produce pure P‐waves without any shear‐wave contaminations in both elliptical and anelliptical TI media. In this paper, we propose a new pure acoustic transversely isotropic wave equation, which can be conveniently solved using the pseudospectral method. Like most other pure acoustic anisotropic wave equations, our equation involves complicated pseudo‐differential operators in space which are difficult to handle using the finite difference method. The advantage of our equation is that all of its model parameters are separable from the spatial differential and pseudo‐differential operators; therefore, the pseudospectral method can be directly applied. We use phase velocity analysis to show that our equation, expressed in a summation form, can be properly truncated to achieve the desired accuracy according to anisotropy strength. This flexibility allows us to save computational time by choosing the right number of summation terms for a given model. We use numerical examples to demonstrate that this new pure acoustic wave equation can produce highly accurate results, completely free from shear‐wave artefacts. This equation can be straightforwardly generalized to tilted TI media. 相似文献
15.
Nonuniqueness in geophysical inverse problems is naturally resolved by incorporating prior information about unknown models into observed data. In practical estimation procedures, the prior information must be quantitatively expressed. We represent the prior information in the same form as observational equations, nonlinear equations with random errors in general, and treat as data. Then we may define a posterior probability density function of model parameters for given observed data and prior data, and use the maximum likelihood criterion to solve the problem. Supposing Gaussian errors both in observed data and prior data, we obtain a simple algorithm for iterative search to find the maximum likelihood estimates. We also obtain an asymptotic expression of covariance for estimation errors, which gives a good approximation to exact covariance when the estimated model is linearly close to a true model. We demonstrate that our approach is a general extension of various inverse methods dealing with Gaussian data. By way of example, we apply the new approach to a problem of inferring the final rupture state of the 1943 Tottori earthquake (M = 7.4) from coseismic geodetic data. The example shows that the use of sufficient prior information effectively suppresses both the nonuniqueness and the nonlinearity of the problem. 相似文献
16.
Theoretical equations of state for Fe at high pressures and temperatures are derived from the expression of the free energy written as a sum of the static energy, the harmonic, the anharmonic and the electronic contributions. All the calculations have been performed for the various crystal structures of Fe using different intermolecular potentials, and namely: Lennard-Jones, Morse and Rydberg functions. The available experimental data do not allow a definite choice between them thus leaving open the problem of the composition of the earth's core. 相似文献
17.
The accuracy of the new equations for long frame-type structures, derived by Kerr and Zarembski, and extended recently by Kerr and Accorsi also to dynamic analyses, is investigated. At first, the natural frequencies of a tall frame-type structure, with twenty floors, are determined using the new equations, and then they are compared with the corresponding results calculated using the finite element method and the common shear-building analysis. The found agreement with the finite element analysis is close. To check further the accuracy of the new equations, small scale models are tested on a shake table. The natural frequencies recorded in the tests are then compared with those determined from the new equations, for a wide range of geometrical parameters. The agreement with the test data is good. The presented study indicates that the new equations are well suited for analyzing the dynamic response of tall frame-type structures. 相似文献
18.
Rainfall replenishes surface and subsurface water but is partially intercepted by a canopy. However, it is challenging to quantify the rainfall passing through the canopy (i.e. throughfall). This study derives simple‐to‐use empirical equations relating throughfall to canopy and rainfall characteristics. Monthly throughfall is calculated by applying a mass balance model on weather data from Singapore; Vancouver, Canada; and Stanford, USA. Regression analysis is then performed on the calculated throughfall with three dependent variables (i.e. maximum canopy storage, average rainfall depth and time interval between two consecutive rainfall in a month) to derive the empirical equations. One local equation is derived for each location using data from that particular location, and one global equation is derived using data from all three locations. The equations are further verified with calculated monthly throughfall from other weather data and actual throughfall field measurements, giving an accuracy of about 80–90%. The global equation is relatively less accurate but is applicable worldwide. Overall, this study provides a global equation through which one can quickly estimate throughfall with only information on the three variables. When additional weather data are available, one can follow the proposed methodology to derive their own equations for better estimates. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
19.
Jianghai Xia Yixian Xu Chao Chen Ronald D. Kaufmann Yinhe Luo 《Soil Dynamics and Earthquake Engineering》2006,26(5):395-403
We discuss five useful equations related to high-frequency surface-wave techniques and their implications in practice. These equations are theoretical results from published literature regarding source selection, data-acquisition parameters, resolution of a dispersion curve image in the frequency–velocity domain, and the cut-off frequency of high modes. The first equation suggests Rayleigh waves appear in the shortest offset when a source is located on the ground surface, which supports our observations that surface impact sources are the best source for surface-wave techniques. The second and third equations, based on the layered earth model, reveal a relationship between the optimal nearest offset in Rayleigh-wave data acquisition and seismic setting—the observed maximum and minimum phase velocities, and the maximum wavelength. Comparison among data acquired with different offsets at one test site confirms the better data were acquired with the suggested optimal nearest offset. The fourth equation illustrates that resolution of a dispersion curve image at a given frequency is directly proportional to the product of a length of a geophone array and the frequency. We used real-world data to verify the fourth equation. The last equation shows that the cut-off frequency of high modes of Love waves for a two-layer model is determined by shear-wave velocities and the thickness of the top layer. We applied this equation to Rayleigh waves and multi-layer models with the average velocity and obtained encouraging results. This equation not only endows with a criterion to distinguish high modes from numerical artifacts but also provides a straightforward means to resolve the depth to the half space of a layered earth model. 相似文献
20.
Constraining the focal mechanism of the Lushan earthquake with observations of the Earth’s free oscillations 总被引：1，自引：0，他引：1
The amplitudes of the Earth's free oscillations have a close relationship to earthquake focal mechanisms. Focal mechanisms of large earthquakes can be well analyzed and constrained with observations of long period free oscillations. Although the 2013 Lushan earthquake was only moderately sized, observable spherical normal modes were excited and clearly observed by a superconductive gravimeter and a broadband seismometer. We compare observed free oscillations with synthetic normal modes corresponding to four different focal mechanisms for the Lushan earthquake. The results show that source parameters can be analyzed and constrained by spherical normal modes in a 2.3–5 mHz frequency band. The scalar seismic moment M0 has a major influence on the amplitudes of free oscillations; additionally, the strike, dip, rake and depth of the hypocenter have minor influences. We found that the synthetic modes corresponding to the focal mechanism determined by the Global Centroid Moment Tensor show agreement to the observed modes, suggesting that earthquake magnitudes predicted in this way can readily reflect the total energy released by the earthquake. The scalar seismic moment obtained by far-field body wave inversion is significantly underestimated. Focal mechanism solutions can be improved by joint inversion of far- and near-field data. 相似文献