CN103106341B - 一种直流保护事故智能分析方法及*** - Google Patents

Abstract

本发明公开了一种直流保护事故智能分析方法及***，其方法包括以下步骤：步骤1：导入故障时刻的信号文件和录波文件；步骤2：将录波文件中的数据解析为计算机可识别的格式；步骤3：对所有动作信号进行筛选，选取动作的保护类型；步骤4：进行智能分析，按照预先给定的判据，自动提取相关录波数据，并分析保护动作过程是否符合设计预期，对于满足判据的数据段，***采用自动截图的方式，生产事故分析报告。采用本发明可大大增强了一线人员故障特征辨识能力，有效提高现场人员直流保护事故处理效率。

Description

一种直流保护事故智能分析方法及***

技术领域

本发明涉及高压直流输电***，尤其涉及的是一种直流保护事故智能分析方法及***。

背景技术

面对技术原理复杂，覆盖面广，出现频率高的继电保护事故，如何快速辨识事故特征，及时定位故障点，尽快恢复故障设备正常运行，是当前直流输电一线生产工作迫切需要解决的难题。

通过文献检索发现，目前国内外对直流保护的研究多数是单个保护算法、原理的研究，几乎没有涉及高压直流保护事故自动分析方法的相关研究，而能同时对西门子TDC、南瑞PCS-9550、西门子SIMADYND三种直流保护***下的事故处理技术进行完整研究的就更是凤毛麟角。

对于现场出现的直流保护事故，伴随着一系列的事故告警信号（简称SER）和故障录波信号（简称TFR），这些信号包含的信息量非常丰富。一般一次跳闸事故，告警信号多达几百条，但真正对我们有用的信号，只有10余条；而故障录波信息就更多了，如云广特高压的一个故障录波文件中，就包含了上万条电流、电压、极控、保护开关量等故障数据。在传统的事故处理过程中，只能依靠人工进行筛选，找出与故障分析有关的SER和TFR，费时费力，而且容易因为分析者的经验差异导致分析结果的不准确。

因此，现有技术有待于进一步发展和改进。

发明内容

本发明的目的在于提供一种直流保护事故智能分析方法及***，旨在解决现有的直流保护事故没有能够实现智能处理的问题。

本发明的技术方案如下：

一种直流保护事故智能分析方法，其包括以下步骤：

步骤1：导入故障时刻的信号文件和录波文件；

步骤2：将录波文件中的数据解析为计算机可识别的格式；

步骤3：对所有动作信号进行筛选，选取动作的保护类型；

步骤4：进行智能分析，按照预先给定的判据，自动提取相关录波数据，并分析保护动作过程是否符合设计预期，对于满足判据的数据段，***采用自动截图的方式，生产事故分析报告，指出可能的故障点。

所述的直流保护事故智能分析方法，其中，将录波文件中的数据解析为计算机可识别的格式首先将数据转化为电力***暂态数据交换公用格式。

所述的直流保护事故智能分析方法，其中，将录波文件中的数据解析为计算机可识别的格式首先将数据转化为Excel格式。

所述的直流保护事故智能分析方法，其中，对所有动作信号进行筛选采用了多重筛选方法，具体包括以下步骤：

步骤31：保护名称筛选；

步骤32：动作后果关键字搜索；

步骤33：根据保护判据锁定录波数据。

所述的直流保护事故智能分析方法，其中，***进行智能分析的具体步骤包括：

步骤41：***根据得出保护动作信息，找出对应的判据方程，将之转化为计算机能识别的格式；

步骤42：***对一些有特殊运算的数据量进行处理，将单一量转化为判据中需要的多元量；

步骤43：***开始录波扫描，根据筛选的结果，对用到的录波数据进行扫描，以找到符合动作判据的录波段；

步骤44：***将扫描出来的录波数据与判据进行对比，做出保护是否正确动作的结论。

所述的直流保护事故智能分析方法，其中，对步骤41中所述的将判据方程转化为计算机能识别的格式的方法分为：方程规范化和方程解析两个过程。

所述的直流保护事故智能分析方法，其中，***对一些有特殊运算的数据量进行处理具体包括：进行有效值计算、进行电压突变量计算、进行交流电流的有效值计算。

