Determining the effective elastic thickness through cross-correlation between isostatic disturbances

The elastic thickness parameter was estimated using the mobile correlation technique between the observed isostatic disturbance and the gravity disturbance calculated through direct gravimetric modeling. We computed the vertical flexure value of the crust for a specific elastic thickness using a given topographic dataset. The gravity disturbance due to the topography was determined after the calculation. A grid of values for the elastic thickness parameter was generated. Then, a moving correlation was performed between the observed gravity data(representing actual surface data) and the calculated data from the forward modeling. The optimum elastic thickness of the particular point corresponded to the highest correlation coefficient. The methodology was tested on synthetic data and showed that the synthetic depth closely matched the original depth, including the elastic thickness value. To validate the results, the described procedure was applied to a real dataset from the Barreirinhas Basin, situated in the northeastern region of Brazil. The results show that the obtained crustal depth is highly correlated with the depth from known models. Additionally, we noted that the elastic thickness behaves as expected, decreasing from the continent towards the ocean. Based on the results, this method has the potential to be employed as a direct estimate of crustal depth and elastic thickness for any region.

An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

People Recognition by RGB and NIR Analysis from Digital Image Database Using Cross-Correlation and Wavelets

This document presents a framework for recognizing people by palm vein distribution analysis using cross-correlation based signatures to obtain descriptors. Haar wavelets are useful in reducing the number of features while maintaining high recognition rates. This experiment achieved 97.5% of individuals classified correctly with two levels of Haar wavelets. This study used twelve-version of RGB and NIR (near infrared) wavelength images per individual. One hundred people were studied;therefore 4,800 instances compose the complete database. A Multilayer Perceptron (MLP) was trained to improve the recognition rate in a k-fold cross-validation test with k = 10. Classification results using MLP neural network were obtained using Weka (open source machine learning software).

Cross-Correlation between Global Temperature and Atmospheric CO2 with a Temperature-Leading Time Lag

The temperature change and rate of CO2 change are correlated with a time lag, as reported in a previous paper. The correlation was investigated by calculating a correlation coefficient r of these changes for selected ENSO events in this study. Annual periodical increases and decreases in the CO2 concentration were considered, with a regular pattern of minimum values in August and maximum values in May each year. An increased deviation in CO2 and temperature was found in response to the occurrence of El Niño, but the increase in CO2 lagged behind the change in temperature by 5 months. This pattern was not observed for La Niña events. An increase in global CO2 emissions and a subsequent increase in global temperature proposed by IPCC were not observed, but an increase in global temperature, an increase in soil respiration, and a subsequent increase in global CO2 emissions were noticed. This natural process can be clearly detected during periods of increasing temperature specifically during El Niño events. The results cast strong doubts that anthropogenic CO2 is the cause of global warming.

The influence of sea water velocity variation on seismic traveltimes,ray paths,and amplitude

The main factors affecting seismic exploration is the propagation velocity of seismic waves in the medium. In the past, during marine seismic data processing, the propagation velocity of sea water was generally taken as a constant 1500 m/s. However, for deep water exploration, the sound velocity varies with the season, time, location, water depth, ocean currents, and etc.. It also results in a layered velocity distribution, so there is a difference of seismic traveltime, ray paths, and amplitude, which affect the migration imaging results if sea water propagation velocity is still taken as constant for the propagation wavefield. In this paper, we will start from an empirical equation of seismic wave velocity in seawater with changes of temperature, salinity, and depth, consider the variation of their values, build a seawater velocity model, and quantitatively analyze the impact of seawater velocity variation on seismic traveltime, ray paths, and amplitude in the seawater velocity model.

Joint 3D traveltime calculation based on fast marching method and wavefront construction

3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method （FMM） is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method （WFC）, which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models.

Finite-difference calculation of traveltimes based on rectangular grid

To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.

Temporal-spatial cross-correlation analysis of non-stationary near-surface wind speed time series

Temporal-spatial cross-correlation analysis of non-stationary wind speed time series plays a crucial role in wind field reconstruction as well as in wind pattern recognition.Firstly,the near-surface wind speed time series recorded at different locations are studied using the detrended fluctuation analysis(DFA),and the corresponding scaling exponents are larger than 1.This indicates that all these wind speed time series have non-stationary characteristics.Secondly,concerning this special feature( i.e.,non-stationarity)of wind signals,a cross-correlation analysis method,namely detrended cross-correlation analysis(DCCA) coefficient,is employed to evaluate the temporal-spatial cross-correlations between non-stationary time series of different anemometer pairs.Finally,experiments on ten wind speed data synchronously collected by the ten anemometers with equidistant arrangement illustrate that the method of DCCA cross-correlation coefficient can accurately analyze full-scale temporal-spatial cross-correlation between non-stationary time series and also can easily identify the seasonal component,while three traditional cross-correlation techniques(i.e.,Pearson coefficient,cross-correlation function,and DCCA method) cannot give us these information directly.

