Tetrolet域利用时空相关性的云图压缩感知重建

针对卫星云图数据量大,但传输通道和存储空间相对狭小的问题,本文基于Tetrolet变换,利用相邻时次云图的时空相关性,实现了一种高重构质量的卫星云图压缩感知方法。该方法将善于表达图像方向纹理及边缘信息的Tetrolet变换引入压缩感知的稀疏表示环节,从而很好地体现了卫星云图细节丰富、纹理结构复杂的特性;同时,考虑到卫星云图序列间的相关性,将时间相邻的卫星云图组成图像组,以中间时刻云图作为参考图像,计算其与相邻时次云图的差异,通过在参考图片及序列差异图片间合理分配采样率,获取测量数据,在压缩感知框架下,采用带平滑处理的投影Landweber算法,重构出相邻时次的图像组。实验结果表明,Tetrolet变换适用于卫星云图的稀疏表示,而且图像组时空相关性的利用可显著提高重构云图的视觉效果及客观评价指标;在采样率低于0.2时,红外1、水汽和可见光3个通道重构云图的峰值信噪比(PSNR)较传统方法平均提高了7.48 dB,13.51 dB和6.15 dB。由此可见,本文方法可以通过获取少数随机测量值,重构出高质量的卫星云图,不仅为云图数据的低比特率压缩提供了一种可行的解决方案,而且对于其他序列图像的压缩采样具有借鉴意义。     

过完备字典稀疏表示的云图超分辨率算法

提出一种基于过完备字典稀疏表示的云图超分辨率算法。首先,联合训练针对低分辨率与高分辨率云图块的两个字典Dl和Dh,保证对应的低分辨率与高分辨率云图块关于各自的字典具有相似的稀疏表示;其次,通过求解优化问题,获得待处理云图每个低分辨率云图块关于Dl的稀疏表示,并将表示系数用于Dh以生成对应的高分辨率云图块;最后,运用最速下降算法,得到满足重构约束的高分辨率云图。红外与可见光云图的数值实验验证了本文算法的有效性,表明本文算法在视觉效果及PSNR指标上均优于插值方法。     

基于多分辨率分析的GF-5和GF-1遥感影像空—谱融合

针对空间分辨率比率较大尺度差异下的高分五号(GF-5)与高分一号(GF-1)卫星影像的空—谱融合问题,提出多传感器影像融合策略:一方面,通过现有空—谱融合方法的分步融合得到融合影像;另一方面,在分步融合理论基础上,推导得出一体化融合基础框架,并进一步提出基于多分辨率分析的多传感器一体化融合方法,缓解现有方法因空间分辨率比率过大导致影像空、谱互补信息难以有效集成的问题。其中,提出的一体化融合方法基于调制传递函数MTF (Modulation Transfer Function)滤波对多传感器影像空间(高频)和光谱(低频)分量进行分解提取,并充分考虑多传感器高空间分辨率影像与高光谱分辨率影像之间的关系,以及高光谱分辨率影像波段间关系,设计合理的融合权重,最终可得到具有最高空间分辨率和最高光谱分辨率的融合影像。通过GF-1全色影像、GF-1多光谱影像、GF-5高光谱影像数据对提出方法进行实验验证,结果表明:本文方法可有效集成多传感器影像间的空、谱互补信息,得到较优融合结果。     

Deriving effects of pressure depletion on elastic framework moduli from sonic logs

We have developed a method to determine the effect of pore pressure depletion on elastic framework moduli, by utilizing sonic logs from wells drilled at different locations through a reservoir at varying depletion stages. This is done by first inverting the sonic logs for elastic framework bulk and shear moduli, thus carefully removing pressure dependent fluid effects. By crossplotting these elastic framework moduli against an increase in net stress (which is directly related to depletion), we derive the stress sensitivity of the elastic framework moduli. We found that the observed stress sensitivity was consistent with time-lapse seismic results and that the sensitivity derived from the sonic logs was much smaller than predicted by hydrostatic measurements on core samples. This method is applicable to depletion scenarios in mature fields and can be used both for modelling and inverting time-lapse seismic data for effects of pore pressure depletion on seismic data.     

主题学习和稀疏表示的MODIS图像超分辨率重建

针对MODIS图像分辨率受传感器限制和噪声干扰,且分辨率局限在一定水平等问题,提出一种采用主题学习和稀疏表示的MODIS图像超分辨率重建方法,该方法通过双边滤波将MODIS图像的平滑及纹理部分分离,并将纹理部分看成是由若干"文档"组成的训练样本;运用概率潜在语义分析提取"文档"的潜在语义特征,从而确定"文档"所属的"主题"。在此基础上,针对每个主题所对应的图像块,采用改进的K-SVD方法训练若干适用于不同主题的高低分辨率字典对,从而可以运用这些字典对,通过稀疏编码实现测试图像相应主题块的超分辨率重建。实验结果表明,重建图像在视觉效果和PSNR等指标上均优于传统方法。     

Recruitment of shrimp (Pandalus borealis) in the Barents Sea related to spawning stock and environment