所述的直流保护事故智能分析方法，其中，***进行录波扫描的具体方法是根据SER信号的分析结果，先找出需要计算的判据，再根据该判据找到对应的开关量名称，在录波文件中查询该开关量所在通道，并且找到开关量发生变化的时刻，并以此刻为起点，对相关的模拟量通道进行扫描。

一种直流保护事故智能分析***，其中，包括：

导入故障时刻的信号文件和录波文件的录波数据上传模块；

识别导入的信号文件和录波文件是文本格式还是非文本格式，若是文本格式则将数据输入SER信号分析模块进行筛选，若是非文本格式，则保护***及保护动作名称输入模块被启动，要求主动输入保护***及保护动作的名称的录波数据识别模块；

输入保护***及保护动作的名称的保护***及保护动作名称输入模块；

进行筛选，筛选出保护动作的SER信号分析模块；

从筛选出的保护动作中提取出它们的名称的保护动作名称提取模块；

先将录波数据格式转换为计算机统一识别的数据格式，再根据自动提取出的保护名称或主动输入的保护名称查询出该保护动作的判据的数据解析和筛选模块；

通过录波扫描找到符合动作判据的录波段，并于与查询出的保护动作的判据进行对比，做出保护是否正确动作的结论的录波数据分析模块；

通过鼠标事件机制将方程计算结果的显示在信息栏的故障片段截取模块；

将分析的结果以WORD格式导出的故障报告输出模块，指出可能的故障点。

所述的直流保护事故智能分析***，其中，***的故障片段截取模块会自动截取一段满足判据的波形图，数据中包含瞬时值、有效值和每个公式的计算结果。

本发明的有益效果：本发明通过基于软件编程技术，根据预先判据、信号，自动解析故障录波，形成了智能分析各类直流保护故障的方法，大大增强了现场人员故障特征辨识能力，有效提高了现场人员直流保护事故处理效率。

附图说明

图1是本发明提供的***框图。

图2是本发明提的***的分析方法流程图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。

参见图1，所述直流保护事故智能分析***包括：录波数据上传模块、录波数据识别模块、SER信号分析模块、保护动作名称提取模块、保护***及保护动作名称输入模块、数据解析和筛选模块、录波数据分析模块、故障片段截取模块和故障报告输出模块。上述模块组成的***主要分为四大部分，其中，录波数据上传模块、录波数据识别模块、SER信号分析模块、保护动作名称提取模块和保护***及保护动作名称输入模块用于完成数据导入；数据解析和筛选模块用于完成数据解析和数据筛选；录波数据分析模块、故障片段截取模块和故障报告输出模块用于完成智能分析。

所述录波数据上传模块用于故障时刻的信号文件和录波文件的导入；所述录波数据识别模块用于识别导入的信号文件和录波文件是文本格式还是非文本格式，若是文本格式则将数据输入SER信号分析模块进行筛选，若是非文本格式，例如图片，则保护***及保护动作名称输入模块被启动，要求主动输入保护***及保护动作的名称；所述保护***及保护动作名称输入模块用于输入保护***及保护动作的名称；所述SER信号分析模块用于进行筛选，筛选出有哪些保护动作；所述保护动作名称提取模块用于从筛选出的保护动作中提取出它们的名称；所述数据解析和筛选模块用于先将录波数据格式转换为计算机统一识别的数据格式，再根据自动提取出的保护名称或用户主动输入的保护名称查询出该保护动作的判据；所述录波数据分析模块用于通过录波扫描找到符合动作判据的录波段，并于与查询出的保护动作的判据进行对比，做出保护是否正确动作的结论；所述故障片段截取模块用于通过鼠标事件机制将方程计算结果的显示在信息栏；所述故障报告输出模块用于将分析的结果以WORD格式导出。

参见图2，本发明提供的直流保护事故智能分析***的分析方法包括以下步骤：

步骤1：导入故障时刻的信号文件和录波文件；

步骤2：将录波文件中的数据解析为计算机可识别的格式；

步骤3：对所有动作信号进行筛选，选取动作的保护类型；

步骤4：按照预先给定的判据，自动提取相关录波数据，并分析保护动作过程是否符合设计预期，对于满足判据的数据段，***采用自动截图的方式，生产事故分析报告，指出可能的故障点。