Finite-difference calculation of traveltimes based on rectangular grid

To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.

Reconstruction of sub cross-correlation cancellation technique for unambiguous acquisition of BOC(kn, n) signals

In order to solve the problem of ambiguous acquisition of BOC signals caused by its property of multiple peaks,an unambiguous acquisition algorithm named reconstruction of sub cross-correlation cancellation technique(RSCCT)for BOC(kn,n)signals is proposed.In this paper,the principle of signal decomposition is combined with the traditional acquisition algorithm structure,and then based on the method of reconstructing the correlation function.The method firstly gets the sub-pseudorandom noise(PRN)code by decomposing the local PRN code,then uses BOC(kn,n)and the sub-PRN code cross-correlation to get the sub cross-correlation function.Finally,the correlation peak with a single peak is obtained by reconstructing the sub cross-correlation function so that the ambiguities of BOC acquisition are removed.The simulation shows that RSCCT can completely eliminate the side peaks of BOC(kn,n)group signals while maintaining the narrow correlation of BOC,and its computational complexity is equivalent to sub carrier phase cancellation(SCPC)and autocorrelation side-peak cancellation technique(ASPeCT),and it reduces the computational complexity relative to BPSK-like.For BOC(n,n),the acquisition sensitivity of RSCCT is 3.25 dB,0.81 dB and 0.25 dB higher than binary phase shift keying(BPSK)-like,SCPC and ASPeCT at the acquisition probability of 90%,respectively.The peak to average power ratio is 1.91,3.0 and 3.7 times higher than ASPeCT,SCPC and BPSK-like at SNR=–20 dB,respectively.For BOC(2n,n),the acquisition sensitivity of RSCCT is 5.5 dB,1.25 dB and 2.69 dB higher than BPSK-like,SCPC and ASPeCT at the acquisition probability of 90%,respectively.The peak to average power ratio is 1.02,1.68 and 2.12 times higher than ASPeCT,SCPC and BPSK-like at SNR=–20 dB,respectively.

Estimation of Cloud Motion Using Cross-Correlation

This paper describes the estimation of cloud motion using lag cross-correlation. In order to compute the lag cross correlation, the Bayes Decision method is used first to identify cloud and surface of earth. Then cloud motion vectors are retrieved at a subset of points through multiple applications of a cross-correlation analysis. An objective analysis is used to define displacement at every satellite pixel throughout the domain and smooth the local inconsistencies. Cloud motions are then produced with a backward trajectory technique by using these displacement vectors.

Seismic traveltime inversion of 3D velocity model with triangulated interfaces

Seismic traveltime tomographic inversion has played an important role in detecting the internal structure of the solid earth. We use a set of blocks to approximate geologically complex media that cannot be well described by layered models or cells. The geological body is described as an aggregate of arbitrarily shaped blocks,which are separated by triangulated interfaces. We can describe the media as homogenous or heterogeneous in each block. We define the velocities at the given rectangle grid points for each block,and the heterogeneous velocities in each block can be calculated by a linear interpolation algorithm. The parameters of the velocity grid positions are independent of the model parameterization,which is advantageous in the joint inversion of the velocities and the node depths of an interface. We implement a segmentally iterative ray tracer to calculate traveltimes in the 3D heterogeneous block models.The damped least squares method is employed in seismic traveltime inversion,which includes the partial derivatives of traveltime with respect to the depths of nodes in the triangulated interfaces and velocities defined in rectangular grids. The numerical tests indicate that the node depths of a triangulated interface and homogeneous velocity distributions can be well inverted in a stratified model.

Detrended cross-correlation analysis of electroencephalogram

In the paper we use detrended cross-correlation analysis （DCCA） to study the electroencephalograms of healthy young subjects and healthy old subjects. It is found that the cross-correlation between different leads of a healthy young subject is larger than that of a healthy old subject. It was shown that the cross-correlation relationship decreases with the aging process and the phenomenon can help to diagnose whether the subject＇s brain function is healthy or not.

Reflection-based traveltime and waveform inversion with second-order optimization

Reflection-based inversion that aims to reconstruct the low-to-intermediate wavenumbers of the subsurface model, can be a complementary to refraction-data-driven full-waveform inversion(FWI), especially for the deep target area where diving waves cannot be acquired at the surface. Nevertheless, as a typical nonlinear inverse problem, reflection waveform inversion may easily suffer from the cycleskipping issue and have a slow convergence rate, if gradient-based first-order optimization methods are used. To improve the accuracy and convergence rate, we introduce the Hessian operator into reflection traveltime inversion(RTI) and reflection waveform inversion(RWI) in the framework of second-order optimization. A practical two-stage workflow is proposed to build the velocity model, in which Gauss-Newton RTI is first applied to mitigate the cycle-skipping problem and then Gauss-Newton RWI is employed to enhance the model resolution. To make the Gauss-Newton iterations more efficiently and robustly for large-scale applications, we introduce proper preconditioning for the Hessian matrix and design appropriate strategies to reduce the computational costs. The example of a real dataset from East China Sea demonstrates that the cascaded Hessian-based RTI and RWI have good potential to improve velocity model building and seismic imaging, especially for the deep targets.

