为了通过HTTP数据传输协议下载JOPENS6系统存储的观测数据，需要编写待下载数据的相关信息，并将此信息提交给系统，从而实现数据的下载。

图1 断计统计执行流程Fig.1 Flowchart of gap count

为了下载数据，本文使用了requests库。这是一个Python的第三方库，通过简单的方式即可实现对网络资源的访问。首先建立一个新的会话对象（session）

#导入requests库

import requests

j6Session = requests. Session（）# 建立一个session对象

然后，构建一个字符串，该字符串包含JOPENS服务器的IP地址、端口号、台网代码、台站代码、数据起止时间。通过会话对象（seesion）将其提交给JOPENS6系统，系统会根据提交的信息返回相应观测数据。关键代码如下：

#构建下载信息字符串

cmd_download = ‘http://JOPENS6IP: 8080 /jopens-ws/app/aws/gwfw; net=BJ; sta=DAX； loc=00；cha=HHZ； beg=2021-01-01 01: 00: 00； end=2021-01-01 02:00:00’

resp = j6Session. get（cmd_download）# 获 取JOPENS6返回信息

with open（outfile_name，'wb'）as f:#打开文件，以二进制形式写入

f.write（resp.content）#将数据写入文件

通过上述操作，即可将需要的观测数据下载到本地，并以二进制miniseed文件形式进行保存。

obspy库有一个内置的read函数，可读取多种格式的地震数据，读入的数据构成流（stream）对象。一个流对象可由一个或多个trace组成，每个trace都有起止时间，可以通过trace的函数获得。需要注意的是，trace的起止时间是UTC时间，为了更符合人们的工作习惯，需要把UTC格式的起止时间替换为北京时间。本文即根据trace的起止时间统计断计结果。流、trace起止时间的示意见图2，其中t₁、t₂、t₃为断计时长。关键代码如下：

#导入dateutil模块和obpsy的read函数

import dateutil

from obspy import read

st=read（‘miniseed文件名’）#读取miniseed文件，构成名为st的流对象

trace1=st[0]#从流对象列表中取一个trace

#trace起止时间的获取

trace1_starttime = dateutil. parser. parse （str （trace1.stats.starttime））\ #提取开始时间

. astimezone（pytz. timezone（'Asia/Shanghai'））\#转换为北京时间

.replace（tzinfo=None）#转换为不含时区的类型trace1_endtime=dateutil.parser.parse（str（trace1. stats.endtime））\

.astimezone（pytz.timezone（'Asia/Shanghai'））\.replace（tzinfo=None）

若一个 stream 包含两个 trace，如 trace1 和trace2，根据以上代码分别获取各自的起止时间，而后将trace2的开始时间与trace1的结束时间相减，即可获得两者之间的断计时间，即

timedelta=trace2_starttime–trace1_endtime

timedelta.seconds#断计的秒数

timedelta.microseconds#断计的微秒数

注意，在统计断计时长时，还需要考虑第1个trace与miniseed文件标称开始时间之差，以及最后一个trace与miniseed文件标称结束时间之差。将这两个时间差与各trace间时间差相加起来，即可获得下载数据的总的断计时长。

至于断计次数，分为两种情况：若第1个trace的开始时间、最后一个trace的结束时间与下载数据时设置的起止时间一致，则：断计次数=trace总数-1。否则，断计次数=trace总数-1+（首或尾的断计次数）

通过上述计算，即可获得下载数据的断计时长和断计次数。将统计时段、断计时长、断计次数输出为csv格式的文件，以便于使用excel等工具进一步开展统计工作，如绘制统计图表，进而对断计分布情况做出分析。断计统计见图3。

图2 流（stream）对象与trace关系示意Fig.2 Diagram of relationship between stream and trace

图3 断计统计图Fig.3 Diagram of gap count

本文软件选用企业微信作为信息发布的媒介主要基于如下三个考虑，一是企业微信是腾讯公司针对企业开发的通信和办公工具，具有和微信一致的交流体验；二是企业微信中可以包括若干应用，各应用对谁可见由管理员根据需要进行控制；三是企业微信提供了开发接口，可以实现与开发者的软件对接，从而实现信息的发布、交互等。为了实现将断计统计结果自动发送给运维人员，本文软件的实现思路是：首先通过企业微信ID、应用的secret，获取接口调用凭证；其次，将统计结果文件提交到腾讯临时媒体库，并获得一个媒体ID；最后，向指定应用发送该媒体ID，从而实现了将文件发送到指定应用的目标。以下为关键代码：

corp_id="企业微信ID"

app_secret="应用的secret"

url = 'https://qyapi. weixin. qq. com / cgi-bin /gettoken?corpid={}&corpsecret={}'

#获取接口调用凭证

access_token = requests. get （url=url. format （corp_id，app_secret））.json（）.get （'access_token'）

file_url='https://qyapi. weixin. qq. com/cgi-bin/media / upload? access_token= {}&type=file'. format （access_token）

#上传文件到临时媒体库

media_id=requests. post（file_url， files= {'file':open（‘文件绝对路径’，'rb'）}）.json （）['media_id']

send_file_url='https://qyapi. weixin. qq. com/cgi-bin / message / send? access_token= {}'. format （access_token）

send_values= {'touser':企业微信的用户名， 'msgtype': 'file'， 'agentid': 应 用 的 ID， 'file':{'media_id':media_id}，}

#将临时媒体库的文件发送到企业微信的应用

requests. post （url=send_file_url， data= json. dumps（send_values））

本文软件使用的定时任务模块为Apscheduler，在windows、linux上均可运行。通过设置任务的开始时间，即可实现周期性的定时任务。通过该模块，每日定时从JOPENS6服务器下载前一天的观测数据，并对其进行断计统计，最后将统计结果发送到企业微信。关键代码如下：

#导入Apscheduler相关模块

from apscheduler. schedulers. background import BlockingScheduler

from apscheduler.triggers.cron import CronTrigger scheduler = BlockingScheduler（）# 建立一个定时任务对象

trigger = CronTrigger（hour='1'， minute='00'） #设置定时任务的开始时间

#周期性执行routine函数，该函数是将前述的数据下载、断计统计整合在一起的函数

scheduler.add_job（routine，trigger，max_instances=10）