feat: filter DC lines out of AC transmission network processing within HIFLD construction

[danielolsen]

Pull Request doc

Purpose

Partial fulfilment of #233: this addresses the HVDC transmission lines within the HIFLD dataset, but does not yet address the B2Bs, which will require a different method of splitting substations into sub-substations for each interconnection.

What the code is doing

In const.py, ratings are added for DC lines with known capacities.

In transmission.py:

• A small helper function split_lines_to_ac_and_dc is added to filter out DC lines via some combination of the TYPE attribute of the HIFLD dataset and manually specified indices.
• build_transmission uses this function, then adds impedances and AC ratings to AC branches only (using functionality introduced in feat: add function to create transformers within substations #218 and feat: add functions to estimate branch parameters (impedance and rating) #219), and uses the pre-defined ratings from const.py to add capacities to the DC lines.

Testing

Tested manually.

Usage Example/Visuals

>>> from prereise.gather.griddata.hifld.data_process.transmission import build_transmission
>>> branch, substations, dc_lines = build_transmission()
dropping 6892 substations of 70857 total due to LINES parameter equal to 0
filter substations based on lines
---------------------------------
assigning substations to lines' endpoints by mapping their rounded (3 digits) coordinates
dropping 3314 lines having same endpoints coordinates after rounding
finding closest substation to unmapped lines' endpoint(s)
100%|███████████████████████████████████| 2624/2624 [00:00<00:00, 19820.94it/s]
dropping 0 lines due to island size filtering
dropping 0 substations due to island size filtering
Reading cached minimum spanning tree
2948 line voltages can't be found via neighbor consensus
No more missing voltages remain after neighbor minimum
>>> branch.columns
Index(['TYPE', 'STATUS', 'NAICS_CODE', 'NAICS_DESC', 'SOURCE', 'SOURCEDATE',
       'VAL_METHOD', 'VAL_DATE', 'OWNER', 'VOLTAGE', 'VOLT_CLASS', 'INFERRED',
       'SUB_1', 'SUB_2', 'COORDINATES', 'SUB_1_ID', 'SUB_2_ID', 'type',
       'length', 'from_bus_id', 'to_bus_id', 'x', 'rateA'],
      dtype='object')
>>> dc_lines.columns
Index(['TYPE', 'STATUS', 'NAICS_CODE', 'NAICS_DESC', 'SOURCE', 'SOURCEDATE',
       'VAL_METHOD', 'VAL_DATE', 'OWNER', 'VOLTAGE', 'VOLT_CLASS', 'INFERRED',
       'SUB_1', 'SUB_2', 'COORDINATES', 'SUB_1_ID', 'SUB_2_ID', 'Pmax',
       'Pmin'],
      dtype='object')
>>> dc_lines["Pmax"]
200823    3100
308464    2400
310053     400
311958       5
157627     330
157629     660
158515    2000
108354     500
113313     660
131914    2000
Name: Pmax, dtype: int64

Branch entries get impedance (x) and rating, DC lines get Pmax and Pmin.

Time estimate

15 minutes.

[BainanXia]

Good to go. Thanks!