Inversion of dispersion coeffcient in water quality model using optimal perturbation algorithm

As a primary parameter in the water quality model for shallow bays, the dispersion coefficient is traditionally determined with a trial-and-error method, which is time-consuming and requires much experience. In this paper, based on the measured data of chemical oxygen demand （COD）, the dispersion coefficient is calculated using an inversion method. In the process, the regularization method is applied to treat the ill-posedness, and an operator identity perturbation method is used to obtain the solu- tion. Using the model with an inverted dispersion coefficient, the distributions of COD, inorganic nitrogen （IN）, and inorganic phosphorus （IP） in Bohai Bay are predicted and compared with the measured data. The results indicate that the method is feasible and the inverted dispersion coefficient can be used to predict other pollutant distribution. This method may also be further extended to the inversion of other parameters in the water quality model.

Temperature inversion in the Huanghai Sea bottom cold water in summer

Using conductivity-Temperature-depth data of a recent cruise during July 22-28, 2008 and historical data, it is found that temperature inversions occur from time to time in the Huanghai Sea（Yellow Sea） cold water mass （HSCWM） in summer. The temperature inversions are produced by the movement of the fresh and cold HSCWM masses above the warm and saline Huanghai Sea Warm Current water at the central bottom of the Huanghai Sea Trough. The non-homogeneous profiles of the temperature and the salinity suggest that vertical mixing in the HSCWM, which is of great importance to the circulation in the Huanghai Sea in summer, is weak. Trajectories of satellite-tracked surface drifters suggest that waters in the northern reach of the Huanghai Sea move southward along the 40-50 m isobaths and descend into the southern Huanghai Sea to form the western core of the HSCWM.

Inversion of Evaporation and Water Vapor Transport Using HY-2 Multi-Sensor Data

HY-2 satellite is the first marine dynamic environment satellite of China.In this study,global evaporation and water vapor transport of the global sea surface are calculated on the basis of HY-2 multi-sensor data from April 1 to 30,2014.The algorithm of evaporation and water vapor transport is discussed in detail,and results are compared with other reanalysis data.The sea surface temperature of HY-2 is in good agreement with the ARGO buoy data.Two clusters are shown in the scatter plot of HY-2 and OAFlux evaporation due to the uneven global distribution of evaporation.To improve the calculation accuracy,we compared the different parameterization schemes and adopted the method of calibrating HY-2 precipitation data by SSM/I and Global Precipitation Climatology Project(GPCP)data.In calculating the water vapor transport,the adjustment scheme is proposed to match the balance of the water cycle for data in the low latitudes.

The sensitivity of the array resistivity log to mud filtrate invasion and its primary five-parameter inversion for improved oil water recognition

In order to improve reservoir fluid recognition, the sensitivity of array resistivity response to the difference of the invasion properties in both oil-bearing layers and water layers is analyzed. Then the primary inversion is carried out based on the array resistivity log. The mud invasion process is numerically simulated based on the oil-water flow equation and water convection diffusion equation. The results show that the radial resistivity of a fresh mud-invaded oil-bearing layer presents complex distribution characteristics, such as nonlinear increase, increasing to decreasing and low resistivity annulus, and the resistive invasion profile of a water layer is monotonic. Under specific conditions, array resistivity log can reflect these changes and the array induction log is more sensitive. Nevertheless, due to the effect of factors like large invasion depth, reservoir physical and oil-bearing properties, the measured apparent resistivity may differ greatly from the actual mud filtrate invasion profile in an oil-bearing layer. We proposed a five-parameter formation model to simulate the complex resistivity distribution of fresh mud-invaded formation. Then, based on the principle of non-linear least squares, the measured array resistivity log is used for inversion with the Marquardt method. It is demonstrated that the inverted resistivity is typically non-monotonic in oil-bearing layers and is monotonic in water layers. Processing of some field data shows that this is helpful in achieving efficient reservoir fluid recognition.

Study on Volume Strain Inversion from Water Level Change of Well-aquifer Systems

Based on linear poroelasfic and hydrogeology theory, a mathenmatical expression describing the relationship between water level clmnge and aquifer volume strain is put forward. Combined with earth tidal theory, we analyze the response characteristics from well-aquifer water level change to earth tide of volume strain and present a method of volume strain inversion from water level clmnge. Comparing the results of inversion with real observed data, we found that there is a good consistency. This suggests that the method of volume strain inversion from water level clmnge is proper. It will offer a reference for learning about hydrogeology characteristics, volume strain and searching for precursor anomalies.

