Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

试论我国城市发展方针

在执行重重工业的发展战略和城乡分割体制的计划经济背景下,我国制定了“控制大城市规模,合理发展中等城市,积极发展小城市”的城市建设方针。但是分析表明,这一方针不能适应市场经济和对外开放形势的需要,也不符合我国人多地少、农村剩余劳动力多的国情,必须调整现有的城市发展方针,选择适宜的城市化发展道路。因此,本文就适度发展大中城市,提高城市现代化水平;加速乡村城市化,积极合理发展小城镇;以及根据不同地区的具体情况,因地制宜地指定城市化途径等方面作了探讨。     

Oil and gas resource assessment: The linkage between geological analysis and discovery process model results

Two widely-used techniques to estimate the volume of remaining oil and gas resources are discovery process modeling and geologic assessment. Both were used in a recent national assessment of oil and gas resources of the United States. Parallel estimates were obtained for 27 provinces. Geological-based estimates can typically see into areas not available to discovery process models (that is areas with little or no exploration history) and thus, on average, yield higher estimates. However, a linear relation does exist between the mean estimates obtained from these two methods. In addition, other variables were found in a multiple regression model that explained much of the difference. Thus, it is possible to perform discovery process modeling and adjust the estimates to yield results that might be expected from geological-based assessments.     

The origin of accretionary lapilli

Experimental investigations in a recirculating wind tunnel of the mechanisms of formation of accretionary lapilli have demonstrated that growth is controlled by collision of liquid-coated particles, due to differences in fall velocities, and binding as a result of surface tension forces and secondary mineral growth. The liquids present on particle surfaces in eruption plumes are acid solutions stable at 100% relative humidity, from which secondary minerals, e.g. calcium sulphate and sodium chloride, precipitate prior to impact of accretionary lapilli with the ground. Concentric grain-size zones within accretionary lapilli build up due to differences in the supply of particular particle sizes during aggregate growth. Accretionary lapilli do not evolve by scavenging of particles by liquid drops followed by evaporation — a process which, in wind tunnel experiments, generates horizontally layered hemispherical aggregates. Size analysis of particles in the wind tunnel air stream and particles adhering to growing aggregates demonstrate that the aggregation coefficient is highly grain-size dependent. Theoretical simulation of accretionary lapilli growth in eruption plumes predicts maximum sizes in the range 0.7–20 mm for ash cloud thicknesses of 0.5–10 km respectively.     

地震综合效应场函数及其在地震预报中的应用──（一）地震综合效应场函数及其计算方法

在系统地分析了目前各种测震学地震预报方法科学思路的基础上,认为测震学地震预报方法基本上可以分为两大类。一类是以已经发生的一些地震作为未来可能发生的地震的“因”,即由于已经发生的地震对区域应力场的影响,导致未来发生较强地震。这一类包括的预报方法较多,如空区、条带、ｂ值、地震迁移、相关地震等等及其由此衍生出来的各种方法。另一类是把已经发生的一些地震作为区域应力场增强的“果”,即已经发生的地震是区域应力场增强过程中的一种反映,而未来地震不一定是已经发生的地震所导致的结果。这一类包括“地震窗口”、小震群活动等方法。针对第一类方法,各种预报方法都是力图从地震三要素中提取未来地震的信息,而具体作法又都是利用地震三要素这个多维空间的某个剖面。为了从地震活动诸要素的多维空间提取综合信息,我们对每个地震加入了破裂面方位,构成了地震第四要素,并依据地震４要素建立了地震综合效应场函数。地震综合效应场函数概括了多种测震学地震预报方法的科学思路和预报经验,从而可以形成测震学的综合预报方法。     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

人工神经网络在白云岩类型识别中的应用

利用地球化学数据,运用人工神经网络方法对美国密苏里州东南Ｂｏｎｎｅｔｅｒｒｅ组（寒武纪）滨海相的白云岩进行了分类、识别,判别率达１００％,结果表明,该方法性能良好,可望成为岩石分类、判别的有效手段。     

A dynamic box model of bioactive elements in the southern Taiwan Strait

A dynamic box model was applied to study the characteristics of biogeochemical cycling of PO_4-P,NO_3-N,AOU,POC and PON in the southern Taiwan Strait region based on field data of the"Minnan Taiwan Bank Fishing Ground Upwelling Ecosystem Study" during the period of Dec.1987-Nov.1988.According to the unique hydrological and topographical features of the region,six boxesand three layers were considered in the model.The variation rates and fluxes of elements induced by hor-izontal current,upwelling,by diffusion,sinking of particles and biogeochemical processes were estimatedrespectively.Results further confirmed that upwellings had important effects in this region.Thenearshore upwelling areas had net input fluxes of nutrients brought by upwelling water,also had high de-pletion rates of nutrients and production rates of particulate organic matter and dissolved oxygen.Theabnormal net production of nutrients in the middle layer(10-30m) indicated the important role of bacte-ria in this high production region.Th     

An integrable case of a rotational motion analogous to that of Lagrange and Poisson for a gyrostat in a Newtonian force field

The scope of the present paper is to provide analytic solutions to the problem of the attitude evolution of a symmetric gyrostat about a fixed point in a central Newtonian force field when the potential function isV ⁽²⁾.We assume that the center of mass and the gyrostatic moment are on the axis of symmetry and that the initial conditions are the following: (t ₀)=₀, (t ₀)=₀, (t ₀)=(t ₀)=₀, ₁(t ₀)=0, ₂(t ₀)=0 and ₃(t ₀)= ₃ ⁰ .The problem is integrated when the third component of the total angular momentum is different from zero (B ₁ 0). There now appear equilibrium solutions that did not exist in the caseB ₁=0, which can be determined in function of the value ofl ₃ ^r (the third component of the gyrostatic momentum).The possible types of solutions (elliptic, trigonometric, stationary) depend upon the nature of the roots of the functiong(u). The solutions for Euler angles are given in terms of functions of the timet. If we cancel the third component of the gyrostatic momentum (l ₃ ^r =0), the obtained solutions are valid for rigid bodies.