Seismic data compression based on integer wavelet transform

Due to the particularity of the seismic data, they must be treated by lossless compression algorithm in some cases. In the paper, based on the integer wavelet transform, the lossless compression algorithm is studied. Comparing with the traditional algorithm, it can better improve the compression rate. CDF (2, n) biorthogonal wavelet family can lead to better compression ratio than other CDF family, SWE and CRF, which is owe to its capability in canceling data redundancies and focusing data characteristics. CDF (2, n) family is suitable as the wavelet function of the lossless compression seismic data. Contribution No.04FE1019, Institute of Geophysics, China Earthquake Administration.     

Application of saddlepoint approximation in reliability analysis of dynamic systems

The application of the saddlepoint approximation to reliability analysis of dynamic systems is investigated. The failure event in reliability problems is formulated as the exceedance of a single performance variable over a prescribed threshold level. The saddlepoint approximation technique provides a choice to estimate the cumulative distribution function (CDF) of the performance variable. The failure probability is obtained as the value of the complement CDF at a specif ied threshold. The method requires computing the saddlepoint from a simple algebraic equation that depends on the cumulant generating function (CGF) of the performance variable. A method for calculating the saddlepoint using random samples of the performance variable is presented. The applicable region of the saddlepoint approximation is discussed in detail. A 10-story shear building model with white noise excitation illustrates the accuracy and effi ciency of the proposed methodology.     

Seismic data compression based on integer wavelet transform

Introduction With the development of the seismological observation technique and deep-going of seismicdata application fields, especially the digitization of data in earthquake station networks, theimprovement of the precision, the data quantity increases as geometric order, which bringdifficulty to saving and transfering these data. To keep all information, seismic data, like medicalimages, should be compressed without error in many applications. In generally, traditionalcompression meth…     

Anti-erosive effectiveness of Eucalyptus coppices through the cover management factor estimate

In this paper sediment yield data, measured from 1978 to 1994 in two small Calabrian basins (W2 and W3) reafforested with eucalyptus trees (Eucalyptus occidentalis Engl.), and the Modified Universal Soil Loss Equation (MUSLE) applied in a distributed form are used to evaluate the anti-erosive effects of eucalyptus cover. At first step the sediment yield measurements observed in W2 basin are used to estimate a single cover management factor representative of the eucalyptus coppice and equal to the median value (0·164) of the cover management factor values calculated for each runoff event. Then, the reliability of the selected representative cover management factor is verified on W3 by comparing the cumulative distribution function (CDF) of the measured sediment yield with the CDF of the calculated one. Finally, the temporal analysis of the crop and management factor is developed searching, at monthly and annual scale, the correlation between crop anti-erosive effectiveness and rainfall erosivity index. © 1998 John Wiley & Sons, Ltd.     

Response uncertainty and time-variant reliability analysis for hysteretic MDF structures

Response uncertainty evaluation and dynamic reliability analysis corresponding to classical stochastic dynamic analysis are usually restricted to the uncertainties of the excitation. The inclusion of the parameter uncertainties contained in structural properties and excitation characteristics has become an increasingly important problem in many areas of dynamics. In the present paper, a point estimate procedure is proposed for the evaluation of stochastic response uncertainty, and a response surface approach procedure in standard normal space is proposed for analysis of time-variant reliability analysis for hysteretic MDF structures having parameter uncertainties. Using the proposed procedures, the response uncertainties and time-variant reliability can be easily obtained by several repetitions of stochastic response analysis under given parameters without conducting sensitivity analysis, which is considered to be one of the primary difficulties associated with conventional methods. In the time-variant reliability analysis, the failure probability can be readily obtained by improving the accuracy of the first-order reliability method using the empirical second-order reliability index. The random variables are divided into two groups, those with CDF and those without CDF. The latter are included via the high-order moment standardization technique. A numerical example of a 15F hysteretic MDF structure that takes into account uncertainties of four structural parameters and three excitation characteristics is performed, based on which the proposed procedures are investigated and the effects of parameter uncertainties are discussed. Copyright © 1999 John Wiley & Sons, Ltd.     

Higher-order approximation techniques for estimating stochastic parameter of a sediment transport model

Higher-order approximation techniques for estimating stochastic parameter of the non-homogeneous Poisson (NHP) model are presented. The NHP model is characterized by a two-parameter cumulative probability distribution function (CDF) of sediment displacement. Those two parameters are the temporal and spatial intensity functions, physically representing the inverse of the average rest period and step length of sediment particles, respectively. Difficulty of estimating the parameters has, however, restricted the applications of the NHP model. The approximation techniques are proposed to address such problem. The basic idea of the method is to approximate a model involving stochastic parameters by Taylor series expansion. The expansion preserves certain higher-order terms of interest. Using the experimental (laboratory or field) data, one can determine the model parameters through a system of equations that are simplified by the approximation technique. The parameters so determined are used to predict the cumulative distribution of sediment displacement. The second-order approximation leads to a significant reduction of the CDF error (of the order of 47%) compared to the first-order approximation. Error analysis is performed to evaluate the accuracy of the first- and second-order approximations with respect to the experimental data. The higher-order approximations provide better estimations of the sediment transport and deposition that are critical factors for such environment as spawning gravel-bed.     

