具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明方法包括以下步骤:
1)在接收站内建立一LNG动态现货预测调度系统
动态现货预测调度系统包括:一DCS1、一LNG动态现货预测调度仿真系统平台2和一GMS终端3。其中,DCS1和GMS终端3在接收站为已有设置,仿真系统平台2通过网络连接接收站的DCS1和GMS终端3。仿真系统平台2包括初始参数数据库21、仿真模块22、指标和报表显示模块23和分析数据库24。仿真系统平台2通过网络连接接收站的DCS1和GMS终端3。
仿真模块22内置有综合了蒙特卡洛算法和运筹学排队论的离散仿真推演过程。其中蒙特卡洛算法是以概率和统计理论方法为基础,把求解的问题同恰当的概率模型相联系,用计算机产生大量的随机数或伪随机数进行统计模拟或抽样,以获得近似解;蒙特卡洛求解算法是计算程序中生成随机数的方法,是公开且广泛应用的技术。运筹学排队论是通过对服务对象到来及服务时间的统计研究,得出包括等待时间、排队长度、忙期长短等数量指标的统计规律;运筹学排队论是排队规则和随机分布函数的设置,这种函数本身是公开且广泛应用的技术。离散仿真方法是综合了运筹学排队论和蒙特卡洛求解算法的具体计算机程序,通过该计算机程序,对各个元素进行细致、完整地描述,对各个元素的操作规则采用同步并行和依照条件按时间顺序两种模式,完全按照客观的实际操作处理离散流和连续流的元素活动,离散仿真方法,即按照连续、均匀间隔的时间序列,根据实际的工作流程和活动顺序,在规定时间内顺序地改变上述元素和实体的工作状态,也是属于公开且广泛应用的技术。本发明根据需要解决问题的具体情况、不同的描述对象和相互之间的对应关系,将上述各种方法引入本发明,并将其具体地应用在各种模型的仿真推演过程中。
2)初始参数数据库21中的原始信息
①初始参数数据库21内预先存储有长期ADP,长期供货点至接收站的往返双向航程航道,执行长期ADP的运输船,长期供货点内部各类设施及附属港口、码头,接收站内部各类设施及附属港口、码头,不同类型运输船进出接收站港口的操作和准备时间,长期供货点与接收站港口和航道的恶劣天气、潮汐条件和日间操作时间限制等信息。
②初始参数数据库21中输入有涵盖工业用户、城镇用户和燃气电厂三类气态用户,以及槽车液态用户在内的下游用户数据信息,包括接收站下游各个用户提出的用气需求报表和波动幅度数据,槽车液态外运和天然气外输用户的年、月、日、小时需求预测等信息。
③DCS1通过网络采集接收站内每一座储罐的初始液位和总罐容信息,并将采集的数据输入初始参数数据库21;在初始参数数据库21中还预存有储罐动态监控指标,储罐动态监控指标包括(如图2所示):
a、逼罐比率参数,其取值范围一般为储罐容量的85%~95%。
逼罐比例参数用于表明储罐内LNG充装量占储罐有效容积的高限百分比,当运输船到达接收站的外港时,如果储罐达到逼罐比率的总空间减去现有LNG罐容的剩余空间,小于运输船的有效船容,则不允许运输船进内港锚地和码头泊位,亦即产生“逼罐”问题。
b、高警戒液位参数,其取值范围一般为储罐容量的85%~90%。
高警戒液位参数类似于逼罐比率参数,但不是功能性参数,可以取与逼罐比率相同数值,仅用于储罐操作运行过程中的高液位报警。
c、低警戒液位参数,其取值范围一般为储罐容量的8%~15%。
低警戒液位参数用于储罐操作运行过程中的低液位报警。
d、警戒罐容满足下游需求临界值参数,其取值范围一般为48小时~72小时。
警戒罐容满足下游需求临界值参数=(储罐的当前剩余罐容-低警戒液位参数对应罐容)/当前日平均小时用气量,该参数以小时为单位进行不间断的数据比对。警戒罐容满足下游需求临界值参数与低警戒液位参数配合,对下游用气需求优先级保供问题进行处理。当储罐剩余库存达到低警戒液位参数,又不足以维持警戒罐容满足下游需求临界值的时长时,如果接收站码头、内港泊位没有运输船卸货,则接收站需按照下游用户用气优先级逐级关停外输泵,停止外输供气(以下简称“保供”),关停顺序为:先槽车,其次燃气电厂用户,最后工业用户,仅保供城镇用户用气。
e、请求现货液位参数,其时间取值范围一般为72小时~150小时。
请求现货液位(维持时长)参数=(储罐当前剩余罐容-低警戒液位参数对应罐容)/当前日平均小时用气量,该参数以小时为单位进行不间断的数据比对。请求现货液位参数用于实现动态现货请求功能。当储罐当前剩余罐容与低警戒液位参数对应罐容的差值,按照当前时刻下游用户用气速率所能维持的时长小于请求现货液位参数时,系统发出一次动态现货请求。
上述各种参数的数值范围也可以根据实际运行经验和需要选取其它值。
3)初始参数数据库21中的信息更新
本发明可以根据需要通过仿真系统平台2启动仿真模块22进行自动仿真推演,所谓需要可以是一年,或更长、或更短时间;进行仿真推演之前,需要先进行初始参数数据库21的信息更新,更新内容包括:
①仿真系统平台2通过接收站的DCS1对每一座储罐的液位信息、总罐容信息和期ADP执行信息等进行实时监测,并对初始参数数据库21进行更新,其中,长期ADP的起始时间可以设定为最近一次运输船完成交付LNG的时间;
②港口和航道的恶劣天气信息等根据最新天气预报资料手动更新;
③槽车液态外运和天然气外输用户的年、月、日、小时需求预测信息可以根据近一年期的历史数据报表、当前签订的购气合同、下游用户最新提出的用气需求和以小时为单位的流量波动幅度等数据手动更新。
4)仿真系统平台2启动仿真模块22进行仿真推演
仿真模块22的仿真推演程序包括以下步骤(如图3所示):
①仿真模块22读入初始参数数据库21更新后的数据信息,给每一个元素赋值,扫描指定时间(通常是指未来三个月、一年或更长时间)内接收站的长期ADP,如果长期ADP为空,则终止仿真推演程序;如果有长期ADP,则进入下一步;
②仿真模块22根据储罐动态监控指标,分析除了长期ADP外,是否还启动了动态现货请求:
如果未启动动态现货请求,则进入步骤③;
如果已启动动态现货请求,则进入步骤④;
③执行长期ADP模式的仿真推演,具体的仿真推演过程包括:
a、在长期供货点内部,LNG物料按照从气田、天然气处理装置、天然气液化装置、储罐、装船设施至运输船满船的顺序流动;与此同步,选定的运输船空船按照进入外港、内港泊位、连接装船设施的操作顺序活动。其中,天然气处理装置和天然气液化装置会考虑一定比例的损耗,在进入外港和内港泊位的作业中,会受到恶劣天气和潮汐条件等不确定性元素的影响,LNG装船速率会有一定的波动性。
b、在接收站内部,LNG物料按照从卸船设施、储罐、BOG(蒸发气)处理设施至LNG外输泵,再从LNG外输泵、LNG气化器、通过外输管道至工业用户、城镇用户和燃气电厂三类气态用户,以及LNG外输泵至槽车液态外运的顺序流动。在流动过程中,根据储罐动态监控指标作相应处理,即在卸船设施与储罐之间判断是否“逼罐”,在储罐与LNG外输泵之间判断是否“保供”和是否“空罐”,处理结果包括:
如果运输船卸料时,储罐剩余罐容不足以存储整船的LNG物料,产生“逼罐”现象,则仿真推演结果中输出一“接收站逼罐明细”报表,记录接收站的ID编号、发生时刻、结束时刻、运输船的ID编号、船型等明细数据,同时还可以记录逼罐船数、平均逼罐时长和总逼罐时长,运输船须等待直至“逼罐”现象消除方能开始卸料;
如果储罐已经“空罐”(到达低警戒液位参数),则仿真推演结果中输出一“站点低液位预警”报表,记录发生低警戒液位报警接收站的ID编号、报警开始时刻、结束时刻等明细数据;
如果储罐已经到达“警戒罐容满足下游需求临界值参数”,需要启动不同等级的“保供”,则仿真推演结果中输出一“站点空罐时段记录”和“站点保民用关停报表”,记录报警时间及空罐时间等详细的运行数据。LNG外输泵须按照“分级关停控制速率”的原则,控制LNG外输泵逐次关停与槽车、燃气电厂用户、工业用户、直至城镇用户之间连接的外输管道。
与此同步,运输船满船按照进入外港、内港泊位、连接卸船设施的操作顺序活动。其中,在内港泊位后可以根据码头数量和功能判断是否同时作业:如仅一个码头则单艘运输船卸料;如两个码头则两艘运输船同时卸料。在进入外港和内港泊位的作业中,会受到恶劣天气和潮汐条件等不确定性元素,以及日间操作时间限制的影响。
日间操作时间限制是指运输船满船只能在日间规定时间内从外港进入内港泊位。在进港和靠泊作业的时间顺序中,产生恶劣天气的不确定性影响和日间限制的规则限制,仅在日间限制时间范围内,允许运输船进港靠泊。其中,接收站的槽车装车速率、运输船装卸船速率存在一定的波动性;接收站的外输泵、气化器外输设备可靠性存在不确定性。
c、接收站的下游用户的用气需求,以设定的时间间隔(比如30分钟)不断地发给接收站,由接收站的LNG气化器和LNG外输泵给予匹配的供应;其间如果接收站供应能力达不到,则按照供应设备的瓶颈能力或限定的LNG供给数量供应,即按照储罐动态监控指标中警戒罐容满足下游需求临界值中的描述执行。其中,下游用户用气的小时用气速率会存在一定的波动性。
d、执行长期ADP的运输船,在长期供货点至接收站的双向航道上来回往返。在抵达长期供货点和离开接收站时,运输船的状态为“空船”;在离开长期供货点和抵达接收站时,运输船的状态为“满船”。
e、输出执行长期ADP模式的仿真推演结果,以用于后期长期ADP的分析。按照本发明方法得到的仿真推演结果会与预设的长期ADP计划存在偏差,属于正常情况,仿真推演结果用于纠正预设长期ADP计划的不合理之处。但仿真推演结果中不会出现动态现货指标,仿真程序结束。
④执行长期ADP与动态现货请求交互模式的仿真推演
仿真模块22在模拟长期ADP的同时,还要模拟根据下游用户的实时用气需求,基于储罐罐容发出的动态现货请求;模拟由现货点和现货运输船予以响应,并根据实际可选的现货点条件和运输航道,预排调度动态现货及动态现货运输船去现货点取货的过程;具体的仿真推演过程包括:
首先执行上述步骤③长期ADP模式的仿真推演步骤a~d;然后执行下述步骤:
e、在进行仿真推演的过程中,根据长期ADP的执行情况和下游用户用气的实际需求,仿真模块22按照实时采集的储罐液位信息与储罐动态监控指标的对比结果,发出动态现货请求,动态现货请求的实现和满足,需要一定的过程和时间,由于系统在发出一次动态现货请求直至其满足的过程中,还会不断产生其它的动态现货请求,因此可以设定仿真模块22只有在一次动态现货请求得到满足后,才会根据最新情况启动新一轮的仿真推演,期间对于其它的动态现货请求不做任何响应。
f、动态现货发货点根据动态现货请求,只确定动态现货的交货量和交货时间,而动态现货发货点是不确定位置和航道、也不体现内部流程的。由于现实中动态现货点不定,相应的航道属性也难以设定,且因实际预留动态现货的准备周期多为2~3个月,因此,不再考虑动态现货点与接收站之间的特定航道问题。
g、执行动态现货运输任务的运输船在动态现货点与接收站之间来回往返;在抵达现货点和离开接收站时,现货运输船的状态为“空船”;在离开现货点和抵达接收站时,现货运输船的状态为“满船”。
在每一次启动仿真模块22进行自动仿真推演时,为了提高仿真推演的准确性,可以设定重复推演多次,比如50次,或更多、或更少次,多次计算后给出均一化的平均结果。
在上述步骤③、步骤④的仿真推演过程中,本发明优选均匀随机分布函数模拟槽车装车速率、运输船装卸船速率,矩形随机分布函数模拟港口恶劣天气影响,正态随机分布函数模拟LNG外输泵和LNG气化器外输设备的可靠性,负指数随机分布函数模拟运输船到港时间延误和槽车灌装到达时间,负指数随机分布函数模拟运输船和槽车的到达时间间隔,泊松随机分布函数模拟运输船和槽车在指定时间内的到达数量,均匀随机分布函数模拟下游用户用气的小时用气速率波动性;同时还采用等待制的排队规则,先到先服务。上述各种函数的使用方法均是公开的,本发明仅是将其引入本发明各种仿真模式的推演过程中。另外在上述各步骤中也可以采用上述指定函数方法之外的其它方法进行模拟推演,在此不再赘述。
5)指标和报表显示模块23显示输出仿真推演结果
仿真推演结束后,在指标和报表显示模块23上显示输出执行长期ADP与动态现货请求交互模式的仿真推演结果,仿真推演结果包括:运输船信息,长期供货点与接收站之间的对应关系,储罐高-低液位与逼罐-空罐报警的时空属性统计和费用统计报表,以及指定时间(未来三个月至一年,或更长时间,取决于给定的仿真时间长度)内的动态现货批次、时间窗口分布等动态现货指标。
6)GMS终端3生成LNG动态现货的调度计划
仿真系统平台2将仿真推演结果及其对应的初始参数数据库21信息,存入分析数据库24,并将指定时间内动态现货需求的数量和时间窗口,通过网络上传至GMS终端3,GMS终端3据此生成接收站内的动态现货的调度计划。这种经过仿真推演预测的动态现货调度计划与原有的长期ADP一起,共同构成可靠、完整的ADP信息,从而能够很好地满足下游用气需求。
上述实施例中,在步骤5)与步骤6)之间可以增加一验证仿真推演结果的步骤,如果选择此步骤,则采用以下方法:
将步骤5)仿真推演结果中的动态现货时间窗口,分布转化成现货ADP,按照时间分布插入原有的长期ADP,并据此更新初始参数数据库21的信息,然后进行以验证为目的的仿真推演;
如果仿真推演结果中不再产生动态现货指标,则可以确认步骤5)的仿真推演结果正确,进入步骤6);
否则返回步骤2),调整初始参数数据库21的信息(主要是检查一下是否有输入性错误或进行一些有关功能性的修改,比如动态现货用运输船的船容等)后,启动仿真模块22进行新一轮的仿真推演。
上述分析数据库24中包含每次仿真推演结果,特别是经过确认后的初始参数列表和对应的长期供货点、接收站、现货点和运输船的联运联供指标,可以用于人工回顾和历史数据统计分析,形成人工分析经验。
上述实施例,仅仅是举例描述一个接收站的情况,在实际操作中,可以是多个接收站共同使用一个仿真系统平台2,即将仿真系统平台2设置在其中一个接收站内,通过网络分别连接多个接收站的DCS1和GMS终端3,对其中任意一接收站进行动态现货的仿真模拟。
上述各实施例仅用于说明本发明,其中一些参数和步骤的设置是可以有所变化的,特别是某些数据的设置内容、设置范围、采用的方法等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。