Underground water quality model inversion of genetic algorithm

The underground water quality model with non-linear inversion problem is ill-posed, and boils down to solving the minimum of nonlinear function. Genetic algorithms are adopted in a number of individuals of groups by iterative search to find the optimal solution of the problem, the encoding strings as its operational objective, and achieving the iterative calculations by the genetic operators. It is an effective method of inverse problems of groundwater, with incomparable advantages and practical significances.

Sea Water Intrusion Modeling in Rashid Area of Nile Delta (Egypt) via the Inversion of DC Resistivity Data

The current research focuses on the detection of sea water intrusion in Rashid area which is located about 75 km east to Alexandria, Egypt. For this purpose, geoelectrical survey was carried out using the Schlumberger Vertical Electric Sounding (VES) to identify freshwater thickness, sea water intrusion and estimate subsurface lithology. Seventeen VES stations were measured with current electrode separation (AB/2) ranging from 1.5 m to 100 m. Then, the VES data was interpreted using 1-D and 2-D inversion schemes of DC resistivity data based on least squares method with smoothness constrains. The inverted resistivity distribution at relatively shallow depth shows an important low resistivity zone that probably reflects salt water alteration zone (northern parts). Depth to the freshwater bearing layer reaches its maximum at the south and decreases towards the north. From quantitative interpretation, invasion of salt water started at depth about 10 m at north in the thickness of freshwater bearing layer ranging from 15 to 25 m, while at depth of about 120 m all the layers were saturated with salt water.

Comparative Study on Coastal Depth Inversion Based on Multi-source Remote Sensing Data