The shrimp spawn in autumn, and the females carry their eggs as out roe until spring when the larvae hatch. Within a period of 2 months the shrimp larvae settle to the bottom. It has been claimed that the year-class strength probably is determined during the larval phase. Today's assessment and forecast of the shrimp stock productivity and potential fishing yields are weak. This is partly due to poor knowledge on population dynamics from hatching until the shrimp are caught in the fishery at the age of 3 or 4 years. We, therefore, here identify the most important abiotic and biotic factors that affect recruitment in addition to spawning stock biomass. Since 1995, a net attached to the underbelly of the survey trawl used at the annual cruise in the Barents Sea has caught juvenile shrimp. The abundance of settled shrimp larvae varies in time and space. The recruitment to the fishery has been quite stable with the exception of the 1996 year-class, which was observed as 1-year-olds but has not been registered since. The temporal pattern of the three youngest year-classes is studied in relation to abiotic factors such as sea temperature, ice index and North Atlantic Oscillation, as well as biotic factors such as spawning stock biomass and presence of copepods, euphausiids and predating cod. Recruitment indices and factors identified by the Spearmann correlation to be significantly correlated with recruitment were used as input in a principal component analysis (PCA) and a generalized additive model (GAM) was applied. Abundance of 1-year-old shrimp is positively correlated to spawning stock biomass the previous year and to temperature of the previous winter, and negatively correlated with the number of 1-year-old cod. Two-year-old shrimp show significant correlation with temperature, whereas there is a strong negative correlation with euphausiids. Three-year-old shrimp are significantly correlated with the number of 2-year-old shrimp the previous year but negatively correlated to temperature at sampling time. This is probably due to less overlap with the main predator cod when cold. Ricker functions indicate an increased density-dependent mortality with age. When predicting the recruitment of shrimp to the fishery, the spawning stock biomass, the abundance of cod and euphausiids, as well as the temperature should be included.     

On the possibility of a kinetic fractionation of nitrogen stable isotopes during natural diamond growth

In a diamond from New South Wales (Australia), cubic and octahedral growth sectors, as identified by cathodoluminescence (CL), show slight differences in N-contents of 29 and 42 ppm respectively but no significant differences in either δ¹³C, δ¹⁵N and nitrogen aggregation state with values at +1.96‰, +19.4‰, and 25% Type IaAB aggregation, respectively.Two gem cubes from the Orapa kimberlite (Botswana) were studied by CL revealing a nonfaceted cubic growth. Accordingly, nine other gem cubes were combusted and yielded δ¹³C-values from -5.33‰ to -6.63‰, δ¹⁵N from -1.0‰ to -5.5‰, and nitrogen contents from 914 to 1168 ppm, with nitrogen aggregation state being only Type IaA (zero % B). The gem cubes show striking similarities to fibrous/coated diamonds, not only in both δ¹³C ranges (less than 3‰ from -5 to -8‰), but also in the high levels of nitrogen (≈ 1000 ppm), suggesting that the two diamond types are related. Additionally, no δ¹⁵N variation was detected between the cube and octahedral growth sectors of the Australian diamond, in the cube sectors of the nine gem cubes from Botswana, nor in fibrous/coated diamonds previously studied. These analyses contrast with an earlier study on a synthetic diamond, which reported a strong kinetic fractionation of N-isotopes of about 40‰ between cube and octahedral growth. The present evidence, therefore, suggests that kinetic fractionation of N-isotopes does not operate during natural diamond formation.     

Object Models with Vector Steering

Object models are widely used to model the distribution of facies in a reservoir. Several computer programs exist for modelling fluvial channels or more general facies objects. This paper focuses on a marked point model with objects that are able to orient locally according to a vector field. In this way, objects with locally varying curvature are created. With this kind of objects it is possible to model complex depositional basins, that are not easily modelled with conventional methods. The new object type is called Backbone objects. The objects have a piecewise linear centerline and are able to follow the direction of a three-dimensional vector field locally in lateral and vertical direction. How well the objects follow the vector field is determined by three parameters. Use of different coordinate systems and mapping between the systems make it possible to generate Gaussian random fields that follow the shape and direction of the objects. The Gaussian fields can be used to model petrophysical variables, which is important for fluid flow modelling.     

Fault displacement modelling using 3D vector fields

In history matching and sensitivity analysis, flexibility in the structural modelling is of great importance. The ability to easily generate multiple realizations of the model will have impact both on the updating workflow in history matching and uncertainty studies based on Monte Carlo simulations. The main contribution to fault modelling by the work presented in this paper is a new algorithm for calculating a 3D displacement field applicable to a wide range of faults due to a flexible representation. This gives the possibility to apply this field to change the displacement and thereby moving horizons and fault lines. The fault is modelled by a parametric format where the fault has a reference plane defined by a centre point, dip and strike angles. The fault itself is represented as a surface defined by a function z = f(x,y), where x, y and z are coordinates local to the reference plane, with the z-axis being normal to the plane. The displacement associated with the fault outside the fault surface is described by a 3D vector field. The displacement on the fault surface can be found by identifying the intersection lines between horizons and the fault surface (fault lines), and using kriging techniques to fill in values at all points on the surface. Away from the fault surface the displacement field is defined by a monotonic decreasing function which ensures zero displacement at a specified distance from the fault. An algorithm is developed where the displacement can be increased or decreased according to user-defined parameters. This means that the whole displacement field is changed and points on horizons around the fault can be moved accordingly by applying the modified displacement field on them. The interaction between several faults influencing the same points is taken care of by truncation rules and the ordering of the faults. The method is demonstrated on a realistic synthetic case based on a real reservoir.