Higher-order approximation techniques for estimating stochastic parameter of a sediment transport model

Higher-order approximation techniques for estimating stochastic parameter of the non-homogeneous Poisson (NHP) model are presented. The NHP model is characterized by a two-parameter cumulative probability distribution function (CDF) of sediment displacement. Those two parameters are the temporal and spatial intensity functions, physically representing the inverse of the average rest period and step length of sediment particles, respectively. Difficulty of estimating the parameters has, however, restricted the applications of the NHP model. The approximation techniques are proposed to address such problem. The basic idea of the method is to approximate a model involving stochastic parameters by Taylor series expansion. The expansion preserves certain higher-order terms of interest. Using the experimental (laboratory or field) data, one can determine the model parameters through a system of equations that are simplified by the approximation technique. The parameters so determined are used to predict the cumulative distribution of sediment displacement. The second-order approximation leads to a significant reduction of the CDF error (of the order of 47%) compared to the first-order approximation. Error analysis is performed to evaluate the accuracy of the first- and second-order approximations with respect to the experimental data. The higher-order approximations provide better estimations of the sediment transport and deposition that are critical factors for such environment as spawning gravel-bed.     

Model simulation of the large-scale high-latitude F-layer modification by powerful HF waves with different modulation

The mathematical model of the high-latitude ionosphere, developed earlier, is applied to investigate how the modulation regime of the ionospheric HF heating facility near Tromso, Scandinavia, affects the large-scale high-latitude F-layer modification. Simulations are made for distinct cases, in which high-power waves have different modulations, namely, for continuous wave transmission and for pulse operation, with the amplitude of the HF wave being square wave-modulated. The calculations are performed for different lengths of pulses and various time intervals between successive pulses. The frequency of HF waves is chosen to be close to the most effective frequency for the large-scale F2-layer modification. Simulations are made for autumn and low geomagnetic activity conditions both for nocturnal and for daytime conditions. The results of modeling indicate that the most considerable decrease in the F-region electron concentration may be achieved when the heater is operated continuously. Moreover, the perceptible decrease in the F-region electron concentration may take place when the heating facility is operated pulsatily. For the pulse operation, the amplitude of the electron concentration variations depends on the ratio of the length of pulses to the time interval between successive pulses. The higher the latter ratio is, the more the electron concentration variation amplitude ought to be.     

基于谱修正方法的非高斯风场模拟

为克服基于H erm ite谱修正方法的缺点,减少该方法中计算H erm ite多项式系数所需耗费的大量机时,提出了一种模拟非高斯风压场的新方法,采用非高斯累积分布函数（CDF）映射技术来代替基于H erm ite的概率密度函数（PDF）修正。选择任意边缘PDF模型作为概率目标模型,采用目标功率谱密度（PSD）作为样本函数,通过迭代修正该样本函数,使其收敛于目标概率密度函数和目标功率谱密度。将该方法应用于实际结构的非高斯风场模拟,模拟结果与目标谱符合良好,表明本文方法模拟非高斯风场具有较高的精确度和计算效率。     

基于提升算法和百分位数软阈值的小波去噪技术

在地震勘探领域,随机噪声一直是影响地震信号信噪比的主要因素之一,如何从被干扰的地震信号中有效去除随机噪声并保护有用信号具有重要的意义.针对经典小波变换在计算效率方面的缺陷,本文推荐应用提升算法实现第二代小波变换的构建,分析和对比了提升算法(Lifting Scheme)下不同小波变换方法的特性,选取更加符合小波域去噪原理的CDF 9/7双正交小波变换作为基本算法,同时应用了简单、有效的百分位数(Percentiles)软阈值进行信噪分离.通过理论模型处理,本方法可以在去噪能力和保护有用信号之间找到很好的平衡点.实际剖面的处理效果表明,此方法不仅能有效的滤除随机噪声,而且很好地保护有用信号,提高地震数据分析的精确性.     

Parameter Uncertainty Propagation in a Rainfall–Runoff Model; Case Study: Karoon-III River Basin