Coastal depth is an important research focus of coastal waters and is also a key factor in coastal environment. Dongluo Island in South China Sea was taken as a typical study area. The band ratio model was established by using measured points and three multispectral images of Landsat-8, SPOT-6(Systeme Probatoire d'Observation de la Terre, No.6) and WorldView-2. The band ratio model with the highest accuracy is selected for the depth inversion respectively. The results show that the accuracy of SPOT-6 image is the highest in the inversion of coastal depth. Meanwhile, analyzing the error of inversion from different depth ranges, the accuracy of the inversion is lower in the range of 0–5 m because of the influence of human activities. The inversion accuracy of 5–10 m is the highest, and the inversion error increases with the increase of water depth in the range of 5–20 m for the three kinds of satellite images. There is no linear relationship between the accuracy of remote sensing water depth inversion and spatial resolution of remote sensing data, and it is affected by performance and parameters of sensor. It is necessary to strengthen the research of remote sensor in order to further improve the accuracy of inversion.

Spatiotemporal variation and mechanisms of temperature inversion in the Bay of Bengal and the eastern equatorial Indian Ocean

In the northern Bay of Bengal,the existence of intense temperature inversion during winter is a widely accepted phenomenon.However,occurrences of temperature inversion during other seasons and the spatial distribution within and adjacent to the Bay of Bengal are not well understood.In this study,a higher resolution spatiotemporal variation of temperature inversion and its mechanisms are examined with mixed layer heat and salt budget analysis utilizing long-term Argo(2004 to 2020)and RAMA(2007 to 2020)profiles data in the Bay of Bengal and eastern equatorial Indian Ocean(EEIO).Temperature inversion exists(17.5%of the total 39293 Argo and 51.6%of the 28894 RAMA profiles)throughout the year in the entire study area.It shows strong seasonal variation,with the highest occurrences in winter and the lowest in spring.Besides winter inversion in the northern Bay of Bengal,two other regions with frequent temperature inversion are identified in this study for the first time:the northeastern part of the Bay of Bengal and the eastern part of the EEIO during summer and autumn.Driving processes of temperature inversion for different subregions are revealed in the current study.Penetration of heat(mean~25 W/m;)below the haline-stratified shallow mixed layer leads to a relatively warmer subsurface layer along with the simultaneous cooling tendency in mixed layer,which controls more occurrence of temperature inversion in the northern Bay of Bengal throughout the year.Comparatively lower cooling tendency due to net surface heat loss and higher mixed layer salinity leaves the southern part of the bay less supportive to the formation of temperature inversion than the northern bay.In the EEIO,slightly cooling tendency in the mixed layer along with the subduction of warm-salty Arabian Sea water beneath the cold-fresher Bay of Bengal water,and downwelling of thermocline creates a favorable environment for forming temperature inversion mainly during summer and autumn.Deeper isothermal layer depth,and thicker barrier layer thickness intensify the temperature inversion both in the Bay of Bengal and EEIO.

Influence of Suspension Parameters on Water Absorbency of Starch-g-poly( sodium acrylate) Synthesized by Inverse Suspension Polymerization

Superabsorbents starch grafted sodium polyacrylate was prepared by inverse suspension polymerization, using toluene as the continuous phase, potassium persulfate as the initiator. The effect of suspension parameters, such as volume ratio of continuous phase and dispersed phase, type and dosage of suspending agents, on water absorbency of the starch grafted polymer was studied. Different starch derivatives were also investigated. Superabsorbents made of cationic starch has higher water absorbency than that made of native corn starch.

A fast inversion method for ocean parameters based on dispersion curves with a single hydrophone

The dispersion characteristics of shallow water can be described by the dispersion curves,which contain substantial ocean parameter information.A fast ocean parameter inversion method based on dispersion curves with a single hydrophone is presented in this paper.The method is achieved through Bayesian theory.Several sets of dispersion curves extracted from measured data are used as the input function.The inversion is performed by matching a replica calculated with a dispersion formula.The bottom characteristics can be described by the bottom reflection phase shift parameter P.The propagation range and the depth can be inverted quickly when the seabed parameters are represented by on parameter P.The inversion results improve the inversion efficiency of the seabed parameters.Consequently,the inversion efficiency and accuracy are improved while the number of inversion parameters is decreased and the computational speed of replica is increased.The inversion results have lower error than the reference values,and the dispersion curves calculated with inversion parameters are also in good agreement with extracted curves from measured data;thus,the effectiveness of the inversion method is demonstrated.

Application of homotopy parameter inversion method in Miyun Reservoir

The large-scale convergence of homotopy parametric inversion method on the water quality model parameters calculated was used,with application in parametric inversion calculation of total phosphorus of Beijing Miyun Reservoir.Through calculated and compared the error of sedimentation rate by homotopy parametric inversion method and genetic inversion calculation method,the results indicate that homotopy parametric inversion method has good stability,calculating speed,and even if the initial selection away from the objective function,the solution still has a good convergence.

Water Quality Evaluation of Chapurson Valley in Hunza Nagar, Gilgit Baltistan, Pakistan, Based on Statistical Analysis and Water Quality Index

Water borne ailments are of serious public health concern in Gilgit Baltistan’s (GB) region of Pakistan. The pollution load on the glacio-fluvial streams and surface water resources of the Chapurson Valley in the Hunza Nagar area of the GB is increasing as a result of anthropogenic activities and tourism. The present study focuses on the public health quality of drinking water of Chapurson valley. The study addressed the fundamental drinking water quality criteria in order to understand the state of the public health in the valley. To ascertain the current status of physico-chemical, metals, and bacteriological parameters, 25 water samples were collected through deterministic sampling strategy and examined accordingly. The physico-chemical parameters of the water samples collected from the valley were found to meet the World Health Organization (WHO) guidelines of drinking water. The water samples showed a pattern of mean metal concentrations in order of Arsenic (As) > Lead (Pb) > Iron (Fe) > Zinc (Zn) > Copper (Cu) > Magnesium (Mg) > Calcium (Ca). As, Cu, Zn, Ca and Mg concentration were under the WHO guidelines range. However, results showed that Pb and Fe are present at much higher concentrations than recommended WHO guidelines. Similarly, the results of the bacteriological analysis indicate that the water samples are heavily contaminated with the organisms of public health importance (including total coliforms (TCC), total faecal coliforms (TFC) and total fecal streptococci (TFS) are more than 3 MPN/100mL). Three principal components, accounting for 48.44% of the total variance, were revealed using principal component analysis (PCA). Bacteriological parameters were shown to be the main determinants of the water quality as depicted by the PCA analysis. The dendrogram of Cluster analysis using the Ward’s method validated the same traits of the sampling locations that were found to be contaminated during geospatial analysis using the Inverse Distance Weight (IDW) method. Based on these findings, it is most likely that those anthropogenic activities and essentially the tourism results in pollution load from upstream channels. Metals may be released into surface and groundwater from a few underlying sources as a result of weathering and erosion. This study suggests that the valley water resources are more susceptible to bacteriological contamination and as such no water treatment facilities or protective measure have been taken to encounter the pollution load. People are drinking the contaminated water without questioning about the quality. It is recommended that the water resources of the valley should be monitored using standard protocol so as to protect not only the public health but to safe guard sustainable tourism in the valley.

兰州市南北两山土壤水分遥感反演及植被需水量估算

探究西北干旱区土壤水分和植被需水量动态变化特征,可为生态恢复不同阶段所需水资源量及水资源优化配置提供科学依据。以兰州市南北两山为研究区,基于Sentinel-2 L2A和Landsat 8 OLI遥感影像,结合实测土壤0~10 cm的111个数据,分别构建垂直干旱指数(Perpendicular Drought Index,PDI)、改进型垂直干旱指数(Modified Perpendicular Drought Index,MPDI)和植被调整垂直干旱指数(Vegetation-adjusted Perpendicular Drought Index,VAPDI)土壤水分反演模型,并采用4种模型指标定量决定系数(R^(2))、平均绝对误差(MAE)、平均相对误差(MRE)、均方根误差(RMSE)对模型反演的效果进行精度评价,选出最优的土壤水分反演模型并结合土壤水分限制系数,与研究区2019年林地、草地和耕地植被面积的空间数据、各站点生长季内的参考作物蒸散量,构建植被生态需水量模型,厘清研究区内土壤水分、植被需水量时空变化特征。结果表明:(1)2种数据源下的PDI、MPDI、VAPDI和实测数据之间均有着不同程度的线性负相关性,其中R^(2)分别为0.37、0.64和0.59,从评价指标的结果来看,MPDI的土壤水分回归模型的拟合决定系数最高,2种遥感数据反演的土壤水分空间分布格局具有一致性。(2)分辨率高的Sentinel-2L2A土壤水分反演更加精细,土壤水分整体呈波动增长趋势,多时段土壤水分的平均值为23.27%,呈现出降低再增加然后下降,总体增幅为74.07%。(3)兰州市南北两山4—10月植被需水量月均值也呈现先增加后下降的趋势,与土壤水分含量变化具有一致性,4—10月中7月植被需水量最大,为3.98×10^(7)m^(3),10月植被生态需水量最小,为0.97×10^(7)m^(3)。随着环境绿化工程的实施,兰州市南北两山从只能生长耐旱草本和低矮灌木的植物,逐步形成以多种类结合的群落结构。本研究可为兰州市南北两山土壤水资源合理利用及植被恢复提供参考。

基于Sentinel-2多光谱遥感影像的小浪底水质反演

多光谱遥感技术可根据遥感波段信息反演水质参数,降低监测成本,提高监测速度和质量,为大范围水环境监测提供了一种新的方法。通过分析小浪底水库的Sentinel-2多光谱影像以及采样点实测水质数据,建立了最佳光谱波段的水质参数反演模型,对小浪底水库的化学需氧量(COD)、总磷(TP)、总氮(TN)和氨氮(NH_3-N)进行了遥感反演,验证了反演模型的精确度和稳定性,并反演了各水质参数的空间分布规律。结果表明:在4种水质参数反演模型中,COD模型精确度和稳定性最高,其次是TP、TN,最低的是NH_3-N,水库出水口和部分边缘COD质量浓度较高,水库中心TN、TP和NH_3-N质量浓度高于边缘处。

敏感属性波形反演技术在大庆长垣水敏储层中的应用

针对大庆长垣油田萨零组储层水敏强、含油性差、厚度薄,以及河道窄小等问题,以长垣油田杏北开发区杏3区块为例,结合测井资料、室内实验及取心井综合解释资料,通过两步法协同敏感属性分析,判别与含油性及水敏强弱关联性最高的属性参数,然后结合属性参数范围,利用基于波形特征指示的井震联合反演技术,精准预测萨零组弱水敏高含油优质储层,从而指导井位部署。结果表明:声波时差与电阻率组合能够划分油层、差油层和非储层,自然伽马可以判别水敏强弱;利用反演预测结果部署1口新井,岩心分析数据表明该井水敏指数为0.62,为中等偏强水敏,含油级别以油浸和含油为主,含油饱和度平均为61.4%,反演结果精准可靠;该新井投产1 a,平均日产油3.5 t,累计年产油1090.6 t,投产效果较好。研究成果为萨零组开发奠定基础,对大庆油田的高质量稳产意义重大。

降雨和库水位变动联合作用下滑坡位移加权反分析及监测预警方法

降雨和库水位变动是诱发大型库岸滑坡体持续变形甚至失稳的主要外在因素。针对以上两种因素联合作用下滑坡体的变形特征,本文提出一种基于多源数据的“反演-评估-融合-预警”四步式滑坡监测预警方法。该方法采用趋势项位移进行力学参数反演,有效避免了降雨、库水位变动等外在因素对反演结果的影响。同时,采用熵权值构建了反演目标函数,借此考虑了不同位置监测点对反演结果的影响程度。此外,考虑了滑坡体自身的地质条件和坡体结构对监测预警的影响,综合运用安全系数和综合变形速率构建了具有地质和力学基础的预警模型。采用本文方法对贵州三板溪水电站东岭信滑坡体进行了力学参数反演和监测预警,结果表明:该滑坡体整个监测期间的安全系数在1.25~1.70之间波动,整体状态呈下降趋势,基于安全系数的角度判断其尚处于稳定状态;综合变形速率呈现上、下波动趋势,与安全系数在变化规律上具有高度的相关性,两者可以从整体上表征该滑坡体在降雨和库水位变动联合作用下稳定性的演进过程;当前监测时刻点的综合变形速率是0.0599,根据分级预警模型判断此滑坡体当前处于稳定变形阶段(Ⅰ级预警),与安全系数的判断结果具有一致性。

基于时空融合算法的水体叶绿素a反演研究

为了准确反演水体中叶绿素a浓度,以黄柏河东支流域为例,采用STNLFFM时空融合算法,对2017年GF-4和Sentinel-2反射率数据进行融合,以重构Sentinel-2影像的时间序列数据,并对应用算法前后获取的水质参数-光谱特征响应关系建立多元线性回归模型,比较模型对叶绿素a的预测效果以验证时空融合算法的可行性,利用重构后影像光谱特征与水质参数的响应关系建立人工神经网络模型,反演2017年黄柏河东支流域各水库水体叶绿素a浓度。结果表明:利用时空融合算法生成的影像接近真实影像,提高了多元线性回归模型预测叶绿素a的效果,R2从融合前0.659提高至融合后0.844,且基于时空融合算法获取的水质参数-光谱关系建立的人工神经网络模型模拟精度较好,R2和MRE达到0.925和9.461%,反演的叶绿素a浓度空间差异性明显。证明了时空融合算法在水质参数反演过程中具有较好的应用前景。

基于多源遥感数据与模型对比的冬小麦土壤含水量区域监测研究

实时、精准的土壤水分含量监测是农业用水管理的基础,探究冬小麦土壤水分反演的最优模型对于提高农业用水效率和可持续发展均具有重要的意义。本研究以河南省鹤壁市浚县冬小麦种植区域的土壤水分含量为研究对象,采用无人机遥感数据、卫星遥感数据、田间采样数据,分别运用温度植被干旱指数模型、水云模型和改进的水云模型3种方法,进行土壤含水量反演对比分析与最优模型选择。结果表明,3种方法中10 cm深度的反演精度均高于20 cm,且R^(2)均大于0.4。其中采用改进的水云模型方法在10 cm深度的R^(2)为0.7055、RMSE为0.0209,20 cm深度的R^(2)为0.5069、RMSE为0.0271,优于水云模型和温度植被干旱指数的反演效果。因此,改进的水云模型是一种适合用于麦田土壤水分反演的方法,它能够提供较高的反演精度。

双线大直径输水隧道建设对既有桥梁及桩基的影响分析

为准确分析双线大直径输水隧道下穿对既有桥梁结构的影响,构建桥梁-地层-桩基-盾构隧道高精度数值仿真模型,对桥梁桩基结构及盾构过程进行精细化模拟,结合基于智能算法的土层参数反演模型,获取更接近真实情况的土层参数,以提高数值模拟精度。结果表明:输水隧道盾构开挖引起的地表最大沉降量为15.01 mm,基于反演的地表变形计算结果与实际监测值的误差减少约70%;盾构隧道双线开挖会引起地表二次沉降变形,两条盾构隧道中间区域的地表二次沉降量最大,约占总沉降量的60%,最大二次沉降量为5.06 mm;隧道两侧的桩基受到影响较大,最大位移为10.03 mm,位于桩基顶部,而隧道开挖对桩基内应力分布影响较小,可以认为无影响。