密度和压缩系数的散射层析成像法

本文在速度成像的基础上研究了同时对密度和压缩系数成像的散射波层析成像法.对不同散射角度的计算可以得到一系列反演图像,拟合这些图像,从而可以有效地达到对密度和压缩系数（或速度）成像的目的.与单纯的速度成像相比,增加了反演的难度.首先是对资料的方位性要求增加;其次是对资料的利用率下降.即便如此,从对较少量的炮点和检波点资料的数值计算来看,仍取得了满意的成像结果.我们对组成字母“A”的散射体结构进行了成像计算,结果能够同时再现密度和压缩系数,成像清晰,表明了方法的可行性,并能应用于复杂结构的成像问题.     

Acoustic plane-wave scattering from a rough surface over a random fluid medium

Acoustic plane-wave scattering from a rough surface overlying a fluid half-space with a sound-speed distribution subject to a small random variation is considered. Under the assumption that the surface roughness and medium randomness are statistically independent, the scattered field may be derived by first solving for the mean field in the medium, and then incorporating with boundary-perturbation method to obtain the total mean field and the power spectral density of the scattered field. The employed algorithm is compatible to the analysis available in the existing literature so that the formulations are conveniently integrated. The results for the power spectral density have shown that the effects of medium inhomogeneities on the rough surface scattering are limited in a spectral regime where the scattered components have shallow grazing angles. The distribution of the power spectral density over the space is primarily governed by the power spectrum and correlation lengths of the rough surface.     

Unified approach to amplitude attenuation and coda excitation in the randomly inhomogeneous lithosphere

A unified model is proposed for explaining the frequency dependent amplitude attenuation and the coda wave excitation on the basis of the single scattering process in the randomly inhomogeneous lithosphere. Adopting Birch's law and a direct proportion between density and wave velocity, we statistically describe the inhomogeneous medium by one random function characterized by the von Karman autocorrelation function. We calculate the amplitude attenuation from the solid angle integral of scattered wave energy on the basis of the Born approxiimation after subtracting the travel-time fluctuation effect caused by slowly varying velocity inhomogeneities. This subtraction is equivalent to neglect energy loss by scattering within a cone around the forward direction. The random inhomogeneity of the von Karman autocorrelation function of order 0.35 with the mean square fractional fluctuation of 7.2×10^–3 1.3×10^–2 and the correlation distance of 2.15.1 km well explains observed backward scattering coefficientg and the ratioQ _P ^–1 /Q _S ^–1 , and observed and partially conjecturedQ _S ^–1 for frequencies between 0.5 Hz and 30 Hz.     

Probing mid-mantle heterogeneity using PKP coda waves

We investigate the utility of PKP coda waves for studying weak scattering from small-scale heterogeneity in the mid-mantle. Coda waves are potentially a useful probe of heterogeneity in the mid-mantle because they are not preferentially scattered near the CMB, as PKP precursors are, but are sensitive to scattering at all depths. PKP coda waves have not been used for this purpose historically because of interference with other late-arriving energy due to near-surface resonance and scattering. Any study of deep mantle scattering using coda waves requires the removal of near-surface effects from the data. We have analyzed 3624 recordings of PKP precursors and coda made by stations in the Incorporated Research Institutions for Seismology (IRIS) Global Seismographic Network (GSN). To study the range and time dependence of the scattered waves, we binned and stacked envelopes of the recordings. We have considered precursors that arrive within a 20 s window before PKP and coda waves in a 60 s window after PKP. The PKP scattered waves increase in amplitude rapidly with range as predicted by scattering theory. At ranges below 125°, we predict and observe essentially no scattered energy preceding PKP. Coda amplitudes at these ranges are independent of range and provide an estimate of energy due to near-surface effects that we can expect at all ranges. We use the average coda amplitude at ranges from 120 to 125° to correct coda amplitudes at other ranges. PKP coda waves show a strong dependence on time and range and are clearly influenced by scattering in the lower mantle. PKP coda waves, however, do not provide a tighter constraint on the vertical distribution of mantle heterogeneity than is provided by precursors. This is due, in part, to relatively large scatter in coda amplitudes as revealed by a resampling analysis. Modeling using Rayleigh–Born scattering theory and an exponential autocorrelation function shows that PKP coda amplitudes are not highly sensitive to the vertical distribution of heterogeneity in the mantle. To illustrate this we consider single-scattering in two extreme models of mantle heterogeneity. One allows heterogeneity just at the CMB; the other includes heterogeneity throughout the mantle. The amplitudes of precursors are tightly constrained by our stack and support our earlier conclusion that small-scale heterogeneity is uniformly distributed throughout the lower mantle. The best-fit model includes 8 km scale length heterogeneity with an rms velocity contrast throughout the mantle of 1%.     

Interaction of water waves with vertical cylinders using null-field integral equations

