Comparison of Seismic Dispersion and Attenuation Models

The frequency-dependent attenuation of seismic waves causes decreased resolution of seismic images with depth, and the difference in transmission losses induces amplitude variations with offset. Transmission losses may occur due to friction or fluid movement, or may result from scattering in thin-layer. Whatever the physical mechanism, they can often be conveniently described using an empirical formulation wherein the elastic moduli and propagation velocity are complex functions of frequency.We have compiled and compared algebraically and numerically eight different models involving complex velocity: the Kolsky-Futterman model, the power-law model, Kjartansson's model, Müller's model, Azimi's second and third model, the Cole-Cole model, and the standard linear-solid model.For two different parameter sets, the attenuation and phase velocity are computed in the seismic frequency band, and the plane-wave propagation of a Ricker wavelet for the other models is compared with that for the Kolsky-Futterman model. The first parameter set consists of parameters for each of the models calculated from expressions given in the appendix. These expressions make the different models behave similarly to the KF model. The second parameter set consists of model parameters that are numerically adapted to the KF model.By selecting proper parameters, all models, except the standard linear-solid model, show behavior similar to that of the Kolsky-Futterman model. The SLS model behaves differently from the other models as the frequency goes to zero or infinity. Broadband measurement data is needed to select a specific model for a given seismic experiment.     

Size Scaling of Cross Correlation Between Multiple Variables

Reservoir models have large uncertainty because of spatial variability and limited sample data. The ultimate aim is to use simultaneously all available data sources to reduce uncertainty and provide reliable reservoir models for resource assessment and flow simulation. Seismic impedance or some other attribute provides a key source of data for reservoir modeling. These seismic data are at a coarser scale than the hard well data and it not an exact measurement of facies proportions or porosity. A requirement for data integration is the cross-covariance between the well and seismic data.The size-scaling behavior of the cross correlation for different measurement scales was nvestigated. The size-scaling relationship is derived theoretically and validated by numerical studies (including an example with real data). The limit properties of the cross-correlation coefficient when the averaging volume becomes large is shown. After some averaging volume, the volume-dependent cross-correlation coefficient reaches a limit value. This plateau value is controlled mainly by the large-scale behavior of the cross and direct variograms.The cross correlation can increase or decrease with volume support depending on the relative importance of long- and short-scale covariance structures. If the direct and cross variograms are proportional, there is no change in the cross correlation as the averaging volume changes. Our study shows that the volume-dependent cross-correlation coefficient is sensitive to the shape of the cross variogram and differences between the direct variograms of the well data and seismic data.     

银川市初夏气溶胶粒子谱分布特征

利用美制APS—3310A型激光粒谱仪，在中国西北城市银川对大气气溶胶粒子进行现场观测，结合气象条件，并通过理论计算分析，给出了银川市初夏大气气溶胶的日变化特征、3种不同粒径粒子分布特点及飘尘谱分布特征，以期为银川市大气污染总量控制方案及对策提供科学依据。     

Statistic forecast of typhoons going over the Chinese coasts

According to the measured data of typhoons going over the Chinese coasts in 1949-2002, a statistic relative equation showing the relation between the central atmospheric pressure of typhoons in a certain region at a certain period of time and their accumulation of frequency is established, and the concept of recurrence interval of typhoons is put forward, which is of actual significance for typhoon disaster reduction along the coastal area.     

桩网二元组合在加固软基中的应用研究

工程实践表明对软土地基工后沉降和不均匀沉降有特殊要求的构筑物,采用桩网加固地基是一种经济有效的处理方法。由于其作用机理、荷载传递方式是一个复杂的相互作用问题,本文运用应力扩散原理估算筋材承受的荷载,并由此计算筋材的拉力和割线模量,为桩网法设计提供理论依据。     

浅水方程数值计算方法的研究

求解以水位为变量的连续方程,并根据Navier-Stokes方程压力修正算法的基本思想,建立了浅水方程的水位修正算法,放宽了对离散时间步长的限制.通过对离散方程系数矩阵的重新构造,建立了高分辨率有限元格式,该格式既具有较高的离散精度又避免了数值解的伪振荡.对动量方程的阻力项做负坡线性化处理,提高了露滩计算的稳定性.数值模拟结果与解析解吻合良好,表明所建立的数值计算方法是正确的和可靠的.     

