High-resolution seismic inversion method based on joint data-driven in the time-frequency domain

Seismic inversion can be divided into time-domain inversion and frequency-domain inversion based on different transform domains.Time-domain inversion has stronger stability and noise resistance compared to frequencydomain inversion.Frequency domain inversion has stronger ability to identify small-scale bodies and higher inversion resolution.Therefore,the research on the joint inversion method in the time-frequency domain is of great significance for improving the inversion resolution,stability,and noise resistance.The introduction of prior information constraints can effectively reduce ambiguity in the inversion process.However,the existing modeldriven time-frequency joint inversion assumes a specific prior distribution of the reservoir.These methods do not consider the original features of the data and are difficult to describe the relationship between time-domain features and frequency-domain features.Therefore,this paper proposes a high-resolution seismic inversion method based on joint data-driven in the time-frequency domain.The method is based on the impedance and reflectivity samples from logging,using joint dictionary learning to obtain adaptive feature information of the reservoir,and using sparse coefficients to capture the intrinsic relationship between impedance and reflectivity.The optimization result of the inversion is achieved through the regularization term of the joint dictionary sparse representation.We have finally achieved an inversion method that combines constraints on time-domain features and frequency features.By testing the model data and field data,the method has higher resolution in the inversion results and good noise resistance.

A Method in Mixed Frequency-time domain for Fatigue Life Estimation of Steel Girders of Cable-stayed Bridges Due to Buffeting

发展了斜拉桥钢主梁抖振疲劳寿命的混合时－频域方法．在本方法中，首先建立了桥面处的风速风向概率密度函数，导出的应力谱因其包括背景分量而不再是窄带谱．据应力谱模拟出了疲劳应力的时间历程，并用雨流计数法统计了其幅值特性．基于修正的Ｍｉｎｅｒ定理导出了考虑风速风向作用的桥梁主梁抖振疲劳寿命的估算公式，并用此分析了杨浦大桥钢主梁的抖振疲劳寿命．结果表明，风向对抖振疲劳寿命有很大影响；杨浦大桥主梁的疲劳寿命远大于设计寿命．

Three-dimensional numerical modeling of gravity anomalies based on Poisson equation in spacewavenumber mixed domain

In gravity-anomaly-based prospecting, the computational and memory requirements for practical numerical modeling are potentially enormous. Achieving an efficient and precise inversion for gravity anomaly imaging over large-scale and complex terrain requires additional methods. To this end, we have proposed a new topography-capable By performing a two-dimensional Fourier transform in the horizontal directions, threedimensional partial differential equations in the spatial domain were transformed into a group of independent, one-dimensional differential equations engaged with different wave numbers. These independent differential equations are highly parallel across different wave numbers. differential equations with different wave numbers, and the efficiency of solving fixedbandwidth linear equations was further improved by a chasing method. In a synthetic test, a prism model was used to verify the accuracy and reliability of the proposed algorithm by comparing the numerical solution with the analytical solution. We studied the computational precision and efficiency with and without topography using different Fourier transform methods. The results showed that the Guass-FFT method has higher numerical precision, while the standard FFT method is superior, in terms of computation time, for inversion and quantitative interpretation under complicated terrain.

Underdetermined DOA estimation and blind separation of non-disjoint sources in time-frequency domain based on sparse representation method

This paper deals with the blind separation of nonstation-ary sources and direction-of-arrival （DOA） estimation in the under-determined case, when there are more sources than sensors. We assume the sources to be time-frequency （TF） disjoint to a certain extent. In particular, the number of sources presented at any TF neighborhood is strictly less than that of sensors. We can identify the real number of active sources and achieve separation in any TF neighborhood by the sparse representation method. Compared with the subspace-based algorithm under the same sparseness assumption, which suffers from the extra noise effect since it can-not estimate the true number of active sources, the proposed algorithm can estimate the number of active sources and their cor-responding TF values in any TF neighborhood simultaneously. An-other contribution of this paper is a new estimation procedure for the DOA of sources in the underdetermined case, which combines the TF sparseness of sources and the clustering technique. Sim-ulation results demonstrate the validity and high performance of the proposed algorithm in both blind source separation （BSS） and DOA estimation.

Surface wave attenuation based polarization attributes in time-frequency domain for multicomponent seismic data

In the paper, we propose a surface wave suppression method in time-frequency domain based on the wavelet transform, considering the characteristic difference of polarization attributes, amplitude energy and apparent velocity between the effective signals and strong surface waves. First, we use the proposed method to obtain time-frequency spectra of seismic signals by using the wavelet transform and calculate the instantaneous polarizability at each point based on instantaneous polarization analysis. Then, we separate the surface wave area from the signal area based on the surface-wave apparent velocity and the average energy of the signal. Finally, we combine the polarizability, energy, and frequency characteristic to identify and suppress the signal noise. Model and field data are used to test the proposed filtering method.

Numerical modeling of 3D DC resistivity method in the mixed space-wavenumber domain

Forward modeling is the basis of inversion imaging and quantitative interpretation for DC resistivity exploration.Currently,a numerical model of the DC resistivity method must be finely divided to obtain a highly accurate solution under complex conditions,resulting in a long calculation time and large storage.Therefore,we propose a 3D numerical simulation method in a mixed space-wavenumber domain to overcome this challenge.The partial differential equation about abnormal potential is transformed into many independent ordinary differential equations with different wavenumbers using a 2D Fourier transform along the x axis and y axis direction.In this way,a large-scale 3D numerical simulation problem is decomposed into several 1D numerical simulation problems,which significantly reduces the computational and storage requirements.In addition,these ordinary 1D differential equations with different wavenumbers are independent of each other and high parallelelism of the algorithm.They are solved using a finite-element algorithm combined with a chasing method,and the obtained solution is modified using a contraction operator.In this method,the vertical direction is reserved as the spatial domain,then grid size can be determined flexibly based on the underground current density distribution,which considers the solution accuracy and calculation efficiency.In addition,for the first time,we use the contraction operator in the integral equation method to iterate the algorithm.The algorithm takes advantage of the high efficiency of the standard Fourier transform and chasing method,as well as the fast convergence of the contraction operator.We verified the accuracy of the algorithm and the convergence of the contraction operator.Compared with a volume integral method and goal-oriented adaptive finite-element method,the proposed algorithm has lower memory requirements and high computational efficiency,making it suitable for calculating a model with large-scale nodes.Moreover,different examples are used to verify the high adaptability and parallelism of the proposed algorithm.The findings show that the 3D numerical simulation method of DC resistivity method in a mixed space-wavenumber domain is highly efficient,precise,and parallel.

Bayesian seismic multi-scale inversion in complex Laplace mixed domains

Seismic inversion performed in the time or frequency domain cannot always recover the long-wavelength background of subsurface parameters due to the lack of low-frequency seismic records. Since the low-frequency response becomes much richer in the Laplace mixed domains, one novel Bayesian impedance inversion approach in the complex Laplace mixed domains is established in this study to solve the model dependency problem. The derivation of a Laplace mixed-domain formula of the Robinson convolution is the first step in our work. With this formula, the Laplace seismic spectrum, the wavelet spectrum and time-domain reflectivity are joined together. Next, to improve inversion stability, the object inversion function accompanied by the initial constraint of the linear increment model is launched under a Bayesian framework. The likelihood function and prior probability distribution can be combined together by Bayesian formula to calculate the posterior probability distribution of subsurface parameters. By achieving the optimal solution corresponding to maximum posterior probability distribution, the low-frequency background of subsurface parameters can be obtained successfully. Then, with the regularization constraint of estimated low frequency in the Laplace mixed domains, multi-scale Bayesian inversion inthe pure frequency domain is exploited to obtain the absolute model parameters. The effectiveness, anti-noise capability and lateral continuity of Laplace mixed-domain inversion are illustrated by synthetic tests. Furthermore,one field case in the east of China is discussed carefully with different input frequency components and different inversion algorithms. This provides adequate proof to illustrate the reliability improvement in low-frequency estimation and resolution enhancement of subsurface parameters, in comparison with conventional Bayesian inversion in the frequency domain.