由于故障录波数据量庞大，且不同的直流工程的录波数据格式各不相同，所以在步骤S2中***首先需采用数据识别技术，将录波数据解析为可识别的格式，然后再进行筛选。

本发明提供的第一种解析方法为：首先将数据转化为电力***暂态数据交换公用格式，其解析步骤为：

步骤A21：读取文件流到内存里，若文件不存在即报错；

步骤A22：新建一个ComtradeImpl类型的录波数据对象（ci）；

步骤A23：以行为单位来读取文件中的数据，读取单行数据后，以逗号为分隔符分割数据，获得一个字符串数组arr；

步骤A24：依据配置文件的格式定义，第一行存储的为站名、特征和修改年份，第二行存储的为通道的数量和类型，接下来的依次为模拟通道信息列表、状态通道信息列表、线路频率、变速采样、采样率信息、第一数据点时间、触发点时间、数据文件类型和时间标记乘数系数；

步骤A25：定义一个整型（int）对象count来记录当前读取到第几行；

步骤A26：使用BufferedReader以行为单位循环读取配置文件流，判断count，若为0，即arr数组存储的是站名、特征和修改年份；若为1，即arr数组存储的是通道的数量和类型，遍历arr数组，分别存储到对象ci里相应的数据结构，依次类推，直到读完文件流，即录波的配置文件读取完毕；

步骤A27：读取配置文件完，开始读取数据文件流到内存里，若文件不存在即报错；

步骤A28：由于数据文件是以二进制来存储数据的，所以需要以字节为单位来读取文件中的数据，存储到字节数组fileByte对象中；

步骤A29：根据通道的数量计算出一条记录所占用的字节数circulate；

步骤A210：新建一个类型为SampleData的记录数据对象（sampleData）；

步骤A211：循环字节数组fileByte，以circulate为单位，读取对应的数据到sampleData，当读取完一条记录后，存储sampleData到ci里的sampleDataList中，依次类推，直至文件流读完，即录波的数据文件读取完毕。至此，ComtradeImpl对象ci里面存储着一个录波文件的完整数据了。

本发明提供的第二种解析方法为：首先将数据转化为Excel格式，相比第一种方法Excel文件格式的解析要简单许多，其主要格式是一个二维表结构，每一列代表一个通道，每一行代表一个时间点，因此在程序解析时通过Java的POI组件读入该Excel表格，按照列代表通道，行代表时间点的格式就可以解析出整个录波数据。

其解析步骤为：

步骤B21：***读取xls文件流到内存里，若文件不存在即报错；

步骤B22：新建一个ComtradeImpl类型的录波数据对象（ci）；

步骤B23：根据Excel文件流创建Workbook对象wb，因为有的***的数据只存储到第一个工作表里，所以这里获取第一个工作表对象sheet，以行（row）为单位来读取文件中的数据，每一个单元格为一个数据点；

步骤B24：工作表sheet前八行存储的为录波的配置信息，依次读取前八行存储到ci里对应的数据结构；

步骤B25：新建一个类型为SampleData的记录数据对象（sampleData）；

步骤B26：从第九行开始循环，每行为一条记录，存储到sampleData，读取完一条记录后，存储sampleData到ci里的sampleDataList，依次类推，直至文件流读完，即录波的数据文件读取完毕。至此，ComtradeImpl对象ci里面存储着一个录波文件的完整数据了。

进行步骤3时，***面对数以万计的故障信号和故障录波数据流，如何快速、准确的找出我们需要的信息，是影响现场事故处理快速性和正确性关键的一步。本发明在此环节采用了多重筛选，具体包括以下步骤：

步骤31：保护名称筛选。每次保护动作过程中，主事件一般都是某项保护动作的报文，如突变量保护，行波保护等；确定了主线后，可以保证后续筛选的正确性。

步骤32：动作后果关键字搜索。由于保护动作后，必然是伴随着一定动作后果的，如跳闸、降功率、双极平衡运行等；可进一步缩小检索范围；

步骤33：根据保护判据锁定录波数据。每项保护都有自身固定的动作判据，如：极差动保护(87DCM)，动作判据为|IdH-IdN|>150A，延时5ms。当***经过第一、第二层筛选，检测到有87DCM保护动作时，则自动关联到录波数据中的IdH（高压母线直流电流）、IdN（中性母线直流电流）两个录波通道，大大提高了后面录波扫描、方程解析和录波分析的效率。