Modeling of Rayleigh waves dispersion in the Sannio region (Southern Italy) from seismic active refraction data

Short period surface waves, recorded during a seismic refractionsurvey in the Sannio region (Southern Italy), have been modeled to infera shallow velocity model for the area. Based on the decrease of resolutionwith depth, due to the bias on group velocity estimates arising frominterference of the Rayleigh waves with higher modes, we carried out aprocedure of fitting, with synthetic seismograms, of selected filtered traceswith a gaussian filter, having a width at half height equal to 1 Hz and acentral frequency lying in the range [1,4] Hz. We estimated the likelihoodbetween synthetic and observed seismograms by measuring their semblance.In this way we were able to infer a more refined local velocity modelcharacterized by a high Vp and Vs vertical gradient in the sedimentarycover. Two ad hoc resolution studies, based on group velocity andamplitude data respectively, indicate that the local velocity model is a goodvelocity model also for the entire studied area. The increase in the numberof available data when using amplitude information allows us to make amore selective choice in the model parameter space (Vp and Vs of eachlayer) and to solve for the Vp/Vs ratio. The inferred Vp velocity in thehalf-space is equal to 2.8 km/s. This value is in excellent agreement withthat inferred by other authors (3 km/s) by modeling P-wave travel timevs. distance. The best-fit model furnish low Vp/Vs for the sedimentarycover so indicating a high degree of the sediment's compaction in thestudied area. The inferred shallow high-velocity gradient indicates thatthe shallow sedimentary layer in the area could trap and focus the energytraveling into it.     

Atmospheric nitrogen inputs into the North Sea: effect on productivity

The ANICE (Atmospheric Nitrogen Inputs into the Coastal Ecosystem) project addressed the atmospheric deposition of nitrogen to the North Sea, with emphasis on coastal effects. ANICE focused on quantifying the deposition of inorganic nitrogen compounds to the North Sea and the governing processes. An overview of the results from modelling and experimental efforts is presented. They serve to identify the role of the atmosphere as a source of biologically essential chemical species to the marine biota. Data from the Weybourne Atmospheric Observatory (UK) are used to evaluate the effect of short episodes with very high atmospheric nitrogen concentrations. One such episode resulted in an average deposition of 0.8 mmol N m⁻² day⁻¹, which has the potential to promote primary productivity of 5.3 mmol C m⁻² day⁻¹. This value is compared to long-term effects determined from model results. The total calculated atmospheric deposition to the North Sea in 1999 is 948 kg N km⁻¹, i.e. 0.19 mmol N m⁻² day⁻¹ which has the potential to promote primary productivity of 1.2 mmol C m⁻² day⁻¹. Detailed results for August 1999 show strong gradients across the North Sea due to adjacent areas where emissions of NO_x and NH₃ are among the highest in Europe. The average atmospheric deposition to the southern part of the North Sea in August 1999 could potentially promote primary production of 2.0 mmol C m⁻² day⁻¹, i.e. 5.5% of the total production at this time of the year in this area of the North Sea. For the entire study area the atmospheric contribution to the primary production per m² is about two-third of this value. Most of the deposition occurs during short periods with high atmospheric concentrations. This atmospheric nitrogen is almost entirely anthropogenic in origin and thus represents a human-induced perturbation of the ecosystem.     

Estimating the orientation error of the Dahan Downhole Accelerometer using the maximum cross-correlation coefficient between the observed and synthetic waves

A new algorithm to correct the orientation error of the accelerometerat the Dahan Downhole Array, Hualien, Taiwan is developed. This algorithmconsists of three stages: (1) rotating two horizontal ground motions on thefree surface to the SH-SVdirection and SH axis offers a reference direction.(2) computing the synthetic downhole SH waves at a downhole station and (3)searching a rotation angle for downhole observation that yields a best waveformmatch between the synthetic and observed downhole seismograms. At this point, the rotated angle corresponding to the best waveform match can be considered as the orientation error. We selected five earthquakes with good data qualityfor analysis. Results show that this algorithm gives a more stable estimationthan a conventional method because it allows the selection of data from a wider time window for analyses. The estimated orientation error of the accelerometers at the Dahan Downhole Array after the 1999 reinstallation are40°, 114° and 285° at depths of 50, 100 and 200 m, respectively.