地震时间序列多重分形计算精度估计

目前计算地震时间分形的方法很多,但具体各种计算方法的误差（或精度）是多少,还没有一个定量的结果。本文针对这一问题,利用Ｃａｎｔｏｒ多分形集理论模型数值模拟来估算其精度（误差）,进而解决了限、离散、点集的多重分形维数谱的计算精度（误差）估计问题。通过一维（时间）理论模型数值模拟得出如下结论：（１）存在尺度问题、样本容量问题、基准点选取问题、不同计算方法差异等;（２）固定质量（ＭＡＳ）法和最小生成树（     

地震多重分形计算的最小生成树法

目前地震分形计算的方法很多,但各方法或多或少都存在一些问题．本文引入一种地震分形计算的新方法——最小生成树法（minimal spanning tree）,对其算法作了重要的修改,并通过对乌什、乌苏、唐山和海城4个区域震例的研究,得出大震前后地震时空分布多重分形谱由集中变疏张,反映了分形结构复杂性的增加、地震时空分布不均匀性的增加．由点集的数值模拟表明:强震前后多重分形谱变化的物理本质是地震分布由均一到不均一,由简单到复杂的变化过程．      

基于显式垂直比例因子的显式分形插值地震数据重建

在前人工作的基础上对分形插值方法作了详细的探讨,给出了分形插值函数的显式表达方式.在量纲分析的基础上给出了垂直比例因子的局部显式表达式,旨在提高地震道插值重建的精度及突出局部信息,并从单道地震图的角度分析其在地震道插值重建中的应用效果.研究了垂直比例因子的变化对分形插值精度的影响.数值实验表明,随着垂直比例因子的增大,分形垂直的误差逐渐增大,二者之间呈显出指数增长的趋势.该法克服了随机分形插值方法必须进行多步迭代的弱点,提高了计算效率.通过对理论地震道插值重建的分析,说明了本文分形插值方法的高精度和高效率.本文提出的显式分形插值方法既能够突出地震道数据的局部信息,又较好地保持了地震道数据的总体变化趋势.     

地震能量空间分布局域标度特性的初步研究

本文利用计算机数值模拟数据及大同地区三次中强地震前后的地震记录资料，在空间能量分布“漫化”的基础上，通过多重分形方法着重对能量空间分布的局域标度性质及可能的震兆信息进行研究。结果表明能量分布普遍地表现出多标度分形特性，但相同ｑ值Ｄｑ的时间变化或Ｄｑ￣ｑ谱形态变化不明显。中强地震前标度指数谱ｆ（α）￣α曲线由两端大体对称或右端略低变为左端明显下降，同时α及ｆ（α）数值分布范围扩展、顶点右移。计算机模     

分形插值地震数据重建方法研究

对分形插值方法作了较详细的探讨,给出了分形插值函数的显式表达方式,同时给出了垂直比例因子的局部显式表达式,旨在提高地震道插值重建的精度及突出局部信息,并从单道地震图的角度分析其在地震道插值重建中的应用效果.利用该方法对理论模型和济阳坳陷实际地震台站资料进行了重建处理,结果表明,分形插值重建的地震道是原始地震道的良好近似,缺失道的振幅和相位都得到了很好的恢复.该法克服了随机分形插值方法必须进行多步迭代的弱点,提高了计算效率.通过对单道地震图插值重建结果的分析,说明了本文分形插值方法具有较高的精度和较高的效率,有深入研究的潜力.本文提出的显式分形插值方法既能够突出地震道数据的局部信息,又较好地保持了地震道数据的总体变化趋势.     

自适应时间步长法及三维堤坝地震液化数值模拟

三维大模型数值计算因巨大的单元和结点数目而非常耗时,在地震响应分析中受计算时间步长的限值则更加耗时。在饱和砂土动力液化计算平台上开发时域离散误差评估方法和时间步长自适应调整的计算程序,并成功应用于三维堤坝地震液化响应分析。时域离散误差包括土骨架的位移误差和单元孔压误差,通过定义孔压误差影响系数计算出混合误差,根据混合误差和设定的误差允许值进行计算步长的自适应调整。在三维堤坝地震液化数值模拟中,采用自适应时间步长法有效避免小步长精确但耗时、大步长省时而不精确的缺点。在大模型和超大模型计算中,最优调整每一步的计算时间步长,完美实现既节省时间又不失精度的时域离散策略。     

基于改进反应位移法的盾构隧道横向抗震规律研究

