专利名称:全球尺度强地震趋势预警系统的制作方法
技术领域:
本发明涉及一种全球尺度强地震趋势预警系统,也就是利用本世纪全球历史强地震数据通过统计相关性原理对全球未来强地震三要素做地震短临预测的方法。本方法可利用全球范围近日已经发生的数次6级以上地震的时空强三要素,得出未来一个月内在全球范围可能发生6级以上地震的区域。预测准确率在50%以上。本发明属统计地学领域。
对于强地震的发生,能否预测、是否存在关联性等问题专家们各持己见,颇有争论而无共识。本发明人是倾向于持客观看法的,在权威人士中例如中国科学家陈绍绪1986年提出中国发生的多数大地震都有很长的余震活动;美国加州大学地球物理学教授B.罗曼诺维奇1993年提出一次大地震能在20-30年后引发地球上其他地方发生一系列大地震;地震能量在全球转移,使一次地震与远距离的另一次地震连起来。中国物理学教授郑联达1996年证明了地震在全球尺度范围存在相互作用力;中国科学家强祖基教授1989年10月开始成功地利用卫星热红外遥感技术观测到强地震发生前地壳运动及其迁移直至发震的现象。
本发明人自1997年10月开始以全球尺度对1900年以来的强地震进行长时间的研究之后,从观察到的重要现象,认为地震有其自身的孕震和发震条件和发生的规律。同时,利用国际互联网及时对全球强地震三要素进行跟踪检索,发现强地震存在关联性。本发明人认为研究和预测地震的发生不应限于孤立的时、空、强概念,不仅以局部震情为分析样本对某一个孤立区域进行分析,而应以全球尺度整体地去观察强地震的孕震、发展与发震以及它们之间的一系列相关现象,从动态的分析角度,引入时间变量,即从地震动力学特征去分析和跟踪强地震。上述重要现象是指1、全球强地震的发生是按时间顺序依次进行的。
全球范围的强地震(特指破坏性大的等于和高于6级地震)每月一般发生10-20次,各次强地震的发生都存在时间间隔,平均在2-3天,至少1小时以上或者说不在同一时刻。这说明全球范围发生的强地震是遵从单一的时间坐标依次发生的自然现象(灾害),即强地震是先后相继发生。2、强地震之间存在相关性。
通过观察分析研究发现全球强地震的发生和走向存在某种相关现象,即强地震之间能量和力的转移能导致未来地震震中沿地震带的位移或转移(1)本人通过长期纵观全球强地震发生的现象之后,发现全球范围的地震其主震源区有时在太平洋中部,即印尼东部至斐济一带。相对而言,每次强地震的发生大多开始于此区域,然后向东西两方向扩展,特别是向西部扩展,影响到印尼西部及菲律宾并引发一系列的地震,然后由菲律宾北上,进而引发中国台湾、琉球群岛和日本的地震,最后延伸和触发日本北部乃至俄罗斯远东堪察加一带的地震。
(2)太平洋板块与欧亚大陆东部板块的相互挤压导致上述地震的相继发生,地震的发生又进一步触发和加速新的地震的形成,进而引发中国大陆及阿富汗、巴基斯坦以及伊朗的强地震。
(3)伊朗的强地震发生向周边释放和转移能量,引发地中海东岸及北岸的地震,包括土耳其、希腊和意大利的地震。
(4)地中海一带的强地震、欧亚大陆的强地震的发生,其能量转移可能以较快的速度穿过地心诱发地心对面的处于中、南美洲的强地震。例如1998年1月10日中国张北6.2级地震之后仅1小时,地心对面的秘鲁发生一次强地震。
(5)整个美洲的地震大多是在欧亚大陆包括地中海强地震发生之后通过地心传递而来的能量所激发,由南向北引发一系列强地震,即由南美洲发展到中美洲,直至通过美国、加拿大到达美国的阿拉斯加。南北美洲的强地震还存在相互呼应发生的形象。
本发明的目的在于提出一种全球尺度强地震趋势预警系统,该系统是指根据观察到的现象,根据各个强地震之间存在相关性,特别是每个强地震发生之前所受其他强地震的影响和它发生之后引发其他强地震特征,如近期发生的几个强地震之间的相关性及其与未来地震之间的关联性特征,利用本世纪以来的强地震目录(1900-1999),应用统计归纳比较的方法找出全球强地震发生的规律,并通过计算方法得出未来强地震的预警系统。
实施本发明全球尺度强地震趋势预警系统所述方法采取的具体步骤如下所述由1900年以来全球大于等于6.0级地震历史目录,选取其中地震三要素,即发震时间、震中位置和震级。
根据美国P.K.Dunbar等人编制的2150B.C.-1991A.D.全球强地震目录中给出的按全球强地震震中分布情况划分的25个地震区,并将其修改29个区域,并对每个区域编号。
设A表示地震区域,这29个地震区域可表达为一个序列A(i),i=1,2……29同时假设每个地震与其发生之前和之后发生的其他地震相关联。
按区域序号,依次记录每个区域从1900年以来所发生的每个地震的前后地震事件。
统计每个区域发生地震之前和之后发生地震事件的概率。
地震预测设近日全球范围发生的6次地震依次为A6,A5,A4,A3,A2,A1,其中A1是刚发生的地震所对应的区域编号,A6是最先发生的。将A1之后可能发生地震的区域及发震概率表示为一个序列A1(i,jh),jh=1,2,3,……20,i为已知定值。
在A1之后所有可能发震的区域序列中确认出未来地震发生的区域,有必要找出A1(i,jh)序列中每个区域发震前发生地震的区域编号及发震概率,可将其表示为A1[A1(i,jh),jq],i为已知定值;jh=jq=1,2,3,……20。
将上式中每个序列的前几项分别与A1至A6序列比较,如果其中某个序列的区域编号与A1至A6序列的区域编号相同的情况越多,该区域序列A1[A1(i,jh),jq]所对应的区域A1(i,jh)即为未来地震发震概率最高的区域,此时可将发震概率最高的几个区域确定为未来地震的发震区域。
为提高预测准确率由统计引入强地震区域的月份(季节)与发震概率的参数。按本发明预警系统预测强地震三要素的范围为震级大于等于6级地震;发震时间1-30日之内;震中某个强地震区域。
现结合附图对本发明全球尺度强地震趋势预警系统所述的方法作进一步的说明。
图1为按全球强地震震中分布情况划分为29个地震区的示意图。未来强地震的预警系统预测基本步骤1、由1900年以来(近百年)全球大于等于6.0级(Ms)地震历史目录,选取其中地震三要素,即发震时间、震中位置和震级。
2、根据美国P.K.Dunbar等人编制的2150B.C.-1991A.D.全球强地震目录中给出的按全球强地震震中分布情况划分的25个地震区,并将其修改为29个区域(见图1),并对每个区域编号,设A表示地震区域,这29个地震区域可表达为一个序列A(i),i=1,2,3,…,29(1)同时,假设每个地震与其发生之前和之后发生的其他地震相关联。
3、按区域序号,依次记录每个区域从1900年以来所发生的每个地震的前后地震事件。
3-1、首先是依次记录被记录区域每个地震发生之前的地震事件。
将该区域每次地震发生之前所发生的五个地震按其对应的区域编号记录下来,这五个地震所对应的区域可能是其他地震区域,也可能是该区域所发生的地震。例如所记录的区域其编号为29,它包括中国东部地区,则对应于1976年7月28日中国唐山7.8级地震,按上述记录法可得19,27,5,26,21,29。这里的19,27,5,26,21五个区域编号标志被记录的29号区域在1976年7月28日地震发生之前发生的其他地震所对应的区域编号。
3-2、按相同方法记录被记录区域每个地震发生之后的地震事件。
将该区域每次地震发生之后所发生的五个地震按其对应的区域编号记录下来,这五个地震所对应的区域可能是其他地震区域,也可能是被统计的区域所发震的地震。例如所记录的区域其编号为29,则对应1976年中国唐山7.8级地震,按记录法可得出29,20,16,11,26,18。这里的20,16,11,26,18五个区域编号标志被记录的29号区域在1976年唐山地震发生之后发生的其他地震所对应的区域编号。
4、统计每个区域发生地震之前和之后发生地震事件的概率。
4-1、每个区域地震之前已发生的其他地震出现的概率。
利用对某个区域从1900年以来所发生的每个地震之前的地震事件(利用3-1节中记录的结果),统计该区域发生地震之前所发生的各地震的概率(以区域编号出现的概率大小排序,概率最大者在前,由左向右递减)。
例如29号区域,其统计结果为jq=1,2,3,…,20。该区域发生地震之前的统计排序(q表示之前)统计序号jq=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20 (2)区域编号23,3,26,6,12,21,14,27,5,16,1,2,11,18,8,22,7,19,28,20这项统计结果说明29号区域发生地震之前所发生的其他地震,其发震区域和发震概率大小的关系。如近百年来的统计显示,在29号区域发震之前,23号区域发生地震的概率最大,其次是3号和26号区域,依次不等差递减。
此时,根据(1)式可将某区域地震前已发生的其他地震的区域编号和发震概率表示为A(i,jq)i=1,2,3,…,29;jq=1,2,3,…,20 (3)4-2、每个区域地震之后将发生的其他地震出现的概率。
利用对某个区域从1900年以来所发生的每个地震之后的地震事件(利用3-2节中记录的结果),统计该区域发生地震之后所发生的各地震的概率(以区域编号出现的概率大小排序,概率最大者在前,由左向右递减)。
例如包括中国东部地区的29号区域,其统计结果为jh=1,2,3,…,20。该区域发生地震之后的统计排序(h表示之后)统计序号jh=1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20 (4)区域编号3,11,6,26,27,18,23,21,16,14,12,5,7,19,20,28,25,10,2,8
这项统计结果说明29号区域发生地震之后所发生的其他地震所在区域和发震概率的关系。如近百年来的统计显示,在29号区域发震之后,3号区域发生地震的概率最大,其次是11号和6号区域。
此时,根据(1)式可将某区域地震后发生其他地震的区域编号和发震概率表示为A(i,jh),i=1,2,3,…,29;jh=1,2,3,…,20(5)5、地震预测5-1、设近日全球范围发生的6次地震依次为A6,A5,A4,A3,A2,A1(当然也可选择更多的区域,要根据预测准确率而定)其中A1是刚发生的地震所对应的区域编号,A2是在A1之前发生的,…,相比之下A6是最先发生的。这里A=A(i),i=1,2,3,…,29。即A1至A6是全球29个强地震区域中随机发生的地震所对应的区域编号。
根据前面的假设,A1至A6序列与未来地震及其发生的区域相关联。
将A1之后可能发生地震的区域及发震概率表示为一个序列A1(i,jh),jh=1,2,3,…,20(6)这里i是已知的定值,即A1表示为29个区域中的某个区域。假设此时i=29,则在A1之后可能发生地震的区域及发震概率构成一个序列,与(4)式相同。
5-2、这一步是要在A1之后所有可能发生地震的区域序列中确认出未来地震发生的区域,因未来地震区域可能发生在A1(i,jh)序列之中。按照这个推理,有必要找出A1(i,jh)序列中每个区域发震前发生地震的区域编号及发震概率,根据(5)式可将其表示为A1[A1(i,jh),jq],i为已知定值;jh=jq=1,2,3,…,20(7)5-3、(7)式反映了A1区域发震之后可能发生的其他地震所对应的区域序列以及该区域序列中各区域发震前发生其他地震的区域序列及其发震概率,或者说(7)式由左向右描述了未来地震发生之前发生地震可能性或发震概率大小的区域序列。这时可以看出,(7)式中的前几项与A1至A6序列存在相似性,即(7)式中的前几项也反映出与未来地震相关联。
将(7)式中的每个序列的前几项分别与A1至A6序列比较,即区域编号之间进行比较,如果(7)式中某一个序列的区域编号与A1至A6序列的区域编号相同的情况越多,则说明该区域序列A1[A1(i,jh),jq]所对应的区域A1(i,jh)即为未来地震发震概率最高的区域。这时可将发震概率最高的几个区域确定为未来地震的发震区域。
本发明人按照上述步骤对1998年以来的地震进行跟踪式(滚动式)预测分析,即设每次刚发生过的地震为A1,A1以前发生的几个地震依次排序得到A1至A6序列,然后进行分析,分析结果给出了未来强地震可能发生的区域及其发震概率。将发震概率最高的区域记录下来并与A1以后实际发生的地震进行统计和比较,发现实发地震所在区域与预测结果给出的发震概率最高的区域存在对应关系,其吻合程度在50%以上。
为了提高预测准确率,可以考虑区域性地震与季节的关系,因为每个强地震区域都有其月份(季节)与发震概率大小的关系。将本世纪6级以上地震按29个区域分别统计每年地震发生的月份,即可得出每个区域每年按月份的发震概率。按照上述方法得出的预测结果,一般是给出发震概率最高的地震发生区域,为在这些区域中进一步确定是哪一个或哪几个区域可能发生地震,可将季节(月份)参数引入,用以参考判断提高预测准确率。
以上预测方法可能给出未来地震发震区域,这是地震预测三要素之一。因为本发明所讨论的样本空间限定在震级大于等于6级地震的范围,所以预测结果反映出未来地震的震级也是在6级以上。又因为本预测方法是根据刚发生过的地震和近期发生的地震作为计算依据,加之全球每个月6级以上地震一般要发生10至20次左右,即平均每两至三天将发生一次6级以上地震,因此按这个方法预测的未来地震发震时段一般在几日至一个月之内。因此可以给出该预测方法预测强地震三要素的范围是震级(Ms) 大于等于6级地震发震时间 1-30日之内震中 某个强地震区域尽管本发明给出的地震预测三要素不能达到更精确的数值,例如没能得出具体的震级、发震时间和震中经纬度。但是,这样的地震预测结果仍具有地震短临预测的性质。这种方法实际上是一个全球尺度强地震趋势性短临预警专家智能系统,通过它不仅可能及时反映全球强地震的走向和发展趋势,还可能减少在应用其他预测方法时所发生的虚报漏报现象。专家可参考此预测结果,再结合相关的地球物理参数,可能作出更为准确的预测判断。如果与卫星热红外地震短临预测技术、地震能量积累法等预测技术相结合,有可能提高地震短临预测三要素的准确度。本系统的特点1、经过实发地震对预测结果的检验可以看出预测结果与实发地震存在对应关系。预测应验结果的统计分析表明,预测准确率在50%以上。
2、该方法为滚动式预测方法每当在全球范围内新发生一个强地震,则可将这个地震作为A1代入预测系统计算,并可立刻得出若干个未来强地震的相关信息。该方法能及时预测出在全球29个强地震区域中将要发生强地震的区域,从而显示全球强地震的发震趋势和走向。
3、预测要素和范围强地震震级大于等于6级;发震时间在1-30日之内;给出若干个发震概率最大的未来强地震区域。
4、该方法计算过程中始终贯穿和保留强地震样本的随机性特征和非线性内涵,包括预测结果给出的并非一个定值,而是发震概率最大的若干个区域,然后进行综合判断。
5、计算方法简单,计算速度快捷(通常一次预测仅需要十分钟时间,包括输入数据时间)。
本发明将这项预测技术作为一种新的分析手段应用到全球强地震趋势预测研究和防震减灾之中。这对于解开当今世界所面临的地震预测问题,特别是短临预测难题具有重要的社会意义和实际应用价值。预测结果与实发地震对照以下是应用新发明的全球尺度强地震专家智能系统对将要发生的强地震的预测和应验结果,其中区域编号一栏中的数字表示被预测出的未来强地震发震区域编号(该编号可与图1对照得出发震区域在全球中的位置),这些区域是在预测计算中由系统给出的发震概率最大的区域,而且其排列顺序也显示出其发震概率的大小程度,即越靠左边的编号,其发震概率越大。因为可能发生强地震的区域编号大都集中在编号系列的左边,故靠右一些的编号被省略掉,用″…″表示。数字下面的横线表示该次强地震预测之后一个月之内得到应验的区域编号,即表示该区域在被预测之后确实发生了强地震。A1下面对应的数字表示刚发生的强地震区域编号,编号右侧括号里的数字表示该区域在本月发生的第几次强地震。例如A1如果对应的区域编号是29(2),则表示出29号区域在当月发生的第2次强地震。预测计算结果中的序号从小到大的顺序表示A1随时间坐标更临近本例讨论的强地震发生的时间。计算使用的全球历史强地震数据采用日本东京大学地震研究所通过互联网在全球全天候公布的强地震目录,该目录在国际上具有权威性。例1、利用历史强地震数据验证本方法的效果(之一)对1997年11月9日中国西藏玛尼7.9级地震(35.0N,87.3E),利用地震发生之前的强地震数据代入本方法可得出以下结果。在区域编号一栏给出的区域有些在计算之后的几日至一个月内发生了强地震,这些实际发震区域如编号下的横线所示。其中27号区域覆盖中国西南部包括西藏地区(如图所示),对应1997.11.9.西藏玛尼7.9级强烈地震。序号 年、月A1 区域编号11997.10. 7(1)20.2.1.3.27.19.12.22.8.…21997.11. 8(1)7.20.21.2.3.23.12.19.27.…31997.11. 12(1) 7.2.3.23.27.8.20.21.…例2、利用历史强地震数据验证本方法的效果(之二)将1998年1月10日中国河北省张北6.2级地震之前的强地震数据代入本方法可得出如下结果。在区域编号一栏给出的区域有些在计算之后的几日至一个月内发生了强地震,这些发震区域如编号下的横线所示。其中29号区域覆盖中国东部地区(如图所示),对应1998.1.10.中国张北6.2级地震(41.3N,114.7E)。序号 年、月A1 区域编号11997.12. 6(2)14.12.20.11.21.29.6.…21998.1. 14(1) 1.21.12.11.29.…31998.1. 22(1) 1.11.12.21.7.27.29.…例3、对1998年7月17日中国台湾嘉义市6.2级地震的预测以下是对台湾嘉义市地震的预测结果。在地震发生前,利用全球已发生的强地震数据代人本方法计算可得出以下结果。在区域编号一栏给出的区域有些在计算之后的几日至一个月内发生了强地震,这些发震区域如编号下的横线所示。其中28号区域覆盖中国东南部包括台湾地区(如图所示),对应1998.7.17.台湾嘉义6.2级地震(23.5N,120.7E)。序号 年、月 A1 区域编号11998.7. 25(1)7.26.23.22.27.28.18.…21998.7. 3(2) 11.18.10.23.16.28.…31998.7. 24(1)11.16.18.23.8.21.7.28.22.…例4、对1998年8月27日中国新疆伽师6.6级地震的预测以下是对1998年8月27日中国新疆伽师6.6级地震的预测结果。在地震发生前,利用全球已发生的强地震数据代人本方法计算可得出以下结果。在区域编号一栏给出的区域有些在计算之后的几日至一个月内发生了强地震,这些发震区域如编号下的横线所示。其中19号区域覆盖中国西北部地区(如图所示),对应1998.8.27.中国新疆伽师6.6级地震(39.7N,77.4E)。序号 年、月 A1区域编号11998.8. 22(1)11.14.12.6.18.19.…21998.8. 6(2) 11.14.12.29.19.…例5、对1999年8月1日美国加州东北部6.0级地震的预测以下是对1999年8月1日美国加州东北部6.0级地震的预测结果。在地震发生前,利用全球已发生的强地震数据代入本方法计算可得出以下结果。在区域编号一栏给出的区域有些在计算之后的几日至一个月内发生了强地震,这些发震区域如编号下的横线所示。其中2号区域覆盖美国西部地区(如图所示),对应1999.8.1.美国加州西部发生的6.0级地震(37.4N,117.1W)。序号 年、月 A1 区域编号11999.7. 6(1) 21.12.27.11.19.2.7.14.…21999.7. 27(1)12.21.11.2.23.…31999.7. 24(2)21.12.1.23.14.11.2.…例6、对1999年8月17日土耳其7.8级地震的预测这是在土耳其大地震之前应用该方法对全球强地震进行滚动式跟踪预测的过程和结果。结果表明,自21(1)号区域发震之后,16号区域始终被显示,其位置在预测结果的顺序中偏左,直至该地震发生。16号区域对应的强地震即为1999.8.17.土耳其7.8级地震(40.6N,29.8E)。同时,从区域编号一栏中可以看出,除了16号区域以外,其他一些区域在预测之后也发生了强地震。序号 年、月 A1 区域编号11999.7. 21(1)11.29.23.16.18.12.14.…21999.7. 3(2) 11.16.18.10.29.…31999.7. 24(3)11.16.18.23.29.…41999.7. 3(3) 11.16.10.18.23.19.14.…51999.8. 1(1) 11.16.18..23.…61999.8. 2(1) 18.1 1.16.23.21.…71999.8. 1(2) 18.11.16.21.…81999.8. 26(1)18.19.1 1.16.5.8.20.…91999.8. 21(1)19.11.16.18.8.5.20.…例7、对1999年9月7日希腊5.9级地震的预测这是在希腊强地震之前应用该方法对全球强地震进行滚动式跟踪预测的过程和结果。结果表明,自6(2)号区域发震之后,12号区域始终被显示,其位置在预测结果的顺序中偏左,直至该地震发生。12号区域对应的强地震即为1999.9.7.希腊5.9级地震(38.1N,23.5E)。同时,从区域编号一栏中可以看出,除了12号区域以外,其他一些区域在预测之后也发生了强地震。序号 年、月 A1区域编号11999.8. 6(2) 14.26.7.16.21.12.…21999.8. 27(1)26.12.22.25.3.…根据对百年来全球强地震进行数以百计的预测试验结果可以认为该方法与自然强地震之间存在一定的对应关系。然而,地震的发生富有极强的随机属性、非线性和不确定性,本方法只能给出将要发生强地震区域的概率。而地震的发生属小概率事件,即发震概率小的不等于不发生。因此,尽管本方法给出的结果不是一个唯一的预测值,或者说本方法在整个计算过程中始终保留自然现象的随机、相关的属性,但是其预测结果与实发地震还是存在某种程度的不一致性,这也许是大自然固有的属性使然,尚须人们进一步探索。
权利要求
1.全球尺度强地震趋势预警系统,其特征在于所述方法的具体步骤如下由1900年以来全球大于等于6.0级地震历史目录,选取其中地震三要素,即发震时间、震中位置和震级;根据美国P.K.Dunbar等人编制的2150B.C.-1991A.D.全球强地震目录中给出的按全球强地震震中分布情况划分的25个地震区,并将其修改29个区域,并对每个区域编号;设A表示地震区域,这29个地震区可表达为一个序列A(i),i=1,2……29同时假设每个地震与其发生之前和之后发生的其他地震相关联;按区域序号,依次记录每个区域从1900年以来所发生的每个地震的前后地震事件统计每个区域发生地震之前和之后发生地震事件的概率;地震预测设近日全球范围发生的6次地震依次为A6,A5,A4,A3,A2,A1,其中A1是刚发生的地震所对应的区域编号A6是最先发生的;将A1之后可能发生地震的区域及发震概率表示为一个序列A1(i,jh),jh=1,2,3,……20,i为已知定值在A1之后所有可能发震的区域序列中确认出未来地震发生的区域,有必要找出A1(i,jh)序列中每个区域发震前发生地震的区域编号及发震概率,可将其表示为A1[A1(i,jh),jq],i为已知定值;jh=jq=1,2,3,……20将上式中每个序列的前几项分别与A1至A6序列比较,如果其中某个序列的区域编号与Al至A6序列的区域编号相同的情况越多,该区域序列A1[A1(i,jh),jq]所对应的区域A1(i,jh)即为未来地震发震概率最高的区域,此时可将发震概率最高的几个区域确定为未来地震的发震区域。
2.如权利要求1所述的全球尺度强地震趋势预警系统,其特征在于为提高预测准确率由统计引入强地震区域的月份(季节)与发震概率的参数。
3.如权利要求1所述的全球尺度强地震趋势预警系统,其特征在于该预警系统预测强地震三要素的范围为震级大于等于6级地震;发震时间1-30日之内;震中某个强地震区域。
全文摘要
本发明全球尺度强地震趋势预警系统所述方法的步骤如下:由本世纪以全球大于等于6.0级地震历史目录选取地震三要素,按全球强地震震中分布划分29个区域,并对其编号以序列表示,同时假设每个地震与其发生前后发生的其他地震相关联,依次记录每个区域从本世纪发生每个地震的前后地震,统计其概率,将刚发生地震之后可能发生和发震前发生地震区域编号和发震概率以序列表示,某区域序列与图示地震区域编号相同情况越多,对应区域即为未来地震发震概率最高区域。
文档编号G01V1/00GK1304049SQ9912622
公开日2001年7月18日 申请日期1999年12月16日 优先权日1999年12月16日
发明者宋期, 郑联达, 宋越 申请人:宋期