CN110635471A - 基于同时段功率和电量的台区电网理论线损计算方法 - Google Patents
基于同时段功率和电量的台区电网理论线损计算方法 Download PDFInfo
- Publication number
- CN110635471A CN110635471A CN201910711520.6A CN201910711520A CN110635471A CN 110635471 A CN110635471 A CN 110635471A CN 201910711520 A CN201910711520 A CN 201910711520A CN 110635471 A CN110635471 A CN 110635471A
- Authority
- CN
- China
- Prior art keywords
- power
- loss
- node
- phase
- endpoint
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Power Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了基于同时段功率和电量的台区电网理论线损计算方法,主要步骤为:1)获取台区电网基本数据。2)基于台区电网基本数据,建立台区电网端点导纳矩阵Y。3)计算平均有功功率和平均无功功率4)对台区电网进行三相潮流计算,得到电网首断面有功损耗Ploss1、末断面有功损耗Ploss2和平均有功损耗5)分别计算首断面有功损耗Ploss1、末断面有功损耗Ploss2和平均有功损耗对应的电量损耗ΔWloss1、电量损耗ΔWloss2、电量损耗和台区电网理论线损ΔWloss。本发明同时应用电量数据和功率数据,计算时间段time内的理论线损,能够更加精确地计算台区电网理论线损,为降损措施的制定等提供科学依据与理论支撑。
Description
技术领域
本发明涉及电力系统调度自动化领域,具体是基于同时段功率和电量的台区电网理论线损计算方法。
背景技术
台区电网是线损管理的“重灾区”,是导致电网总体线损率高、达标率低的关键因素之一,其不仅具备巨大的降损潜力,更是电网线损管理与降损工作的重中之重。线损可分为理论线损和管理线损两部分,其中,理论线损是由相应时段内运行参数和设备参数决定,是电网实际运行中线路和变压器所消耗的理论电量值。管理线损主要由窃电、台区隶属关系混乱或计量设备故障等原因引起,可通过采取必要的组织管理措施予以避免或减少。
380/220V台区电网典型日理论线损的精确计算不仅需要电网的拓扑信息及其阻抗导纳参数,而且还需要电源与终端负荷全天24小时的三相功率与电压数据。实际电网中,电源与终端负荷的三相电力数据以及台区电网的结构参数都很难同时获取,因此,电力部门的典型日理论线损计算一般都采用近似方法,且台区电网的有功功率数据大多基于电量数据进行近似,无功数据凭经验估计取值。理论线损的计算常采用最大负荷损耗时间法、代表日均方根电流法和等效电阻法等给出代表日电网的理论线损指标。但明显地,以上方法所得结果不够精确,不能真实反应负荷和电网运行方式变化对理论线损的影响。
近年来,随着智能配电网建设的不断发展和持续推进,具有三相功率、电流与电压及电量数据采集和双向通讯功能的智能电表在我国配电网中得到初步的推广应用,这为台区电网的理论线损计算提供了必要的实时数据。
发明内容
本发明的目的是解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,基于同时段功率和电量的台区电网理论线损计算方法,主要包括以下步骤:
1)获取台区电网基本数据。
进一步,所述台区电网基本数据主要包括台区电网参数、台区电网结构和用户智能电表的量测数据。
所述台区电网参数主要包括元件电阻、元件电抗、元件电纳和元件额定电压。
所述台区电网结构为网络拓扑连接关系。
2)基于台区电网基本数据,建立台区电网端点导纳矩阵Y。
式中,为端点导纳矩阵中的元素。自导纳为与端点i的节点d直接相联的所有支路导纳之和。互导纳为端点i的节点d和节点t之间支路导纳的相反数与跟端点i相连的端点的节点d和节点t之间所有支路导纳之和。互导纳为端点i和端点k中节点d和t之间支路导纳的相反数。φi为不包括端点i且与端点i直接相连端点的集合。为端点i节点d、节点X间并联支路导纳。为端点i节点d、节点t间并联支路导纳。为端点i至端点j支路节点d、节点t间的支路互导纳元素,且 为端点i至端点k支路节点d、t间并联支路导纳元素,且B1={a,b,c,n},表示三相与中性点的集合。
首末断面总时间段长度time内平均功率因素如下所示:
进一步,对台区电网进行三相潮流计算的主要步骤如下:
4.1)将配电变压器低压侧端点作为平衡端点S。其中,平衡端点S的中性点作为零电位参考点,平衡端点S的A相节点作为零相位参考点。设定平衡端点三相电压幅值为额定相电压幅值,三相相角分别相差120°。
4.2)设定所有端点的三相电压幅值初始值为台区电网的额定相电压。任意端点中,3个三相电压相角两两相差120°,中性点电压为0。
4.3)初始化最大迭代次数Tmax,并设置收敛精度ε。设定初始迭代次数iter=1。
4.4)计算三相潮流计算中任意量测变量的残差,即:
r=z-h(x)。 (7)
式中,z为m维的量测变量。m为量测变量个数。x为n维的状态变量。n为状态变量个数。m=n。h(x)为量测方程。量测方程h(x)为包含零注入端点的节点注入功率对应的等效注入电流量测方程。r为量测残差。
其中,包含零注入端点的节点注入功率对应的等效注入电流量测方程如下所示:
式中,为d相节点的等效注入电流,d∈Bp,Bp={a,b,c}表示三相节点的集合。为中性点的等效注入电流。为包括端点i且与端点i直接相连端点的集合。B1={a,b,c,n}为三相节点和中性点的集合。为端点导纳矩阵中端点i中d相节点和端点k中t相节点的导纳元素。为端点导纳矩阵中端点i中中性点和端点k中t相节点的导纳元素。为端点k中t相节点电压相量。为端点i处d相对于中性点的注入功率,且 为端点i处d相对中性点的发电机功率。为端点i处d相对中性点的负荷功率。为端点i的d相节点电压相量。为端点i的中性点电压相量。
4.5)建立用于三相潮流计算的雅克比矩阵H,主要分为以下两种情况:
式中,Hik为雅可比矩阵的子矩阵,且k≠i。为端点i中节点t和端点k中节点δ之间支路电导的相反数,δ∈B1,t∈B1,B1={a,b,c,n}。为端点i中节点t和端点k中节点δ之间支路电纳的相反数,δ∈B1,t∈B1,B1={a,b,c,n}。为电压实部,为电压虚部。
式中,Hii为雅可比子矩阵矩阵。为端点i中节点t和δ之间支路电导的相反数与跟端点i相连的端点中节点t和δ之间所有支路电导之和,δ∈B1,t∈B1,B1={a,b,c,n}。为端点i中节点t和δ之间支路电纳的相反数与跟端点i相连的端点中节点t和δ之间所有支路电纳之和。
其中,矩阵ΔHii(1:8,1:8)如下所示:
式中,为端点i中d相节点电压的实部,d∈Bp,Bp={a,b,c}。fi d端点i中d相节点电压的虚部。为端点i中中性点电压的实部。fi n端点i中中性点电压的虚部。Pi dn为端点i中d相对于中性点的有功功率。为端点i中d相对于中性点的无功功率。
4.6)计算状态变量的修正量Δx(iter),即:
Δx(iter)=inv(H(iter))*r(iter)。 (14)
式中,iter为计算迭代次数。x(iter)为第iter次迭代时的状态变量。为第iter次迭代时量测变量的雅可比矩阵。inv为求逆运算。r(iter)=z-h(x(iter))为迭代值为x(iter)时的残差。
4.7)更新状态变量,得到状态变量迭代值x(iter+1)=x(iter)+Δx(iter),iter=iter+1,直至满足收敛标准max(Δx(iter))<ε为止。满足收敛标准的状态变量迭代值记为x'。
4.8)根据状态变量x',计算全网潮流和有功损耗Ploss,有功损耗Ploss计算公式如下:
5)分别计算首断面有功损耗Ploss1、末断面有功损耗Ploss2和平均有功损耗对应的电量损耗ΔWloss1、电量损耗ΔWloss2和电量损耗并根据电量损耗ΔWloss1、电量损耗ΔWloss2和电量损耗计算得到台区电网理论线损ΔWloss。
进一步,计算台区电网理论线损ΔWloss的主要步骤如下:
ΔWloss1=Ploss1*1/3time。 (16)
ΔWloss2=Ploss2*1/3time。 (17)
5.2)计算台区电网理论线损ΔWloss,即:
式中,ΔWloss为总时间段长度time内基于同时段功率和电量计算的理论线损。
本发明的技术效果是毋庸置疑的。本发明计算理论线损时,同时应用到了时间段time内的电量数据和功率数据。本发明虽然针对的是时间段time,但类似的对每个时间段time进行计算,即可累计计算得到日理论线损和月度理论线损,故而本发明也适用于日理论线损计算和月度理论线损计算。
本发明同时应用电量数据和功率数据,计算时间段time内的理论线损,能够更加精确地计算台区电网理论线损,为降损措施的制定等提供科学依据与理论支撑。
附图说明
图1为基于同时段功率和电量的台区电网理论线损计算方法的流程示意图;
图2为同时段功率和电量关系示意图;
图3为台区电网三相潮流计算流程示意图;
图4为IEEE-13节点配电系统的接线图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
参见图1,基于同时段功率和电量的台区电网理论线损计算方法,主要包括以下步骤:
1)获取台区电网基本数据。
进一步,所述台区电网基本数据主要包括台区电网参数、台区电网结构和用户智能电表的量测数据。
所述台区电网参数主要包括元件电阻、元件电抗、元件电纳和元件额定电压。
所述台区电网结构为网络拓扑连接关系。
2)基于台区电网基本数据,建立台区电网端点导纳矩阵Y。
式中,为端点导纳矩阵中的元素。自导纳为与端点i的节点d直接相联的所有支路导纳之和。互导纳为端点i的节点d和节点t之间支路导纳的相反数与跟端点i相连的端点的节点d和节点t之间所有支路导纳之和。互导纳为端点i和端点k中节点d和t之间支路导纳的相反数。φi为不包括端点i且与端点i直接相连端点的集合。为端点i节点d、节点X间并联支路导纳。为端点i节点d、节点t间并联支路导纳。为端点i至端点j支路节点d、节点t间的支路互导纳元素,且 为端点i至端点k支路节点d、t间并联支路导纳元素,且B1={a,b,c,n},表示三相与中性点的集合。
3)计算平均有功功率和平均无功功率d∈Bp,Bp={a,b,c}。a、b、c分别表示a相节点、b相节点和c相节点。
进一步,对台区电网进行三相潮流计算的主要步骤如下:
4.1)将配电变压器低压侧端点作为平衡端点S。其中,平衡端点S的中性点作为零电位参考点,平衡端点S的A相节点作为零相位参考点。设定平衡端点三相电压幅值为额定相电压幅值,三相相角分别相差120°。
4.2)设定所有端点的三相电压幅值初始值为台区电网的额定相电压。任意端点中,3个三相电压相角两两相差120°,中性点电压为0。
4.3)初始化最大迭代次数Tmax,并设置收敛精度ε。设定初始迭代次数iter=1。
4.4)计算三相潮流计算中任意量测变量的残差,即:
r=z-h(x)。 (7)
式中,z为m维的量测变量。m为量测变量个数。x为n维的状态变量。n为状态变量个数。m=n。h(x)为量测方程。量测方程h(x)为包含零注入端点的节点注入功率对应的等效注入电流量测方程。r为量测残差。
其中,包含零注入端点的节点注入功率对应的等效注入电流量测方程如下所示:
式中,为d相节点的等效注入电流,d∈Bp,Bp={a,b,c}表示三相节点的集合。为中性点的等效注入电流。为包括端点i且与端点i直接相连端点的集合。B1={a,b,c,n}为三相节点和中性点的集合。为端点导纳矩阵中端点i中d相节点和端点k中t相节点的导纳元素。为端点导纳矩阵中端点i中中性点和端点k中t相节点的导纳元素。为端点k中t相节点电压相量。为端点i处d相对于中性点的注入功率,且 为端点i处d相对中性点的发电机功率。为端点i处d相对中性点的负荷功率。为端点i的d相节点电压相量。为端点i的中性点电压相量。
4.5)建立用于三相潮流计算的雅克比矩阵H,主要分为以下两种情况:
I)对台区电网中的任意端点i,当状态变量为其它端点三相和中性点电压实部和虚部且k≠i时,雅可比矩阵的子矩阵Hik(1:8,1:8)如下所示:
式中,Hik为雅可比矩阵的子矩阵,且k≠i。为端点i中节点t和端点k中节点δ之间支路电导的相反数,δ∈B1,t∈B1,B1={a,b,c,n}。为端点i中节点t和端点k中节点δ之间支路电纳的相反数,δ∈B1,t∈B1,B1={a,b,c,n}。
式中,Hii为雅可比子矩阵矩阵。为端点i中节点t和δ之间支路电导的相反数与跟端点i相连的端点中节点t和δ之间所有支路电导之和,δ∈B1,t∈B1,B1={a,b,c,n}。为端点i中节点t和δ之间支路电纳的相反数与跟端点i相连的端点中节点t和δ之间所有支路电纳之和。
其中,矩阵ΔHii(1:8,1:8)如下所示:
式中,为端点i中d相节点电压的实部,d∈Bp,Bp={a,b,c}。fi d端点i中d相节点电压的虚部。为端点i中中性点电压的实部。fi n端点i中中性点电压的虚部。Pi dn为端点i中d相对于中性点的有功功率。为端点i中d相对于中性点的无功功率。
4.6)计算状态变量的修正量Δx(iter),即:
Δx(iter)=inv(H(iter))*r(iter)。 (14)
式中,iter为计算迭代次数。x(iter)为第iter次迭代时的状态变量。为第iter次迭代时量测变量的雅可比矩阵。inv为求逆运算。r(iter)=z-h(x(iter))为迭代值为x(iter)时的残差。
4.7)更新状态变量,得到状态变量迭代值x(iter+1)=x(iter)+Δx(iter),iter=iter+1,直至满足收敛标准max(Δx(iter))<ε为止。满足收敛标准的状态变量迭代值记为x'。
4.8)根据状态变量x',计算全网潮流和有功损耗Ploss,有功损耗Ploss计算公式如下:
5)分别计算首断面有功损耗Ploss1、末断面有功损耗Ploss2和平均有功损耗对应的电量损耗ΔWloss1、电量损耗ΔWloss2和电量损耗并根据电量损耗ΔWloss1、电量损耗ΔWloss2和电量损耗计算得到台区电网理论线损ΔWloss。
进一步,计算台区电网理论线损ΔWloss的主要步骤如下:
ΔWloss1=Ploss1*1/3time。 (16)
ΔWloss2=Ploss2*1/3time。 (17)
5.2)计算台区电网理论线损ΔWloss,即:
式中,ΔWloss为总时间段长度time内基于同时段功率和电量计算的理论线损。
实施例2:
基于同时段功率和电量的台区电网理论线损计算方法,主要包括以下步骤:
1)获取台区电网基本数据。
2)基于台区电网基本数据,建立台区电网端点导纳矩阵Y。
5)分别计算首断面有功损耗Ploss1、末断面有功损耗Ploss2和平均有功损耗对应的电量损耗ΔWloss1、电量损耗ΔWloss2和电量损耗并根据电量损耗ΔWloss1、电量损耗ΔWloss2和电量损耗计算得到台区电网理论线损ΔWloss。
实施例3:
式中,为端点导纳矩阵中的元素。自导纳为与端点i的节点d直接相联的所有支路导纳之和。互导纳为端点i的节点d和节点t之间支路导纳的相反数与跟端点i相连的端点的节点d和节点t之间所有支路导纳之和。互导纳为端点i和端点k中节点d和t之间支路导纳的相反数。φi为不包括端点i且与端点i直接相连端点的集合。为端点i节点d、节点X间并联支路导纳。为端点i节点d、节点t间并联支路导纳。为端点i至端点j支路节点d、节点t间的支路互导纳元素,且 为端点i至端点k支路节点d、t间并联支路导纳元素,且B1={a,b,c,n},表示三相与中性点的集合。
实施例4:
1)根据电量消耗数据和总时间段长度time,计算平均有功功率即:
式中,表示用户智能电表测得的总时间段长度time内d相有功电量,d∈Bp,Bp={a,b,c}表示三相节点的集合。time表示总时间段长度。
实施例5:
基于同时段功率和电量的台区电网理论线损计算方法,主要步骤见实施例2,其中,对台区电网进行三相潮流计算的主要步骤如下:
1)将配电变压器低压侧端点作为平衡端点S。其中,平衡端点S的中性点作为零电位参考点,平衡端点S的A相节点作为零相位参考点。设定平衡端点三相电压幅值为额定相电压幅值,三相相角分别相差120°。
2)设定所有端点的三相电压幅值初始值为台区电网的额定相电压。任意端点中,3个三相电压相角两两相差120°,中性点电压为0。
3)初始化最大迭代次数Tmax,并设置收敛精度ε。设定初始迭代次数iter=1。
4)计算三相潮流计算中任意量测变量的残差,即:
r=z-h(x)。 (1)
式中,z为m维的量测变量。m为量测变量个数。x为n维的状态变量。n为状态变量个数。m=n。h(x)为量测方程。量测方程h(x)为包含零注入端点的节点注入功率对应的等效注入电流量测方程。r为量测残差。
其中,包含零注入端点的节点注入功率对应的等效注入电流量测方程如下所示:
式中,为d相节点的等效注入电流,d∈Bp,Bp={a,b,c}表示三相节点的集合。为中性点的等效注入电流。为包括端点i且与端点i直接相连端点的集合。B1={a,b,c,n}为三相节点和中性点的集合。为端点导纳矩阵中端点i中d相节点和端点k中t相节点的导纳元素。为端点导纳矩阵中端点i中中性点和端点k中t相节点的导纳元素。为端点k中t相节点电压相量。为端点i处d相对于中性点的注入功率,且 为端点i处d相对中性点的发电机功率。为端点i处d相对中性点的负荷功率。为端点i的d相节点电压相量。为端点i的中性点电压相量。
5)建立用于三相潮流计算的雅克比矩阵H,主要分为以下两种情况:
式中,Hik为雅可比矩阵的子矩阵,且k≠i。为端点i中节点t和端点k中节点δ之间支路电导的相反数,δ∈B1,t∈B1,B1={a,b,c,n}。为端点i中节点t和端点k中节点δ之间支路电纳的相反数,δ∈B1,t∈B1,B1={a,b,c,n}。
式中,Hii为雅可比子矩阵矩阵。为端点i中节点t和δ之间支路电导的相反数与跟端点i相连的端点中节点t和δ之间所有支路电导之和,δ∈B1,t∈B1,B1={a,b,c,n}。为端点i中节点t和δ之间支路电纳的相反数与跟端点i相连的端点中节点t和δ之间所有支路电纳之和。
其中,矩阵ΔHii(1:8,1:8)如下所示:
其中,矩阵元素和矩阵元素分别如下所示:
式中,为端点i中d相节点电压的实部,d∈Bp,Bp={a,b,c}。fi d端点i中d相节点电压的虚部。为端点i中中性点电压的实部。fi n端点i中中性点电压的虚部。Pi dn为端点i中d相对于中性点的有功功率。为端点i中d相对于中性点的无功功率。
6)计算状态变量的修正量Δx(iter),即:
Δx(iter)=inv(H(iter))*r(iter)。 (8)
式中,iter为计算迭代次数。x(iter)为第iter次迭代时的状态变量。为第iter次迭代时量测变量的雅可比矩阵。inv为求逆运算。r(iter)=z-h(x(iter))为迭代值为x(iter)时的残差。
7)更新状态变量,得到状态变量迭代值x(iter+1)=x(iter)+Δx(iter),iter=iter+1,直至满足收敛标准max(Δx(iter))<ε为止。满足收敛标准的状态变量迭代值记为x'。
8)根据状态变量x',计算全网潮流和有功损耗Ploss,有功损耗Ploss计算公式如下:
实施例6:
基于同时段功率和电量的台区电网理论线损计算方法,主要步骤见实施例2,其中,计算台区电网理论线损ΔWloss的主要步骤如下:
ΔWloss1=Ploss1*1/3time。 (1)
ΔWloss2=Ploss2*1/3time。 (2)
2)计算台区电网理论线损ΔWloss,即:
式中,ΔWloss为总时间段长度time内基于同时段功率和电量计算的理论线损。
实施例7:
参见图2至图3,一种验证基于同时段功率和电量的台区电网理论线损计算方法的实验,主要包括以下步骤:
1)获取台区电网的基础数据,主要步骤如下:
获取台区电网的参数、网络结构和用户智能电表的量测数据;
所述电网的参数主要包括电网中元件的电阻、电抗、电纳和额定电压;
所述网络结构信息为网络拓扑连接关系;
2)计算台区电网端点导纳矩阵Y,主要步骤如下:
式中,为端点导纳矩阵中的元素;自导纳为与端点i的节点d直接相联的所有支路导纳之和;互导纳为端点i的节点d和节点t之间支路导纳的相反数与跟端点i相连的端点的节点d和节点t之间所有支路导纳之和;互导纳为端点i和端点k中节点d和t之间支路导纳的相反数;φi为不包括端点i且与端点i直接相连端点的集合;为端点i节点d、节点X间并联支路导纳;为端点i节点d、节点t间并联支路导纳;为端点i至端点j支路节点d、节点t间的支路互导纳元素,且 为端点i至端点k支路节点d、t间并联支路导纳元素,且B1={a,b,c,n},表示三相与中性点的集合。
式中,表示首断面的采集的d相有功功率;表示首断面的采集的d相无功功率;表示末断面的采集的d相有功功率;表示末断面的采集的d相无功功率;表示首断面d相功率因素;表示末断面d相功率因素;d∈Bp,Bp={a,b,c}表示三相节点的集合。
计算首末断面总时间段长度time内平均功率因素
4)利用总时间段长度time内的首断面有功和无功末断面有功和无功以及总时间段长度time内平均有功和平均无功d∈Bp,Bp={a,b,c},依据电流注入型牛顿法,分别对台区电网做三次潮流计算,台区电网电流注入型牛顿法潮流计算主要步骤如下:
4.1)将配电变压器低压侧端点作为平衡端点S;其中,平衡端点S的中性点作为零电位参考点,平衡端点S的A相节点作为零相位参考点;给定平衡端点三相电压幅值,取为额定相电压幅值,三相相角分别相差120°。
4.2)在台区电网进行三相潮流计算时,所有端点的三相电压幅值初始值均为所述台区电网的额定相电压;任意端点中,3个三相电压相角两两相差120°,中性点电压为0。
4.3)初始化最大迭代次数Tmax,并设置收敛精度ε;设定初始迭代次数iter=1。
4.4)计算三相潮流计算中各量测变量的残差,即。
r=z-h(x) (7)
式中,z为m维的量测变量;m为量测变量个数;x为n维的状态变量;n为状态变量个数;m=n;h(x)为量测方程;量测方程h(x)主要包括包含零注入端点的节点注入功率对应的等效注入电流量测方程;r为量测残差。
其中,包含零注入端点的节点注入功率对应的等效注入电流量测方程如下所示:
式中,为d相节点的等效注入电流,d∈Bp,Bp={a,b,c}表示三相节点的集合;为中性点的等效注入电流;为包括端点i且与端点i直接相连端点的集合;B1={a,b,c,n}为三相节点和中性点的集合;为端点导纳矩阵中端点i中d相节点和端点k中t相节点的导纳元素;为端点导纳矩阵中端点i中中性点和端点k中t相节点的导纳元素;为端点k中t相节点电压相量;为端点i处d相对于中性点的注入功率,且 为端点i处d相相对中性点的发电机功率;为端点i处d相相对中性点的负荷功率;为端点i的d相节点电压相量;为端点i的中性点电压相量。
4.5)计算三相潮流计算中雅克比矩阵H:
式中,Hik为雅可比矩阵的子矩阵,且k≠i;为端点i中节点t和端点k中节点δ之间支路电导的相反数,δ∈B1,t∈B1,B1={a,b,c,n};为端点i中节点t和端点k中节点δ之间支路电纳的相反数,δ∈B1,t∈B1,B1={a,b,c,n}。
式中,Hii为雅可比子矩阵矩阵;为端点i中节点t和δ之间支路电导的相反数与跟端点i相连的端点中节点t和δ之间所有支路电导之和,δ∈B1,t∈B1,B1={a,b,c,n};为端点i中节点t和δ之间支路电纳的相反数与跟端点i相连的端点中节点t和δ之间所有支路电纳之和。
其中,矩阵ΔHii(1:8,1:8)如下所示:
式中,为端点i中d相节点电压的实部,d∈Bp,Bp={a,b,c};fi d端点i中d相节点电压的虚部;为端点i中中性点电压的实部。fi n端点i中中性点电压的虚部;Pi dn为端点i中d相对于中性点的有功功率;为端点i中d相对于中性点的无功功率。
4.6)计算状态变量的修正量Δx(iter)
Δx(iter)=inv(H(iter))*r(iter) (14)
式中,iter为计算迭代次数;x(iter)为第iter次迭代时的状态变量;为第iter次迭代时量测变量的雅可比矩阵;inv为求逆运算;r(iter)=z-h(x(iter))为迭代值为x(iter)时的残差。
4.7)更新状态变量,得到状态变量新值x(iter+1)=x(iter)+Δx(iter),iter=iter+1,直至满足收敛标准max(Δx(iter))<ε为止。
4.8)依据迭代收敛的状态变量x,计算全网潮流和有功损耗Ploss,有功损耗Ploss计算公式如下:
5)利用三次潮流计算的有功损耗Ploss1、Ploss2和每次潮流计算结果分别占总时间段长度time内的1/3,计算1/3总时间段长度time内的电量损耗ΔWloss1、ΔWloss2和将三次电量损耗相加将即可求得总时间段长度time内的理论线损ΔWloss,主要步骤如下:
ΔWloss1=Ploss1*1/3time (16)
ΔWloss2=Ploss2*1/3time (17)
式中,Ploss1、Ploss2和分别为三次潮流计算的有功损耗;ΔWloss1、ΔWloss2和分别为三次潮流计算的理论线损;time为总时间段长度;ΔWloss为总时间段长度time内基于同时段功率和电量计算的理论线损。
实施案例9:
参见图4,基于同时段功率和电量的台区电网理论线损计算方法的对比实验,主要步骤如下:
1)在IEEE-13节点标准配电(12.47kV/4.16kV)系统基础上,构建一个修正系统。IEEE-13节点修正系统的构建过程如下:
1.1)将各支路的型号均设置为501,且相间距相同均设置为1.1292m;
1.2)线路7-11的a相和b相单位长度电阻设为0;
1.3)所有的负荷都是恒定功率Y型接线的负荷;
1.4)忽略了系统中的并联电容器、电压调节器和配电变压器。系统中的5、7和8号端点为零注入端点,1号端点为平衡端点,平衡端点中性点接大地作为零电位参考。
2)在IEEE-13节点配电修正系统的基础上,基于某台区配变低压侧总表一个月采集的电量和功率数据,将总表电量和功率数据按照IEEE-13节点配电修正系统的各节点与根节点的功率和电量比例关系,分配到各个负荷节点,并将标准系统的线电压修改为380V,模拟实际低压配电网,进行月度理论线损的模拟计算。分三种方案计算:
(1)常规的基于每日冻结电量的台区电网理论线损计算方法。用日电量的平均功率来做潮流计算,计算理论线损。将总表的日消耗电量数据除以24算得总表一天三相总的有功功率,进一步利用总表的功率因数,假设为0.9,计算出总表的总无功功率,并将该功率按照IEEE-13节点修正系统上各节点各相的功率相对于配变低压侧的总功率的比例进行分配,并做潮流计算,计算月度理论线损。
(2)基于智能电表小时级电量数据的台区电网理论线损计算方法。用小时级的电量来进行潮流计算,模拟方法同(1),计算月度理论线损。
(3)基于本发明的一种基于同时段功率和电量的台区电网理论线损计算方法。在(2)的基础上结合每个时间段内首末两个断面功率进行3次潮流计算,每次占20min时间,计算月度理论线损。
以上三种方案的月度理论线损的计算结果如表1所示。
表1三种模拟方案计算的月度理论线损
由表1可知,经三种模拟方法计算后,方案(1)计算的月度理论线损为0.143174kWh,方案(2)计算的月度理论线损为0.153038kWh,方案(3)计算的月度理论线损为0.155743kWh。由于方案(1)和方案(2)只用到了电量数据,且方案(1)用的是日冻结电量数据,故而理论线损计算的误差会较大,方案(2)用的是小时级电量数据,计算精度会比方案(1)高。在方案(2)的基础上,方案(3)采用本发明的一种基于同时段功率和电量的台区电网理论线损计算方法,同时利用了功率和电量数据来计算理论线损,计算精度最高。
Claims (6)
1.基于同时段功率和电量的台区电网理论线损计算方法,其特征在于,主要包括以下步骤:
1)获取所述台区电网基本数据。
2)基于台区电网基本数据,建立台区电网端点导纳矩阵Y;
式中,表示用户智能电表测得的总时间段长度time内d相有功电量,d∈Bp,Bp={a,b,c}表示三相节点的集合;time表示总时间段长度;
5.根据权利要求1所述的一种基于同时段功率和电量的台区电网理论线损计算方法,其特征在于,对台区电网进行三相潮流计算的主要步骤如下:
1)将配电变压器低压侧端点作为平衡端点S;其中,平衡端点S的中性点作为零电位参考点,平衡端点S的A相节点作为零相位参考点;设定平衡端点三相电压幅值为额定相电压幅值,三相相角分别相差120°;
2)设定所有端点的三相电压幅值初始值为台区电网的额定相电压;任意端点中,3个三相电压相角两两相差120°,中性点电压为0;
3)初始化最大迭代次数Tmax,并设置收敛精度ε;设定初始迭代次数iter=1。
4)计算三相潮流计算中任意量测变量的残差,即:
r=z-h(x); (7)
式中,z为m维的量测变量;m为量测变量个数;x为n维的状态变量;n为状态变量个数;m=n;h(x)为量测方程;量测方程h(x)为包含零注入端点的节点注入功率对应的等效注入电流量测方程;r为量测残差;
其中,包含零注入端点的节点注入功率对应的等效注入电流量测方程如下所示:
式中,为d相节点的等效注入电流,d∈Bp,Bp={a,b,c}表示三相节点的集合;为中性点的等效注入电流;为包括端点i且与端点i直接相连端点的集合;B1={a,b,c,n}为三相节点和中性点的集合;为端点导纳矩阵中端点i中d相节点和端点k中t相节点的导纳元素;为端点导纳矩阵中端点i中中性点和端点k中t相节点的导纳元素;为端点k中t相节点电压相量;为端点i处d相对于中性点的注入功率,且 为端点i处d相相对中性点的发电机功率;为端点i处d相相对中性点的负荷功率;为端点i的d相节点电压相量;为端点i的中性点电压相量。
5)建立用于三相潮流计算的雅克比矩阵H,主要分为以下两种情况:
式中,Hik为雅可比矩阵的子矩阵,且k≠i;为端点i中节点t和端点k中节点δ之间支路电导的相反数,δ∈B1,t∈B1,B1={a,b,c,n};为端点i中节点t和端点k中节点δ之间支路电纳的相反数,δ∈B1,t∈B1,B1={a,b,c,n};
式中,Hii为雅可比子矩阵矩阵;为端点i中节点t和δ之间支路电导的相反数与跟端点i相连的端点中节点t和δ之间所有支路电导之和,δ∈B1,t∈B1,B1={a,b,c,n};为端点i中节点t和δ之间支路电纳的相反数与跟端点i相连的端点中节点t和δ之间所有支路电纳之和;
其中,矩阵ΔHii(1:8,1:8)如下所示:
式中,为端点i中d相节点电压的实部,d∈Bp,Bp={a,b,c};fi d端点i中d相节点电压的虚部;为端点i中中性点电压的实部。fi n端点i中中性点电压的虚部;Pi dn为端点i中d相对于中性点的有功功率;为端点i中d相对于中性点的无功功率;
6)计算状态变量的修正量Δx(iter),即:
Δx(iter)=inv(H(iter))*r(iter); (14)
式中,iter为计算迭代次数;x(iter)为第iter次迭代时的状态变量;为第iter次迭代时量测变量的雅可比矩阵;inv为求逆运算;r(iter)=z-h(x(iter))为迭代值为x(iter)时的残差;
7)更新状态变量,得到状态变量迭代值x(iter+1)=x(iter)+Δx(iter),iter=iter+1,直至满足收敛标准max(Δx(iter))<ε为止;满足收敛标准的状态变量迭代值记为x';
8)根据状态变量x',计算全网潮流和有功损耗Ploss;有功损耗Ploss计算公式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910711520.6A CN110635471B (zh) | 2019-08-02 | 2019-08-02 | 基于同时段功率和电量的台区电网理论线损计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910711520.6A CN110635471B (zh) | 2019-08-02 | 2019-08-02 | 基于同时段功率和电量的台区电网理论线损计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110635471A true CN110635471A (zh) | 2019-12-31 |
CN110635471B CN110635471B (zh) | 2023-06-16 |
Family
ID=68969247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910711520.6A Active CN110635471B (zh) | 2019-08-02 | 2019-08-02 | 基于同时段功率和电量的台区电网理论线损计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110635471B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112418701A (zh) * | 2020-12-02 | 2021-02-26 | 广东电网有限责任公司广州供电局 | 配电网线路等效电阻计算方法、装置、终端及存储介质 |
CN113884820A (zh) * | 2021-11-19 | 2022-01-04 | 安徽南瑞中天电力电子有限公司 | 基于电能表的线路阻抗测量方法及线路故障类型分析方法 |
CN115201742A (zh) * | 2022-09-16 | 2022-10-18 | 北京志翔科技股份有限公司 | 低压台区计量点误差的确定方法 |
CN115276004A (zh) * | 2022-09-21 | 2022-11-01 | 石家庄科林电气股份有限公司 | 电网区块线损计算方法、装置及终端设备 |
WO2024037351A1 (zh) * | 2022-08-19 | 2024-02-22 | 西门子(中国)有限公司 | 配电系统的非技术损失检测方法、电子设备及存储介质 |
CN118068082A (zh) * | 2024-04-17 | 2024-05-24 | 威海双城电气有限公司 | 一种基于物联网的电线用电量监测系统及方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100872725B1 (ko) * | 2008-11-10 | 2008-12-05 | (주)제일전기기술단 | 배전계통의 배전손실 산정 방법 |
CN107046286A (zh) * | 2017-05-12 | 2017-08-15 | 国网上海市电力公司 | 基于前推回代法的低压配电台区理论线损计算方法 |
CN109188204A (zh) * | 2018-09-07 | 2019-01-11 | 济南大学 | 计及小电源的配电网理论线损计算方法 |
CN109214684A (zh) * | 2018-09-06 | 2019-01-15 | 国网福建省电力有限公司电力科学研究院 | 一种配网线路理论线损在线计算方法 |
CN109217303A (zh) * | 2018-10-15 | 2019-01-15 | 国网重庆市电力公司长寿供电分公司 | 基于配电变压器二次侧实抄电量的10kV配电台区理论线损计算方法 |
CN109449920A (zh) * | 2018-09-27 | 2019-03-08 | 中国电力科学研究院有限公司 | 一种低压配电网理论线损计算方法 |
CN109560540A (zh) * | 2018-11-05 | 2019-04-02 | 广西大学 | 基于线损负荷特性曲线的配电网理论线损计算方法 |
CN109861230A (zh) * | 2018-12-28 | 2019-06-07 | 重庆大学 | 一种含光伏逆变电源的三相四线制低压配电网改进潮流计算方法 |
-
2019
- 2019-08-02 CN CN201910711520.6A patent/CN110635471B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100872725B1 (ko) * | 2008-11-10 | 2008-12-05 | (주)제일전기기술단 | 배전계통의 배전손실 산정 방법 |
CN107046286A (zh) * | 2017-05-12 | 2017-08-15 | 国网上海市电力公司 | 基于前推回代法的低压配电台区理论线损计算方法 |
CN109214684A (zh) * | 2018-09-06 | 2019-01-15 | 国网福建省电力有限公司电力科学研究院 | 一种配网线路理论线损在线计算方法 |
CN109188204A (zh) * | 2018-09-07 | 2019-01-11 | 济南大学 | 计及小电源的配电网理论线损计算方法 |
CN109449920A (zh) * | 2018-09-27 | 2019-03-08 | 中国电力科学研究院有限公司 | 一种低压配电网理论线损计算方法 |
CN109217303A (zh) * | 2018-10-15 | 2019-01-15 | 国网重庆市电力公司长寿供电分公司 | 基于配电变压器二次侧实抄电量的10kV配电台区理论线损计算方法 |
CN109560540A (zh) * | 2018-11-05 | 2019-04-02 | 广西大学 | 基于线损负荷特性曲线的配电网理论线损计算方法 |
CN109861230A (zh) * | 2018-12-28 | 2019-06-07 | 重庆大学 | 一种含光伏逆变电源的三相四线制低压配电网改进潮流计算方法 |
Non-Patent Citations (2)
Title |
---|
王晓明 等: ""采用改进牛顿法计算配电网理论线损"" * |
陈潇一 等: ""基于智能电表和PMU混合量测的低压配电网 三相状态估计"" * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112418701A (zh) * | 2020-12-02 | 2021-02-26 | 广东电网有限责任公司广州供电局 | 配电网线路等效电阻计算方法、装置、终端及存储介质 |
CN112418701B (zh) * | 2020-12-02 | 2022-10-25 | 广东电网有限责任公司广州供电局 | 配电网线路等效电阻计算方法、装置、终端及存储介质 |
CN113884820A (zh) * | 2021-11-19 | 2022-01-04 | 安徽南瑞中天电力电子有限公司 | 基于电能表的线路阻抗测量方法及线路故障类型分析方法 |
WO2024037351A1 (zh) * | 2022-08-19 | 2024-02-22 | 西门子(中国)有限公司 | 配电系统的非技术损失检测方法、电子设备及存储介质 |
CN115201742A (zh) * | 2022-09-16 | 2022-10-18 | 北京志翔科技股份有限公司 | 低压台区计量点误差的确定方法 |
CN115276004A (zh) * | 2022-09-21 | 2022-11-01 | 石家庄科林电气股份有限公司 | 电网区块线损计算方法、装置及终端设备 |
CN115276004B (zh) * | 2022-09-21 | 2023-02-03 | 石家庄科林电气股份有限公司 | 电网区块线损计算方法、装置及终端设备 |
CN118068082A (zh) * | 2024-04-17 | 2024-05-24 | 威海双城电气有限公司 | 一种基于物联网的电线用电量监测系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110635471B (zh) | 2023-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110635471A (zh) | 基于同时段功率和电量的台区电网理论线损计算方法 | |
Džafić et al. | Real time estimation of loads in radial and unsymmetrical three-phase distribution networks | |
CN103532137B (zh) | 一种三相四线低压配电网的状态估计方法 | |
CN108306284B (zh) | 一种基于局部智能量测的在线负荷建模方法 | |
CN109217303B (zh) | 基于配电变压器二次侧实抄电量的10kV配电台区理论线损计算方法 | |
CN105826918A (zh) | 电网理论线损的计算方法和系统 | |
CN102509159A (zh) | 基于用电信息采集系统的配电网理论线损评估方法 | |
CN103236691B (zh) | 基于复仿射数学理论的三相不平衡潮流计算方法 | |
CN106026086B (zh) | 一种电网中运行状态的动态估计方法 | |
CN102623993A (zh) | 一种分布式电力系统状态估计方法 | |
CN103208797A (zh) | 一种基于智能优化技术的含新能源配电网状态估计方法 | |
CN106655152B (zh) | 一种基于ami量测特性的配电网状态估计方法 | |
CN100563078C (zh) | 一种准确计算电力线路损耗的方法 | |
CN105891612B (zh) | 一种三相不平衡治理效果量化评估方法 | |
CN111861136A (zh) | 一种台区线损率评价标杆值计算方法 | |
CN109861230B (zh) | 一种含光伏逆变电源的三相四线制低压配电网改进潮流计算方法 | |
CN113076618B (zh) | 一种配电网单线图校正方法、系统、终端设备和存储介质 | |
CN103972889B (zh) | 一种配电线路阻抗在线辨识方法 | |
Gupta et al. | Power system network equivalents: Key issues and challenges | |
CN105977963B (zh) | 一种基于两阶段法的主配网一体化分布式状态估计方法 | |
Sharma et al. | PyPSA: Open Source Python Tool for Load Flow Study | |
CN103838962A (zh) | 一种计及pmu量测的分步线性状态估计方法 | |
CN116819230A (zh) | 高密度分布式光伏电源的低压台区线损计算方法及系统 | |
CN103986231A (zh) | 一种网络化变电站智能计量系统 | |
CN116599035A (zh) | 低压配网日理论线损概率分析方法、装置、电子设备及介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |