CN103595039B - The distributed cyclization risk analysis method of network system - Google Patents

The distributed cyclization risk analysis method of network system Download PDF

Info

Publication number
CN103595039B
CN103595039B CN201310488081.XA CN201310488081A CN103595039B CN 103595039 B CN103595039 B CN 103595039B CN 201310488081 A CN201310488081 A CN 201310488081A CN 103595039 B CN103595039 B CN 103595039B
Authority
CN
China
Prior art keywords
bus
cyclization
node
network
risk analysis
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.)
Active
Application number
CN201310488081.XA
Other languages
Chinese (zh)
Other versions
CN103595039A (en
Inventor
潘凯岩
王清玲
李吉勇
刘仲尧
刘尚伟
赵勇
李文朝
胡铁斌
叶云琴
宋学清
董文杰
刘全昌
谭大帅
何群玮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Maoming Power Supply Bureau of Guangdong Power Grid Co Ltd
Dongfang Electronics Co Ltd
Original Assignee
Maoming Power Supply Bureau of Guangdong Power Grid Co Ltd
Dongfang Electronics Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Maoming Power Supply Bureau of Guangdong Power Grid Co Ltd, Dongfang Electronics Co Ltd filed Critical Maoming Power Supply Bureau of Guangdong Power Grid Co Ltd
Priority to CN201310488081.XA priority Critical patent/CN103595039B/en
Publication of CN103595039A publication Critical patent/CN103595039A/en
Application granted granted Critical
Publication of CN103595039B publication Critical patent/CN103595039B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The object of the present invention is to provide the distributed cyclization risk analysis method of a kind of network system, comprise and carry out dynamic partition process centered by 220kV bus, Tidal Simulation analysis more than 220kV electric pressure and following cyclization risk analysis analog computation, 220kV voltage.This method is run by ring network power supply mode for equipment more than 220kV electric pressure, the equipment of below 220kV electric pressure is pressed in the environment that radiation network runs, and the actual motion feature invention for area power grid provides a kind of distributed cyclization risk analysis processing method based on interrupter technique and Sub-region and hierarchical.The present invention aims to provide a general cyclization risk analysis method, and the cyclization simulated operation computing speed of 110kV and following electric pressure is improved greatly, and the stability simultaneously restrained obtains very large improvement.

Description

The distributed cyclization risk analysis method of network system
Technical field
The present invention relates to power system automatic field, relate in particular to a kind of distributed cyclization risk analysis processing method based on interrupter technique and Sub-region and hierarchical.
Background technology
Along with the development of electric power system and carrying out in a deep going way of Power grid commercial operation, the scheduling and controlling of electric power system becomes and becomes increasingly complex.The data acquisition that dispatch automated system is the most basic and monitoring function SCADA can not meet the requirement of traffic control, for improving the safety of electrical network, stable, high-quality and economic operation level, the electrical network advanced applied software PAS that can analyze electric network composition and real time data is more and more subject to the favor of management and running personnel, in the simulated operation of senior application, cyclization risk analysis is very important functional module, dispatcher in any real-time operation as cyclization, whether all will carry out simulated operation in advance before unlinking etc. with verification operation can cause the out-of-limit of busbar voltage, whether can cause the overload of circuit or transformer, loop power flow changing situation etc. after operation.But based on the Load flow calculation of the whole network often sometimes because the model in other sheet net or region is incomplete or the inaccurate cyclization analytical calculation to one's respective area impacts and maybe can not use, sometimes convergence also can receive impact, based on this, propose the distributed cyclization risk analysis method based on interrupter technique and Sub-region and hierarchical herein, make the cyclization simulated operation of local area not by the impact of other sheet net of system, stability is greatly improved.
Summary of the invention
The object of the present invention is to provide the distributed cyclization risk analysis method of a kind of network system, the method provides more general, stable computational methods for the cyclization simulated operation of area power grid 110kV and following electric pressure, the method is not subject to the impact of other sheet net by the subnet that interrupter technique and Dynamic partitioning techniques make plan operate, thus substantially increases the stability of calculating and improve computational speed.
Object of the present invention realizes by following technical measures:
The distributed cyclization risk analysis method of a kind of network system, comprises the whole network dynamic partition, subnet distributed cyclization risk analysis calculating, major network Load flow calculation three part computational process, wherein:
Described the whole network dynamic partition carries out dynamic partition process centered by 220kV bus, carries out topology process according to the ground state topological relation of electrical network and the real time position of switch, during process, 220kV main transformer uprised switch simulation and disconnect, the subregion that final formation is different.
The detailed process of described the whole network dynamic partition is: the dynamic partition determining the subregion of 500kV electric pressure major network and the subnet centered by 220kV bus according to the real-time topology relation of electrical network and the on off state of correspondence, need 220kV transformer to be uprised switch to disconnect when determining the subregion of 500kV electric pressure, need 220kV line switching to disconnect when determining the dynamic partition centered by 220kV bus, then 220kV bus and 500kV bus are set to balance node and carry out dynamic partition, final result is the corresponding subregion of each 220kV bus, the corresponding dynamic partition of 500kV bus.
Wherein dynamic partition with the 220kV bus in 220kV transformer station for root node, 220kV line switching is disconnected and blocks, search for downwards by depth-first search rule, search for all 110kV transformer stations of this 220kV transformer station band, using this 220kV transformer station and the 110kV transformer station that searches as a region.In general, how many 220kV transformer stations are had just to have how many districts.Carry out major network subregion after having determined subnet subregion, 220kV transformer is uprised switch simulation and disconnect, search element with any one 500kV bus for initial point carries out depth-first downwards, finally determine the major network dynamic partition comprising 500kV bus and 220kV bus.
Then, according to the information of the switch of the simulated operation provided, determine the dynamic partition belonging to it, in corresponding subregion, carry out the operation of the Tidal Simulation after cyclization.
The detailed process that the distributed cyclization risk analysis of described subnet calculates is: with the 220kV bus in the subregion of current sub network as balance node, the cyclization switch of intending carrying out simulated operation in Home Network is closed, start Load flow calculation, the power flow changing of the trend distribution situation obtaining equipment in this subregion especially on looped network.Owing to calculating the subregion only related among a small circle, be not subject to the influence on tidal flow of major network and other subregion, thus computational speed improves greatly, and computing time obviously shortens, and convergence is greatly improved.
First the switch of intending carrying out Loop Closing Operation is closed, and determine the subregion at switch place, computing node is obtained by topology process by switch junctions model, final electrical island information is determined again by node branch model, 220kV bus in this electrical island is set to balance node, because algorithm quicksort is faster than Newton method in arithmetic speed, PQ decomposition method is adopted for the tidal current analysis after the cyclization of subnet, finally obtain voltage phase angle and the Branch Power Flow of each node as calculated, the Power Flow Information of cyclization switch and the active reactive injection rate of balance node, thus can judge whether Loop Closing Operation can cause apparatus overload.
Described PQ method carries out corrected Calculation by following iteration form:
V -1ΔP=B'VΔδ
V -1ΔQ=B''ΔV
Wherein Δ P, Δ Q are respectively the meritorious of each iteration and idle correction, and B' is n-1 rank admittance battle array, and B'' is n-m rank admittance battle array (m is PV node number), and Δ δ, Δ V are respectively the phase angle of each iteration and the correction of voltage magnitude.V -1for the inverse matrix of node voltage diagonal matrix.
Described major network Load flow calculation adopts the power after cyclization risk analysis simulated operation to the power on the 220kV border bus of the subregion carrying out cyclization simulated operation, to the power that the power of the 220kV border bus do not carried out in the subregion of simulated operation adopts bus current by the mode of equivalence, carry out the Load flow calculation of major network on this basis.
Finally, with the 5000kV bus in the subregion of current major network as balance node, adopt equivalence method to process to the injecting power on the bus of 220kV border and obtain the meritorious of border bus and idle injection rate, in major network, the node type of other node is constant, start Load flow calculation, thus obtain the trend distribution situation of main interior equipment.Only relate to the subregion of major network owing to calculating, process border bus by equivalence, thus achieve the decoupling zero process of major network and border netting, thus computational speed improves greatly, and computing time obviously shortens, and convergence is greatly improved.
The detailed process of described major network Load flow calculation is: the node injecting power first obtaining 220kV border bus, choose bus in 500kV to connect and occur maximum being set to balance node, computing node is obtained by topology process by switch junctions model, final island information is determined again by node branch model, because Newton method stability is strong, major network adopts Newton method to process, and finally obtains the voltage phase angle of each node and Branch Power Flow, Branch Power Flow thus can judge whether the equipment of major network can transship as calculated.
Described Newton method carries out corrected Calculation by following iteration form
ΔP ΔQ = H N J L Δδ V - 1 ΔV
Wherein Δ P, Δ Q are respectively the meritorious of each iteration and idle correction, and H, N, J, L are the Jacobian matrix of object, Δ δ, and Δ V is respectively the phase angle of each iteration and the correction of voltage magnitude.V -1for the inverse matrix of node voltage diagonal matrix.
Accompanying drawing explanation
Fig. 1 is dynamic partition schematic diagram of the present invention;
Fig. 2 is that subnet cyclization sunykatuib analysis of the present invention calculates flow chart;
Fig. 3 is major network Load flow calculation flow chart of the present invention.
Embodiment
Below in conjunction with accompanying drawing and concrete enforcement, the present invention will be further described.
Suppose that the switch that will carry out simulating Loop Closing Operation is in the dynamic partition centered by 220kVC station.Computational process is as follows:
As Fig. 1, carry out the whole network dynamic partition: dynamic partition with the 220kV bus in 220kV transformer station for root node, 220kV line switching is disconnected and blocks, search for downwards by depth-first search rule, search for all 110kV transformer stations of this 220kV transformer station band, using this 220kV transformer station and the 110kV transformer station that searches as a region.As shown in Figure 1,220kVC station is a district with transformer station 1,2,3,4.In general, how many 220kV transformer stations are had just to have how many districts.Carry out major network subregion after having determined subnet subregion, 220kV transformer is uprised switch simulation and disconnect, search element with any one 500kV bus for initial point carries out depth-first downwards, finally determine the major network dynamic partition comprising 500kV bus and 220kV bus.
As Fig. 2, subnet cyclization sunykatuib analysis calculates: first closed by the switch of intending carrying out Loop Closing Operation, and determine the subregion at switch place, computing node is obtained by topology process by switch junctions model, final island information is determined again by node branch model, 220kV bus in the island proper is set to balance node, because PQ decomposition method is faster than Newton method in arithmetic speed, PQ decomposition method is adopted for the tidal current analysis after the cyclization of subnet, finally obtain voltage phase angle and the Branch Power Flow of each node as calculated, Branch Power Flow, the Power Flow Information of cyclization switch and the active reactive injection rate of balance node, thus can judge whether Loop Closing Operation can cause apparatus overload.
As Fig. 3, major network Load flow calculation: final major network subregion is the main partition comprising 220kV circuit and 220kV bus centered by 500kVA station, first the node injecting power of 220kV border bus is obtained, choose bus in 500kV to connect and occur maximum being set to balance node, computing node is obtained by topology process by switch junctions model, final island information is determined again by node branch model, because Newton method stability is strong, major network adopts Newton method to process, finally obtain voltage phase angle and the Branch Power Flow of each node as calculated, Branch Power Flow, thus can judge whether the equipment of major network can transship.
Embodiments of the present invention are not limited thereto; under stating basic fundamental thought prerequisite on the invention; according to the ordinary technical knowledge of this area and customary means to content of the present invention make the amendment of other various ways, replacement or change, all drop within rights protection scope of the present invention.

Claims (6)

1. the distributed cyclization risk analysis method of network system, is characterized in that: comprise the whole network dynamic partition, subnet distributed cyclization risk analysis calculating, major network Load flow calculation three part computational process, wherein:
Described the whole network dynamic partition carries out dynamic partition process centered by 220kV bus, carries out topology process according to the ground state topological relation of electrical network and the real time position of switch, during process, 220kV main transformer uprised switch simulation and disconnect, the subregion that final formation is different;
It is using the 220kV bus in the subregion of current sub network as balance node that the distributed cyclization risk analysis of described subnet calculates, the cyclization switch of intending carrying out simulated operation in Home Network is closed, start Load flow calculation, the power flow changing of the trend distribution situation obtaining equipment in this subregion especially on looped network;
Described major network Load flow calculation adopts the power after cyclization risk analysis simulated operation to the power on the 220kV border bus of the subregion carrying out cyclization simulated operation, to the power that the power of the 220kV border bus do not carried out in the subregion of simulated operation adopts bus current by the mode of equivalence, carry out the Load flow calculation of major network on this basis;
The detailed process of described the whole network dynamic partition is: the dynamic partition determining the subregion of 500kV electric pressure major network and the subnet centered by 220kV bus according to the real-time topology relation of electrical network and the on off state of correspondence, need 220kV transformer to be uprised switch to disconnect when determining the subregion of 500kV electric pressure, need 220kV line switching to disconnect when determining the dynamic partition centered by 220kV bus, then 220kV bus and 500kV bus are set to balance node and carry out dynamic partition, final result is the corresponding subregion of each 220kV bus, the corresponding dynamic partition of 500kV bus.
2. the distributed cyclization risk analysis method of network system according to claim 1, it is characterized in that: the detailed process of described the whole network dynamic partition is: wherein dynamic partition with the 220kV bus in 220kV transformer station for root node, 220kV line switching is disconnected and blocks, search for downwards by depth-first search rule, search for all 110kV transformer stations of this 220kV transformer station band, using this 220kV transformer station and the 110kV transformer station that searches as a region; Carry out major network subregion after having determined subnet subregion, 220kV transformer is uprised switch simulation and disconnect, search element with any one 500kV bus for initial point carries out depth-first downwards, finally determine the major network dynamic partition comprising 500kV bus and 220kV bus.
3. the distributed cyclization risk analysis method of network system according to claim 1, it is characterized in that: the detailed process that the distributed cyclization risk analysis of described subnet calculates is: first the switch of intending carrying out Loop Closing Operation is closed, and determine the subregion at switch place, computing node is obtained by topology process by switch junctions model, final electrical island information is determined again by node branch model, 220kV bus in this electrical island is set to balance node, PQ decomposition method is adopted for the tidal current analysis after the cyclization of subnet, finally obtain voltage phase angle and the Branch Power Flow of each node as calculated, the Power Flow Information of cyclization switch and the active reactive injection rate of balance node, thus can judge whether Loop Closing Operation can cause apparatus overload.
4. the distributed cyclization risk analysis method of network system according to claim 3, is characterized in that: described PQ decomposition method carries out corrected Calculation by following iteration form:
V -1ΔP=B'VΔδ
V -1ΔQ=B″ΔV
Wherein Δ P, Δ Q are respectively the meritorious of each iteration and idle correction, and B' is n-1 rank admittance battle array, B, and " be n-m rank admittance battle array, m is PV node number, and Δ δ, Δ V are respectively the phase angle of each iteration and the correction of voltage magnitude, V -1for the inverse matrix of node voltage diagonal matrix.
5. the distributed cyclization risk analysis method of network system according to claim 1, is characterized in that:
The detailed process of described major network Load flow calculation is: the node injecting power first obtaining 220kV border bus, choose bus in 500kV to connect and occur maximum being set to balance node, computing node is obtained by topology process by switch junctions model, final island information is determined again by node branch model, because Newton method stability is strong, major network adopts Newton method to process, and finally obtains the voltage phase angle of each node and Branch Power Flow, Branch Power Flow thus can judge whether the equipment of major network can transship as calculated.
6. the distributed cyclization risk analysis method of network system according to claim 5, is characterized in that: described Newton method carries out corrected Calculation by following iteration form
ΔP ΔQ = H N J L Δδ V - 1 ΔV
Wherein Δ P, Δ Q are respectively the meritorious of each iteration and idle correction, and H, N, J, L are the Jacobian matrix of object, Δ δ, and Δ V is respectively the phase angle of each iteration and the correction of voltage magnitude; V -1for the inverse matrix of node voltage diagonal matrix.
CN201310488081.XA 2013-10-17 2013-10-17 The distributed cyclization risk analysis method of network system Active CN103595039B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310488081.XA CN103595039B (en) 2013-10-17 2013-10-17 The distributed cyclization risk analysis method of network system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310488081.XA CN103595039B (en) 2013-10-17 2013-10-17 The distributed cyclization risk analysis method of network system

Publications (2)

Publication Number Publication Date
CN103595039A CN103595039A (en) 2014-02-19
CN103595039B true CN103595039B (en) 2015-08-05

Family

ID=50085074

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310488081.XA Active CN103595039B (en) 2013-10-17 2013-10-17 The distributed cyclization risk analysis method of network system

Country Status (1)

Country Link
CN (1) CN103595039B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103872681A (en) * 2014-03-25 2014-06-18 国家电网公司 Online real-time loop closing method based on integration of major network and distribution network
CN104008462A (en) * 2014-06-16 2014-08-27 深圳供电局有限公司 Reliability management information system and method covering high-voltage, medium-voltage and low-voltage customers
CN105207199B (en) * 2015-09-08 2017-05-10 江苏省电力公司 Partition power grid main supply element overload quantity filling method
CN106300337A (en) * 2016-08-09 2017-01-04 国网浙江德清县供电公司 Power distribution network partition method based on dynamic Decomposition method
CN107480109B (en) * 2017-10-23 2021-06-01 积成电子股份有限公司 Incomplete parallel state estimation realization method based on network structure dynamic partitioning

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101431255A (en) * 2008-11-17 2009-05-13 国电南瑞科技股份有限公司 Closed loop operation risk analysis method for power system
CN102904346B (en) * 2012-11-01 2015-05-20 宁波电业局 On-off control method and device of loop closing switch of 10kV power distribution network

Also Published As

Publication number Publication date
CN103595039A (en) 2014-02-19

Similar Documents

Publication Publication Date Title
CN103595039B (en) The distributed cyclization risk analysis method of network system
CN102193504B (en) Safety and stability control system modeling method in power system dynamic simulation
CN103825267B (en) A kind of computational methods of MMC-MTDC dc-side short-circuit electric current
CN101750562B (en) Non-PMU measure point dynamic process estimation method based on flow equation sensitiveness analysis
CN103973203B (en) A kind of online equivalent modeling method in large-sized photovoltaic power station suitable in security and stability analysis
CN103872681A (en) Online real-time loop closing method based on integration of major network and distribution network
CN102403720B (en) Hyper-real-time setting method for superposed time sequences based on transient voltage safety margin
CN100466412C (en) United accident-prevention rehearse method for provincial and city power system
CN107453357A (en) A kind of State Estimation for Distribution Network based on hierarchical solving
CN103475010B (en) Transient stabilization control method for interconnection electrical network after impact disturbance
CN103001214B (en) A kind of power distribution network Three Phase Power Flow based on neutral point excursion
CN103065016A (en) Method and system for determining alternating-current system fault range resulting in direct-current commutation failure
CN107516896B (en) Load recovery scheme generation method based on decision tree pruning algorithm
CN102832618A (en) Practical method for calculation of closed-loop impact current of power distribution network based on combination of power transmission network and power distribution network
CN107392498A (en) A kind of power network topology analysis system
CN103020726A (en) All-PMU (phase measurement unit)-orientated robust state estimation method
CN103400303A (en) Static equivalence method for external network based on internal network routine Ward equivalent circuit
CN103199528A (en) Status estimating and coordinating method of wide-area power system
CN104113061A (en) Three-phase load flow calculation method of power distribution network with distributed power supply
CN105893714B (en) Based on the detection of automatic voltage control system closed loop and appraisal procedure under Bulk power system simulation system
CN103475003A (en) Automatic regional power grid voltage reactive power control method based on regional strategy optimization
CN104156531A (en) Fault current limiter distribution point optimization and capacity selection method taking operation loss into consideration
CN103632031B (en) A kind of rural area based on load curve decomposition load type load modeling method
CN103810915B (en) Regulation and control integrative simulation method and system
Fang et al. Realization of electromechanical transient and electromagnetic transient real time hybrid simulation in power system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: DONGFANG ELECTRONICS CORP.

Effective date: 20141021

C41 Transfer of patent application or patent right or utility model
C53 Correction of patent of invention or patent application
CB03 Change of inventor or designer information

Inventor after: Pan Kaiyan

Inventor after: Song Xueqing

Inventor after: Dong Wenjie

Inventor after: Liu Quanchang

Inventor after: Tan Dashuai

Inventor after: He Qunwei

Inventor after: Wang Qingling

Inventor after: Li Jiyong

Inventor after: Liu Zhongyao

Inventor after: Liu Shangwei

Inventor after: Zhao Yong

Inventor after: Li Wenchao

Inventor after: Hu Tiebin

Inventor after: Ye Yunqin

Inventor before: Wang Qingling

Inventor before: Song Xueqing

Inventor before: Zhao Yong

Inventor before: Dong Wenjie

Inventor before: Liu Quanchang

Inventor before: Li Wenchao

Inventor before: Hu Tiebin

Inventor before: Ye Yunqin

Inventor before: He Qunwei

Inventor before: Pan Kaiyan

Inventor before: Li Jiyong

Inventor before: Liu Zhongyao

Inventor before: Liu Shangwei

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: WANG QINGLING LI WENCHAO HU TIEBIN YE YUNQIN HE QUNWEI PAN KAIYAN LI JIYONG LIU ZHONGYAO LIU SHANGWEI SONG XUEQING ZHAO YONG DONG WENJIE LIU QUANCHANG TO: PAN KAIYAN WANG QINGLING LI JIYONG LIU ZHONGYAO LIU SHANGWEI ZHAO YONG LI WENCHAO HU TIEBIN YE YUNQIN SONG XUEQING DONG WENJIE LIU QUANCHANG TAN DASHUAI HE QUNWEI

TA01 Transfer of patent application right

Effective date of registration: 20141021

Address after: 525000 Guangdong Province, Maoming City, Guan Shan Road, No. four, No. 2

Applicant after: MAOMING POWER SUPPLY BUREAU OF GUANGDONG POWER GRID Corp.

Applicant after: DONGFANG ELECTRONICS Co.,Ltd.

Address before: 525000 Guangdong Province, Maoming City, Guan Shan Road, No. four, No. 2

Applicant before: MAOMING POWER SUPPLY BUREAU OF GUANGDONG POWER GRID Corp.

C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 525000 No.2, Guanshan 4th Road, Maoming City, Guangdong Province

Patentee after: MAOMING POWER SUPPLY BUREAU, GUANGDONG POWER GRID Co.,Ltd.

Patentee after: DONGFANG ELECTRONICS Co.,Ltd.

Address before: 525000 No.2, Guanshan 4th Road, Maoming City, Guangdong Province

Patentee before: MAOMING POWER SUPPLY BUREAU OF GUANGDONG POWER GRID Corp.

Patentee before: DONGFANG ELECTRONICS Co.,Ltd.

TR01 Transfer of patent right

Effective date of registration: 20210224

Address after: 525000 No.2, Guanshan 4th Road, Maoming City, Guangdong Province

Patentee after: MAOMING POWER SUPPLY BUREAU, GUANGDONG POWER GRID Co.,Ltd.

Address before: 525000 No.2, Guanshan 4th Road, Maoming City, Guangdong Province

Patentee before: MAOMING POWER SUPPLY BUREAU, GUANGDONG POWER GRID Co.,Ltd.

Patentee before: DONGFANG ELECTRONICS Co.,Ltd.

TR01 Transfer of patent right