Chromodomain-helicase-DNA binding protein 5, 7 and pronecrotic mixed lineage kinase domain-like protein serve as potential prognostic biomarkers in patients with resected pancreatic adenocarcinomas

Pancreatic cancer is one of the deadliest cancers with a very poor prognosis. Recently, there has been a significant increase in research directed towards identifying potential biomarkers that can be used to diagnose and provide prognostic information for pancreatic cancer. These markers can be used clinically to optimize and personalize therapy for individual patients. In this review, we focused on 3 biomarkers involved in the DNA damage response pathway and the necroptosis pathway: Chromodomainhelicase-DNA binding protein 5, chromodomain-helicaseDNA binding protein 7, and mixed lineage kinase domain-like protein. The aim of this article is to review present literature provided for these biomarkers and current studies in which their effectiveness as prognostic biomarkers are analyzed in order to determine their future use as biomarkers in clinical medicine. Based on the data presented, these biomarkers warrant further investigation,and should be validated in future studies.

ANALYSIS OF AUGMENTED THREE-FIELD MACRO-HYBRID MIXED FINITE ELEMENT SCHEMES

On the basis of composition duality principles, augmented three-field macrohybrid mixed variational problems and finite element schemes are analyzed. The compatibility condition adopted here, for compositional dualization, is the coupling operator surjectivity, property that expresses in a general operator sense the Ladysenskaja-Babulka-Brezzi inf-sup condition. Variational macro-hybridization is performed under the assumption of decomposable primal and dual spaces relative to nonoverlapping domain decompositions. Then, through compositional dualization macro-hybrid mixed problems are obtained, with internal boundary dual traces as Lagrange multipliers. Also, ＂mass＂ preconditioned aug- mentation of three-field formulations are derived, stabilizing macro-hybrid mixed finite element schemes and rendering possible speed up of rates of convergence. Dual mixed incompressible Darcy flow problems illustrate the theory throughout the paper.

More general results on mixed extreme value distributions

The sequences {Zi,n, 1≤i≤n}, n≥1 are multi-nomial distribution among i.i.d, random variables {X1,i, i≥1}, {X2,i, i≥1 } {Xm,i, i≥1 }. The extreme value distribution Gz（x） of this particular triangular array of i.i,d, random variables Z1,n, Z2 n,...,Zn,n is discussed. A new type of not max-stable extreme value distributions which are Fréchet mixture, Gumbel mixture and Weibull mixture has been found if Fj,…… Fm belong to the same MDA. Whether mixtures of different types of extreme value distributions exist or not and the more general case are discussed in this paper. We found that Gz（x） does not exist as mixture forms of the different types of extreme value distributions after we investigated all cases.

Rapid identification of time-frequency domain gravitational wave signals from binary black holes using deep learning

Recent developments in deep learning techniques have provided alternative and complementary approaches to the traditional matched-filtering methods for identifying gravitational wave(GW)signals.The rapid and accurate identification of GW signals is crucial to the advancement of GW physics and multi-messenger astronomy,particularly considering the upcoming fourth and fifth observing runs of LIGO-Virgo-KAGRA.In this study,we used the 2D U-Net algorithm to identify time-frequency domain GW signals from stellar-mass binary black hole(BBH)mergers.We simulated BBH mergers with component masses ranging from 7 to 50 M_(⊙)and accounted for the LIGO detector noise.We found that the GW events in the first and second observation runs could all be clearly and rapidly identified.For the third observing run,approximately 80% of the GW events could be identified.In contrast to traditional convolutional neural networks,the U-Net algorithm can output time-frequency domain signal images corresponding to probabilities,providing a more intuitive analysis.In conclusion,the U-Net algorithm can rapidly identify the time-frequency domain GW signals from BBH mergers.

Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients

The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.

Extreme value distributions of mixing two sequences with different MDA＇s

Suppose {Xi, i≥1} and {Yi, i≥1} are two independent sequences with distribution functions FX(x) and FY(x), respectively. Zi is the combination of Xi and Yi with a probability pn for each i with 1≤i≤n. The extreme value distribution ,n GZ(x) of this particular triangular array of the i.i.d. random variables Z1, , Z2, ,…, Zn n n ,nis discussed. We found a new form of the extreme value distribution ΛA(ρx)Λ(x)(0<ρ <1), which is not max-stable. It occurs if FX(x) and FY(x) belong to the same MDA(Λ). GZ(x) does not exist as mixture forms of the different types of extreme value distributions.

关于混合Domain的几个问题

研究了Lawson紧的连续domainD的混合幂domainM(D),给出了M(D)的极大元的刻画.证明了M(D)是其极大点空间Max(M(D))的domain环境,同时证明了FB domain关于混合幂domain构造是封闭的.

Extreme value distributions of mixing two sequences with the same MDA

Suppose {Xi, i 1} and {Yi, i 1} are two independent sequences with distribution functions ()XFx and ()YFx, respectively. Zi,n is the combination of Xi and Yi with a probability np for each i with 1 in. The extreme value distribution GZ(x) of this particular triangular array of the i.i.d. random variables Z1,n, Z2,n, ,L Zn,n is discussed. We found a new form of the extreme value distributions i) 12()()AxxaaFF and ii) 12()()AxxaaYY (a1<a2), which are not max-stable. It occurs if FX and FY belong to the same MDA(? or MDA(?.

A weak invariance principle for self-normalized products of sums of mixing sequences

Let variables in the {X, Xn, n ≥ 1} be a sequence of strictly stationary φ-mixing positive random domain of attraction of the normal law. Under some suitable conditions the principle for self-normalized products of partial sums is obtained.

基于主动尾流控制的风电机群协同优化调度

针对机组间尾流效应严重影响风电机组发电效率的问题,提出了风电机组安全偏航约束计算方法、尾流特性混合半机理建模方法以及风电机群多目标协同优化调度方法。基于FAST.FARM平台完善了多自由度可控机组与尾流的动态交互集成仿真环境,对比分析了2台机组串列式排布以及华东地区某海上风电场7台机组实际排布下的协同运行优化性能。结果表明:所建立的集成仿真模型能够合理表征风电机群与空气流场的多领域动态交互特性,所提方法能够有效提升风电机群发电效能,促进经济效益、资源利用和成本控制的均衡优化。

外周血受体相互作用蛋白激酶3、混合系列蛋白激酶样结构域水平与新生儿坏死性小肠结肠炎病情严重程度的关系

目的分析新生儿坏死性小肠结肠炎(NEC)患儿外周血受体相互作用蛋白激酶3(RIPK3)、混合系列蛋白激酶样结构域(MLKL)的表达情况及其与病情严重程度的关系。方法选取92例NEC患儿纳入NEC组,并根据病情严重程度进一步分为轻度NEC组(Ⅰ级)60例和重度NEC组(Ⅱ~Ⅲ级)32例,另选取同期诊治的60例腹股沟斜疝患儿纳入对照组。采用实时荧光定量聚合酶链反应检测外周血RIPK3 mRNA、MLKL mRNA表达;采用Pearson相关分析法明确NEC组外周血RIPK3 mRNA与MLKL mRNA表达的相关性;采用免疫印迹法检测NEC回肠组织和正常回肠组织中RIPK3、MLKL蛋白表达;采用多因素Logistic回归分析明确重度NEC发生的独立危险因素。绘制受试者工作特征曲线,分析外周血RIPK3 mRNA、MLKL mRNA单独及联合预测重度NEC的价值。结果NEC组外周血RIPK3 mRNA、MLKL mRNA相对表达量分别为(2.41±0.52)、(3.03±0.64),高于对照组的(1.02±0.21)、(0.93±0.20),差异有统计学意义(P<0.001)。NEC回肠组织中RIPK3、MLKL蛋白相对灰度值分别为(1.20±0.21)、(1.13±0.24),高于正常回肠组织的(0.34±0.12)、(0.32±0.11),差异有统计学意义(P<0.05)。NEC组患儿外周血RIPK3 mRNA与MLKL mRNA相对表达量呈正相关(r=0.623,P<0.001)。重度NEC组合并气腹征、多器官功能障碍综合征、败血症者占比和RIPK3 mRNA、MLKL mRNA相对表达量均高于轻度NEC组,差异有统计学意义(P<0.05);RIPK3 mRNA、MLKL mRNA相对表达量升高是重度NEC发生的独立危险因素(P<0.05)。外周血RIPK3 mRNA、MLKL mRNA联合预测重度NEC的曲线下面积大于RIPK3 mRNA、MLKL mRNA单独预测(Z=4.127、4.261,P<0.05)。结论RIPK3 mRNA、MLKL mRNA在NEC患儿外周血中表达升高,两者均与NEC病情严重程度有关,且两者联合检测对重度NEC具有较高的预测价值。

任意起伏地形下重力异常三维正演及并行计算

为了进一步提高空间-波数域三维重力异常正演算法的适用范围和计算效率,本文采用任意傅里叶变换算法实现了空间-波数域三维重力异常正演,且在NVIDIA CUDA平台上进行CPU-GPU并行加速.任意傅里叶变换算法的基本思想是将二维傅里叶变换转化为两个一维傅里叶变换,一维傅里叶变换积分离散为多个单元积分累加和,离散单元中原函数采用二次插值形函数拟合,求出单元积分的解析表达式.相比现有的傅里叶变换算法,新方法具有采样灵活、积分精度高、计算速度快和傅里叶变换的截断效应小等优势.利用空间-波数域算法的高度并行性,采用CPU并行求解常微分方程,GPU并行计算任意傅里叶变换,实现了CPU-GPU并行加速方案,进一步提升了本文算法效率.利用常密度模型,对比数值解和解析解,结果表明本文算法正确;利用变密度模型对比了任意傅里叶变换算法与高斯快速傅里叶变换算法的计算效率与精度,在相近的数值精度下,本文算法波数选取少,效率高;测试CPU-GPU并行效果,结果表明相比CPU串行算法,CPU-GPU并行算法的计算效率大大提升,千万数量级节点数模型正演仅耗时数秒.最后利用实际地形数据进行三维重力异常场数值模拟,证明了新方法的高效性与实用性,对实现大规模复杂条件下重力异常精细化反演成像与综合解释有重要意义.

基于界面电荷极化特性的变压器油纸绝缘寿命预测方法研究

准确预测变压器油纸绝缘寿命对于保障电力系统安全稳定运行具有重要意义。根据微观介电响应机制剖析界面极化与油纸绝缘老化的内在机理联系,并基于混合极化模型——一种在扩展德拜模型基础上引入界面极化支路的等效电路拓扑,探究界面电荷极化特性在油纸绝缘寿命预测中的应用。首先,分析温升过程中频谱特性变化背后蕴含的普适弛豫极化规律,基于混合极化模型重构回复电压极化谱。结合温度对主时间常数和模型参数的影响规律,提出一种单一温度测试下的油纸绝缘活化能计算方法,通过界面极化电容构建出频温平移因子。其次,剖析介质损耗因数随老化加深的变化规律,探究损耗峰偏移现象与界面极化弛豫时间的内在联系,凝练全新表征油纸绝缘聚合度含量的频域特征参量——界面极化支路极点Ph2。最后,联合界面极化支路极点Ph2、频温平移因子和损失累积动力学方程构建油纸绝缘寿命预测模型,并将不同运行年限变压器现场实测数据代入模型验证实效性,为油纸绝缘寿命预测提供了一种新思路。