共查询到20条相似文献，搜索用时 15 毫秒
1.
A technique for reconstruction of the 2d surface velocity field from radar observations is proposed. The method consecutively employs two processing techniques: At the first stage raw radial velocity data are subject to EOF analysis, which enables to fill gaps in observations and provides estimates of the noise level and integral parameters characterizing small-scale variability of the sea surface circulation. These parameters are utilized at the second stage, when the cost function for variational interpolation is constructed, and the updated radial velocities are interpolated on the regular grid.Experiments with simulated and real data are used to assess the method's skill and compare it with the conventional 2d variational (2dVar) approach. It is shown that the proposed technique consistently improves performance of the 2dVar algorithm and becomes particularly effective when a radar stops operating for 1–2 days and/or a persistent gap emerges in spatial coverage of a basin by the HFR network. 相似文献
2.
Data interpolation is an important step for seismic data analysis because many processing tasks, such as multiple attenuation and migration, are based on regularly sampled seismic data. Failed interpolations may introduce artifacts and eventually lead to inaccurate final processing results. In this paper, we generalised seismic data interpolation as a basis pursuit problem and proposed an iteration framework for recovering missing data. The method is based on non‐linear iteration and sparse transform. A modified Bregman iteration is used for solving the constrained minimisation problem based on compressed sensing. The new iterative strategy guarantees fast convergence by using a fixed threshold value. We also propose a generalised velocity‐dependent formulation of the seislet transform as an effective sparse transform, in which the non‐hyperbolic normal moveout equation serves as a bridge between local slope patterns and moveout parametres in the common‐midpoint domain. It can also be reduced to the traditional velocity‐dependent seislet if special heterogeneity parametre is selected. The generalised velocity‐dependent seislet transform predicts prestack reflection data in offset coordinates, which provides a high compression of reflection events. The method was applied to synthetic and field data examples, and the results show that the generalised velocity‐dependent seislet transform can reconstruct missing data with the help of the modified Bregman iteration even for non‐hyperbolic reflections under complex conditions, such as vertical transverse isotropic (VTI) media or aliasing. 相似文献
3.
A short-term predictive system for surface currents from a rapidly deployed coastal HF radar network
Donald Barrick Vicente Fernandez Maria I. Ferrer Chad Whelan Øyvind Breivik 《Ocean Dynamics》2012,62(5):725-740
In order to address the need for surface trajectory forecasts following deployment of coastal HF radar systems during emergency-response situations (e.g., search and rescue, oil spill), a short-term predictive system (STPS) based on only a few hours data background is presented. First, open-modal analysis (OMA) coefficients are fitted to 1-D surface currents from all available radar stations at each time interval. OMA has the effect of applying a spatial low-pass filter to the data, fills gaps, and can extend coverage to areas where radial vectors are available from a single radar only. Then, a set of temporal modes is fitted to the time series of OMA coefficients, typically over a short 12-h trailing period. These modes include tidal and inertial harmonics, as well as constant and linear trends. This temporal model is the STPS basis for producing up to a 12-h current vector forecast from which a trajectory forecast can be derived. We show results of this method applied to data gathered during the September 2010 rapid-response demonstration in northern Norway. Forecasted coefficients, currents, and trajectories are compared with the same measured quantities, and statistics of skill are assessed employing 16 24-h data sets. Forecasted and measured kinetic variances of the OMA coefficients typically agreed to within 10–15%. In one case where errors were larger, strong wind changes are suspected and examined as the cause. Sudden wind variability is not included properly within the STPS attack we presently employ and will be a subject for future improvement. 相似文献
4.
描述地震波衰减特征的品质因子Q对地震数据处理和油藏描述非常重要,在地震勘探领域,Q值一般通过垂直地震剖面(VSP)数据或地面地震数据得到.由于叠前地面地震数据具有复杂的射线路径且存在噪声、调谐干涉效应等影响,从叠前地震数据中准确估计Q值相对困难.本文以地震波射线传播为基础,根据同相轴局部斜率和射线参数的映射关系,将多射线波形频谱同时带入谱比法联合反演估计Q值,提出了基于多射线联合反演的速度无关叠前Q值估计方法.该方法通过局部斜率属性避开了速度对Q值估计的影响,局部斜率携带地震波传播的速度信息,具有相同局部斜率的地震反射波具有相同的传播射线参数.同相轴局部斜率是地震数据域的属性,而速度是模型域的参数,在估计Q值中采用数据域的属性参数可以直接应用于数据的联合反演,而不需要通过速度对其做进一步的转化,从而提高了Q值估计的精度.同时,本方法采用预测映射(predictive mapping)技术将非零炮检距反射信息映射到零炮检距处,从而获得零偏移距走时对应的Q值.模拟和实际算例验证了本文方法的有效性. 相似文献
5.
Near‐surface problem is a common challenge faced by land seismic data processing, where often, due to near‐surface anomalies, events of interest are obscured. One method to handle this challenge is near‐surface layer replacement, which is a wavefield reconstruction process based on downward wavefield extrapolation with the near‐surface velocity model and upward wavefield extrapolation with a replacement velocity model. This requires, in theory, that the original wavefield should be densely sampled. In reality, data acquisition is always sparse due to economic reasons, and as a result in the near‐surface layer replacement data interpolation should be resorted to. For datasets with near‐surface challenges, because of the complex event behaviour, a suitable interpolation scheme by itself is a challenging problem, and this, in turn, makes it difficult to carry out the near‐surface layer replacement. In this research note, we first point out that the final objective of the near‐surface layer replacement is not to obtain a newly reconstructed wavefield but to obtain a better final image. Next, based upon this finding, we propose a new thinking, interpolation‐free near‐surface layer replacement, which can handle complex datasets without any interpolation. Data volume expansion is the key idea, and with its help, the interpolation‐free near‐surface layer replacement is capable of preserving the valuable information of areas of interest in the original dataset. Two datasets, i.e., a two‐dimensional synthetic dataset and a three‐dimensional field dataset, are used to demonstrate this idea. One conclusion that can be drawn is that an attempt to interpolate data before layer replacement may deteriorate the final image after layer replacement, whereas interpolation‐free near‐surface layer replacement preserves all image details in the subsurface. 相似文献
6.
Mapping the mesoscale surface velocity stream function by combining estimates of surface height from satellite altimetry
and surface currents from sequential infrared (sea-surface temperature) imagery using optimal interpolation is described.
Surface currents are computed from infrared images by the method of maximum cross-correlations (MCC) and are combined with
altimeter sea-level anomaly data from the TOPEX/Poseidon and ERS satellites. The analysis method was applied to 6 years of
data from the East Australian Current region. The covariance of velocity and sea-level data is consistent with the statistical
assumptions of homogeneous, isotropic turbulence, with typical length scales of order 220 km and time scales of 10 days in
this region. Augmenting the analysis of altimeter data with MCC velocity observations improves the resolution of the surface
currents, especially near the Australian coast, and demonstrates that the two data sources provide consistent and complementary
observations of the surface mesoscale circulation. The volume of MCC data is comparable to that from a satellite altimeter,
but with a more variable distribution of spatial and temporal resolution. In concert with altimetry, satellite radiometer
velocimetry represents a technique useful for retrospective analysis of currents from high-resolution satellite radiometer
data-sets.
Received: 3 July 2001 / Accepted: 16 November 2001 相似文献
7.
Local seismic event slopes contain subsurface velocity information and can be used to estimate seismic stacking velocity. In this paper, we propose a novel approach to estimate the stacking velocity automatically from seismic reflection data using similarity‐weighted k‐means clustering, in which the weights are local similarity between each trace in common midpoint gather and a reference trace. Local similarity reflects the local signal‐to‐noise ratio in common midpoint gather. We select the data points with high signal‐to‐noise ratio to be used in the velocity estimation with large weights in mapped traveltime and velocity domain by similarity‐weighted k‐means clustering with thresholding. By using weighted k‐means clustering, we make clustering centroids closer to those data points with large weights, which are more reliable and have higher signal‐to‐noise ratio. The interpolation is used to obtain the whole velocity volume after we have got velocity points calculated by weighted k‐means clustering. Using the proposed method, one obtains a more accurate estimate of the stacking velocity because the similarity‐based weighting in clustering takes into account the signal‐to‐noise ratio and reliability of different data points in mapped traveltime and velocity domain. In order to demonstrate that, we apply the proposed method to synthetic and field data examples, and the resulting images are of higher quality when compared with the ones obtained using existing methods. 相似文献
8.
9.
This work introduces a new variational Bayes data assimilation method for the stochastic estimation of precipitation dynamics using radar observations for short term probabilistic forecasting (nowcasting). A previously developed spatial rainfall model based on the decomposition of the observed precipitation field using a basis function expansion captures the precipitation intensity from radar images as a set of ‘rain cells’. The prior distributions for the basis function parameters are carefully chosen to have a conjugate structure for the precipitation field model to allow a novel variational Bayes method to be applied to estimate the posterior distributions in closed form, based on solving an optimisation problem, in a spirit similar to 3D VAR analysis, but seeking approximations to the posterior distribution rather than simply the most probable state. A hierarchical Kalman filter is used to estimate the advection field based on the assimilated precipitation fields at two times. The model is applied to tracking precipitation dynamics in a realistic setting, using UK Met Office radar data from both a summer convective event and a winter frontal event. The performance of the model is assessed both traditionally and using probabilistic measures of fit based on ROC curves. The model is shown to provide very good assimilation characteristics, and promising forecast skill. Improvements to the forecasting scheme are discussed. 相似文献
10.
S. Abhilash A. K. Sahai K. Mohankumar John P. George Someshwar Das 《Pure and Applied Geophysics》2012,169(11):2047-2070
In this paper the impact of Doppler weather radar (DWR) reflectivity and radial velocity observations for the short range forecasting of a tropical storm and associated rainfall event have been examined. Doppler radar observations of a tropical storm case that occurred during 29–30 October 2006 from SHARDWR (13.6° N, 80.2° E) are assimilated in the WRF 3DVAR system. The observation operator for radar reflectivity and radial velocity is included within latest version of WRF 3DVAR system. Keeping all model physics the same, three experiments were conducted at a horizontal resolution of 30?km. In the control experiment (CTRL), NCEP Final Analysis (FNL) interpolated to the model grid was used as the initial condition for 48-h free forecast. In the second experiment (NODWR), 6-h assimilation cycles have been carried out using all conventional (radiosonde and surface data) and non-conventional (satellite) observations from the Global Telecommunication System (GTS). The third experiment (DWR) is the same as the second, except Doppler radar radial velocity and reflectivity observations are also used in the assimilation cycle. Continuous 6-h assimilation cycle employed in the WRF-3DVAR system shows positive impact on the rainfall forecast. Assimilation of DWR data creates several small scale features near the storm centre. Additional sensitivity experiments were conducted to study the individual impact of reflectivity and radial velocity in the assimilation cycle. Radar data assimilation with reflectivity alone produced large analysis response on both thermodynamical and dynamical fields. However, radial velocity assimilation impacted only on dynamical fields. Analysis increments with radar reflectivity and radial velocity produce adjustments in both dynamical and thermodynamical fields. Verification of QPF skill shows that radar data assimilation has a considerable impact on the short range precipitation forecast. Improvement of the QPF skill with radar data assimilation is more clearly seen in the heavy rainfall (for thresholds >7?mm) event than light rainfall (for thresholds of 1 and 3?mm). The spatial pattern of rainfall is well simulated by the DWR experiment and is comparable to TRMM observations. 相似文献
11.
J. A. CORBYN 《Geophysical Prospecting》1971,19(4):698-717
The application of automatic methods to the processing of potential field survey results has led to the need for the development of special procedures for the interpolation and extrapolation of such data. A method is proposed which utilizes the mathematical properties of potential field data to achieve a minimum distortion solution to the interpolation/extrapolation problem. The proposed method is suitable for use with any data point distribution. Digital computer programs have been prepared using the method; results obtained through the use of the programs are presented. 相似文献
12.
《中国科学:地球科学(英文版)》2017,(3)
The Proper Orthogonal Decomposition(POD)-based ensemble four-dimensional variational(4DVar) assimilation method(POD4DEnVar) was proposed to combine the strengths of EnKF(i.e.,the ensemble Kalman filter) and 4DVar assimilation methods.Recently,a POD4DEnVar-based radar data assimilation scheme(PRAS) was built and its effectiveness was demonstrated.POD4 DEnVar is based on the assumption of a linear relationship between the model perturbations(MPs)and the observation perturbations(OPs);however,this assumption is likely to be destroyed by the highly non-linear forecast model or observation operator.To address this issue,using the Gauss-Newton iterative method,the nonlinear least squares enhanced POD4 DEnVar algorithm(referred to as NLS-4DVar) was proposed.Naturally,the PRAS was upgraded to form the NLS-4DVar-based radar data assimilation scheme(NRAS).To evaluate the performance of NRAS against PRAS,observing system simulation experiments(OSSEs) were conducted to assimilate reflectivity and radial velocity individually,with one,two,and three iterations.The results demonstrated that the NRAS outperformed PRAS in improving the initial condition and the forecasting of model variables and rainfall.The NRAS,with a smaller number of iterations,can yield a convergent result.In contrast to the situation when assimilating radial velocity,the advantages of NRAS over PRAS were more obvious when assimilating reflectivity. 相似文献
13.
Two very high-frequency radars (VHFR) operating on the Opal coast of eastern English Channel provided a nearly continuous 35-day long dataset of surface currents over a 500 km2 area at 0.6–1.8 km resolution. Argo drifter tracking and CTD soundings complemented the VHFR observations, which extended approximately 25 km offshore. The radar data resolve three basic modes of the surface velocity variation in the area, that are driven by tides, winds and freshwater fluxes associated with seasonal river discharge. The first mode, accounting for 90% of variability, is characterized by an along-shore flow pattern, whereas the second and third modes exhibit cross-shore, and eddy-like structures in the current velocity field. All the three modes show the dominant semi-diurnal variability and low-frequency modulation by the neap-spring tidal cycle. Although tidal forcing provides the major contribution to variability of local currents, baroclinicity plays an important role in shaping the 3D velocity field averaged over the tidal cycle and may strongly affect tracer dynamics on larger time scales. An empirical orthogonal function (EOF) decomposition and a spectral rotary analysis of the VHFR data reveal a discontinuity in the velocity field occurring approximately 10 km offshore which was caused by the reversal in the sign of rotation of the current vector. This feature of local circulation is responsible for surface current convergence on ebb, divergence on flood and strong oscillatory vertical motion. Spectral analysis of the observed currents and the results of the Agro drifter tracking indicate that the line of convergence approximately follows the 30-m isobath. The most pronounced feature of the radar-derived residual circulation is the along-coast intensification of surface currents with velocity magnitude of 0.25 m/s typical for the Regions of Freshwater Influence (ROFI). The analysis has provided a useful, exploratory examination of surface currents, suggesting that the circulation off the Opal coast is governed by ROFI dynamics on the hypertidal background. 相似文献
14.
Tomographic inversions for velocity variations in western Washington indicate a high correlation with surface geology and geophysical measurements, including gravity observations. By assuming a simple linear relationship between density and velocity (Birch's law) it is possible to calculate the gravity field predicted from the velocity perturbations obtained by local tomographic inversion. While the predicted gravity matches observations in parts of the model, the overall correlation is not satisfactory. In this paper we suggest a method of constraining the tomographic inversion to fit the gravity observations simultaneously with the seismic travel time data. The method is shown to work well with synthetic data in 3 dimensions where the assumption of Birch's law holds strictly. If the sources of the gravity anomalies are assumed to be spatially localized, integration can be carried out over a relatively small volume below the observation points and sparse matrix techniques can be applied. We have applied the constrained inversion method to western Washington using 4,387 shallow earthquakes, to depths of 40.0 km, (36,865 raypaths) convering a 150×250 km region and found that the gravitational constraints may be satisfied with minor effect on the degree of misfit to the seismic data. 相似文献
15.
人工地震方法由于受到野外观测系统和经济因素等的限制,采集的数据在空间方向总是不规则分布.但是,许多地震数据处理技术的应用（如：多次波衰减,偏移和时移地震）都基于空间规则分布条件下的地震数据体.因此,数据插值技术是地震数据处理流程中关键环节之一.失败的插值方法往往会引入虚假信息,给后续处理环节带来严重的影响.迭代插值方法是目前广泛应用的地震数据重建思路,但是常规的迭代插值方法往往很难保证插值精度,并且迭代收敛速度较慢,尤其存在随机噪声的情况下,插值地震道与原始地震道之间存在较大的信噪比差异.因此开发快速的、有效的迭代数据插值方法具有重要的工业价值.本文将地震数据插值归纳为数学基追踪问题,在压缩感知理论框架下,提出新的非线性Bregman整形迭代算法来求解约束最小化问题,同时在迭代过程中提出两种匹配的迭代控制准则,通过有效的稀疏变换对缺失数据进行重建.通过理论模型和实际数据测试本文方法,并且与常规迭代插值算法进行比较,结果表明Bregman整形迭代插值方法能够更加有效地恢复含有随机噪声的缺失地震信息. 相似文献
16.
本文利用检层法、表面波法实测的波速资料。通过对实测数据的处理和回归分析,求出了本溪市区不同岩土波速与埋深的回归关系。同时收集了市区已有的钻孔资料,经分析计算,给出了本溪市区平均剪切波速的分布。 相似文献
17.
本文介绍了利用地脉动中的面波信息探测地壳浅部横波速度结构的观测方法、仪器设备、野外数据采集和用空间自相关方法测定相速度及相关的数据处理方法,并给出了在两个场地的探测实例. 通过与钻探结果对比表明,由该方法推测的地壳浅部1～2 km 深度范围的浅层横波速度结构,对于确定物性差异大的地层界面是有效和可靠的. 由于该方法使用的是地脉动噪声中的面波信息,不需要专门的震源, 因此具有经济快捷、应用范围广泛等优点,将可能成为一种新的具有很好应用前景的地球物理勘探方法. 相似文献
18.
三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先,在常规双线性插值公式推导过程中,充分利用平面波双线性假设的结论,获得了二元极小值超越方程的解析解,进而推导出了准确的局部走时计算公式,同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次,引入迎风差分的基本思想,提出迎风双线性插值的局部走时计算策略,该计算策略能简化算法、提高效率且保证无条件稳定性;然后,将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合,获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后,通过精度和效率分析检验了新算法的精度、效率和正确性,并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性. 相似文献
19.
许多地震资料处理方法需要完整的数据信息,但是受野外施工条件等因素的影响,观测系统很难记录完整的地震波场,如空间采样率不足和地震道缺失等现象,尤其是缺失的叠前地震数据时常产生空间假频现象,给后续处理流程中很多重要环节带来严重的影响.传统数据插值方法通常很难同时解决数据缺失和空间假频问题,因此开发有效的反空间假频数据插值方法具有重要的意义.本文通过同时改变时间和空间方向采样比例,利用预测误差滤波器的尺度缩放不变性,计算反空间假频地震倾角模式,构建可有效压缩含空间假频不完整地震数据的反假频seislet变换方法,通过压缩感知Bregman迭代算法,对缺失地震数据进行反假频插值.理论模型和实际数据的处理结果验证了基于seislet变换的迭代插值方法可以有效地恢复含有假频的缺失地震信息. 相似文献
20.
Three choices of control variables for meteorological variational analysis (3DVAR or 4DVAR) are associated with horizontal
wind: (1) streamfunction and velocity potential, (2) eastward and northward velocity, and (3) vorticity and divergence. This
study shows theoretical and numerical differences of these variables in practical 3DVAR data assimilation through statistical
analysis and numerical experiments. This paper demonstrates that (a) streamfunction and velocity potential could potentially
introduce analysis errors; (b) A 3DVAR using velocity or vorticity and divergence provides a natural scale dependent influence
radius in addition to the covariance; (c) for a regional analysis, streamfunction and velocity potential are retrieved from
the background velocity field with Neumann boundary condition. Improper boundary conditions could result in further analysis
errors; (d) a variational data assimilation or an inverse problem using derivatives as control variables yields smoother analyses,
for example, a 3DVAR using vorticity and divergence as controls yields smoother wind analyses than those analyses obtained
by a 3DVAR using either velocity or streamfunction/velocity potential as control variables; and (e) statistical errors of
higher order derivatives of variables are more independent, e.g., the statistical correlation between U and V is smaller than the one between streamfunction and velocity potential, and thus the variables in higher derivatives are more
appropriate for a variational system when a cross-correlation between variables is neglected for efficiency or other reasons.
In summary, eastward and northward velocity, or vorticity and divergence are preferable control variables for variational
systems and the former is more attractive because of its numerical efficiency. Numerical experiments are presented using analytic
functions and real atmospheric observations. 相似文献