CN105468874A - Field-circuit coupling modeling method for alternating current power grid - Google Patents
Field-circuit coupling modeling method for alternating current power grid Download PDFInfo
- Publication number
- CN105468874A CN105468874A CN201510989387.2A CN201510989387A CN105468874A CN 105468874 A CN105468874 A CN 105468874A CN 201510989387 A CN201510989387 A CN 201510989387A CN 105468874 A CN105468874 A CN 105468874A
- Authority
- CN
- China
- Prior art keywords
- transformer station
- node
- transformer
- neutral point
- formula
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
- G06F30/398—Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention discloses a field-circuit coupling modeling method for an alternating current power grid. Data modeling is performed on the alternating current power grid, so as to analyze parameters in the alternating current power grid and achieve a purpose of accurately judging the risk of direct current magnetic bias.
Description
Technical field
The present invention relates to power domain, particularly, relate to a kind of field circuit method modeling method of AC network.
Background technology
DC magnetic biasing is a kind of abnormal operating state of transformer, and because the former limit equiva lent impedance of transformer only presents resistance characteristic to DC component, and resistance is very little.Therefore, very little DC component will form very large DC magnetizing magnetic potential in the windings, and this direct current magnetic potential acts on transformer primary side together with interchange magnetic potential, causes the work magnetization curve of transformer core to offset, occur asymmetric about initial point, i.e. transformer bias phenomenon.Cause the reason of transformer DC magnetic bias different, mainly comprise direct current monopolar operation and earth induction.
The exciting current of large-scale power transformer is smaller, and a small amount of DC current flowing through transformer just may cause DC magnetic biasing, causes iron core cutter, causes current waveform to distort, and produces higher hamonic wave, the safe operation of harm transformer and electric system.The D.C. high voltage transmission monopolar ground return method of operation easily causes AC Substation transformer around to occur DC magnetic bias phenomena, and converter power transformer is multiple DC magnetic biasing harm also.The earth current of current China extra-high voltage direct-current transmission is larger than common DC transmission engineering: now generally ± and 500kV direct current transportation earth current is 3000A, the wide extra-high voltage of cloud is 3125A, Burner zone-Shanghai and Xi Luodu-West Zhejiang Province extra-high voltage are 4000A, and silk screen-southern Jiangsu is 4500A.A large amount of earth currents will cause more serious DC magnetic biasing harm, jeopardize AC network safe operation.
Summary of the invention
The object of the invention is to, for the problems referred to above, propose a kind of field circuit method modeling method of AC network, to realize the advantage accurately judging DC magnetic biasing risk.
For achieving the above object, the technical solution used in the present invention is:
A field circuit method modeling method for AC network, is characterized in that, comprise the following steps:
If AC network is total m transformer station always, b bus nodes, n independent neutral point, then had by the nodal method of analysis:
YV=J(1)
In formula: V represents grid nodes column voltage vector, V=[V
s; V
n; V
b], V
s, V
n, V
brepresent transformer station's node voltage, independent neutral point voltage and busbar voltage column vector respectively;
Y represents grid nodes conductance matrix, Y=H
tg+Q; H is transformer station's node and all internodal incidence matrix, H
tfor the transposed matrix of H, H
m × (m+n+b)=[E
m0
m × n0
m × b], E
mfor m rank unit matrix; G is substation grounding conductance battle array, G=R
-1, R=diag (R
g1, R
g2..., R
gm), R
gibe i-th transforming plant DC stake resistance; Q is AC network ground upper network node conductance matrix;
J represents grid nodes Injection Current column vector:
J=[J
S;J
N;J
B]=[GP;0;0]=H
TGP(2)
In formula: J
s, J
n, J
brepresent transformer station's node, independent neutral point and bus nodes Injection Current column vector respectively;
P represents the induced potential column vector of transformer station, is had by ground connection theory:
P=MI
D+NI
A(3)
In formula: I
drepresent direct current pole earth current, I
arepresent the DC current injecting grounding net of transformer substation, M represents mutual resistance matrix between direct current pole and transformer station, and N represents the mutual resistance matrix not comprising self-acting between transformer station; P represents transformer station's induced potential, refers to the inlet potential between neutral point and zero point;
The DC current now injecting grounding net of transformer substation is:
I
A=G(V
A-P)(4)
In formula: V
arepresent transformer station's node voltage, it is defined as:
V
A=HV(5)
Simultaneous formula 1 to formula 5 draws:
(R-ZN)I
A=ZMI
D(6)
In formula (6), Z=HY
-1h
tg – E, HY
-1h
ttransformer station's node-resistance matrix.
Preferably, if independent neutral point and transformer station's node short circuit, then leave out the independent neutral point at this station, only retain transformer station's node.
Technical scheme of the present invention has following beneficial effect:
Technical scheme of the present invention, by carrying out data modeling to AC network, thus analyzes the parameter in AC network, reaches the object accurately judging DC magnetic biasing risk.
Below by drawings and Examples, technical scheme of the present invention is described in further detail.
Accompanying drawing explanation
The 750kV transformer station theory diagram that Fig. 1 runs for the two self couplings change medium voltage side bus sections described in the embodiment of the present invention;
The Equivalent DC modular concept block diagram that Fig. 2 is the transformer station described in the embodiment of the present invention;
Fig. 3 a is autotransformer winding type schematic diagram;
Fig. 3 b is other Transformer Winding type schematic diagram except autotransformer.
By reference to the accompanying drawings, in the embodiment of the present invention, Reference numeral is as follows:
1-1 transformer 220kV bus; 2-2 transformer 220kV bus; 3-750kV bus; 4-1 transformer; 5-2 transformer; 6-neutral point series devices; 7-grounding net of transformer substation; 8-circuit; 9-750kV bus nodes; 10-750kV winding; 11-220kV bus nodes; 12-220kV winding; 13-neutral point; 14-transformer station node; 15-stake resistance; 16-induced potential.
Embodiment
Below in conjunction with accompanying drawing, the preferred embodiments of the present invention are described, should be appreciated that preferred embodiment described herein is only for instruction and explanation of the present invention, is not intended to limit the present invention.
For direct current transportation earth current, modern AC electrical network is a huge DC network, and the field circuit method model that the technical program proposes needs to set up corresponding realistic model to circuit, transformer station and other grid parts.
Circuit model: circuit model is exactly three resistance branch connecting transformer station's three-phase bus in fact, and its key parameter has two: direct current resistance and the mode of connection.The estimation equation of circuit direct current resistance is:
In formula, R
lcircuit one phase resistance (Ω), the sectional area (m2) of ρ to be line material resistivity coefficient (Ω m), S be a phase circuit.Direct current resistance as the aluminum stranded conductor of the 2 × 400mm2 of long 1km is:
Can find out, the direct current resistance of 1km circuit is very little, and the equivalent d.c. resistance of three-phase line is only 0.02125 Ω, therefore low-resistance circuit is the important factor in order of DC current distribution.
In power system component, the mode of connection of circuit often changes, and circuit can divide multiple electric pressure.The circuit of 110kV and above thereof is only considered when the calculating of DC current distribution.In addition, three-phase dc resistance and the mode of connection of circuit are consistent, can consider by parallel way.
Transformer station model: in DC network model, transformer station model is the most complicated, can be subdivided into aerial part and under ground portion.Transformer station's aerial part refers to more than the ground shown in Fig. 1 not containing the part of circuit, and under ground portion refers to the below ground part shown in Fig. 1.Fig. 2 is shown in by the Equivalent DC model of transformer station shown in Fig. 1.
(1) transformer station's aerial part model: transformer station's aerial part model belongs to pure circuit problem, model mainly comprises transformer station's node, transformer model, neutral point series devices model.Transformer station's nodal analysis method is see Fig. 3 a and Fig. 3 b.
Transformer model comprises Transformer Winding type, winding D.C. resistance, transformer bus.
During DC current distribution calculates, differentiating transformer is only needed whether to be autotransformer.Transformer Winding type schematic diagram is see Fig. 3 a and Fig. 3 b.
In Fig. 3 a, autotransformer series connection winding switching high voltage bus and middle pressure bus, middle pressure bus connects neutral point or neutral point series devices (as low reactance, neutral point suppress DC current equipment) through public winding again, in Fig. 3 b, other form transformers (as three-winding transformer) are then that high voltage bus and middle pressure bus connect neutral point or neutral point series devices through respective winding.
Transformer bus is the end points of circuit, as the voltage-dropping type autotransformer used a large amount of in current 500kV transformer station, mostly is neutral-point solid ground, and medium voltage side bus section often will be adopted in reality to run restriction single-phase short circuit current.As shown in Figure 2, two autotransformers show as a 500kV node and connect two 220kV bus nodes through series connection winding, and two 220kV bus nodes connect respective neutral point or neutral point series devices through respective public winding.
The core content of DC current distribution is the DC current distribution on winding, and transformer model has significant impact to it.
Neutral point series devices model is relatively simple, its equivalent resistance branch road connection transformer neutral point and transformer station's node.If neutral point series devices is electric capacity, then disconnect the connection of neutral point and transformer station's node.
(2) transformer station's under ground portion model: the under ground portion model of transformer station is field parameters, model carries out equivalence according to Dai Weinan principle: under whole neutral point disconnection, stake resistance is the equivalent resistance between transformer station's node and zero point, and induced potential is the entrance electromotive force between transformer station's node and zero point.
(3) the asymmetric disposal route of triphase parameter: if AC network three-phase dc parameter is symmetrical, the single-phase circuit that can be reduced to as Fig. 2 is analyzed; If asymmetrical three-phase, can when earth's surface be constant with drag, with the above model in the earth's surface of Fig. 2 to wherein one carrying out modeling mutually, then transformer station's node of short circuit three-phase.
The technical program is symmetrical by DC network parameter, and the model that three-phase is merged into a phase carries out.Meanwhile, the technical program derivation result has versatility, also comparatively flexible during practical application, if triphase parameter asymmetric (as inconsistent in Three-Phase Transformer winding D.C. resistance), as long as the correlation parameter of amendment model becomes triphase parameter.
If AC network is total m transformer station always, b bus nodes, n independent neutral point, then had by the nodal method of analysis:
YV=J(2)
In formula (2):
V---grid nodes column voltage vector, V=[V
s; V
n; V
b], V
s, V
n, V
bfor transformer station's node voltage (V), independent neutral point voltage (V), busbar voltage column vector (V).If neutral point and transformer station's node short circuit, then leave out the neutral point at this station in this paper model, only retain transformer station's node.
Y---grid nodes conductance matrix, Y=H
tg+Q.H is transformer station's node and all internodal incidence matrix, H
tfor the transposition of H, H
m × (m+n+b)=[E
m0
m × n0
m × b], E
mfor m rank unit matrix; G is substation grounding conductance battle array, G=R
-1, R=diag (R
g1, R
g2..., R
gm), R
gibe i-th transforming plant DC stake resistance (Ω).Q is AC network ground upper network node conductance matrix.
J---grid nodes Injection Current column vector:
J=[J
S;J
N;J
B]=[GP;0;0]=H
TGP(3)
In formula (3):
J
s, J
n, J
b---transformer station's node, independent neutral point, bus nodes Injection Current column vector (A);
The induced potential column vector (V) of P---transformer station, is had by ground connection theory:
P=MI
D+NI
A(4)
In formula (4):
I
d---direct current pole earth current (A).
I
a---inject the DC current (A) of grounding net of transformer substation.
M---mutual resistance matrix (Ω) between direct current pole and transformer station.
N---the mutual resistance matrix (Ω) (not comprising self-acting) between transformer station.
P---transformer station's induced potential, refers to the inlet potential (V) between neutral point and zero point.
The definition of substation grounding resistance and induced potential can be called " the Dai Weinan equivalence " carried out at neutral point.The DC current now injecting grounding net of transformer substation is:
I
A=G(V
A-P)(5)
In formula (5):
V
a---transformer station's node voltage (V), it is defined as:
V
A=HV(6)
Simultaneous formula (2)-(6) have:
(R-ZN)I
A=ZMI
D(7)
In formula (7), Z=HY
-1h
tg – E.HY
-1h
tbe equivalent to take out the part of node-resistance matrix about transformer station's node.As can be seen from formula (2.24), removing magnitude of current I
aand I
d, R, N and M are pure field amount, and Z matrix is the amount of a field circuit method.Make transformer station node-resistance matrix F=HY
-1h
t, therefore Z (i, i)=F (i, i)/R
gi-1, Z (i, j)=F (i, j)/R
gi, i ≠ j.In the model of direct current distribution, if transformer station's neutral point is in the state of cut-offfing, then there is Z (i, i)=0, I
a(i)=0; If substation grounding resistance F (i, i)/R
gi(this situation is limiting case to ≈ 1, and namely soil is good conductor, R
gimuch smaller than circuit and transformer resistance; Or circuit and transformer resistance much larger than stake resistance time), then have Z (i, i) ≈ 0, I
a(i) ≈ 0.In fact soil is not good conductor, circuit and transformer resistance less, these two kinds of factors cause AC network direct current distribution.From another perspective, if the line that this also illustrates a transformer station is more, the risk of this station DC magnetic biasing is also larger.
Last it is noted that the foregoing is only the preferred embodiments of the present invention, be not limited to the present invention, although with reference to previous embodiment to invention has been detailed description, for a person skilled in the art, it still can be modified to the technical scheme described in foregoing embodiments, or carries out equivalent replacement to wherein portion of techniques feature.Within the spirit and principles in the present invention all, any amendment done, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (2)
1. a field circuit method modeling method for AC network, is characterized in that, comprise the following steps:
If AC network is total m transformer station always, b bus nodes, n independent neutral point, then had by the nodal method of analysis:
YV=J(1)
In formula: V represents grid nodes column voltage vector, V=[V
s; V
n; V
b], V
s, V
n, V
brepresent transformer station's node voltage, independent neutral point voltage and busbar voltage column vector respectively;
Y represents grid nodes conductance matrix, Y=H
tg+Q; H is transformer station's node and all internodal incidence matrix, H
tfor the transposed matrix of H, H
m × (m+n+b)=[E
m0
m × n0
m × b], E
mfor m rank unit matrix; G is substation grounding conductance battle array, G=R
-1, R=diag (R
g1, R
g2..., R
gm), R
gibe i-th transforming plant DC stake resistance; Q is AC network ground upper network node conductance matrix;
J represents grid nodes Injection Current column vector:
J=[J
S;J
N;J
B]=[GP;0;0]=H
TGP(2)
In formula: J
s, J
n, J
brepresent transformer station's node, independent neutral point and bus nodes Injection Current column vector respectively;
P represents the induced potential column vector of transformer station, is had by ground connection theory:
P=MI
D+NI
A(3)
In formula: I
drepresent direct current pole earth current, I
arepresent the DC current injecting grounding net of transformer substation, M represents mutual resistance matrix between direct current pole and transformer station, and N represents the mutual resistance matrix not comprising self-acting between transformer station; P represents transformer station's induced potential, refers to the inlet potential between neutral point and zero point;
The DC current now injecting grounding net of transformer substation is:
I
A=G(V
A-P)(4)
In formula: V
arepresent transformer station's node voltage, it is defined as:
V
A=HV(5)
Simultaneous formula 1 to formula 5 draws:
(R-ZN)I
A=ZMI
D(6)
In formula (6), Z=HY
-1h
tg – E, HY
-1h
ttransformer station's node-resistance matrix.
2. the field circuit method modeling method of AC network according to claim 1, is characterized in that, if independent neutral point and transformer station's node short circuit, then leaves out the independent neutral point at this station, only retains transformer station's node.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510989387.2A CN105468874A (en) | 2015-12-25 | 2015-12-25 | Field-circuit coupling modeling method for alternating current power grid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510989387.2A CN105468874A (en) | 2015-12-25 | 2015-12-25 | Field-circuit coupling modeling method for alternating current power grid |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105468874A true CN105468874A (en) | 2016-04-06 |
Family
ID=55606569
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510989387.2A Pending CN105468874A (en) | 2015-12-25 | 2015-12-25 | Field-circuit coupling modeling method for alternating current power grid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105468874A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112308395A (en) * | 2020-10-27 | 2021-02-02 | 国网江西省电力有限公司电力科学研究院 | Screening method and device for important power transmission channel |
CN117113733A (en) * | 2023-10-24 | 2023-11-24 | 国家电网有限公司西北分部 | Method and device for acquiring bias current of direct current near zone of power grid |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592062A (en) * | 2012-03-14 | 2012-07-18 | 东北电力大学 | Method for calculating DC (direct current) magnetic bias dynamic leakage inductance of transformer |
CN103439565A (en) * | 2013-09-04 | 2013-12-11 | 博爱县电业公司 | Calculation method for DC magnetic bias current flowing into transformer of high-voltage DC power transmission system |
CN104316810A (en) * | 2014-11-21 | 2015-01-28 | 武汉大学 | Direct-current magnetic bias multistage progressive comprehensive treatment assessment method |
CN104777384A (en) * | 2015-04-17 | 2015-07-15 | 华北电力大学(保定) | Determination method for magnetic hysteresis characteristics and depletion characteristics of transformer in DC magnetic bias state |
-
2015
- 2015-12-25 CN CN201510989387.2A patent/CN105468874A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592062A (en) * | 2012-03-14 | 2012-07-18 | 东北电力大学 | Method for calculating DC (direct current) magnetic bias dynamic leakage inductance of transformer |
CN103439565A (en) * | 2013-09-04 | 2013-12-11 | 博爱县电业公司 | Calculation method for DC magnetic bias current flowing into transformer of high-voltage DC power transmission system |
CN104316810A (en) * | 2014-11-21 | 2015-01-28 | 武汉大学 | Direct-current magnetic bias multistage progressive comprehensive treatment assessment method |
CN104777384A (en) * | 2015-04-17 | 2015-07-15 | 华北电力大学(保定) | Determination method for magnetic hysteresis characteristics and depletion characteristics of transformer in DC magnetic bias state |
Non-Patent Citations (2)
Title |
---|
于永军等: "天中直流工程入地电流对新疆哈密地区交流电网的影响", 《电网技术》 * |
王建等: "±800kV天—中直流对哈密电网变压器直流偏磁的影响", 《四川电力技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112308395A (en) * | 2020-10-27 | 2021-02-02 | 国网江西省电力有限公司电力科学研究院 | Screening method and device for important power transmission channel |
CN112308395B (en) * | 2020-10-27 | 2022-06-14 | 国网江西省电力有限公司电力科学研究院 | Screening method and device for important power transmission channel |
CN117113733A (en) * | 2023-10-24 | 2023-11-24 | 国家电网有限公司西北分部 | Method and device for acquiring bias current of direct current near zone of power grid |
CN117113733B (en) * | 2023-10-24 | 2024-02-02 | 国家电网有限公司西北分部 | Method and device for acquiring bias current of direct current near zone of power grid |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Boteler et al. | Comparison of methods for modelling geomagnetically induced currents | |
Haginomori et al. | Power system transient analysis: theory and practice using simulation programs (ATP-EMTP) | |
Abu-Hashim et al. | Test systems for harmonics modeling and simulation | |
Andersson | Power system analysis | |
CA2567519A1 (en) | Method and equipment for the protection of power systems against geomagnetically induced currents | |
Benato et al. | Zero sequence behaviour of a double-circuit overhead line | |
Berge et al. | Laboratory validation of the relationship between Geomagnetically Induced Current (GIC) and transformer absorbed reactive power | |
Chiesa et al. | Calculation of inrush currents–Benchmarking of transformer models | |
Zadeh et al. | Investigation of neutral reactor performance in reducing secondary arc current | |
CN103412199A (en) | Computing method of unbalancedness of multi-circuit power transmission lines on same tower | |
Yalçın et al. | A study of symmetrical and unsymmetrical short circuit fault analyses in power systems | |
Piccinelli et al. | North Europe power transmission system vulnerability during extreme space weather | |
CN105468874A (en) | Field-circuit coupling modeling method for alternating current power grid | |
Bu et al. | Unified magnetic equivalent circuit-based model of two-core Sen transformer for low-frequency electromagnetic transients and its numerical calculation | |
CN104483639B (en) | Residual magnetism estimation method for non-fault tripping of YNd11 type three-phase combined transformer | |
Roy | Generalised polyphase fault-analysis program: calculation of cross-country fault | |
Abdullah et al. | The study of triplen harmonics currents produced by salient pole synchronous generator | |
CN105631119A (en) | Simulation modeling method for grounding grid | |
Li et al. | Research on effects of transformer DC bias on negative sequence protection | |
Bellan et al. | Circuit representation of voltage unbalance emission due to line asymmetry | |
Lambert et al. | Topological transformer leakage modeling with losses | |
Novitskiy et al. | Some aspects of steady state simulations of AC/DC hybrid transmission lines | |
Kim et al. | Bus-voltage sag suppressing and fault current limiting characteristics of the SFCL due to its application location in a power distribution system | |
Lopez-Fernandez et al. | RNM2D 0 Fast Stray Losses Hazard Evaluation on Transformer Tank Wall & Cover due to Zero Sequence | |
Bernardić et al. | NEUTRAL POINT CONNECTIONS IN MV POWER NETWORKS WITH GROUNDING ZIGZAG TRANSFORMERS–ANALYSIS AND SIMULATIONS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160406 |
|
RJ01 | Rejection of invention patent application after publication |