RESEARCH ON THE ALGORITHM OF BINARY IMAGE CROSS-CORRELATION FOR UNSTEADY FLOW FIELD MEASUREMENT

The algorithm of Binary Image Cross-Correlation (BICC) was developed to measure the unsteady flow field. A vortex flow field was used to test the algorithm by numerical simulation. The results show that BICC is an effective algorithm for particle identification from consecutive images, the accurate velocity vector field can be obtained. The real velocity field in a valve chamber was measured by BICC in this study. From the full-field velocity information, the pressure and vorticity fields were also extracted by post-processing. (Edited author abstract) 6 Refs.

2-D crustal structure of Changbaishan-Tianchi volcanic region determined by seismic traveltime inversion

2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.

Model-data-driven P-wave impedance inversion using ResNets and the normalized zero-lag cross-correlation objective function

Model-driven and data-driven inversions are two prominent methods for obtaining P-wave impedance,which is significant in reservoir description and identification.Based on proper initial models,most model-driven methods primarily use the limited frequency bandwidth information of seismic data and can invert P-wave impedance with high accuracy,but not high resolution.Conventional data-driven methods mainly employ the information from well-log data and can provide high-accuracy and highresolution P-wave impedance owing to the superior nonlinear curve fitting capacity of neural networks.However,these methods require a significant number of training samples,which are frequently insufficient.To obtain P-wave impedance with both high accuracy and high resolution,we propose a model-data-driven inversion method using Res Nets and the normalized zero-lag cross-correlation objective function which is effective for avoiding local minima and suppressing random noise.By using initial models and training samples,the proposed model-data-driven method can invert P-wave impedance with satisfactory accuracy and resolution.Tests on synthetic and field data demonstrate the proposed method’s efficacy and practicability.

Joint tomographic inversion of first-arrival and reflection traveltimes for recovering 2-D seismic velocity structure with an irregular free surface

Irregular surface flattening,which is based on a boundary conforming grid and the transformation between curvilinear and Cartesian coordinate systems,is a mathematical method that can elegantly handle irregular surfaces,but has been limited to obtaining first arrivals only.By combining a multistage scheme with the fast-sweeping method(FSM,the method to obtain first-arrival traveltime in curvilinear coordinates),the reflected waves from a crustal interface can be traced in a topographic model,in which the reflected wavefront is obtained by reinitializing traveltimes in the interface for upwind branches.A local triangulation is applied to make a connection between velocity and interface nodes.Then a joint inversion of first-arrival and reflection traveltimes for imaging seismic velocity structures in complex terrains is presented.Numerical examples all perform well with different seismic velocity models.The increasing topographic complexity and even use of a high curvature reflector in these models demonstrate the reliability,accuracy and robustness of the new working scheme;checkerboard testing illustrates the method's high resolution.Noise tolerance testing indicates the method's ability to yield practical traveltime tomography.Further development of the multistage scheme will allow other later arrivals to be traced and used in the traveltime inversion.

Public cloud computing for seismological research:Calculating large-scale noise cross-correlations using ALIYUN

The amount of seismological data is rapidly increasing with accumulating observational time and increasing number of stations, requiring modern technique to provide adequate computing power. In present study, we proposed a framework to calculate large-scale noise crosscorrelation functions(NCFs) using public cloud service from ALIYUN. The entire computation is factorized into small pieces which are performed parallelly on specified number of virtual servers provided by the cloud. Using data from most seismic stations in China, five NCF databases are built. The results show that, comparing to the time cost using a single server, the entire time can be reduced over two orders of magnitude depending number of evoked virtual servers. This could reduce computation time from months to less than 12 hours. Based on obtained massive NCFs, the global body waves are retrieved through array interferometry and agree well with those from earthquakes. This leads to a solution to process massive seismic dataset within an affordable time and is applicable to other large-scale computing in seismological researches.

DETECTING TECHNIQUE OF WEAK PERIODIC PULSE SIGNAL VIA SYNTHESIS OF CROSS-CORRELATION AND CHAOTIC SYSTEM

In this letter, with the synthesis of usual cross-correlation detecting method andchaotic detecting method, a new detecting system for the weak periodic pulse signal is constituted,in which the two methods can play respective preponderance. Theoretical analyses and simulationstudies have shown that the detecting system is very sensitive to the periodic pulse signal understrong noise background and has exceedingly powerful capability of suppressing complex noise.