所属的技术人员能够理解，本发明的各个方面可以实现为系统、方法或程序产品。因此，本发明的各个方面可以具体实现为以下形式，即：完全的硬件实施方式、完全的软件实施方式(包括固件、微代码等)，或硬件和软件方面结合的实施方式，这里可以统称为“电路”、“模块”或“平台”。本发明再一个实施例中，提供一种考虑电流限幅的100％可再生能源系统稳定域估计系统，该系统能够用于实现上述考虑电流限幅的100％可再生能源系统稳定域估计方法，具体的，该考虑电流限幅的100％可再生能源系统稳定域估计系统包括构建模块、函数模块以及估计模块。其中，构建模块，建立考虑电流限幅的100％可再生能源系统模型；函数模块，基于虑电流限幅的100％可再生能源系统模型构建暂态能量函数；估计模块，基于暂态能量函数，根据区域寻找定理确定系统保守稳定域，确保电力系统不进入限流控制模式的最大区域。本发明再一个实施例中，提供了一种终端设备，该终端设备包括处理器以及存储器，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器用于执行所述计算机存储介质存储的程序指令。处理器可能是中央处理单元(central processingunit，cpu)，还可以是其他通用处理器、数字信号处理器(digital signal processor、dsp)、专用集成电路(application specific integrated circuit，asic)、现成可编程门阵列(field-programmable gate array，fpga)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等，其是终端的计算核心以及控制核心，其适于实现一条或一条以上指令，具体适于加载并执行一条或一条以上指令从而实现相应方法流程或相应功能；本发明实施例所述的处理器可以用于考虑电流限幅的100％可再生能源系统稳定域估计方法的操作，包括：建立考虑电流限幅的100％可再生能源系统模型；基于虑电流限幅的100％可再生能源系统模型构建暂态能量函数；基于暂态能量函数，根据区域寻找定理确定系统保守稳定域，确保电力系统不进入限流控制模式的最大区域。本发明再一个实施例中，本发明还提供了一种存储介质，具体为计算机可读存储介质(memory)，所述计算机可读存储介质是终端设备中的记忆设备，用于存放程序和数据。可以理解的是，此处的计算机可读存储介质既可以包括终端设备中的内置存储介质，当然也可以包括终端设备所支持的扩展存储介质，可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。计算机可读存储介质提供存储空间，该存储空间存储了终端的操作系统。并且，在该存储空间中还存放了适于被处理器加载并执行的一条或一条以上的指令，这些指令可以是一个或一个以上的计算机程序(包括程序代码)。需要说明的是，此处的计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、光纤、便携式紧凑盘只读存储器(cd-rom)、光存储器件、磁存储器件、或者上述的任一合适的组合。计算机可读存储介质还包括在基带中或者作为载波一部分传播的数据信号，其中承载了可读程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任一合适的组合。可读存储介质还可以是可读存储介质以外的任何可读介质，该可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。可读存储介质上包含的程序代码可以用任何适当的介质传输，包括但不限于无线、有线、光缆、rf等等，或者上述的任一合适的组合。可以以一种或多种程序设计语言的任一组合来编写用于执行本发明操作的程序代码，程序设计语言包括面向对象的程序设计语言—诸如java、c++等，还包括常规的过程式程序设计语言—诸如“c”语言或类似的程序设计语言。程序代码可以完全地在用户计算设备上执行、部分地在用户设备上执行、作为一个独立的软件包执行、部分在用户计算设备上部分在远程计算设备上执行、或者完全在远程计算设备或服务器上执行。在涉及远程计算设备的情形中，远程计算设备可以通过任一种类的网络，包括局域网(lan)或广域网(wan)，连接到用户计算设备，或者，可以连接到外部计算设备(例如利用因特网服务提供商来通过因特网连接)。可由处理器加载并执行计算机可读存储介质中存放的一条或一条以上指令，以实现上述实施例中有关考虑电流限幅的100％可再生能源系统稳定域估计方法的相应步骤；计算机可读存储介质中的一条或一条以上指令由处理器加载并执行如下步骤：建立考虑电流限幅的100％可再生能源系统模型；基于虑电流限幅的100％可再生能源系统模型构建暂态能量函数；基于暂态能量函数，根据区域寻找定理确定系统保守稳定域，确保电力系统不进入限流控制模式的最大区域。请参阅图5，终端设备为计算机设备，该实施例的计算机设备60包括：处理器61、存储器62以及存储在存储器62中并可在处理器61上运行的计算机程序63，该计算机程序63被处理器61执行时实现实施例中的储层改造井筒中流体组成计算方法，为避免重复，此处不一一赘述。或者，该计算机程序63被处理器61执行时实现实施例储层改造井筒中流体组成计算系统中各模型/单元的功能，为避免重复，此处不一一赘述。计算机设备60可以是桌上型计算机、笔记本、掌上电脑及云端服务器等计算设备。计算机设备60可包括，但不仅限于，处理器61、存储器62。本领域技术人员可以理解，图5仅仅是计算机设备60的示例，并不构成对计算机设备60的限定，可以包括比图示更多或更少的部件，或者组合某些部件，或者不同的部件，例如计算机设备还可以包括输入输出设备、网络接入设备、总线等。所称处理器61可以是中央处理单元(central processing unit，cpu)，还可以是其它通用处理器、中央处理器、图形处理器、数字信号处理器(digital signal processor，dsp)、专用集成电路(application specific integrated circuit，asic)、现场可编程门阵列(field-programmable gate array，fpga)或者其它可编程逻辑器件、分立门或者晶体管逻辑器件、基于量子计算的数据处理逻辑器、分立硬件组件等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。存储器62可以是计算机设备60的内部存储单元，例如计算机设备60的硬盘或内存。存储器62也可以是计算机设备60的外部存储设备，例如计算机设备60上配备的插接式硬盘，智能存储卡(smart media card，smc)，安全数字(secure digital,sd)卡，闪存卡(flash card)等。进一步地，存储器62还可以既包括计算机设备60的内部存储单元也包括外部存储设备。存储器62用于存储计算机程序以及计算机设备所需的其它程序和数据。存储器62还可以用于暂时地存储已经输出或者将要输出的数据。本技术所提供的各实施例中所使用的对存储器、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的至少一种。非易失性存储器可包括只读存储器(read-only memory，rom)、磁带、软盘、闪存、光存储器、高密度嵌入式非易失性存储器、阻变存储器(reram)、磁变存储器(magnetoresistive random access memory，mram)、铁电存储器(ferroelectric random access memory，fram)、相变存储器(phase change memory，pcm)、石墨烯存储器等。易失性存储器可包括随机存取存储器(random access memory，ram)或外部高速缓冲存储器等。作为说明而非局限，ram可以是多种形式，比如静态随机存取存储器(static random access memory，sram)或动态随机存取存储器(dynamic randomaccess memory，dram)等。本技术所提供的各实施例中所涉及的数据库可包括关系型数据库和非关系型数据库中至少一种。非关系型数据库可包括基于区块链的分布式数据库等，不限于此。本技术所提供的各实施例中所涉及的处理器可为通用处理器、中央处理器、图形处理器、数字信号处理器、可编程逻辑器、基于量子计算的数据处理逻辑器等，不限于此。请参阅图6，终端设备600为电子设备，电子设备以通用计算设备的形式表现。电子设备的组件可以包括但不限于：至少一个处理单元610、至少一个存储单元620、连接不同平台组件(包括存储单元620和处理单元610)的总线630、显示单元640等。其中，存储单元存储有程序代码，程序代码可以被处理单元610执行，使得处理单元610执行本说明书上述方法部分中描述的根据本发明各种示例性实施方式的步骤。例如，处理单元610可以执行如图1中所示的步骤。存储单元620可以包括易失性存储单元形式的可读介质，例如随机存取存储单元(ram)6201和/或高速缓存存储单元6202，还可以进一步包括只读存储单元(rom)6203。存储单元620还可以包括具有一组(至少一个)程序模块6205的程序/实用工具6204，这样的程序模块6205包括但不限于：操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。总线630可以为表示几类总线结构中的一种或多种，包括存储单元总线或者存储单元控制器、外围总线、图形加速端口、处理单元或者使用多种总线结构中的任一总线结构的局域总线。电子设备600也可以与一个或多个外部设备700(例如键盘、指向设备、蓝牙设备等)通信，还可与一个或者多个使得用户能与该电子设备600交互的设备通信，和/或与使得该电子设备600能与一个或多个其它计算设备进行通信的任何设备(例如路由器、调制解调器等等)通信。这种通信可以通过输入/输出(i/o)接口650进行。并且，电子设备600还可以通过网络适配器660与一个或者多个网络(例如局域网(lan)，广域网(wan)和/或公共网络，例如因特网)通信。网络适配器660可以通过总线630与电子设备600的其它模块通信。应当明白，尽管图中未示出，可以结合电子设备600使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、raid系统、磁带驱动器以及数据备份存储平台等。为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。通常在此处附图中的描述和所示的本发明实施例的组件可以通过各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。请参阅图2，图2为测试系统的拓扑图。请参阅图3，图3为系统故障后暂态能量小于本发明所得的临界能量的系统动态过程。在0.1秒时，母线9三相接地短路，故障持续35ms后被清除，故障清除瞬间系统的暂态能量为0.0619，小于求解优化问题后所得的临界能量0.0645，最后系统保持稳定。因此，采用本发明提出的稳定域寻找方法可以保证系统的暂态稳定性。请参阅图4，图4为系统故障后暂态能量大于本发明所得的临界能量的系统动态过程。在0.1秒时，母线9三相接地短路，故障持续43ms后被清除，故障清除瞬间系统的暂态能量为0.0893，大于求解优化问题后所得的临界能量0.0645，最后系统失稳。请参阅表1，表1为验证所提控制策略所采用的系统参数。表1综上所述，本发明一种考虑电流限幅的100％可再生能源系统稳定域估计方法及系统，通过建立100％可再生能源的系统模型，构造暂态能量函数，最后求解优化问题得到一个保证系统暂态稳定的稳定域。所属领域的技术人员可以清楚地了解到，为了描述的方便和简洁，仅以上述各功能单元、模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能单元、模块完成，即将所述装置的内部结构划分成不同的功能单元或模块，以完成以上描述的全部或者部分功能。实施例中的各功能单元、模块可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中，上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。另外，各功能单元、模块的具体名称也只是为了便于相互区分，并不用于限制本技术的保护范围。上述系统中单元、模块的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述或记载的部分，可以参见其它实施例的相关描述。本领域普通技术人员可以意识到，结合本发明中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。在本发明所提供的实施例中，应该理解到，所揭露的装置/终端和方法，可以通过其它的方式实现。例如，以上所描述的装置/终端实施例仅仅是示意性的，例如，所述模块或单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通讯连接可以是通过一些接口，装置或单元的间接耦合或通讯连接，可以是电性，机械或其它的形式。所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。所述集成的模块/单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，所述计算机程序包括计算机程序代码，所述计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。所述计算机可读介质可以包括：能够携带所述计算机程序代码的任何实体或装置、记录介质、u盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(read-only memory，rom)、随机存取存储器(randomaccess memory，ram)、电载波信号、电信信号以及软件分发介质等，需要说明的是，所述计算机可读介质包含的内容可以根据司法管辖区内立法和专利实践的要求进行适当的增减，例如在某些司法管辖区，根据立法和专利实践，计算机可读介质不包括是电载波信号和电信信号。本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。以上内容仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明权利要求书的保护范围之内。
背景技术：
1、近年来，由传统发电技术带来的环境问题严峻，推动了能源结构改革；以风电和光伏为代表的可再生能源在电力系统中的比例不断增加，预计到2030年将达到56.1％。新型电力系统的建设是以非同步电源为特征的新能源比重不断上升的过程，大量的变流器电源将打破传统同步发电机的主导地位，深刻影响电力系统的动态特性和稳定机理。电压源型变流器在大规模可再生能源并网方面具有独特的优势。
2、目前，随着可再生能源的大量接入，电力电子换流器型电源占比不断提升。构网型换流器比跟网型转换器越来越受欢迎，因其可在电力系统需要时自发提供惯性或频率支撑。但是构网型换流器耐受电流能力有限，需要采用额外的限流措施来保障换流器安全运行。
3、在100％可再生能源电力系统中，换流器的限流控制会使大干扰下的系统成为一个切换系统，其动态特性与传统同步机主导的系统有很大的不同，这导致系统稳定域更加难以得到。
技术实现思路
1、本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种考虑电流限幅的100％可再生能源系统稳定域估计方法及系统，用于解决寻找考虑电流限幅的100％可再生能源系统稳定域的技术问题。
2、本发明采用以下技术方案：
3、一种考虑电流限幅的100％可再生能源系统稳定域估计方法，包括以下步骤：
4、建立考虑电流限幅的100％可再生能源系统模型；
5、基于虑电流限幅的100％可再生能源系统模型构建暂态能量函数；
6、根据区域寻找定理确定系统保守稳定域，确保电力系统不进入限流控制模式的最大区域。
7、优选地，考虑电流限幅的100％可再生能源系统模型具体为：
8、
9、其中，ωi为第i台vsg的转速；ω0为coi的转速，pi为第i个vsg输出的有功功率，ωci为第i台vsg的实际转速，tji为第i台vsg的惯量，pisp为第i台vsg的参考功率，pcoi为coi的功率。
10、优选地，第i个vsg输出的有功功率pi推导如下：
11、
12、其中，|ici|≤imax，gii为第i个vsg的自电导，ui为第i个vsg的电压，cij为第i个vsg与第j个vsg的互电导，θij为第i个vsg与第j个vsg的相角差，dij为第i个vsg与第j个vsg的互电纳，n为vsg的总数。
13、优选地，暂态能量函数具体为：
14、
15、其中，为暂态能量函数，tji为第i台vsg的惯量，pi0为第i台vsg的参考功率，θci为第i台vsg的相角，θsi为稳定平衡点，cij为第i个vsg与第j个vsg的互电导，为coi下第i个vsg与第j个vsg的相角差，为第i个vsg与第j个vsg的相角差，为常数，dij为第i个vsg与第j个vsg的互电纳，n为vsg的总数，θ为vsg的相角，为vsg的角速度。
16、优选地，第i台vsg的参考功率pi0为：
17、
18、其中，pisp为第i台vsg的参考功率，ui为第i个vsg的电压，gii为第i个vsg的自电导。
19、优选地，常数为：
20、
21、其中，θci为第i台vsg的相角，θcj为第j台vsg的相角，θsi为第i台vsg的实际相角，θsj为第j台vsg的实际相角，为coi下第i个vsg与第j个vsg的相角差。
22、优选地，根据区域寻找定理确定系统保守稳定域具体为：
23、当所有等式约束的拉格朗日乘子小于零或关于等式约束的目标函数值的求导大于零时，对于的任意点，存在一个上界vcr，如果则具有并且vcr＝min{vcri,i∈n}，d为系统保守稳定域。
24、优选地，通过求解以下n个子问题得到vcr如下：
25、vcri(θ)＝minvp(θ)
26、
27、其中，vcri(θ)为最优暂态能量，vp(θ)为暂态能量，ici为第i台vsg的电流，imax为电流极限值，ick为第k台vsg的电流，ykj为第i个vsg与第k个vsg的互导纳，uj为第j台vsg的电压，θj为第j台vsg的相角，k为第k台vsg，n为vsg的总数。
28、优选地，电力系统包括n台虚拟同步发电机，负荷采用恒阻抗模型，简化网络中消除所有负载总线。
29、第二方面，本发明实施例提供了一种考虑电流限幅的100％可再生能源系统稳定域估计系统，包括：
30、构建模块，建立考虑电流限幅的100％可再生能源系统模型；
31、函数模块，基于虑电流限幅的100％可再生能源系统模型构建暂态能量函数；
32、估计模块，基于暂态能量函数，根据区域寻找定理确定系统保守稳定域，确保电力系统不进入限流控制模式的最大区域。
33、第三方面，一种计算机设备，包括存储器、处理器以及存储在所述存储器中并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述考虑电流限幅的100％可再生能源系统稳定域估计方法的步骤。
34、第四方面，本发明实施例提供了一种计算机可读存储介质，包括计算机程序，所述计算机程序被处理器执行时实现上述考虑电流限幅的100％可再生能源系统稳定域估计方法的步骤。
35、与现有技术相比，本发明至少具有以下有益效果：
36、一种考虑电流限幅的100％可再生能源系统稳定域估计方法，建立了一个考虑电流限幅的100％可再生能源系统模型。进而在系统模型的基础上构建暂态能量函数和优化模型。最后根据区域寻找定理提出了一种考虑电流限幅的100％可再生能源系统保守稳定域估计方法以找到一个使系统不进入限流控制模式的最大区域；通过求解优化模型，找到了一个使系统不进入限流控制模式的最大区域，有利于提升100％可再生能源系统安全稳定运行能力。
37、进一步的，通过设置考虑电流限幅的100％可再生能源系统模型，有利于后续暂态能量函数的建立。
38、进一步的，设置第i个vsg输出的有功功率pi为后续暂态稳定性分析提供理论基础。
39、进一步的，设置暂态能量函数为后续构造优化问题进而求解稳定域奠定理论基础。
40、进一步的，根据区域寻找定理确定系统保守稳定域，能够实现系统在所提稳定域内不进入限流模式，有利于提升系统的暂态稳定性。
41、可以理解的是，上述第二方面的有益效果可以参见上述第一方面中的相关描述，在此不再赘述。
42、综上所述，相对与传统同步发电机主导的电力系统，本发明理论清晰，有利于寻找考虑电流限幅的100％可再生能源电力系统的保守稳定域，对系统安全稳定运行提供了理论帮助。
43、下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
1.一种考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，考虑电流限幅的100％可再生能源系统模型具体为：
3.根据权利要求2所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，第i个vsg输出的有功功率pi推导如下：
4.根据权利要求1所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，暂态能量函数具体为：
5.根据权利要求4所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，第i台vsg的参考功率pi0为：
6.根据权利要求4所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，常数为：
7.根据权利要求1所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，根据区域寻找定理确定系统保守稳定域具体为：
8.根据权利要求7所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，通过求解以下n个子问题得到vcr如下：
9.根据权利要求1所述的考虑电流限幅的100％可再生能源系统稳定域估计方法，其特征在于，电力系统包括n台虚拟同步发电机，负荷采用恒阻抗模型，简化网络中消除所有负载总线。
10.一种考虑电流限幅的100％可再生能源系统稳定域估计系统，其特征在于，包括：