在地震动作用下,小间距盾构隧道夹土层区域的动力反应影响结构的内力与变形,而反应位移法未考虑夹土层弹簧参数的变化,本文针对这一误差来源,提出了改进反应位移法。一方面在反应位移法合理内容的基础上完善计算模型,另一方面通过数值计算将改进反应位移法与动力时程法进行对比分析,并拟合理论模型对改进的方法进行修正。结果表明,两种算法下结构的地震响应规律一致;间距d≥8m时,结构地震动反应变化不明显,可按反应位移法计算;理论模型提高了改进反应位移法计算结果的精度。     

基于半精度浮点数优化与OpenMP的地震数值模拟

研究地震波场的传播特征是地震勘探的基础,利用有限差分法求解波动方程进行地震正演,其优势在于占用内存低、计算速度快、易于实现。但在三维情况下,计算量和数据量会急剧增加,传统串行地震数值模拟将无法满足计算效率的需求。因此,本文提出一种基于半精度浮点数优化与OpenMP的三维波动方程地震数值模拟方法,该方法首先利用半精度浮点数对地震常用的浮点型数据进行优化;其次利用应用程序接口OpenMP在多核CPU下通过以分割波场计算区域的方式实现并行计算;在保证计算结果满足精度需求的同时,能有效提高三维地震数值模拟的计算效率,并减少近一半的内存需求。通过数值试验证明该方法的有效性和实用性。     

复杂介质地震波传播的褶积微分算子数值模拟h

将佛尔塞（Forsyte）广义正交多项式微分算子地震波场正演模拟算法运用于复杂非均匀介质模型的波场数值模拟中,并比较了该方法在计算效率和计算精度方面与有限差分方法和伪谱法的差异. 数值结果表明,这种广义正交多项式微分算子法计算速度快、精度高,对计算资源需求低,是一种颇具潜力的数值模拟方法.      

地震序列自仿射分形及自相似分形的计算与比较

本用标准偏差法和标度变换法，计算了云南省五个地震区带上１９６５－１９９５年间地震序列的自仿射分形的分维值Ｄ△ｔ（地震时间间隔序列自仿射分形）ＤＭ－ｔ（地震震级一时间序列自仿射分形）及自相似分形的分维值Ｄ０（地震时间分布自相似分形）得到了相应的时间变化曲线，ＤＭ－ｔ～ｔ，Ｄ△ｔ～ｔ和Ｄ０～ｔ。通过相互比较发现，其中ＤＭ－ｔ的随时间变化在中强地震前的反应最为敏感。     

The Precision Estimate of Multi-Fractal Methods of Time Sequence of Earthquakes

At present,there are many methods of calculating seismic time-fractal.However,there isn’t aquantitative result about the precision of every method.So,in this Paper,we use the digitalimitation of theoretic model to solve precision estimate problems of calculating the precision ofone dimension distribution of theoretic models with Cantor multi-fractal set,we obtained someresults as follows:(1)There exists many problems such as rules,numbers of samples,basicpoint selection,the diffence resulted from different methods and so on.(2)The fixed-massmethod(MAS)and the minimal spanning tree method(MST)can give good structure charac-teristics with different q value,while the counting-boxes method can’t.And the error of thefixed-radius method(RAD)in the range of-q is too big.(3)There are scale problems of rulesfor multi-fractal,it is objective reflection for non-rule area.(4)MST has the boundary prob-lem,while MAS and RAD don’t.(5)With increasing sample number,the precision of allfractal-dimension values becomes     

基于空间-光谱特征和稀疏表达的高光谱图像分类算法（英文）

针对传统的高光谱数据分类方法分类精度不高、没有充分地利用空间信息等缺陷,提出一种基于Gabor空间纹理特征（Gabor spatial texture features）及无参数加权光谱特征（Nonparametric weighted spectral features）和稀疏表达分类（Sparse representation classification）的高光谱图像分类算法,可以简写为Gabor-NW SF和SRC,即GNWSF-SRC。所提出的GNWSF-SRC分类方法首先通过融合高光谱的Gabor空间特征和无参数加权光谱特征来更好地描述高光谱图像,然后通过其进行稀疏表达,最终通过对比其重构误差获得分类结果。在训练集比例不同的情况下,用所提出的方法对两组典型的高光谱数据进行处理,理论研究和仿真结果表明：与传统的分类方法相比,所提出算法能够提高分类精度、Kappa系数等,取得了较好的分类效果。     

模拟震源动力学破裂过程的数值方法研究

为模拟震源的动力学破裂过程,本文讨论一种求解剪切裂纹动力学扩展问题的有限差分法。在研究二维反平面破裂的基础上,我们把所用方法推广到三维问题。研究了许多简单震源模型,通过对比同一问题的解析解和其它数值解,对方法的正确性和计算精度进行了检验,结果表明我们的数值方法是可行的。最后,作为初步应用的例子,我们研究了非均匀断层的自发破裂问题。     

南天山西段数字化形变观测的潮汐参数分析