The scattering of water waves by bottom-mounted vertical circular cylinders is solved by using the null-field integral equations in conjunction with degenerate kernels and Fourier series to avoid calculating the Cauchy and Hadamard principal values. In the implementation, the null-field point can be exactly located on the real boundary owing to the introduction of degenerate kernels for fundamental solutions. An adaptive observer system of polar coordinates is considered to fully employ the properties of degenerate kernels. For the hypersingular equation, vector decomposition for the radial and tangential gradients is carefully considered. This method can be seen as a semi-analytical approach since errors attribute from the truncation of Fourier series. Neither hypersingularity in the Burton and Miller approach nor the CHIEF concept was required to deal with the problem of irregular frequencies. Five advantages of free of calculating principal value, well-posed algebraic system, convergence rate of exponential order, meshfree and elimination of boundary-layer effect, are achieved by using the present approach. Numerical results are given for the forces and free-surface elevation around the circular boundaries. Also, the near-trapped behavior arisen from the physical resonance is detected. A general-purpose program for water wave impinging several circular cylinders with arbitrary number, radii, and positions was developed. Several examples of water wave structure interaction by vertical circular cylinders were demonstrated to see the validity of the present formulation.     

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.     

植被电磁散射的半空间模型研究

本文给出了地面植被电磁散射的半空间模型研究方法.在以往的相关文献中,均采用自由空间格林函数求解单个叶片散射体的散射场,本文利用半空间并矢格林函数求出了单个散射体的散射场,然后运用Monte Carlo方法模拟生成地面植被层,得到了在半空间下植被层的电磁散射特性,计算出了单、双站雷达散射截面,并与自由空间下的散射场做了相互对比.结果表明,在半空间格林函数下得到的散射场更为适用于描述地面植被的电磁散射特性.     

Meteorology and haze structure during AGASP-II,Part 1: Alaskan Arctic flights, 2–10 April 1986

The second Arctic Gas and Aerosol Sampling Program (AGASP-II) was conducted across the Alaskan and Canadian Arctic in April 1986, to study the in situ aerosol, and the chemical and optical properties of Arctic haze. The NOAA WP-3D aircraft, with special instrumentation added, made six flights during AGASP-II. Measurements of wind, pressure, temperature, ozone, water vapor, condensation nuclei (CN) concentration, and aerosol scattering extinction (b_sp) were used to determine the location of significant haze layers. The measurements made on the first three flights, over the Arctic Ocean north of Barrow and over the Beaufort Sea north of Barter Island, Alaska are discussed in detail in this report of the first phase of AGASP II. In the Alaskan Arctic the WP-3D detected a large and persistent region of haze between 960 and 750 mb, in a thermally stable layer, on 2, 8, and 9 April 1986. At its most dense, the haze contained CN concentrations >10,000 cm^–3 and b_sp of 80×10^–6 m^–1 suggesting active SO₂ to H₂SO₄ gas-to-particle conversion. Calculations based upon observed SO₂ concentrations and ambient relative humidities suggest that 10⁴–10⁵ small H₂SO₄ droplets could have been produced in the haze layers. High concentrations of sub-micron H₂SO₄ droplets were collected in haze. Ozone concentrations were 5–10 ppb higher in the haze layers than in the surrounding troposphere. Outside the regions of haze, CN concentrations ranged from 100 to 400 cm^–3 and b_sp values were about (20–40)×10^–6 m^–1. Air mass trajectories were computed to depict the air flow upwind of regions in which haze was observed. In two cases the back trajectories and ground measurements suggested the source to be in central Europe.     

Multiple scattering of seismic waves in fractured media: Velocity and effective attenuation of the coherent components ofP waves andS waves

We show that the multiple scattering by small fractures of seismic waves with wavelengths long compared to the fracture size and fracture spacing is indistinguishable from multiple-scattering effects produced by regular porosity, except for an orientation factor due to fracture alignment. The fractures reduce theP-wave andS-wave velocities and produce an effective attenuation of the coherent component of the seismic waves. The attenuation corresponds to 1000/Q of about unity for a Gaussian spectrum of fractures, and it varies with frequencyf asf ³. For a Kolmogorov spectrum of fractures of spectral index the attenuation is an order of magnitude or so larger and varies with frequency asf ^3-v The precise degree of attenuation depends upon the matrix properties, the fracture porosity, the degree of fracture anisotropy, the type of fluid filling the fractures, and the incidence angle of the wave.For fracture porosities less than about 15% theP-wave andS-wave velocities are decreased by the order of 5–10% with a lesser dependence on the type of fluid filling the fractures (gas, oil, or brine) and with a dependence on both the degree of anisotropy and the incident angle made by the wave. The tendency of fractures to occur perpendicularly to bedding suggests that the best way to measure seismically fractured rock behavior in situ is by using the travel-time delay and reflection amplitude. As both the offset and the azimuth of receivers vary from a shot, the travel-time delay and reflection amplitude should both show an elliptical pattern of behavior—the travel-time delay in response to the varying seismic speed, and the reflection amplitude in response to angular variations in the multiple scattering. Observations of attenuation at several frequencies should permit (a) determination of the spectrum of fractures (Gaussian versus Kolmogorovian) and (b) determination of the contribution of viscous damping to the effective attenuation.     

圆弧型沉积盆地对平面SH波的散射

采用波函数展开方法给出圆弧型沉积盆地对面ＳＨ波二维散射问题的封闭级数解答。利用外域Ｇｒａｆ加法公式将解答归结为无穷代数方程组的求解。通过截断计算得到解答的数值结果，并通过连续条件的满足程度检验了截断计算的精度。将本文结果已有的近似解析解进行比较，指出了近似解答的误差来源和适用范围。给出一些典型算例说明盆地深度比对地震动的复杂作用。