第一章 引论 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）加工原始地震目录的程序DATPC 78
（2）合并地震目录的程序LSAP 79
（3）子程序LSIN（K） 80
（4）固定区域的N值、b值时间和个数扫描程序TNFR 81
（5）按凸多边形选地震目录子程序SOA 82
（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）打印扫描区地震简目的子程序SOC 91
（14）儒略日转换子程序DATE（NN，KM，KD） 91
（15）固定区域数据块TID 92
（16）固定时期N值、b值空间扫描主程序FPSC 92
（17）扫描子程序SCANP（IY，IX，IM） 94
（18）子程序SOQUD（D，XC，YC） 95
（19）空间扫描拟合子程序REGRN 95
附录C 唐山震源区地震简目 97
参考资料 101