利用Venedikov调和分析方法,计算了南天山西段8套数字化形变观测资料的潮汐参数：O1波和M2波的潮汐因子、中误差、观测振幅以及M2波的相位滞后与中误差,分析并评价了资料的精度和稳定性,为今后的资料使用和异常判定提供科学的参考依据。结果分析表明：（1）阿合奇、哈拉峻和乌什3个台地倾斜资料、乌什洞体应变仪资料和乌什重力仪资料潮汐分析结果较好,中误差小,精度均达到了中国地震局形变I类台的标准;（2）马场台井下倾斜仪EW向γM2偏小,受格值或者仪器常数影响较大;（3）喀什栏杆乡和乌恰井下倾斜仪两分向的定位存在较大偏差。     

基于条分模式的边坡可靠度近似计算方法

通过对边坡稳定分析方法中的条分理论和响应面法的研究,针对边坡可靠性计算往往没有明确的解析表达式,以及稳定性系数计算方法和响应面法（RSM）的特点,将响应面法中的有限元数值模拟以条分模式中的稳定性系数隐式方程的迭代计算方法代替,建立了条分模式下的边坡可靠性计算的极限状态方程,从而形成了一种新的边坡稳定可靠性响应面分析方法。本文提出的改进的响应面法原理简单,计算效率较高并具有一定的精度,适用于对边坡可靠度的近似计算。     

SNREI地球模型负荷勒夫数数值计算的新进展

基于Ｌｏｎｇｍａｎ和Ｆａｒｒｅｌｌ的表面负荷形变理论，深入研究计算ＳＮＲＥＩ弹性地球模型表面负荷勒夫数的数值计算方法，取得新的进展：采用了一种新的代换方法─ｒ￣ｎ因子法，比传统的无量纲化法简单易行，且能达到很高的计算精度；确立了检验表面边界条件、微分方程积分线性性和计算结果的收敛性，用于评价数值计算结果精度的有效准则；给出了ＰＲＥＭ地球模型的一组高精度计算结果．     

怀来地震台伸缩仪EW向改进试验

对怀来台伸缩仪EW向基线固定端、传感器进行布设改进试验,将原有28m长基线改为两段短基线（东段：14．6m,西段：13．4m）。改进后,两条短基线取得了较好的观测成果,其中东段基线改造效果非常明显,观测曲线干扰减少,光滑度显著提高,观测数据精度明显优于改造前的观测数据。这次改造试验解决了困扰怀来台多年的难题。     

信号的瞬时参数与正交小波基

对于可表成ｕ（ｔ）＝ａ（ｔ）ｃｏｓθ（ｔ）的信号,称ａ（ｔ）为其瞬时振幅,θ（ｔ）为瞬时位相,为瞬时频率．在许多实际问题中,瞬时振幅、位相、频率等瞬时参数的提取有重要意义．但已知的一些瞬时参数计算方法均各有其局限性．本文构造了一种频谱支集紧且有解析表达式的小波,利用它得到一种计算信号瞬时参数的能抗噪声干扰的快速、局部算法．理论分析和数值模拟的结果表明,该方法是实用而有效的,可望能满足许多实际问题的需要．     

基于OpenSees的空心钢管混凝土柱恢复力计算公式研究

采用OpenSees建立基于纤维截面的空心钢管混凝土柱有限元模型,对空心钢管混凝土柱的滞回性能进行计算。在该有限元模型基础上,考虑轴压比、长细比、截面空心率、截面含钢率、混凝土强度和钢材强度等因数对骨架曲线恢复力模型的影响,结合正交试验设计方法,对骨架曲线恢复力计算公式进行回归分析。研究结果表明：1数值计算结果与试验结果吻合良好;2得到的空心钢管混凝土柱骨架曲线恢复力计算公式在一定设计参数范围内具有较高的计算精度。     

Modifying the intensity distribution by assessing the reliability

This article presents an application of a procedure to modify the intensity distribution by assessing the reliability. There are two potential possibilities that may influence the intensity distribution: (1) For the interpolation error, we generate a measured grid across the calculation region. When the point to station spacing is <5 km, we consider the results precise; however, some points have less precision because these are farther from the corresponding stations. When the spacing is between 5 and 50 km, we consider the results imprecise and define a reliability factor that correlates with the distance. (2) Some records may have errors that result from local site conditions, equipment problems, or some disturbance such as lightning stroke, which will lead to some grid points having an incorrect intensity. We regress the attenuation relation for sites with abnormal intensities and consider the results to be accurate when the standard deviation (STD) is <σ and inaccurate when the STD is > 2σ. We then define a reliability factor to correlate with STD between σ and 2σ, such that the intensity distribution is in accord with both wave propagation theory and the investigation intensity.