俗话说，人无远虑必有近忧。天灾如洪水、地震造成的灾难和损失不可估量，如果能够提前预知，做些防范就能尽量减轻灾害。
天灾是自然灾害，既然是自然灾害，那就把研究对象锁定在自然上。自然中对地球影响最大的就是太阳和月亮。
自然灾害是灾害应力积累到某个极限，最后能量释放的结果。要积累就会有一个运动过程，过程就意味着时间，有时间就有周期。
万事万物逃不过数，只要运动发展的事物，就有定数和周期。
定数包括事物发生的时间，变化的时间，消失的时间。周期，事有往复循环，因为自然中也不过就那么几件事，发生的还会再次发生，重复往来。洪水、地震就是几千年来重复发生的事情，重复发生的事情必然有数，有定数，有周期在内。
中国术数早已经验证了事情是可以预测的，有象必有数，有数就能测。不管是定数还是周期，都受到自然界宏观运动的支配。既然自然中对地球影响最大的就是太阳和月亮。而干支就包含太阳和月亮运动的周期规律，那么不排出天灾的周期与干支有关，因此可以运用干支方法作为天灾预测的工具之一。
1、古代天文预测的两个范例
中国古代天文很早就注意到木星的恒星周期差不多等于12年（实际是11.86年），人们猜想，它和12支以及回归年的丨2.37个朔望月有联系。
《越绝• 计倪内经》一书在介绍公元前4世纪南方的博物知识时，有下列一段话：
“太阴三岁处金则穰，三岁处水则毁，三岁处木则康，三岁处火则旱。故散有时，积有时。领则决万物， 不 过 三 岁 而 发 矣 .以智论之，以决断之。以道佐之，断长续短。一岁再倍，其次一倍，其次而反。水则资车，旱则资舟，物之理也。天下六岁一穰，六岁一康，凡十二岁一饥，是以民相离也。故圣人早知天地之反，为之预备”
大意是：
最初三年，太阴的位置在金（西的方向，各地都得到丰收；当在水（北）的方向时，就有三年歉收。当它在木（东）的方向时，有三年富足，当它在火(南）的方向时，则有三年旱灾。因此，有时适于囤积农产品，有时却要把米粮分散出去。囤积品不必超过三年的需要，只要明智的考虑并进行决断，人们就可以靠自然界规律的帮助，用盈余来补救不足，第一年可以有两倍的丰产，第二年娃正常收成，第三年歉收。水灾时期就应该想到制造车子，干旱时期也要想到准备舟船。天下每六年有一次大丰收，每六年有一次正常收入，十二年有一次饥荒。所以 圣人既能顶见自然界的反复，也能对未来的灾变早作准备。
这 真 是 一 段 绝 妙 文 ， 它 不 仅 为 我 们 传 递 了 祖先 的十 二 支 与 天 气 灾 样 关 系 的 信 息 .更 可 贵 的 是 向 我 们 提 出了 天灾是可以预知的 忠 告 。时 隔 数 千 年， 今 口 读 来， 犹似 昨 天 一 样 的 殷 切 叮 嘱。这 是 在 古 代 处 于 雏 科 学 阶 段 的事 ，斗转星移 ， 进 人 现 代 科 学 的 2 1 世 纪 ，我们更应该继续 继 承 、维 护 、发 展 、弘 扬祖先的文化。
从 另 一 份 资 料 获 悉 ， 我 国 最 早 预测 天 气 的 文 献 是 :
湖 南 一 带 民 间 流 行 的 《 娄景书 》 ，成 书 时 间 大 槪 在 公元 前 2 0 6 年 （ 2 0 6B.C., 汉 高 祖 刘 邦 元 年 ） 前 后 ， 内 容以 预 測 农 业 气 象 为 主 。 这 本 书 历 经 2 0 0 0 多 年 的 考 验 , 仍 存 生 命 力 ， 说 明 其 有 – 定 的 客 观 根 据 。但 自 然 现 象 ,不 可 能 对 任 何 地 区 都 有 严 格 而 简 单 的 单 一 周 期 。
本 书下 面 论 述 的 目 的 ， 就 是 要 探 索 干 支 6 0 年 和 6 0 周 期 的根据和事实。顺 便 ， 以 近 儿 个 长 江 水 系 洪 水 年 份 ， 与 2 0 0 0 多 年 前《娄景书》的预測作一比较（见图)，可见该书有一定的可信度。
本文中有时引用了中国古代延用下来的“相位”名称来描写60 (年或日）周期的相位，但不考虑它们的含 义 。
近年据有的研究者调查：湖南安乡县气象站曾访问过数百位有看天经验的老农民、老船民，发现他们往往 用10年、60年周期变化来预测水旱趋势。
1968年有许多老农民说：“明年是乙西年，老乙酉年（公元1849 年）大水，前乙西年（公元1909年）也大水，明年又遇上了 60年大水周期，该站根据民间经验较正确地估计预报了 年的大水。查阅《娄景书》，果然有“乙酉年来雨水倾，夏秋流郎（即指雨水）略有增”之 句。
由 此 可 见 ， 早 在秦汉时代，就有运用干支周期进行 天时预测的确实纪录。虽然《娄镜书》说的是湖南一隅 的天时变化，难免其有局部性，但是它的原则却有普遍 意义。
2、日干支预测
日干 支 周 期 即 从 甲 子 到 下 一 个 甲 子 的 6 0 日 周 期 ^。
根据我 们的研 究 ， 对 于地震来说， 有 两 种 运 动 规 律 ， 可以导出两个经验公式。
(1) 日干支第一式
中国有几次地震时间有一个特殊的60 日周期序列，它的经验公式如下：
y = 1923.2269+0.1642746×í
式 中 ， y是 这 个 特 殊 时 间 序 列 内 ， 在 中 国 大 陆 发 生 的 7
级 或 更大地震的 时 间 ：1923为 年， 0.2269为 10进制年内的天 数 （ 月、日），0.1642746为10进制年内 60天的数值；í是整数。
在 1 9 2 3 年 3 月 2 2 日 〜 1 9 8 5 年 8 月 2 5 日 期 间 ， 6 2年 内 共 发 生 7 级 或 更大地震 3 9 次 ， 现将其中 9 次根据上式计算的 y(í ) 值 和实际地震要素对比列干 表内。
1923年3月22日〜1985年8月25日期间计算地震与实际地震比较表
计算值是指用公式（1）推算的数值，实际要素是指实际发生地震的时间、地点、震级， x-y为误差资。
实际计算时可以每年365.2422日/ 12月=30.44日 (月平均日数，避免查对闰月有关的差异）计算。
如1923年3月22 日的10进制月，日为：
30.44 x 2 = 60.88 日（1、2两月日数）
60.88+22 (3 月份 22 号前的日数） =82.88 日
82.88/365.2422 = 0.2269
以上表中第 2 列í=8为例计算如下 ：
1923.2269+0.1642746 x 8
= 1923.2269+1.3142
= 1924.5411
=1924 年+0.5411
因为 0.5411 x 365.2422 = 197.63 天
197.63/30.44=6.49月
30.44×0.49=14.9，约为15日
故y= 1 9 2 4 年 7 月 1 5 日
上表中所 列 的 9 次 计 算地震时间与实际地地震时间基本相符。
计算过程是 ： 已 知 在 1 9 2 3 年 〜 1 9 8 5 年 6 2 年中发生3 9 次 7 级 或 更 大 的地震， 其 地 震 是 已 知 的 。
按公式 （ 1 ) 以i=0,1，2, 3 , … 380 代入式中计算， 求出y。
在这些计算结 果 中 有 9 次地震时间与实际地震时间x接 近 ，x-y离 差 值 范 围 在负 3 天 （ 即 提 前 3 天 ）〜+2天（即晚2天）区间内。
即-3、-2、-1、0、 + 1、+ 2 共 6 天 ， 以 符合 无 偏 条 件 ，有时可加1天则为 7
天，总地震次数 3 9 除 以 60 ，即 3 9 次 / 6 0 = 0 .6 5 次 ， 再乘以离差区间7 , 即 0 .6 5×7 = 4 .5 5 次 。
换句话说 ，假如地震均匀分 布， 则应发生符合于公式（1） 条件的地 震 4 .5 5 次 ， 而 实 际 上 是 9 次 .根 据 一 种 检 查 方 法，可以一定可信度否定了均匀分布假说。
从表11上看，还可以发现一个有趣而值得里重视的现象。那就是地震发生地点的纬度都在北纬31 °—39.2 °之间，其平均值为北纬35.1度，这个数值恰与地球体积不变而扁度变化的中和线吻合。
从以上的讨论中可以看到，1923年3月22 日在60日周期序列中的日干支是“甲午”，符合干支纳音歌诀“甲午乙未沙中金”，因而可把这种地震序列称为“沙中金序列”。
(2) 日干支第二式
另一个地震时间的经验公式，可写成下列形式：
1966.2396+0.164275×i ⑵
式中 ， y是地震时间 ， （i 是 整 数 ， 在 19 6 6 年 3 月2 2 日 〜 1976 年7月28 日的 1 0 年 时 期 内 ， 近 7 级 或更大地震的时间的计算值和实际值比较如表13。
用这一经验公式进行远程外推，可得1990年，北京4级或更大级别的时间。如下表。
注：远程外推是根据公式采用适当i 值以计算预测未来地震。
1966年3月22日在60日周期序列中的干支是“庚辰”符合干支讷音歌诀“庚辰辛巳白腊金”，可称为白腊金序列。
两个序列的地理分布：式（1）的沙中金序列和式(2)的白腊金序列分别代表两类序列，这两列地震序列的震中分布显然不同的.如下图12所示。
从图12中可以明显看出：式（1)沙中金序列地震的震中，用空心圆圈表示，它们分布在北纬35度两边；而式（2）白腊金序列地震的震中，用黑色方块表示，
它们分布在北北东方向（图的左上角，有一点除外).
以 上分布 的两 个方向， 为我们 提供 了有关地震机制的猜想。
北纬 3 5 度线是地球体积不变的条件下 ， 扁率变化的中和线 。
式 （1） 序 列 地 震 发 生 的 机 制 ， 可 能 是 由于地球扁率的变化。
北北东走向： 是环太平洋的走向， 式 （2） 序列地震发生的机制，可能与环太平洋运动的构造运动有关(尤指太平洋西岸）
（3）北京暴雨
北京暴雨的日子，也反映出沙中金序列的特征，如用上列公式（1）计算，即能表示这一序列，在1893〜1980半 的87年 内 ， 每 年日最大降水量大于100mm的
共有20次，其中7次暴雨日和公式（1）的计算值，差别不大于正负4天，如下表15所示。
因为公式（1）不是针对北京暴雨日期模拟的，所以符合区间取（-2， 4 ) 或 7 日，在均匀分布假设下，7 天 的平均符合 频 数 ， 只有 2 .3 3 ( 算法同前 2 0 / 6 0×7
=2.33。用一种假设检查方法，可得否定均匀分布假设的置信水平大大于90%。
进一步将公式 （1）用于研究1976〜1991年间降水量大于30mm/d的雨天，可得表16。
从上表可以祈出：公式（1）适用于两个时间段: 第一段从1976到1982年，第二段从1985到1991 年 。 在 每一时 间 段 中 ， 降水量大体从大到小， 再转大。
用其它方法预测，1991年北京暴雨日期有：1991 年7月 16日 ，1991年7月23曰 和1991年8月6日 ，
表 中y (416）=1991 年7月24 日， 和第二次预测可能暴雨日期相接近，1991 年7月24日 的 千 支 是 “ 乙未”，是“沙中金”第二天，十分接近。
（4）大江洪峰
某些 大型的 水文站 ，特大洪峰过站的日 期，在 一段时间内 ， 也可用白腊金序列的 式 （2）近 似 地 表 述 ， 下面是几个例子。
1）松 花 江 哈 尔 滨 站 在 18 9 8 〜 1 9 7 9 的 81年间 ， 一 年中最大日流量大于或 等于 4 9 0 0m³/s的 ， 有 1 4年 ， 其 中 6 次 的日子 ， 是式 （ 2 ) 答数日期前 5 ~ 1 4 天，如表17所示。
洪峰日期在y(í) 前 5 〜 1 4 天 ， 或 （ — 5 , — 1 4 ) 天 ， 其间的间隔区间为 1 0 天 ， 假如大洪峰日期为均匀分 布 ， 则在间隔区间10天 内 ， 应 为 2.33天（10 / 6 0 ×I4 = 2.33）现在实际 有 6天 ， 用一种 假设检查方 法 ， 可以得置信水平为 8 5 %，否 定 均 匀 分 布 的 假 设 。
实 际洪峰过站日期 ， 比 式（2） 的 计 算 值y(í) , 早 （或减 ） 9天 ， 令 y =1991.9.13，减9天得1991.9.4，可预测：1991年8月底、9月初，可能出现洪峰.
2 ) 淮河蚌埠站 在 1 9 5 4 〜 1 9 7 9 年 间 ， 日洪峰流量等于或大于 5 0 0 0 m3/ s 的有 1 0 年 ， 其中有 6 年的洪峰过站日期 比 白 腊 金 公 式 （2） 的 答 数 y (í) 晚 8 〜 1 7天，如表18所示。
在均匀分布假设下 ， 比式 （ 2 ) 答案时间 晚 8-17 天区间的平均天数只有 1 .6 7 天 （ 1 0 / 6 0 x1 0 ) , 而实际有6 天 .用 一 种假设检查方 法 ， 可得 ： 否 定 均 匀 分布的置信水平大于 9 0 % ， 洪峰过站日期平均比 式 （2）计算的y (í) ， 晚 （ 或 加 ） 1 3 天 。
令í=1 5 4 ， 可得y ( 154 ) =1991.7.4，加 13 天，得：1991.7.17 日前后几天内 ， 可能有轻度洪峰过蚌埠站 。
3 ) 长江武汉站 在 1 9 7 6 〜 19 8 9 年的 13年间，长江武汉站 ， 最大洪水日流量 ， 等于或大于 58000m3/s的日子共 有 7 次 ， 其中实际日期比式中（2 ) 的 y (í)
值早 7 〜 1 2 天区间之内 ， 即(-7)〜（-12）天区间之内的有 5 次 ， 如表 1 9 所示 。
同样， 根据均匀分布假设， 计算次数为 0 .7 次（6 / 6 0 × 7 ) ,实 际 为 5 次 ， 否 定均匀假设的置信水平为 9 0 % 以 上 ， 实际 洪峰日期比 式 （ 2 ) 计算值 y (í) ，早9天，和松花江、哈尔滨站校正值相同。
