共查询到20条相似文献，搜索用时 562 毫秒
1.
Xing Fang 《Studia Geophysica et Geodaetica》2014,58(1):1-16
A new proof is presented of the desirable property of the weighted total least-squares (WTLS) approach in preserving the structure of the coefficient matrix in terms of the functional independent elements. The WTLS considers the full covariance matrix of observed quantities in the observation vector and in the coefficient matrix; possible correlation between entries in the observation vector and the coefficient matrix are also considered. The WTLS approach is then equipped with constraints in order to produce the constrained structured TLS (CSTLS) solution. The proposed approach considers the correlation between the observation vector and the coefficient matrix of an Error-In-Variables model, which is not considered in other, recently proposed approaches. A rigid transformation problem is done by preservation of the structure and satisfying the constraints simultaneously. 相似文献
2.
地震反演是储层定量描述和地震油气识别的关键技术,反演结果在复杂构造区域的横向连续性和保真性是影响地震资料定量解释精度的重要因素.基于此,本文发展了地震数据互相关驱动的多道反演方法.考虑地层反射系数与地震数据在结构上具有相似性的特点,基于地震数据互相关描述地层反射系数的结构特征,并将其作为多道地震反演的横向约束条件;此外,为改善地震数据本身横向连续性差对反演结果的影响,在目标泛函的惩罚项中引入局部优化算子,构建了一个易于求解的多道地震反演目标泛函.与常规多道地震反演方法相比,本文方法能够设计更合理、更符合实际情况的横向约束算子,提高反演结果的横向连续性,并且能有效降低地震资料质量对反演结果的影响.模型测试和实际应用验证了本方法的可靠性和稳定性. 相似文献
3.
4.
《Advances in water resources》2007,30(4):954-965
This paper exploits the unique feature of the Ant Colony Optimization Algorithm (ACOA), namely incremental solution building mechanism, to develop partially constraint ACO algorithms for the solution of optimization problems with explicit constraints. The method is based on the provision of a tabu list for each ant at each decision point of the problem so that some constraints of the problem are satisfied. The application of the method to the problem of storm water network design is formulated and presented. The network nodes are considered as the decision points and the nodal elevations of the network are used as the decision variables of the optimization problem. Two partially constrained ACO algorithms are formulated and applied to a benchmark example of storm water network design and the results are compared with those of the original unconstrained algorithm and existing methods. In the first algorithm the positive slope constraints are satisfied explicitly and the rest are satisfied by using the penalty method while in the second one the satisfaction of constraints regarding the maximum ratio of flow depth to the diameter are also achieved explicitly via the tabu list. The method is shown to be very effective and efficient in locating the optimal solutions and in terms of the convergence characteristics of the resulting ACO algorithms. The proposed algorithms are also shown to be relatively insensitive to the initial colony used compared to the original algorithm. Furthermore, the method proves itself capable of finding an optimal or near-optimal solution, independent of the discretisation level and the size of the colony used. 相似文献
5.
地震勘探方法在深部固体矿产资源勘探中发展潜力巨大,同时也面临挑战.由于固体矿产资源地下分布呈现陡峭构造、尺度小,物性差异小的特点,常规偏移方法对小尺度矿体成像的分辨率提高有限.本文研究了一种基于稀疏促进约束的最小二乘逆时偏移方法.首先,将非均匀分布的矿体等效为随机介质,建立小尺度扰动的矿体模型;其次,改进现有最小二乘偏移方法,以稀疏模型为先验信息约束成像结果,并通过Curvelet变换压缩成像空间,经过多次迭代计算,可以提高小尺度散射体的成像分辨率;再次,对炮域记录进行随机震源编码,减少成像所需的炮集个数,通过稀疏促进约束条件,降低串扰噪声引起的成像误差.通过庐枞金属矿模型数值计算,验证本文方法可以较好的成像包含小尺度散射体的金属矿地质模型. 相似文献
6.
为了合理计算山区桥梁支座刚度,针对桥墩高度不相同的特点,考虑上部结构对桥墩顶部的转动约束作用,提出在横桥向可将墩顶视为自由约束,而在纵桥向将墩顶视为定向约束。分别按照地震作用下各墩底剪力和弯矩相等的原则,推导桥梁支座纵、横桥向的刚度设计公式,并给出各桥墩支座的设计方法。为验证方法的正确性,以墩底剪力相等的原则为例,利用OpenSees建立一座墩高不等的5跨连续梁桥模型,并依支座刚度取值不同分三种工况:工况一各桥墩支座刚度相同;工况二按墩顶自由计算各支座的纵、横桥向刚度;工况三按墩顶定向约束计算各支座的纵、横桥向刚度。分别对三种工况下的桥梁结构输入三条地震动记录进行时程分析,考察各桥墩的底部剪力。分析结果表明:工况一各桥墩纵、横桥向的底部剪力均不相同;工况二各桥墩横桥向的底部剪力相同而纵桥向的底部剪力不同;工况三各桥墩纵桥向的底部剪力相同而横桥向的底部剪力不同。上述结果表明在桥梁支座设计时,横桥向桥墩的抗推刚度应按墩顶自由计算,而纵桥向桥墩的抗推刚度应按墩顶为定向约束计算。 相似文献
7.
J. J. Gu P. Guo G. H. Huang 《Stochastic Environmental Research and Risk Assessment (SERRA)》2013,27(5):1207-1219
In this study, an interval parameter multistage joint-probability programming (IMJP) approach has been developed to deal with water resources allocation under uncertainty. The IMJP can be used not only to deal with uncertainties in terms of joint-probability and intervals, but also to examine the risk of violating joint probabilistic constraints in the context of multistage. The proposed model can handle the economic expenditure caused by regional water shortage and flood control. The model can also reflect the related dynamic changes in the multi-stage cases and the system safety under uncertainty. The developed method is applied to a case study of water resources allocation in Shandong, China, under multistage, multi-reservoir and multi-industry. The violating reservoir constraints are addressed in terms of joint-probability. Different risk levels of constraint lead to different planning. The obtained results can help water resources managers to identify desired system designs under various economic, environment and system reliability scenarios. 相似文献
8.
全波形反演是一种高精度的反演方法,其目标函数是一个强非线性函数,易受局部极值影响,而且反演过程计算量较大.波场重构反演是近几年提出的一种改进的全波形反演理论.该反演方法通过将波动方程作为惩罚项引入到目标函数中,通过拓宽解的寻找空间减弱了局部极小值的影响,而且反演过程不需要计算伴随波场,提高了计算效率.但该反演方法一直缺少准确的惩罚因子算法,直接影响到该方法的准确度.本文将波场重构反演拓展到时间域并利用梯度法进行波场重构.频率域的惩罚因子用来加强波动方程的约束,而时间域惩罚因子表现为调节模拟波场和实际波场的权重因子.为此,我们根据约束优化理论,在波动方程准确以及重构波场与反演参数解耦的假设下,提出以波动方程为目标函数的新的惩罚因子算法.根据波形反演在应用时普遍存在的噪音干扰、子波错误和低频信息缺失的情况下,应用部分Sigsbee2A模型合成数据对本文提出的算法进行实验.数值实验结果表明：基于新的惩罚因子算法,在其他信息不准确的情况下,波场重构反演可以给出高精度的反演结果. 相似文献
9.
Typical pump-and-treat (PAT) optimization problems involve design of pumping schemes, while minimizing cost and meeting a set of constraints. Due to scarcity of information about the hydrogeological system, stochastic modeling approaches can be used to assess tradeoffs between optimality and reliability. Using a stochastic approach, the constrained, single-objective problem may be turned into a multiobjective problem by substituting constraint inequalities with an additional objective function (OF) that accounts for the reliability of the PAT process. In this work, two approaches are analyzed: in one case, the additional OF consists of the probability of failure of a given remediation policy; in another, the OF additional is represented by the recourse, namely the penalty cost induced by the violation of constraints. In order to overcome the overwhelming computational cost required by stochastic simulation, surrogate forms of the OFs are introduced. In the test case under investigation, such functions are estimated by a kriging interpolation of the OF over a series of data points obtained from stochastic simulations of flow and transport, and calibrated against stochastic optimization solutions. The analysis of the two approaches for addressing the tradeoff of cost vs. reliability indicates that recourse accounts not only for the frequency of constraint violations, as the probability of failure does, but also for the intensity with which these occur. Ultimately, the recourse method allows considering less restrictive policies, although these may be highly sensitive to the choice of penalty functions. 相似文献
10.
重力卫星可以在相同误差尺度下对全球质量变化进行连续重复观测,并在近十余年来取得了巨大成功,探索重力卫星数据精化处理方法和相关应用研究具有重要意义.本文基于三维加速度点质量模型法的基本原理,进一步发展建立了时变重力场模型球谐位系数的变化和地面点质量变化的关系,可有效考虑地表质量变化导致的负荷形变的影响;引入等权形式、线性形式、指数形式和高斯形式的空间约束方法处理南北条带噪声和向下延拓导致的病态问题,并与零阶Tikhonov正则化方法进行对比分析.采用模拟数据和一个月的实测GRACE时变重力场模型计算全球质量变化,对三维加速度点质量模型法和几种空间约束方法进行对比分析验证.计算结果表明,对于3°等面积的全球格网质量点,高斯和指数形式空间约束方法的最优相关距离约为500km,等权和线性形式空间约束方法的最优相关距离约为600km,各方法均可有效处理条带噪声的影响,四种空间约束方法的计算效果优于零阶Tikhonov正则化方法,本文的相关方法为进一步利用三维加速度点质量模型法监测全球质量变化提供了借鉴. 相似文献
11.
To date, an outstanding issue in hydrologic data assimilation is a proper way of dealing with forecast bias. A frequently used method to bypass this problem is to rescale the observations to the model climatology. While this approach improves the variability in the modeled soil wetness and discharge, it is not designed to correct the results for any bias. Alternatively, attempts have been made towards incorporating dynamic bias estimates into the assimilation algorithm. Persistent bias models are most often used to propagate the bias estimate, where the a priori forecast bias error covariance is calculated as a constant fraction of the unbiased a priori state error covariance. The latter approach is a simplification to the explicit propagation of the bias error covariance. The objective of this paper is to examine to which extent the choice for the propagation of the bias estimate and its error covariance influence the filter performance. An Observation System Simulation Experiment (OSSE) has been performed, in which ground water storage observations are assimilated into a biased conceptual hydrologic model. The magnitudes of the forecast bias and state error covariances are calibrated by optimizing the innovation statistics of groundwater storage. The obtained bias propagation models are found to be identical to persistent bias models. After calibration, both approaches for the estimation of the forecast bias error covariance lead to similar results, with a realistic attribution of error variances to the bias and state estimate, and significant reductions of the bias in both the estimates of groundwater storage and discharge. Overall, the results in this paper justify the use of the traditional approach for online bias estimation with a persistent bias model and a simplified forecast bias error covariance estimation. 相似文献
12.
Energy balance closure at ChinaFLUX sites 总被引：1，自引：0，他引：1
Network of eddy covariance observation is measuring long-term carbon and water fluxes in contrasting ecosystems and climates. As one important reference of independently evaluating scalar flux estimates from eddy covariance, energy balance closure is used widely in study of carbon and water fluxes. Energy balance closure in ChinaFLUX was evaluated by statistical regression of turbulent energy fluxes (sensible and latent heat) against available energy (net radiation, soil heat flux, canopy heat storage) and the energy balance ratio (EBR) and the frequency distribution of relative errors of energy balance (δ). The trends of diurnal and seasonal variation of energy balance in ChinaFLUX were analyzed. The results indicated that the imbalance was prevalent in all observation sites, but there were little differences among sites because of the properties variation of sites. The imbalance was greater during nocturnal periods than daytime and closure was improved with friction velocity intensifying. Generally the results suggested that estimates of the scalar turbulent fluxes of sensible and latent heat were underestimated and/or that available energy was overestimated. Finally, we discussed certain factors that are contributed to the imbalance of energy, such as systematic errors associated with the sampling mismatch, systematic instrument bias, neglected energy sinks, low and high frequency loss of turbulent fluxes and advection of heat and water vapor. 相似文献
13.
This study has applied evolutionary algorithm to address the data assimilation problem in a distributed hydrological model. The evolutionary data assimilation (EDA) method uses multi-objective evolutionary strategy to continuously evolve ensemble of model states and parameter sets where it adaptively determines the model error and the penalty function for different assimilation time steps. The assimilation was determined by applying the penalty function to merge background information (i.e., model forecast) with perturbed observation data. The assimilation was based on updated estimates of the model state and its parameterizations, and was complemented by a continuous evolution of competitive solutions.The EDA was illustrated in an integrated assimilation approach to estimate model state using soil moisture, which in turn was incorporated into the soil and water assessment tool (SWAT) to assimilate streamflow. Soil moisture was independently assimilated to allow estimation of its model error, where the estimated model state was integrated into SWAT to determine background streamflow information before they are merged with perturbed observation data. Application of the EDA in Spencer Creek watershed in southern Ontario, Canada generates a time series of soil moisture and streamflow. Evaluation of soil moisture and streamflow assimilation results demonstrates the capability of the EDA to simultaneously estimate model state and parameterizations for real-time forecasting operations. The results show improvement in both streamflow and soil moisture estimates when compared to open-loop simulation, and a close matching between the background and the assimilation illustrates the forecasting performance of the EDA approach. 相似文献
14.
J. Sénégas H. Wackernagel W. Rosenthal T. Wolf 《Stochastic Environmental Research and Risk Assessment (SERRA)》2001,15(1):65-86
The efficiency of a sequential data assimilation scheme relies on the capability to describe the error covariance. This aspect
is all the more relevant if one needs accurate statistics on the estimation error. Frequently an ad hoc function depending
on a few parameters is proposed, and these parameters are tuned, estimated or updated. This usually requires that the covariance
is second-order stationary (i.e. depends only on the distance between two points). In this paper, we discuss this feature
and show that even in simple applications (such as one-dimensional hydrodynamics), this assumption does not hold and may lead
to poorly described estimation errors. We propose a method relying on the analysis of the error term and the use of the hydrodynamical
model to generate one part of the covariance function, the other part being modeled using a second-order stationary approach.
This method is discussed using a twin experiment in the case where a physical parameter is erroneous, and improves significantly
the results: the model bias is strongly reduced and the estimation error is well described. Moreover, it enables a better
adaptation of the Kalman gain to the actual estimation error. 相似文献
15.
AbstractA new methodology is proposed for the calibration of distributed hydrological models at the basin scale by constraining an internal model variable using satellite data of land surface temperature (LST). The model algorithm solves the system of energy and mass balances in terms of a representative equilibrium temperature that governs the fluxes of energy and mass over the basin domain. This equilibrium surface temperature, which is a critical model state variable, is compared to operational satellite LST, while calibrating soil hydraulic parameters and vegetation variables differently in each pixel, minimizing the errors. This procedure is compared to the traditional calibration using only discharge measurements. The distributed energy water balance model, Flash-flood Event-based Spatially-distributed rainfall–runoff Transformation – Energy Water Balance model (FEST-EWB), is used to test this approach. This methodology is applied to the Upper Yangtze River basin (China) using MODIS LST retrieved from satellite data in the framework of the NRSCC-ESA DRAGON-2 Programme. The calibration procedure based on LST seems to outperform the calibration based on discharge, with lower relative error and higher Nash-Sutcliffe efficiency index on cumulated volume.
Editor D. Koutsoyiannis; Associate editor C. Perrin 相似文献
16.
A new variational method is proposed for grid generation in numerical solution of groundwater flow problems requiring imposition of constraints on the space step. The proposed procedure utilizes the penalty function method for the solution of an optimization problem with inequality constraints. The efficiency of the method is illustrated by the generation of solution-dependent grids, which allow such constraints to be satisfied for complex-geometry domains. 相似文献
17.
18.
基于模型的波阻抗(速度)反演中,约束是不可缺少的,它可以有效地减少反演多解性和提高运算效率.国内外大量学者均详细论述了一些反演的约束手段,但对于地震道与地震道之间的约束,除了采用测井外推模式之外,就鲜有介绍,相同的约束条件直接应用于全部的地震道,不再变化,这样势必会增加运算时间和影响反演效果.针对这些问题,本文描述了一种地震相关性分段约束速度反演的方法,其目的在于增加内在的道与道之间的约束,主要采用了地震道的相关特征去优化约束条件与减少反演计算量.本文先分析了地震反演的不同内在约束方法,并利用地震的相关性去分段自动选择约束的强弱与范围大小,而不是在全部地震道反演中采用相同的约束参数,即可大大优化运算效率;其次设计两个数学模型进行反演验证,在反演效果比较符合实际的基础上,运算时间大大减少;最后引入中国东部某地区的实例对该方法进行检验,反演剖面上可以清楚观察出冲积扇的沉积特点,各期砂体清晰可见,便于追踪解释与储层预测,具有较高的分辨率与可信度.结果表明该方法是一种有利于提高反演效果和效率的方法. 相似文献
19.
The development of coastal ocean modeling in the recent years has allowed an improved representation of the associated complex
physics. Such models have become more realistic, to the point that they can now be used to design observation networks in
coastal areas, with the idea that a “good” network is a network that controls model state error. To test this ability without
performing data assimilation, we set up a technique called Representer Matrix Spectra (RMS) technique that combines the model
state and observation error covariance matrices into a single scaled representer matrix. Examination of the spectrum and the
eigenvectors of that matrix informs us on which model state error modes a network can detect and constrain amidst the observation
error background. We applied our technique to a 3D coastal model in the Bay of Biscay, with a focus on mesoscale activity,
and tested the performance of various altimetry networks and an in situ array deployment strategy. It appears that a single
nadir altimeter is not efficient enough at capturing coastal mesoscale physics, while a wide swath altimeter would do a much
better job. Testing various local in situ array configurations confirms that adding a current meter to a vertical temperature
measurement array improves the detection of secondary variability modes, while shifting the array higher on the shelf break
would obviously enhance the model constraint along the coast. The RMS technique is easily set up and used as a “black box,”
but the utility of its results is maximized by previous knowledge of model state error physics. The technique provides both
quantitative (eigenvalues) and qualitative (eigenvectors) tools to study and compare various network options. The qualitative
approach is essential to discard possibly inconsistent modes. 相似文献
20.
This paper addresses the parametric inverse problem of locating the point of release of atmospheric pollution. A finite set
of observed mixing ratios is compared, by use of least squares, with the analogous mixing ratios computed by an adjoint dispersion
model for all possible locations of the release. Classically, the least squares are weighted using the covariance matrix of
the measurement errors. However, in practice, this matrix cannot be determined for the prevailing part of these errors arising
from the limited representativity of the dispersion model. An alternative weighting proposed here is related to a unified
approach of the parametric and assimilative inverse problems corresponding, respectively, to identification of the point of
emission or estimation of the distributed emissions. The proposed weighting is shown to optimize the resolution and numerical
stability of the inversion. The importance of the most common monitoring networks, with point detectors at various locations,
is stressed as a misleading singular case. During the procedure it is also shown that a monitoring network, under given meteorological
conditions, itself contains natural statistics about the emissions, irrespective of prior assumptions. 相似文献