Conceptual hydrological models are popular tools for simulating land phase of hydrological cycle. Uncertainty arises from a variety of sources such as input error, calibration and parameters. Hydrologic modeling researches indicate that parametric uncertainty has been considered as one of the most important source. The objective of this study was to evaluate parameter uncertainty and its propagation in rainfall-runoff modeling. This study tried to model daily flows and calculate uncertainty bounds for Karoon-III basin, Southwest of Iran, using HEC-HMS (SMA). The parameters were represented by probability distribution functions (PDF), and the effect on simulated runoff was investigated using Latin Hypercube Sampling (LHS) on Monte Carlo (MC). Three chosen parameters, based on sensitivity analysis, were saturated-hydraulic-conductivity (Ks), Clark storage coefficient (R) and time of concentration (t _c ). Uncertainty associated with parameters were accounted for, by representing each with a probability distribution. Uncertainty bounds was calculated, using parameter sets captured from LHS on parameters PDF of sub-basins and propagating to the model. Results showed that maximum reliability (11%) resulted from Ks propagating. For three parameters, underestimation was more than overestimation. Maximum sharpness and standard deviation (STD) was resulted from propagating Ks. Cumulative Distribution Function (CDF) of flow and uncertainty bounds showed that as flow increased, the width of uncertainty bounds increased for all parameters.     

中国大陆及邻区强震活动主体地区形成的数值模拟

在一定时期内,中国大陆及邻区强震活动发生在相对集中的地区,本文简称为主体地区;不同时期的主体地区在空间上有较大距离的转移与变换．本文在前人研究基础上,对主体地区的形成机制进行探讨,认为岩石圈韧性层的变形控制了强震活动主体地区的主要位置和形态,而周围板块边界作用是形成主体地区的主要驱动力．\r\n 以中国大陆及邻区韧性层为模拟对象,建立二维粘弹性有限元模型模拟韧性层的变形及应力场分布．考虑到板块边界活动在时空上的不均匀性,将板块边界分段,以不同板块边界段组合作为各时期强震活动主体地区形成的主要约束,模拟大陆内部韧性层应力集中区的形成．模拟结果显示了6种与实际强震活动主体地区有可比性的应力分布,说明在相应的板块边界段组合作用下,在韧性层可以形成与各主体地区的形态和方位基本一致的应力集中区．     

Statistical distributions of phytoplankton concentration in river water

Daily time series of phytoplankton concentration are theoretically analyzed with the use of population dynamics equation with fluctuations in the growth factor and environmental carrying capacity taken into account. The statistical distributions of phytoplankton cell concentration are shown to obey different laws, depending on the season of the year: lognormal in the winter and logarithmic in the vegetation period. The probability of extremely high concentration numbers is described by the normal law. The verification of the obtained relationships against a body of empirical data confirmed the theoretical forecasts. The obtained results make it possible to predict the probabilities of various phytoplankton concentration values within a wide range, including the domain of large values, which are of greatest hazard in terms of water quality, water treatment processes, and aquatic ecosystem well-being.     

Non-extensive statistical physics approach to fault population distribution. A case study from the southern Hellenic arc (Central Crete)

Fault population statistics play a key role in the understanding of any statistical seismicity approach. In the present work a non-extensive statistical physics approach is formulated and tested for the local fault length distribution. The approach is composed of the following parts: (i) Tsallis entropy, S _q , (ii) maximization of the Tsallis entropy under appropriate constrains, and (iii) derivation of the cumulative distribution function (CDF) of the fault length population. This model is tested using fault length data from the Central Crete graben in front of the Hellenic arc and estimated a thermodynamic q parameter equal to 1.16, which supports the conclusion that the fault system in Central Crete graben is a sub-extensive one.     

On the relations between the hydrological dynamical systems of water budget,travel time,response time and tracer concentrations

Using previous results on extended Petri Nets (EPN), we present the relations between various hydrological dynamical systems (HDSys) derived from the water budget. Once the water budget has been implemented, there is a consistent way of getting the equations for backward travel time distributions, for forward response time distributions and for the concentration of a solute or tracer. We show that the water budget has a correspondence of one to many with the backward travel time distributions. In fact, to any one of the water budget equations there correspond as many equations as there are input precipitation events. The backward travel time distributions are related to the response time distributions by Niemi's relationship and, if there are n outputs, by the definition of n − 1 partition functions. These determine what fraction of the water volume injected into the control volume at a specific time t_in goes asymptotically to a specific output. Given the backward travel time distributions, the output concentration of a solute or tracer also depends on the input concentration. The paper clarifies the complicated relations described above by taking [Hydrology and Earth System Sciences, 20, 299–328] as an example from literature. Once the appropriate information is presented, it is shown how these HDSys can be solved simultaneously without duplicating calculations. Then, it is also shown that, under the hypothesis of uniform mixing of water ages within each reservoir, these systems can be solved exactly.     

Effect of channel size on solute residence time distributions in rivers

