Two-Dimensional Spectral Estimation: A Radon Transform Approach

A new technique for two-dimensional (2-D) spectral estimation of a stationary random field (SRF) is investigated in this paper. This is based on the extension of the Radon transform theory to stationary random fields (SRF's), proposed by Jain and Ansari [19]. Using the Radon transform, the 2-D estimation problem is reduced to a set of one-dimensional (1-D) independent problems, which could then be solved using 1-D linear prediction (LP) or by any other high-resolution estimation procedure. This is unlike previous methods which obtain the 2-D power spectral density OPSD) estimate by using 1-D high-resolution techniques in the spirit of a separable estimator [2]. Examples are provided to illustrate the performance of the new technique. Various features of this approach are highlighted.     

A novel technique in analyzing non-linear wave-wave interaction

During wave growth non-linear wave–wave interactions cause transfer of some wave energy from lower to higher wave periods as the spectrum grows. Wavelet bicoherence, which is a new technique in the analysis of wind–wave and wave–wave interactions, is used to analyze non-linear wave–wave interactions. A selected record of wind wave that contains the maximum wave height observed during 6 h of wave generation is divided into five segments and wavelet bicoherence is computed for the whole record, and for all divided segments. The study shows that the non-linear wave–wave interaction occurs at different bicoherence levels and these levels are different from one segment to another due to the non-stationarity feature of the examined data set.     

Wavelet bicoherence analysis of wind–wave interaction

A new technique in the analysis of wind–wave interaction, wavelet bicoherence, will be applied in this article. Wavelet bicohence has the ability to detect phase coupling and nonlinear interactions of the quadratic order with time resolution. It is used in this study to analyze wind–wave interaction during wave growth in a Mistral event. A selected record of simultaneously measured wind and wave data during Mistral is divided into five segments and the computations of the wavelet bicoherence are conducted for the whole record and for all divided segments. The results show that the phase coupling occurs between wind speed and wave height over a certain range of frequencies and that the range is different from one segment to another due to the non-stationary feature of the time series.     

平衡地质剖面技术发展现状与实际应用——以苏北盆地溱潼凹陷为例

平衡剖面技术广泛应用于构造解释合理性检验和构造演化史分析。平衡剖面发展历史与现状表明,二维平衡剖面技术理论与应用已趋于成熟,真三维平衡恢复算法的精确实现和实际应用是今后的主要发展方向。在详细总结二维平衡剖面技术原理和制作过程后,以苏北盆地溱潼凹陷为例,紧密结合区域地质与地球物理资料、选择合理的剖面走向、采用变速时深转换,多次修改地震构造解释方案并最终制作平衡剖面。结果表明,合理运用平衡剖面技术可检验地震构造解释的合理性并分析构造演化历史,指导油气资源勘探。     

Laboratory Study of the Nonlinear Transformation of Irregular Waves over A Mild Slope

This paper considers the nonlinear transformation of irregular waves propagating over a mild slope (1?40). Two cases of irregular waves, which are mechanically generated based on JONSWAP spectra, are used for this purpose. The results indicate that the wave heights obey the Rayleigh distribution at the offshore location; however, in the shoaling region, the heights of the largest waves are underestimated by the theoretical distributions. In the surf zone, the wave heights can be approximated by the composite Weibull distribution. In addition, the nonlinear phase coupling within the irregular waves is investigated by the wavelet-based bicoherence. The bicoherence spectra reflect that the number of frequency modes participating in the phase coupling increases with the decreasing water depth, as does the degree of phase coupling. After the incipient breaking, even though the degree of phase coupling decreases, a great number of higher harmonic wave modes are also involved in nonlinear interactions. Moreover, the summed bicoherence indicates that the frequency mode related to the strongest local nonlinear interactions shifts to higher harmonics with the decreasing water depth.     

Laboratory study of the nonlinear transformation of irregular waves over a mild slope

This paper considers the nonlinear transformation of irregular waves propagating over a mild slope （1：40）. Two cases of irregular waves, which are mechanically generated based on JONSWAP spectra, are used for this purpose. The results indicate that the wave heights obey the Rayleigh distribution at the offshore location; however, in the shoaling region, the heights of the largest waves are underestimated by the theoretical distributions. In the surf zone, the wave heights can be approximated by the composite Weibull distribution. In addition, the nonlinear phase coupling within the irregular waves is investigated by the wavelet-based bicoherence. The bicoherence spectra reflect that the number of frequency modes participating in the phase coupling increases with the decreasing water depth, as does the degree of phase coupling. After the incipient breaking, even though the degree of phase coupling decreases, a great number of higher harmonic wave modes are also involved in nonlinear interactions. Moreover, the summed bicoherence indicates that the frequency mode related to the strongest local nonlinear interactions shifts to higher harmonics with the decreasing water depth.     

裂隙识别技术及其在油气和水合物勘探中的应用

裂隙是油气藏和天然气水合物成藏成矿的有利疏导体系, 但是受地震分辨率的影响, 微裂隙不易用常规手段识别。为研发一套实用的裂隙识别技术, 文章分析了裂隙的地震波运动学和动力学特征, 结合相关技术功能原理, 融合三维可视化技术、相干技术等方法, 在莺歌海盆地的底辟构造中成功刻画出裂隙通道, 并据此技术发现大气田。经实践及分析, 该技术亦可应用于琼东南盆地天然气水合物的疏导体系研究, 是一套行之有效的裂隙识别技术。     

Experimental study of wave–wave nonlinear interactions using the wavelet-based bicoherence

The wavelet-based bicoherence, which is a new and powerful tool in the analysis of nonlinear phase coupling, is used to study the nonlinear wave–wave interactions of breaking and non-breaking gravity waves propagating over a sill. Two cases of mechanically generated random waves based on Jonswap spectra are used for this purpose. Values of relative depth, k_ph (k_p is the wave number of the spectral peak and h is the water depth) for this study range between 0.38 and 1.22. The variations of wavelet-based total bicoherence for the test cases indicate that the degree of quadratic phase coupling increases in the shoaling region consistent with a wave profile that is pitched shoreward, relative to a vertical axis as seen in the experiments, but decreases in the de-shoaling region. For the non-breaking case, the degree of quadratic phase coupling continues to increase until waves reach the top of the sill. Breaking waves, however, achieve their highest level of quadratic phase coupling immediately before incipient breaking and the degree of phase coupling decreases sharply following breaking. In addition the wavelet-based bicoherence spectra provide evidence of the harmonics' growth which is reflected in the energy spectra. The bicoherence spectra also show that quadratic phase coupling between modes within the peak frequency as well as between modes of the peak frequency and its higher harmonics are dominant in the shoaling region, even though there are relatively high levels of quadratic phase coupling occurring between other frequencies. Furthermore, using the temporal resolution property of the wavelet-based bicoherence, we find that the quadratic wave interactions occur more readily during segments of time with large change of wave amplitude, rather than those segments having large wave amplitudes, but small gradients in amplitude.     

An integral eigenvalue approach to three-dimensional field problems--A low-frequency variation of SEM

The three-dimensional (3-D) eddy-current transient field problem is formulated first using theu-vmethod. This method breaks the vector Helmholtz equation into two scalar Helmholtz equations. Null-field integral equations and the appropriate boundary conditions germane to the problem are used to set up an identification matrix which is independent of null-field point locations. Embedded in the identification matrix are the unknown eigenvalues of the problem representing its impulse response in time. These eigenvalues are found by equating the determinant of the identification matrix to zero. The eigenvalues, which can be equated with temporal response, are found to be intimately linked to the initial forcing function which triggers the transient in question. When this initial forcing function is Fourier decomposed into its respective spatial harmonics, it is possible to associate with each Fourier component a unique eigenvalue by this technique. The true transient solution comes through a convolution of the impulse response so obtained with the particular imposed external field governing the problem at hand. The technique is applied to the FELIX medium cylinder (a conducting cylinder placed in a collapsing external field) and compared to data. A pseudoanalytic confirmation of the eigenvalues so obtained is formulated to validate the procedure. The technique proposed is applied in the low-frequency regime where the near-field effects must be considered. Application of the technique to a high frequency follows directly if the Coulomb gauge is adopted to represent the vector potential.     

Planar SAS for Sea Bottom and Subbottom Imaging: Concept Validation in Tank

The paper deals with a new concept of imaging technique based on planar synthetic aperture sonar (P-SAS). P-SAS processing takes advantage of the overlap of both ping-to-ping and track to track footprints to improve the imaging resolution in both directions. After describing the method, its implementation is discussed. The exact implementation [three-dimensional (3-D)] is defined and an approximated and faster approach [2timestwo-dimensional (2-D)] is then proposed. The gain in computation time and the possible loss in performances (resolutions, echo amplitude) are investigated by simulating a simple point target. The point spread function is studied in details. The new concept is finally validated by tank experiments (scale about 100) for both proud and buried target imaging     

Pseudorandom realization of transient acoustic waves scattered fromrandomly rough patches

A simple numerical technique is developed for generating pseudorandom realizations of three-dimensional (3-D) transient acoustic waves that are scattered from two-dimensional (2-D) patches of randomly rough surfaces. The rough surface height of a patch is represented numerically in the 2-D horizontal wavenumber plane by choosing a scheme for interpolation between pseudorandom complex coefficients. Using this approach, the realizations of the patches can be generated from experimentally measured roughness power spectra, and phase information is generated in the frequency domain that leads to time spreads in the time domain. The acoustic scattering is modeled here with first-order perturbation theory. The boundary conditions considered here are pressure-release, rigid, and fluid-fluid. Three different spatial windows are considered for defining the patches. In the time domain, the time spreads of the scattered waveforms agree with predictions. In the frequency domain, the phase is seen as a random walk. The solutions developed here can be used with normal mode propagation models or ray propagation models     

Measurements of short water waves using stereo matched image sequences

Image analysis techniques are used for retrieving water surface elevation fields spatially and temporally from CCD-images and CCD-image-sequences. The technique proposed herein utilizes binocular stereogrammetry to recover topographic information from a sequence of synchronous, overlapping video images. The method used differs from the traditional stereo-photogrammetric analysis of a single stereo-pair because the use of video allows for a continuous sequence of stereo-images to be digitally sampled and analyzed. For data acquisition two synchronized progressive-scan cameras were used.A partially supervised 3-D stereo system (called WASS, Wave Acquisition Stereo System) is shown here. It is used to reconstruct the 3-D shape of water surface waves, acquired at frame rate, with small computational time needed. The stereo method is presented, including the derivation of a relationship relating the geometry of the stereo rig and the expected errors. Finally, the 3-D calculated scattered points give the complete spatio-temporal distribution of the water surface elevations. The measurable length-scales depend on the pixel resolution, the triangulation accuracy, and the acquisition frame rate. Limitations in the stereo measurements are also discussed.Two experiments to test and to demonstrate the system took place: one on the Venice lagoon, north of the city of Venice in September 2004 and the second on the coast of California at San Diego in December 2005. For the second experiment, qualitative and quantitative intercomparisons of the stereo-matching and in-situ sensor measurements are presented. All the measurements of water surface waves indicate that the proposed approach is both accurate and applicable for measuring water surface waves. Moreover, shape estimates are accurately and extremely dense both in space and time, and the remote location of the instruments alleviates some difficulties associated with in situ instrumentation.     

Tomographic inversion of the El Nino profile by using amatched-mode processing (MMP) method

A high-resolution mode-matching (HRMM) estimator is used to invert El Nino profile perturbations. By matching a proper set of modal travel time perturbations in a 2-D parameter space constructed by the first two coefficients of the empirical orthogonal function (EOF), the El Nino profile can be efficiently inverted. As contrasted with the conventional linear inversion scheme, the performance of the matched-mode scheme is significantly improved because the nonlinearity of the modal travel time perturbation has been taken into account     

A fast algorithm for high-accuracy frequency measurement:application to ultrasonic Doppler sonar

In an attempt to investigate the technical feasibility of a CW Doppler sonar, we have examined a method of measuring low velocities with a high-velocity resolution, or frequency resolution, by use of a simple circuit configuration employing digital signal processing technique. The following discussion presents the results of the investigation. In the measuring method described, the fast Fourier transform (FFT) of undersampled data is calculated and the Doppler shift is obtained by searching for a peak frequency of the power spectrum. To achieve the intended frequency resolution of 1 Hz by FFT operation, measurement of data for a minimum measuring period of 1 s is essential. If the sampling frequency is set to 50 kHz, the number of samples obtained during the minimum measuring period of I s would amount to 50000. This is not practical in light of the time required for the FFT operation. To overcome this problem, our new measuring method employs a decimation technique for reducing the number of samples down to 1024 while maintaining a frequency resolution of about 1 Hz. This paper describes how the processing time can be drastically reduced to about 1/300th compared to the conventional technique by a combination of complex exponential functions, filtering and decimation, and thereby indicates the possibility of real-time CW Doppler data processing     

洞庭湖地区水系水动力耦合数值模型

针对洞庭湖地区复杂水系和复杂水情下的实际水流运动特征,建立了完全基于水力学方法的洞庭湖地区水系的一、二维耦合的全局水动力学模型,实现了不同水体数值模拟的自动有机衔接。模型包括湖泊、河网水系两部分。对长江干流荆江段、湘资沅澧四水、三口分流河道,运用能够反映河道主槽和边滩不同行洪特征的扩展一维水量模型进行模拟;对于洞庭湖湖区部分,采用二维非结构的有限体积法建立水动力模型;应用重叠投影法实现模型的耦合。模型具有动边界的自动处理,河道分区计算,分洪、溃堤过程实时模拟等功能。剖分的网格较为细致地刻画了湖区的地形,使模型能够模拟"高水湖相、低水河相"的湖流特征。运用1996年7月的实际洪水过程,进行了模型的验证,较为准确地模拟了实际水流特征。     

Waveform Shaping of Sonar Transducers for Improving the Vertical Resolution in Sub-bottom Sediments Profiling

Vertical resolution is of fundamental importance in sonar exploration and is directly related to the duration of the acoustic pulse generated by the transducer. The shorter the radiated pulse, the higher the vertical resolution. Many sub-bottom profiling sonar systems use piezoelectric transducers because they are reversible and well understood. Piezoelectric projectors are normally resonant transducers, which are intrinsically narrowband. A piezoelectric transducer is usually driven by a tone-burst. However, it is possible to use Fourier techniques to find a pre-compensated electrical driving function so that the transducer radiates a prescribed wider band acoustic waveform. This technique can be applied to synthesize zero-phase cosine-magnitude, Gaussian, and bionic pulses, with a conventional sandwich transducer. Zero-phase cosine-magnitude waveforms provide minimum length pulses (and therefore maximum resolution) within a prescribed frequency band.The aim of this paper is to illustrate the synthesis of wideband acoustic pulses using an underwater piezoelectric projector. The conventional acoustic waveform radiated when a Tonpiltz transducer is transiently excited using a “click” and allows its frequency response function to be measured. This function is used to design the electrical signal which then drives the transducer so that it radiates the shortest pulse compatible with its mechanical response. The significant resolution enhancement of the waveform shaping process is illustrated by its application to a sediment wedge model.     

床面冲淤计算机处理方法

冲淤计算是河流，河口，海岸等水域中分析，预测床面变化的基础，利用非结构网络自动剖分技术形成覆盖计算区域的三角网，能较好地适应边界和特殊地形，基于高程内插建立三角网数字高程模型，采用垂直三棱柱体作为基本单位计算不同时段冲淤量，从而提高计算精度，通过网络剖分，利用在计算区域内任意选择子域进行淤量的计算和分析。同时将生成的非结构网格直接有于二维床面冲淤数学模型，便于对每个节点的计算结果进行检证以提高模型的可靠性。     

基于八向链码的矩的快速计算及2—D目标识别

提出一种直接基于八向链码的矩的快速计算方法，给出了与多边形近似法矩的计算及形状识别时间比较，最后使用不变性矩和人工神经网络完成2—D形状识别。实验结果表明，本文所提出的方法具有更高的计算和识别速度。     

Temporal and Spatial Sampling Constraints in High-Resolution 3-D Seismic Surveys

An experimental, high-resolution 3-D seismic survey was acquired over a 1.5 km2 section of an incised fluvial valley. The data were acquired as a near-zero offset, single-channel survey using a 15 in3 water gun as the source, and differential GPS for navigation and positioning. The objective was to acquire a 3-D seismic data volume suitable for calculating the volume of shallow sand deposits. Horizontal time sections from the 3-D volume clearly show the flanks of the incised valley, as well as high-amplitude reflections interpreted as coarse-grained channel-lag deposits. The volume of this lag deposit can be calculated using the combination of the horizontal and vertical sections from the high-resolution 3-D seismic data set. The results of the experiment also illustrate the importance of spatial sampling in 3-D seismic surveying.