《地震频度-震级关系的时空扫描》

第一章引论1

1 频度—震级关系1

2 b值的物理意义2

3 频度—震级关系的唯象推导4

第二章b值统计分析方法7

1 基本概念7

2最大似然法8

2.1 b的估值8

2.2 最大似然估值的性质10

2.3 误差估计13

3线性最小二乘拟合14

3.1 a和b的估值14

3.2 线性最小二乘估值的性质和误差估计17

4 指数分布的非线性最小二乘拟合23

5 时空扫描的具体作法26

6数据整理28

6.1 关于震级问题的讨论28

6.2 数据的归纳整理31

第三章京津唐渤张地区时空扫描结果37

1唐山、宁河、卢龙地区b值时间扫描37

1.1 唐山震源区b值时间扫描38

1.2 震源附近不同区域的比较41

2全区地震个数N值和b值的时空扫描43

2.1 全区地震个数N值和b值的时间扫描43

2.2 全区地震个数N值的空间扫描43

2.3 全区b值空间扫描46

2.4 全区ML≥4.0的地震震中分布47

2.5 b值的分块时间扫描48

3 其它震例的比较50

4 对唐山地震孕育过程的一点认识50

第四章统计检验和统计方法评述53

1不同方法所得结果的比较53

1.1 结果的一致性53

1.2 误差估计54

1.3 个别大震的影响55

1.4 起算震级和分档方式的影响57

2统计检验58

2.1 相关系数检验59

2.2 x2检验61

2.3 柯尔莫哥洛夫-斯米尔诺夫检验62

2.4 F检验65

2.5 唐山震例统计检验结果66

3 统计分析方法评述68

第五章讨论71

1 b值分析的历史回顾71

2 主要结论和有待进一步研究的问题72

附录A程序和数据的组织74

1 地震目录的穿孔74

2 目录的合并和挑选75

3 固定区域的时间和个数扫描75

4 固定时期的空间扫描76

附录B主要程序模块78

（1）加工原始地震目录的程序DATPC78

（2）合并地震目录的程序LSAP79

（3）子程序LSIN（K）80

（4）固定区域的N值、b值时间和个数扫描程序TNFR81

（5）按凸多边形选地震目录子程序SOA82

（6）判断震中是否在凸多边形内的子程序SOBLK（X，Y，KEY）83

（7）时间或个数扫描控制子程序SOBL（KEY，LFG）84

（8）拟合控制子程序FIT（KB，KE，IY，KERR）85

（9）最大似然法子程序MAXLH（N，X，IY，TN，B，DB）87

（10）线性最小二乘拟合子程序LSQR（KEY，N，X，Y，W，A，B，DB，R，KERR）87

（11）非线性最小二乘拟合子程序TNG2V（N，X，Y，W，A，B，DB，R，EPS，KMAX，LAG，KEY，KERR）88

（12）x2和K-S检验子程序CHECK（N，X，Y，W，B，FKAI，EMAX）90

（13）打印扫描区地震简目的子程序SOC91

（14）儒略日转换子程序DATE（NN，KM，KD）91

（15）固定区域数据块TID92

（16）固定时期N值、b值空间扫描主程序FPSC92

（17）扫描子程序SCANP（IY，IX，IM）94

（18）子程序SOQUD（D，XC，YC）95

（19）空间扫描拟合子程序REGRN95

附录C唐山震源区地震简目97

参考资料101