智能事故处理***通过上述三个步骤的筛选，可以从海量数据中，准确挑选出我们要的事故信号和录波数据，有效排除干扰信息。

经过前面三步筛选的准备后，***开始进行智能分析，其分析方法包括以下步骤：

步骤41：***会根据得出保护动作信息，找出对应的判据方程，将之转化为计算机能识别的格式；

步骤42：***对一些有特殊运算的数据量进行处理，以便于将现场采集到的单一量转化为判据中需要的多元量；

步骤43：***开始录波扫描，根据筛选的结果，对用到的录波数据进行扫描，以找到符合动作判据的录波段；

步骤44：***将扫描出来的录波数据与判据进行对比，做出保护是否正确动作的结论。

对步骤41中所述的判据方程，将之转化为计算机能识别的格式的方法分为两个过程，分为：方程规范化和方程解析。

所谓方程规范化是指，将用户提供的方程整理成计算机***能够理解的规范格式。因为用户书写方式的时候会比较随意，方程中可能会出现任意的空格、中文括号或英文括号、绝对值的表达有时是ABS函数，有时是双竖线的数学表示，这种不统一的公式输入是无法被计算机程序简单识别的，这就需要对方程进行规范化。而方程的解析过程就是将规范化后的数学公式转换为计算机程序，以进行数据处理。

本发明中，方程规范化主要有如下约定：

1、方程中不要出现空格

2、所有函数使用大写字母表示，比如MIN，MAX等

3、绝对值使用双竖线表示

4、定值使用Δ来表示

5、关系“和”使用&&表示，关系“或”使用or来表示

规范化后的方程使用了一套名为IKExpression的表达式解析引擎来进行分析和处理。从表面上来看，只要把公式和公式中的参数传给表达式解析引擎，就可以获得最终的计算结果，大大方便了众多复杂公式的解析过程，也节省了工作时间。

各个直流***的录波文件中，给出的都只是单个数值量的录波，如：直流电压UdL、直流电流IdL等，但某些保护判据中需要用到多元量，如：电压突变量du/dt，阀星侧交流电流有效值IacY等。对于这些特殊量，本发明按照一定运算方法进行处理，如：有效值计算，电压突变量计算，交流电流的有效值计算。

关于有效值计算：对于交流模拟量，例如IacY等，录波数据中记录的是该量的瞬时值，在进行数据计算时，瞬时值没有明显的意义（尤其是在跟直流值作比较的时候），需要计算转换为有效值。虽然我国使用的三相交流电是标准的正余弦波形，但是这只是理论上的结果。实际中，由于测量仪器的精度、外界干扰等问题，不可能严格按照正余弦的波形变换，而且当***发生故障时，波形会出现明显的抖动，种种原因都使得交流模拟量的有效值，不能按照理论公式进行推算。正余弦波的有效值与最大值之间在理论上有如下的关系：

$I_{\mathrm{RMS}}=\sqrt{2}\·I_{\mathrm{MAX}}$ 公式1

式中：I_RMS是有效值，I_MAX是最大值。

任意波形的有效值（即均方根值），可以通过如下的积分获得：

$I_{\mathrm{RMS}}=\sqrt{\frac{1}{T}\underset{0}{\overset{T}{\∫}}{i}^2\left(t\right)\·\mathrm{dt}}$ 公式2

式中：I_RMS是有效值，T是周期，i(t)是t时刻的瞬时值。

但是这样的计算方法也只对连续函数有效。很明显在现实中，通过录波得来的数据必定在时间点上是离散的，因此需要将以上公式离散化为：

$I_{\mathrm{RMS}}=\sqrt{\frac{1}{T}\underset{0}{\overset{T}{\mathrm{\Σ}}}{i}^2\left(n\right)}$ 公式3

式中：I_RMS是有效值，T是周期，i(n)是在n时刻的采样数据且要求n是2的正整数次幂。

只要将一个周波内的录波数据代入公式3，就可以计算得到某个时间点上对应的有效值。

关于电压突变量计算：有些的录波数据文件中，没有录得du/dt的数据，因此需要***自行计算。du/dt是个微分量，因此需要取一个极小时间，根据贵广录波的频率6400Hz，算得其周期为0.156ms，可以符合极小的条件，因此在计算du/dt的时候，***取得两个临近采样点的UdL数值，做差以后再除以0.156ms的周期即得，如下：

