seisjun

各种方法合成地震图的对应单位

摘要：最近尝试了各种理论地震图的合成方法，主要分为以下三种： fk:结合syn一般来说合成的理论地震图为速度，单位是cm/s;其中需要注意的是syn合成理论地震图时的duration_time和rise_time的设置。 tel3:和fk合成理论地震图的单位类似，也是cm/s;但目前还是有一些问题；后面阅读全文

posted @ 2018-01-07 15:21 seisjun 阅读(936) 评论(0) 推荐(2)

matla互相关协方差的计算和理解

摘要：其中不太明白的是为什么有一个偏移，因为互协方差的定义为： Rxy(m)=E[x(n)y(n+m)]=lim1/N(∑x(n)y(m+n)); Cxy(m)=E{[x(n)-μx][y(n+m)-μy]}=Rxy-μxμy; 阅读全文

posted @ 2018-01-02 16:27 seisjun 阅读(1302) 评论(0) 推荐(0)

matlab运行中出现“Caught "std::exception" Exception message is: Message Catalog MATLAB:builtins was not loaded from the file."

摘要：在我运行过程中，经常爆出这一不确定是什么的问题，经排查后发现，原来是fopen 文件后，没有及时fclose导致的。阅读全文

posted @ 2018-01-01 11:59 seisjun 阅读(6087) 评论(0) 推荐(0)

拉普拉斯算子球面坐标式的推导

摘要：阅读全文

posted @ 2017-12-27 21:41 seisjun 阅读(4510) 评论(0) 推荐(0)

对拉普拉斯平滑的认识

摘要：拉普拉斯平滑是一个二阶微分的过程，简单地用图形解释为其实，要把其应用到实际中，其实就是求出的不同位置的S值的权重因子不同的问题；这样，我们来推： (S(i,j)-S(i-1,j)-(S(i+1,j)-S(i,j)))/L是横着的二阶微分形式； (S(i,j+1)-S(i,j)-(S(i,j+1)- 阅读全文

posted @ 2017-12-25 10:33 seisjun 阅读(1489) 评论(0) 推荐(0)

关于小波变换

摘要：小波就像是不同尺度的放大器，检测这一频率的波的突变点，因为它可以反映这一频率波的幅值的变化。所以小波可以反映在某一时间窗内某一频率波的能量变化。阅读全文

posted @ 2017-12-15 21:47 seisjun 阅读(195) 评论(0) 推荐(0)

随笔学习

摘要：目前在做破裂反演，用的远震数据是频率比较低的数据；所以它对模型的分辨率不高，但近震对模型的分辨率很高，近震拟合了绝对位置和到时，对绝对位置有很好的分辨。震级越低，越往高频带移动。还有基于小波基的破裂反演方法，将地震波作小波变换，反演每个基的系数。但这个小波系数反映的是“该频率所占的比重”。阅读全文

posted @ 2017-12-15 20:58 seisjun 阅读(165) 评论(0) 推荐(1)

关于滤波频段的选取（近震波形）

摘要：讨论两个问题：破裂过程用于反演的两个频段：一是近震波形常用的频段范围，二是远震频段常用的范围；做破裂过程，必须满足波长远远大于子断层尺度，台站震中距远远大于子断层尺度，这两个条件缺一不可。例如，假设子断层尺度为0.1°*0.1°，高频部分的截止频率为0.1HZ. 有一个经验的关系如下：而低频阅读全文

posted @ 2017-12-08 09:03 seisjun 阅读(688) 评论(0) 推荐(0)

pssac plot

摘要： for multi-waveforms with different colors: R-Xmin/-Xmax/Ymin/Ymax -Ba/b, a: x delta; b:y delta -MYmax-Ymin for normalize y values *.z X Y plot x posit 阅读全文

posted @ 2017-12-05 20:48 seisjun 阅读(274) 评论(0) 推荐(0)

关于单层地壳介质模型和多层介质模型中计算地震波震相及走时

摘要：首先计算的具体原理要参考万老师书《地震学导论》里第六章6.7中6.7.1和6.7.2中的内容。在此可以总结为：对于单层地壳介质模型，由于离源角可以通过δ=tan-1(震中距/震源深度)直接求出;因此走时相对好计算。而对于多层地壳介质模型，需要将地球分为n层，而假设震源在第j层，震源到第j层顶厚度阅读全文

posted @ 2017-11-24 11:28 seisjun 阅读(829) 评论(0) 推荐(0)