The effect of channel size on residence time distributions (RTDs) of solute in rivers is investigated in this paper using tracer test data and the variable residence time (VART) model. Specifically, the investigation focuses on the influence of shear dispersion and hyporheic exchange on the shape of solute RTD, and how these two transport processes prevail in larger and smaller streams, respectively, leading to distinct tails of RTD. Simulation results show that (1) RTDs are dispersion-dependent and thereby channel-size (scale) dependent. RTDs increasing longitudinal dispersion coefficient. Small streams with negligible dispersion coefficient may display various types of RTD from upward curving patterns to a straight line (power-law distributions) and further to downward curving lognormal distributions when plotted in log–log coordinates. Moderate-sized rivers are transitional in terms of RTDs and commonly exhibit lognormal and power-law RTDs; (2) the incorporation of water and solute losses/gains in the VART model can improve simulation results and make parameter values more reasonable; (3) the ratio of time to peak concentration to the minimum mean residence time is equal to the recovery ratio of tracer. The relation provides a simple method for determining the minimum mean residence time; and (4) the VART model is able to reproduce various RTDs observed in rivers with 3–4 fitting parameters while no user-specified RTD functions are needed.     

Stochastic modeling of solute transport in 3-D heterogeneous porous media with random source condition

During probabilistic analysis of flow and transport in porous media, the uncertainty due to spatial heterogeneity of governing parameters are often taken into account. The randomness in the source conditions also play a major role on the stochastic behavior in distribution of the dependent variable. The present paper is focused on studying the effect of both uncertainty in the governing system parameters as well as the input source conditions. Under such circumstances, a method is proposed which combines with stochastic finite element method (SFEM) and is illustrated for probabilistic analysis of concentration distribution in a 3-D heterogeneous porous media under the influence of random source condition. In the first step SFEM used for probabilistic solution due to spatial heterogeneity of governing parameters for a unit source pulse. Further, the results from the unit source pulse case have been used for the analysis of multiple pulse case using the numerical convolution when the source condition is a random process. The source condition is modeled as a discrete release of random amount of masses at fixed intervals of time. The mean and standard deviation of concentration is compared for the deterministic and the stochastic system scenarios as well as for different values of system parameters. The effect of uncertainty of source condition is also demonstrated in terms of mean and standard deviation of concentration at various locations in the domain.     

Time and Magnitude of Peak Concentration of Reactive Groundwater Contaminants Discharged to a River

An analytical solution for calculating the contaminant discharge rate in an aquifer following an instantaneous release of reactive contaminant mass to groundwater is used to derive relationships for the time and magnitude of peak concentration in a river receiving the transported material. Relationships are developed for the time of peak concentration relative to the time of travel for the contaminant, and the magnitude of peak concentration relative to the concentration calculated at the time of travel. Both quantities are found to be a function of two dimensionless parameters characterizing advective‐dispersive‐reactive transport—the Peclet number and the Damkohler number. It is shown that the time to peak concentration may occur before the time of travel, considering advection and retardation only, depending on the magnitudes of the Peclet and Damkohler numbers. Similarly, the magnitude of peak concentration may exceed the concentration calculated assuming that the time of peak concentration coincides with time of travel for the contaminant. For large Damkohler numbers, equating the time of peak concentration with the time of travel for the contaminant can significantly underestimate peak concentrations.     

Signification de la présence de radon 222 dans les fluides géothermiques

The²²²Rn content of a subterranean fluid depends on the physico-chemical nature of the media through which the fluid passes and on the time spent by the fluid in the various geological formations encountered. Consequently,²²²Rn measurement in a geofluid gives a clue to the subsurface transport conditions of this geofluid. Simple models are proposed for various types of buried structures connected with the surface by a geothermal well; they all lead to the²²²Rn concentration expressed as a function of the flow-rate of the carrier fluid. In real systems, the variations of this concentration with flow-rate would allow the type of structure feading the well to be determined, and, in propitious circumstances, the volume of the possible reservoir to be estimated.     

Equilibration Times of Gas-Filled Diffusion Samplers in Slow-Moving Ground Water Systems

The installation of gas-filled diffusion samplers into small-diameter boreholes results in a significant reduction of the dissolved gas concentration around the sampler. In aquifers where the diffusive flux of solutes outpaces advective transport, the process that governs the equilibration time of a sampler is the resupply of solutes by diffusion from the aquifer. We have derived a solution that can be used to estimate the time required for a diffusion sampler to reach equilibrium with the dissolved gas concentration in the aquifer, where diffusion is the only solute transport mechanism. Thus the solutions provide equilibration times for cells placed in aquifers where diffusion dominates and maximum equilibration times for cells placed in aquifers where advection can also resupply solutes. The solutions are generic and are functions of nondimensionalized variables, therefore providing estimates of equilibration times for any type of solute, sampler volume, bore dimensions, and aquifer porosity. Examples are given for various sized gas-filled helium samplers placed in boreholes of different radii.