$\mathrm{du}/\mathrm{dt}=\frac{\mathrm{UdL}\left(n\right)-\mathrm{UdL}(n+1)}{0.156}$ 公式4

式中：UdL(n)为n时刻的采样数据。

关于交流电流的有效值Iac的计算：***还需要计算的量就是Iac（Uac同理）。在录波文件中，可以录到分项的IacY和IacD，即阀星侧A项电流、B项电流和C项电流，以及阀角侧A项电流、B项电流和C项电流。想要获得Iac必须先算得不分相的阀星侧电流和阀角侧电流。计算公式如下：

$\mathrm{IacY}=\frac{\mathrm{IacY}{A}_{\mathrm{RMS}}+\mathrm{IacY}{B}_{\mathrm{RMS}}+\mathrm{IacY}{C}_{\mathrm{RMS}}}{2}$ 公式5

$\mathrm{IacD}=\frac{\mathrm{IacD}{A}_{\mathrm{RMS}}+\mathrm{IacD}{B}_{\mathrm{RMS}}+\mathrm{IacD}{C}_{\mathrm{RMS}}}{2}$

得到IacY和IacD之后，Iac的计算公式为：

Iac=max(IacY,IacD)公式6

其中，步骤43所述的录波扫描是根据SER信号的分析结果，***会先找出需要计算的判据，根据该判据找到对应的开关量名称，在录波文件中查询该开关量所在通道，并且找到开关量发生变化的时刻，并以此刻为起点，对相关的模拟量通道进行扫描。扫描的长度大约为判据出口时间的长度，这样一方面把数据量限定在判据所需的若干通道上，而且通过开关量的变化，还可以缩短需要计算的录波时间段，显著减少了参与运算的数据量，明显提高计算分析效率。

处理到这一步，已经可以获得判据方程，以及所需的录波数据或根据录波数据计算来的数据。录波分析的过程就变成了对录波数据的扫描。

以时间点为索引，扫描录波数据，获取每个时间点上模拟量的数值，代入公式计算。如果满足判据要求，则检查出口时间，如果出口时间为0，则说明该保护出口动作；如果出口时间不为0，则记录下该点数据，继续计算下一个点，直到计算结果不满足判据，或出口时间已到。

以上是计算保护动作正确的情况。如果是保护不正确动作的情况，***就没有开关量可以参考，只能取全部录波数据进行计算，以判断是否有符合定值和出口时间的录波片段存在。

对于判据中有逻辑关系的情况，会分别计算每个方程的结果，再取相应的逻辑计算即可。对于分段的判据，就将每一段当作一个独立的判据对待，判据分多少段，就对录波数据进行多少次扫描，以判定具体是哪段保护发生了动作。

录波分析完成以后，***的故障片段截取模块会自动截取一段满足判据的波形图，数据中包含瞬时值、有效值（如果存在的话）和每个公式的计算结果。通过Rickshaw提供的鼠标事件机制，可以获得当前鼠标所在位置，通过该位置就可以得到波形图上鼠标所在点的具体坐标。将x坐标所在的位置画上一条竖线以示指针，而y坐标（即各通道的值）显示在图形下方的信息栏中。同时每一组图形的x轴指针是可以跟随移动的，而每一组还有一个用来显示方程计算结果的信息栏，其中的数据也会随光标移动发生变化。

满足判据的故障录波显示出来之后，***的故障报告输出模块还提供导出报告的功能。点击导出报告的按钮，会先弹出对话框让用户确认生成报告的内容，可以根据需要修改，然后点击确定，就可以导出Word格式的录波分析报告了，可以大大节省后续出分析报告的时间。

本发明将传统依靠人工进行的事故分析过程，通过数据解析和软件编程方法，转变为自动分析过程，在国内外事故处理技术中都是一种创新和突破。本高压直流输电保护事故处理***根据现场使用情况来看，项目成果大大增强了一线人员故障特征辨识能力，有效提高现场人员直流保护事故处理效率。

在以往的直流保护事故处理过程中，我们还没有一套***性的分析处理方法，遇到紧急事故后总是在现场临时查找二次回路、软件逻辑，同时还要花费较长时间来分析保护录波和相关信号，认知上的不足，严重影响了我们事故处理的效率，不利于尽快恢复故障设备投入运行。

该***应用到现场以后，用户只需要导入对应的SER和TFR数据，该***可根据事先设定的判据自动分析保护是否正确动作，还可以自动生成事故分析报告。在节省了人工分析时间的同时，为快速定位故障定，恢复设备运行提供了强有力的技术保障。

例如，±500kV广州换流站SER报“换流器保护换相失败（87CFP）被检测到”，继而发生双极换相失败。运行人员将对应TFR和SER导入《直流保护事故智能分析***》后，***自动分析对应的开关量和模拟量，很快就得出了保护正确动作的结论，并生成了固定格式的分析报告，与人工分析结论一致。此外，***根据历史经验数据，直接指出了可能的故障点，帮助现场很快查出了事故原因，原本需要12小时才能完成的事故处理，这次仅用了3个多小时。

应当理解的是，本发明的应用不限于上述的举例，对本领域普通技术人员来说，可以根据上述说明加以改进或变换，所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (6)

1.一种直流保护事故智能分析方法，其特征在于，包括以下步骤：

步骤1：导入故障时刻的信号文件和录波文件；

步骤2：将录波文件中的数据解析为计算机可识别的格式；

步骤3：对所有动作信号进行筛选，选取动作的保护类型；

步骤4：进行智能分析，按照预先给定的判据，自动提取相关录波数据，并分析保护动作过程是否符合设计预期，对于满足判据的数据段，***采用自动截图的方式，生产事故分析报告，指出可能的故障点；其中，***进行智能分析的具体步骤包括：

步骤41：***根据得到的保护动作信息，找出对应的判据方程，将之转化为计算机能识别的格式；

步骤42：***对一些有特殊运算的数据量进行处理，将单一量转化为判据中需要的多元量；

步骤43：***开始录波扫描，根据筛选的结果，对用到的录波数据进行扫描，以找到符合动作判据的录波段；

步骤44：***将扫描出来的录波数据与判据进行对比，做出保护是否正确动作的结论。

2.根据权利要求1所述的直流保护事故智能分析方法，其特征在于，对步骤41中所述的将判据方程转化为计算机能识别的格式的方法分为：方程规范化和方程解析两个过程。

3.根据权利要求1所述的直流保护事故智能分析方法，其特征在于，***对一些有特殊运算的数据量进行处理具体包括：进行有效值计算；进行电压突变量计算以及进行交流电流的有效值计算。

4.根据权利要求1所述的直流保护事故智能分析方法，其特征在于，***进行录波扫描的具体方法是根据SER信号的分析结果，先找出需要计算的判据，再根据该判据找到对应的开关量名称，在录波文件中查询该开关量所在通道，并且找到开关量发生变化的时刻，并以此刻为起点，对相关的模拟量通道进行扫描。

5.一种直流保护事故智能分析***，其特征在于，包括：

导入故障时刻的信号文件和录波文件的录波数据上传模块；

识别导入的信号文件和录波文件是文本格式还是非文本格式，若是文本格式则将数据输入SER信号分析模块进行筛选，若是非文本格式，则保护***及保护动作名称输入模块被启动，要求主动输入保护***及保护动作的名称的录波数据识别模块；

输入保护***及保护动作的名称的保护***及保护动作名称输入模块；

进行筛选，筛选出保护动作的SER信号分析模块；

从筛选出的保护动作中提取出它们的名称的保护动作名称提取模块；

先将录波数据格式转换为计算机统一识别的数据格式，再根据自动提取出的保护名称或主动输入的保护名称查询出该保护动作的判据的数据解析和筛选模块；

通过录波扫描找到符合动作判据的录波段，并与查询出的保护动作的判据进行对比，做出保护是否正确动作的结论的录波数据分析模块；

通过鼠标事件机制将方程计算结果显示在信息栏的故障片段截取模块；

将分析的结果以WORD格式导出的故障报告输出模块，指出可能的故障点。

6.根据权利要求5所述的直流保护事故智能分析***，其特征在于，***的故障片段截取模块会自动截取一段满足判据的波形图，数据中包含瞬时值、有效值和每个公式的计算结果。