Household distributed photovoltaic power statistics and transaction settlement method

Abstract

The present application relates to photovoltaic power metering and settlement technical field, disclose a kind of household distributed photovoltaic power statistical and transaction settlement method, comprising the following steps: obtaining normalized metering data by obtaining original cumulative power reading, extract authoritative cumulative anchor point and calculate authoritative power increment, generate settlement microelement, verify independent time-sharing source to obtain reference power increment and independent source available identifier, using reverse type identification reverse residual pair, form multi-channel residual signature and determine abnormal type, obtain automatic correction result or time mapping version;Determine microelement power and determine identifiable label by constructing constraint feasible solution set, temporarily settlement record is generated when identifiable, determine settlement influence domain according to abnormal type and propagation rule and local re-settlement, generate final settlement record and settlement archiving single after grid single check.This application identifies intermediate cumulative anchor point deviation as a single anomaly and distinguishes abnormal type, settles when power is identifiable in price segment, and local re-settlement according to settlement influence domain.

Description

Technical Field

[0001] This invention relates to the field of photovoltaic power metering and settlement technology, specifically a method for statistical analysis and transaction settlement of household distributed photovoltaic power. Background Technology

[0002] Residential distributed photovoltaic systems are equipped with bidirectional settlement electricity meters and generation metering devices to record cumulative grid-connected electricity, cumulative off-grid electricity, and generated electricity, and to conduct transaction settlement according to time-of-use pricing, contract version, and settlement cycle. Existing methods include differential settlement of cumulative meter readings, time-of-use calculation at multiple metering points, supplementary recording of key nodes, verification, estimation, and editing of general metering data, and settlement evidence storage based on a trusted ledger.

[0003] When a deviation occurs at the intermediate cumulative anchor point, the deviation will simultaneously enter two adjacent anchoring intervals. Existing methods often process them separately by interval, splitting the same anchor point anomaly into two interval anomalies; time alignment deviations and actual changes in electricity consumption also lack distinguishing rules, making it difficult to classify and process them according to anchor point value anomalies, time alignment anomalies, missing time-sharing data, and physical relationship anomalies.

[0004] When time-of-use electricity is written into the settlement, existing methods usually directly provide a supplementary entry or allocation result. After the cumulative total is closed, the electricity in micro-elements will still not be unique; existing methods do not distinguish whether the distribution of micro-elements is unique and whether the electricity in price segments can be uniquely determined, and will include electricity with uncertain price attribution in the settlement.

[0005] After data correction or version change, existing methods often determine the scope of resettlement based on the reference relationship of account items or the overall recalculation of the accounting period, without taking into account the anomaly type and the power propagation relationship of power generation, grid connection, and grid disconnection channels. Irrelevant settlement micro-elements and account items are easily recalculated together. Summary of the Invention

[0006] To address the shortcomings of existing technologies, this invention provides a method for statistical analysis and transaction settlement of residential distributed photovoltaic power generation. This method solves the problems of locating and classifying anomalies at intermediate cumulative anchor points, generating settlements under the condition that the power generation in price segments is identifiable, and local re-settlement constrained by anomaly propagation relationships.

[0007] To achieve the above objectives, the present invention provides the following technical solution:

[0008] This invention provides a method for statistical analysis and transaction settlement of residential distributed photovoltaic power generation, comprising the following steps:

[0009] Obtain the raw cumulative electricity readings, normalize and segment them into events to obtain normalized metering data;

[0010] Authoritative cumulative anchor points are extracted from standardized metering data to form anchoring intervals and calculate authoritative electricity increments. Settlement micro-elements are generated. Using the authoritative cumulative anchor points as cumulative anchor points to be verified, independent time-sharing sources not derived from the cumulative anchor points to be verified are verified, and reference electricity increments and available identifiers of independent sources are obtained.

[0011] By combining authoritative power increment, reference power increment and independent source available identifier to calculate interval residuals, a reverse method is used to identify reverse residual pairs with opposite directions and closed sums on the anchor interval, forming a multi-channel residual signature and determining the anomaly type. When the anomaly type is an anchor value anomaly, the anchor correction amount is calculated to obtain the automatic correction result. When the anomaly type is a time alignment anomaly, a time mapping version is generated.

[0012] Based on the authoritative electricity increment, settlement micro-element and the automatic correction result or time mapping version corresponding to the anomaly type, construct a set of constrained feasible solutions, determine the micro-element electricity and identify the identifiable marker, and generate a temporary settlement record when the identifiable marker is identifiable.

[0013] The settlement impact domain is determined based on the anomaly type and preset propagation rules, and a local re-settlement is performed. Temporary settlement records are then checked against the power grid single checklist to generate the final settlement record and settlement archive sheet.

[0014] Furthermore, the steps for obtaining raw cumulative electricity readings, standardizing and segmenting them by event to obtain standardized metering data include: registering the metering topology and obtaining the metering resolution, collecting raw data and obtaining raw cumulative electricity readings and grid settlement statements, standardizing and converting the raw cumulative electricity readings, dividing the cumulative segments by event and calculating the effective increment to obtain standardized metering data.

[0015] Furthermore, the steps for obtaining the reference power increment and the independent source availability identifier include: extracting authoritative cumulative anchor points from the standardized metering data, forming an anchoring interval and obtaining the authoritative power increment, merging time boundaries to generate settlement micro-elements; when the independent time-of-use source is not derived from the cumulative anchor point to be verified, calculating the reference power increment and the reference micro-element power, and the independent source availability identifier is available; otherwise, the independent source availability identifier is unavailable.

[0016] Furthermore, the steps for generating multi-channel residual signatures and determining anomaly types include:

[0017] When an independent source is marked as available, the interval residual is obtained by subtracting the reference energy increment from the authoritative energy increment, and the allowable deviation of the interval is determined; when the absolute value of the interval residual is not greater than the allowable deviation of the interval, the anchored interval is marked as a closable interval; the authoritative energy increments of the generation channel, grid connection channel and grid disconnection channel are mapped to the generation increment, grid connection increment and grid disconnection increment, respectively.

[0018] The authoritative cumulative anchor point shared by adjacent anchoring intervals is used as the intermediate cumulative anchor point. A reverse method is adopted. When the available identifiers of independent sources in both intervals are available, the residual directions of the intervals are opposite, the absolute values ​​of each interval are greater than the allowable deviation of their respective intervals, and the absolute value of the sum of the two is not greater than the sum of the allowable deviations of the two intervals, a reverse residual pair is formed, candidate abnormal anchor points are determined, and the reverse residual direction status is recorded.

[0019] When the absolute values ​​of the residuals of both intervals are greater than their respective interval allowable deviations but in the same direction, or when the absolute value of the sum of the residuals of both intervals is greater than the sum of the two interval allowable deviations, a time mapping candidate is generated and tested; when the absolute values ​​of the interval residuals of multiple consecutive anchoring intervals are greater than their respective interval allowable deviations, they are merged into a continuous abnormal anchor point group; the physical feasibility state is determined based on the grid connection increment, power generation increment, and grid disconnection increment.

[0020] When the increment of grid connection is not greater than the increment of power generation, the increment of power generation is not greater than the product of the rated AC output power and the anchoring interval duration, and the increment of grid disconnection is not less than zero, the physically feasible state is physically feasible; otherwise, it is physically infeasible. The independent source availability identifiers, reverse residual direction states, and physically feasible states of the power generation channel, grid connection channel, and grid disconnection channel are summarized to form a multi-channel residual signature. When the physically feasible state is physically infeasible and the reverse formula is not valid, the anomaly type is physical relationship anomaly.

[0021] When there is a unique time mapping candidate such that the absolute value of the interval residual of each channel after mapping is not greater than its respective interval allowable deviation, the anomaly type is time alignment anomaly; when a reverse residual pair is formed and it does not belong to time alignment anomaly, the anomaly type is anchor value anomaly; when the independent source availability identifier is unavailable, the anomaly type is time-division data missing; when the absolute value of the interval residual is greater than the interval allowable deviation and does not belong to the above anomaly types, the anomaly type is unlocated anomaly.

[0022] Furthermore, the steps for obtaining the automatic correction results include: determining authoritative cumulative anchor points that are both closable intervals and are not candidate anomaly anchor points as reliable anchor points; when there is only a single candidate anomaly anchor point, calculating the anchor point correction amount of the intermediate cumulative anchor point using a correction formula to form an automatic correction result and write it into the correction view; when a continuous group of anomaly anchor points is formed, using the reliable anchor points on both sides as fixed boundaries, calculating the anchor point correction amount of each intermediate cumulative anchor point using a group correction formula to form an automatic correction result and write it into the correction view.

[0023] Furthermore, the steps for generating a time mapping version include: generating time mapping candidates and recalculating the interval residuals, taking the optimal mapping to obtain the time mapping version; when the optimal mapping is not unique and its impact on the electricity consumption of the price segment is greater than the metering resolution, a non-unique mapping state is formed and marked as pending confirmation.

[0024] Furthermore, the steps for constructing a set of constrained feasible solutions and determining the micro-element power include: constructing a set of constrained feasible solutions, which needs to satisfy preset cumulative total constraints, power generation constraints, directional constraints, or energy storage metering constraints when directional constraints are disabled; when an independent source is marked as available, the reference micro-element power is obtained from the independent time-sharing source; when the set of constrained feasible solutions is not empty, the micro-element power is solved by projection, so that the weighted squared difference between the micro-element power of each channel and the reference micro-element power is minimized.

[0025] Furthermore, the step of generating a temporary settlement record when the identifiable marker is identified includes: obtaining a price segment from the settlement micro-element, calculating the lower and upper bounds of the segment's electricity volume using the identifiable method; when the difference between the upper and lower bounds of the segment's electricity volume is not greater than the metering resolution, the identifiable marker is identified, the identifiable electricity volume is obtained, and a temporary settlement record is generated; when the difference between the upper and lower bounds of the segment's electricity volume is greater than the metering resolution, the settlement micro-element to which the price segment belongs is set to pending settlement, and a pending settlement record is generated.

[0026] Furthermore, the steps of determining the settlement impact domain and performing partial resettlement based on the anomaly type and preset propagation rules include: constructing the settlement impact domain according to the propagation rules, where changes in power generation channels affect power generation and consumption within the boundary, changes in grid connection channels affect grid connection and consumption within the boundary, changes in off-grid channels affect off-grid consumption, changes in time mapping affect the union of settlement micro-elements involving the old and new time mapping versions, and changes in price versions affect amount accounts while the basic metered electricity remains unchanged; when correcting a single anchor point, only the two anchor intervals before and after are processed, and continuous abnormal anchor point groups are processed with the two credible anchor points as boundaries; when the amount of issued accounts decreases, a refund is performed, when it increases, a supplementary collection is performed, and when the entire account is replaced, a reversal is performed, generating a replacement record and obtaining the resettled account.

[0027] Furthermore, the steps for generating the final settlement record and settlement archive include: obtaining the power grid settlement statement from the original data, and comparing the power grid settlement statement with the local verification quantity; when a corrected total amount exists, the local verification quantity is taken as the corrected total amount, otherwise the authoritative electricity increment is taken; when the verification results are consistent, the final settlement record is generated from the temporary settlement record, and a settlement archive is generated in conjunction with the resettled items; when the verification results are inconsistent, an unsettled discrepancy statement is formed, and the resettled items are resettled locally according to the settlement impact domain and then verified again; when the re-verification results are consistent, the final settlement record and settlement archive are generated.

[0028] This invention provides a method for statistical analysis and transaction settlement of residential distributed photovoltaic power generation. It has the following beneficial effects:

[0029] 1. This invention obtains the reference power increment through an independent time-division source, uses a reverse identification method to identify reverse residual pairs with opposite directions and closed summation, and forms a multi-channel residual signature to determine the anomaly type. It treats the deviation between the preceding and following intervals caused by the same intermediate cumulative anchor point as a single anchor point value anomaly, and processes it separately from time alignment anomalies, time-division data missing, and physical relationship anomalies.

[0030] 2. This invention constructs a set of feasible solutions under the constraints of authoritative electricity increment and settlement micro-element, uses projection to determine the micro-element electricity, and uses distinguishability to calculate the lower and upper bounds of the electricity in the price segment. When the difference between the two is not greater than the metering resolution, it is determined to be distinguishable and a temporary settlement record is generated.

[0031] 3. This invention determines the settlement impact domain based on the anomaly type and propagation rules and performs local resettlement. When correcting a single anchor point, only the two anchor intervals before and after are processed. Issued accounts are refunded, supplemented, or reversed to generate alternative records. After verification by the power grid, the final settlement record and settlement archive are generated. Attached Figure Description

[0032] Figure 1 This is a system architecture diagram of the present invention;

[0033] Figure 2 This is a flowchart of the method of the present invention;

[0034] Figure 3 This is a flowchart illustrating the standardization process of metering topology registration, data acquisition, and cumulative meter codes in this invention.

[0035] Figure 4 This is a flowchart illustrating the authoritative cumulative anchor point extraction, settlement micro-element generation, and independent time-sharing source verification process of this invention.

[0036] Figure 5 This is a flowchart of the interval residual calculation, anomaly diagnosis and correction process of the present invention;

[0037] Figure 6 This is a flowchart illustrating the closed projection of electricity, price identification and determination, and transaction settlement process of the present invention.

[0038] Figure 7 This is a flowchart illustrating the settlement dependency establishment, partial resettlement, and power grid verification processes of the present invention.

[0039] Figure 8 This is a box plot of the absolute error of electricity consumption in price ranges in an application embodiment of the present invention;

[0040] Figure 9 This is a box plot of the absolute error of the settlement amount in an application embodiment of the present invention. Detailed Implementation

[0041] The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0042] Please see the appendix Figure 1 This invention provides a residential distributed photovoltaic power generation statistics and transaction settlement system, comprising:

[0043] The metering standardization module is used to register metering topologies, determine capability assessment results and processing modes, collect raw data and save it in layers, standardize the raw cumulative electricity readings and segment events to obtain standardized metering data.

[0044] The anchor point micro-element module is used to extract authoritative cumulative anchor points, calculate authoritative power increments, merge time boundaries to generate settlement micro-elements, and verify independent time-sharing sources to obtain reference power increments and reference micro-element power.

[0045] The anomaly diagnosis module is used to calculate the allowable deviation and residual of the interval. It uses a reverse method to identify reverse residual pairs, forms a multi-channel residual signature and determines the anomaly type. It uses a correction method or a group correction method to correct the authoritative cumulative anchor point, or generates a time-mapped version.

[0046] The projection settlement module is used to construct a set of constrained feasible solutions, determine the infinitesimal energy using projection, determine the identifiable markers and identifiable energy using distinguishability, calculate the energy absorbed within the boundary, and generate temporary settlement records or records to be settled.

[0047] The recalculation and verification module is used to establish settlement dependencies, determine the settlement impact domain, perform local recalculation, and generate final settlement records and settlement archives based on the power grid single verification results.

[0048] See appendix Figure 2 This invention also provides a method for statistical analysis and transaction settlement of residential distributed photovoltaic power generation, comprising the following steps:

[0049] S1, the metering specification module registers the metering topology, determines the capability determination results and processing mode, collects raw data and saves it in layers, standardizes the raw cumulative electricity readings and segments events to obtain standardized metering data.

[0050] S2, the anchor point micro-element module extracts the authoritative cumulative anchor point from the normalized power sequence, calculates the authoritative power increment, merges the time boundaries to generate the settlement micro-element, and verifies the independent time-sharing source to obtain the reference power increment and the reference micro-element power.

[0051] S3, the anomaly diagnosis module calculates the allowable deviation and residual of the interval, uses reverse identification to identify reverse residual pairs, forms a multi-channel residual signature and determines the anomaly type; for anchor value anomalies, it uses correction or group correction to calculate the anchor correction amount; for time alignment anomalies, it generates a time mapping version.

[0052] S4, the projection settlement module constructs a set of feasible solutions based on preset cumulative total constraints, power generation constraints, directional constraints or energy storage metering constraints, uses projection to determine the micro-element power, uses distinguishable methods to determine distinguishable markers and distinguishable power, calculates the power consumed within the boundary, and generates temporary settlement records or records to be settled.

[0053] S5, the recalculation and verification module establishes settlement dependencies, determines the settlement impact domain based on data change events and performs local recalculation, and then performs grid single verification according to preset grid single verification rules to generate final settlement records and settlement archive sheets.

[0054] The technical solutions in the embodiments of the present invention will be further explained below in conjunction with the specific processing procedures.

[0055] See appendix Figure 3 In this embodiment, step S1 registers the metering topology and determines the capability determination result and processing mode, collects raw data and saves it in layers, and then normalizes and segments the raw cumulative power readings to obtain normalized metering data.

[0056] S101, the metering specification module registers the metering topology for each household photovoltaic (PV) project. The registration content includes: household number, PV project identifier, bidirectional settlement energy meter identifier, power generation metering device identifier, inverter identifier, register address, rated AC output power, metering direction, multiplier, meter code resolution, equipment accuracy class, equipment sampling time source, server receiving time source, meter replacement rules, zeroing rules, rollback rules, register full scale, power generation metering point, grid-connected metering point, metering boundary, PV operation contract, authorized entity, auxiliary power ownership, line loss ownership, energy storage status, other power generation status, allowable clock error, communication recovery uncertainty window, metering resolution, commissioning time, and outage time. Each registered value is determined based on the equipment nameplate, metering verification data, register configuration, equipment operation records, and PV operation contract, and is associated with the household number and PV project identifier.

[0057] When energy storage is in use, the energy storage metering point, contract ownership rules, and energy storage metering constraints should be registered together. When the photovoltaic operation contract has defined the measurable consumption, the measurable consumption and its independent metering basis should be registered together. The sampling resolution and integral discrete error of the independent time-sharing source are determined based on the equipment sampling configuration and metering verification data. The time-sharing error variance of each channel is obtained by statistical analysis of the historical time-sharing data of the corresponding independent time-sharing source.

[0058] When calculating the variance of time-of-use errors, historical intervals with the same metering channel, the same sampling period, and cumulatively recorded confirmed results from authorized entities are selected. The sum of the cumulative electricity increment and the time-of-use electricity for each historical interval is calculated. When the reference source is time-of-use power, the time-of-use power is integrated over time according to the sampling time and sampling period. The difference between the cumulative electricity increment and the sum of the time-of-use electricity or the time integration result is divided by the number of valid time-of-use samples to obtain the time-of-use error sample. The arithmetic mean of the time-of-use error samples is calculated, and the sum of the squares of the differences between each time-of-use error sample and the arithmetic mean is divided by the number of samples minus 1 to obtain the time-of-use error variance. The preset number of historical intervals used for calculating the variance of time-of-use errors ranges from 30 to 10,000. When the actual number of usable historical intervals is less than the preset minimum number of historical intervals of 30, the square of the upper limit of time-of-use error specified in the metrological verification data is used. When the calculation result is zero, the square of the error corresponding to the sampling resolution is used to avoid division by zero in the calculation of the micro-element projection weight.

[0059] The meter code resolution equivalent is converted to the meter code resolution using a multiplier. When the multiplier represents the conversion ratio between the primary side electricity and the electricity recorded in the register, the meter code resolution equivalent is the product of the meter code resolution and the multiplier. If the multiplier is already included in the register output value, the conversion ratio is taken as 1. The maximum possible power of the channel is determined according to the rated AC output power or the contractually allowed power. The metering resolution is the meter code resolution; when the metering accuracy is specified on the power grid settlement statement, the metering resolution is taken as that metering accuracy.

[0060] The metering topology also registers preset coverage thresholds, preset rollback criteria, preset propagation rules, and grid single-check rules. The coverage threshold is set according to the sampling period, the allowable degree of data loss, and the settlement period, for example, 90% to 99%; the allowable clock error is set according to the equipment clock accuracy and clock calibration records, for example, 1s to 300s; the uncertainty window for communication recovery is set according to the data buffering and batch upload method, for example, one to six sampling periods.

[0061] The preset rollback criterion is set based on the register's full-scale range, the meter code resolution equivalent, and the channel's maximum possible power. For example, if the meter code is within one to five meter code resolutions below the register's full-scale range before rollback, and within zero to five meter code resolutions above zero after rollback, and the interval increment calculated by rollback is not greater than the product of the channel's maximum possible power and the corresponding duration plus the meter code resolution equivalent, then the preset rollback criterion is valid. The propagation rules are set based on the reference relationship between the metering channel, settlement micro-element, and electricity category; the grid single-checking rules are set based on the metering resolution, or the metering accuracy and rounding rules specified in the grid settlement statement.

[0062] When a channel has a cumulative register sequence that can be normalized according to multiplier, metering direction, and event segmentation information, the authoritative anchor capability is present. When a channel has a time-of-use power sequence or time-of-use electricity sequence that is not derived from the cumulative anchor to be verified, the independent source capability is present. When at least two of the generation channel, grid connection channel, and off-grid channel have authoritative anchor capabilities, and at least one channel has independent source capabilities, the multi-channel joint capability is present. The authoritative anchor capability, independent source capability, and multi-channel joint capability are stored as capability fields for the channel or project, respectively.

[0063] When only cumulative anchor points for the start and end of the settlement period exist and the independent source capabilities of each channel are unavailable, the processing mode is the total volume mode, outputting only the cumulative total volume. When the combined capability of multiple channels is unavailable, and only one channel has both authoritative anchor point capability and independent source capability, the processing mode is the single-channel mode, performing residual analysis and correction only on that channel. When the combined capability of multiple channels is available, the processing mode is the multi-channel mode, performing multi-channel residual signature, physical feasibility determination, and price segment identification determination. If the capability determination result does not meet the entry conditions of the above processing modes, automatic residual analysis and correction are not performed, the already obtained cumulative total volume is retained, and the relevant settlement micro-elements are set to pending settlement, awaiting re-determination after the addition of new cumulative records or independent time-sharing sources.

[0064] When the energy storage status is "no energy storage" and the other power source status is "no other power source," the boundary absorption activation status is calculated based on the difference. When the energy storage status is "with energy storage" or the other power source status is "with other power source," and the measurable absorption capacity has been registered and there is an independent absorption measurement basis, the content in the directional constraint that the on-grid electricity is not greater than the power generation of the corresponding settlement micro-element is disabled, and the boundary absorption activation status is calculated based on the measurable absorption capacity; when no measurable absorption capacity has been registered, this content in the directional constraint is disabled, and the boundary absorption activation status is disabled.

[0065] S102, the metering specification module collects raw data according to the metering topology and processing mode, including raw cumulative power readings, grid connection and grid disconnection meter codes, time-of-use power, time-of-use electricity, device sampling time, server receiving time, data quality code, device event logs, clock calibration records, batch upload identifiers, communication recovery events, communication verification results, sampling period, price version, contract version, settlement period, grid settlement statement, issued accounts, reversal records, data change events, authorized entity confirmation results, and power generation metering confirmation status; when the energy storage status is that energy storage is present, the energy storage charging power and energy storage discharging power are collected simultaneously.

[0066] The equipment's event log includes events such as meter replacement, zeroing, rollback, commissioning, shutdown, and communication recovery; batch upload identifiers include batch identifiers and the sequence number or relative sampling time of packets that can be obtained within the batch; reversal records save the reversal data that has been generated; data change events are generated by data access logs, version change records, and manual confirmation records.

[0067] The original cumulative power reading, time-of-use power, time-of-use electricity, device sampling time, server receiving time, data quality code, device event log, clock calibration record, batch upload identifier, communication recovery event, communication verification result, sampling period, energy storage charging power and energy storage discharging power are recorded as original metering data and are only written to the original metering layer.

[0068] Standardized electricity sequences and event segmentation information are written only to the standardized metering layer, calibration views are written only to the calibration view layer, temporary settlement records are written only to the temporary settlement layer, and pending settlement records are written only to the pending settlement layer. After grid single verification is completed and a settlement archive is generated, the final settlement record and settlement archive are written to the final settlement layer. Records at each layer are associated with each other through the account number, photovoltaic original cumulative electricity reading project identifier, metering channel, and time range. Subsequent processing does not overwrite the original metering records.

[0069] S103, the metering specification module converts the original cumulative electricity reading into a standardized cumulative electricity reading based on the multiplier. When a meter is replaced, reset, or register is rolled back, the cumulative segment is divided by event time according to the device's event record and meter replacement, reset, or rollback rules. The cumulative difference is only calculated within the same cumulative segment, and the effective increments of each cumulative segment are added in chronological order.

[0070] When replacing the meter, the old meter's baseline and the new meter's baseline are respectively assigned to the cumulative segments on both sides of the meter replacement time, and the increments before the old meter is deactivated and after the new meter is activated are calculated respectively; when resetting, the increments corresponding to the valid values ​​before resetting and the increments corresponding to the starting values ​​after resetting are calculated respectively, so that the valid values ​​before resetting and the starting values ​​after resetting are connected in segments.

[0071] If the cumulative value is less than the cumulative value of the previous anchor point, and there is a rollback event or the preset rollback criterion is met, the remaining rollback amount is obtained by subtracting the pre-rollback table code from the full scale of the register. The remaining rollback amount is then added to the starting value after rollback to obtain the corresponding interval increment.

[0072] If the cumulative value is less than the cumulative value of the previous anchor point, but the rollback event does not exist and the preset rollback criterion is not met, the corresponding cumulative segment will be marked as pending confirmation. If the event record of the device does not specify the event time that can be determined, or the valid table code before and after the event is incomplete, the corresponding cumulative segment will also be marked as pending confirmation, and the cumulative difference will not be calculated across the event.

[0073] After completing the calculation of each cumulative segment, a standardized electricity sequence is formed in chronological order, and the event type, event time, cumulative segment boundary and processing result are recorded to form event segment information. The standardized metering data is obtained by combining the standardized electricity sequence and the event segment information.

[0074] See appendix Figure 4 In this embodiment, step S2 extracts the authoritative cumulative anchor point from the normalized power sequence and event segmentation information, calculates the authoritative power increment, merges the time boundary to generate the settlement micro-element, and then verifies the independent time-sharing source to obtain the reference power increment and the reference micro-element power.

[0075] The authoritative cumulative anchor point is used to determine the cumulative total within the anchored interval, while the independent time-division source is used to characterize the time distribution of the cumulative total within the anchored interval. The independent time-division source does not replace the authoritative cumulative anchor point, nor does it change the authoritative power increment; subsequent processing uses the authoritative power increment as the cumulative total constraint, and the reference power increment and reference infinitesimal power as reference data for residual calculation and infinitesimal power closure projection.

[0076] S201, the anchor point micro-element module filters authoritative cumulative anchor points based on the normalized energy sequence and event segmentation information. The cumulative records participating in the filtering come from cumulative segments not marked as pending confirmation and meet the following basic conditions:

[0077] The cumulatively recorded register identifier is consistent with the register address registered in the metering topology; the device identifier is consistent with the identifier of the bidirectional settlement energy meter, the generator metering device, or the inverter; the metering direction is consistent with the metering direction registered in the metering topology; the data quality code is normal or an equivalent usable code; the timestamp is no earlier than the commissioning time and no later than the shutdown time; and the communication verification result is passed.

[0078] If there are no unprocessed table replacement, zeroing, rollback, commissioning, or shutdown events between the starting boundary of the cumulative segment and the cumulative record, the cumulative record will be used as the starting anchor of the authoritative cumulative anchor sequence, and no incremental determination relative to the previous authoritative cumulative anchor will be performed.

[0079] Starting from the second authoritative cumulative anchor point, it must also meet the following conditions: there are no unprocessed table replacement, zeroing, rollback, commissioning, or shutdown events between the cumulative record and the previous authoritative cumulative anchor point; the increment relative to the previous authoritative cumulative anchor point is not less than zero, and not greater than the product of the channel's maximum possible power and the anchoring interval duration plus the table code resolution equivalent.

[0080] The equivalent usable code is determined according to the communication protocol of the metering device, indicating that the cumulative record is complete, the value is valid, and there are no data anomalies that affect the calculation of the cumulative increment. Meter replacement, zeroing, or rollback events that have been segmented and have had their valid increments calculated based on the event segmentation information are not considered unprocessed equipment events; when the event time cannot be determined, the valid meter codes before and after the event are incomplete, or the corresponding cumulative segment is marked as pending confirmation, the relevant cumulative record is not used as the authoritative cumulative anchor point.

[0081] The maximum possible power of the power generation channel is taken as the rated AC output power, and the maximum possible power of the grid connection channel and the downstream channel is taken as the maximum power allowed by the contract or the maximum power allowed by the equipment. When neither the contract nor the equipment specifies the maximum allowable power, the power upper limit is not used to exclude cumulative records; only the determination that the increment is not less than zero is retained, and the equivalent of the meter resolution is used as the numerical comparison tolerance.

[0082] An anchoring interval is formed between two adjacent authoritative cumulative anchor points in the same metering channel. The duration of the anchoring interval is the time difference between the two authoritative cumulative anchor points. When there is no meter replacement, zeroing, or register rollback, the authoritative electricity increment is the normalized cumulative electricity of the later authoritative cumulative anchor point minus the normalized cumulative electricity of the earlier authoritative cumulative anchor point. When there is a meter replacement, zeroing, or register rollback, the effective increment of each cumulative segment is calculated based on the event segmentation information and then added in chronological order. The authoritative electricity increments for the generation channel, grid connection channel, and grid disconnection channel are mapped to the generation increment, grid connection increment, and grid disconnection increment, respectively. All authoritative cumulative anchor points are sorted according to the timestamps in the normalized electricity sequence to form an authoritative cumulative anchor point sequence. The authoritative cumulative anchor points involved in the extraction are simultaneously recorded as the original cumulative anchor points.

[0083] S202, the anchor point micro-element module extracts the authoritative cumulative anchor point time, the effective and ineffective times of the price version, the effective and ineffective times of the contract version, the start and end times of the settlement cycle, and the times of table replacement, zeroing, rollback, commissioning, and shutdown for each channel, and uses these times as time boundaries. All time boundaries are sorted using the same time base; when multiple time boundaries correspond to the same moment, they are merged into one time boundary, while retaining the corresponding authoritative cumulative anchor point, price version, contract version, and equipment event records.

[0084] A settlement micro-element is formed between two adjacent time boundaries, and the duration of the settlement micro-element is the difference between the two time boundaries. Each settlement micro-element falls only within one anchored interval, one price version, one contract version, and one device status of the corresponding channel, and does not cross the change time of the price version, contract version, or device status. When the price version or contract version lacks an executable effective time or invalidation time, the relevant time range cannot form a definite price segment or contract version affiliation, and the relevant settlement micro-element is set to pending settlement.

[0085] For each metering channel, settlement micro-elements within the same anchoring interval are recorded, forming a channel micro-element set; for each price segment, settlement micro-elements covered by that price segment are recorded, forming a segment micro-element set. A correspondence is established between settlement micro-elements and anchoring intervals, price segments, and contract versions through start and end times, forming the settlement micro-element division.

[0086] S203, the anchor point micro-element module verifies independent time-sharing sources based on authoritative cumulative anchor point sequences and anchoring intervals, combined with time-of-use power and time-of-use electricity. The independent time-of-use source for the power generation channel is taken from the inverter's AC power integral, or from the power generation time-of-use metering independent of the settlement cumulative register; the independent time-of-use sources for the grid connection channel and the grid disconnection channel are taken from time-of-use telemetry using different sampling periods or different register systems than the cumulative anchor points. If the time-of-use value is obtained by cutting or differentially dividing the same cumulative register sequence, the independent time-of-use source does not exist.

[0087] When the device's sampling time source is valid, the timestamp of the independent time-division source is taken from the device's sampling time; when the device's sampling time source is invalid, the timestamp of the independent time-division source is taken from the server's receiving time, and the receiving time mapping is listed as the preferred candidate in the subsequent S305 time mapping candidates. The device's sampling time source is valid if the clock calibration record is valid and the clock difference reflected in the clock calibration record is not greater than the allowable clock error; the device's sampling time source is invalid if the clock calibration record is missing, invalid, or the clock difference is greater than the allowable clock error.

[0088] A valid time-sharing sample is one whose data quality code is normal or equivalent, whose timestamp is within the anchor interval, and whose communication verification result is passed. When time-sharing power or time-sharing electricity represents the quantity within a sampling period, the overlapping portion of the corresponding sampling period within the anchor interval is used as the coverage period. When time-sharing power represents instantaneous power, the period between the sampling times of two adjacent valid time-sharing samples is used as the coverage period. If the difference between adjacent sampling times is greater than the sum of the allowable clock error and a sampling period, this period is not included in the valid coverage duration. The effective coverage duration is obtained by merging the coverage periods and deducting the overlapping portion.

[0089] An independent time-division source is considered available if it exists, the ratio of its effective coverage duration to the anchor interval duration is not less than the coverage threshold, and there are no undefinable consecutive missing data points. If any of these conditions are not met, the independent source is marked as unavailable and recorded as unavailable. Undefinable consecutive missing data points refer to situations where the start and end of the missing period cannot be determined based on the data quality code, device event records, or adjacent effective samples, or where the length of the consecutive missing data points exceeds the sum of the allowable clock error and one sampling period. The coverage threshold is set based on the sampling period, settlement period, and the allowable degree of data loss, for example, 90% to 99%.

[0090] When an independent source is marked as available and the reference source is time-of-use electricity, the time-of-use electricity within the anchoring interval is added together to obtain the reference electricity increment. When the time-of-use electricity crosses the boundary of the settlement micro-element, it is allocated according to the ratio of the sampling period to the overlap duration of each settlement micro-element to obtain the reference micro-element electricity.

[0091] When the reference source is time-division power, if the power value represents the average power over the sampling period, multiply the power value by the corresponding effective duration; if the power value represents the instantaneous power, multiply the arithmetic mean of adjacent power values ​​by the difference between adjacent sampling times. Add the calculation results of each effective time period within the anchored interval to obtain the reference power increment; divide the effective time period according to the start and end times of the settlement micro-element and use the same calculation method to obtain the reference micro-element power. If the difference between adjacent sampling times is greater than the sum of the allowable clock error and one sampling period, integration is not performed across the corresponding consecutive missing time periods.

[0092] Automatic anchor point correction is performed only when the independent source availability flag is available. When the independent source availability flag is unavailable, the cumulative total of the corresponding anchoring interval is retained; if the anchoring interval corresponds to only one price version and the status is neither pending confirmation nor physical relationship abnormal, it enters the cumulative total entry of subsequent S401; in other cases, the settlement micro-element to which the anchoring interval belongs is set to pending settlement.

[0093] See appendix Figure 5 In this embodiment, step S3 calculates the allowable deviation and residual of the interval, identifies the reverse residual pairs and forms a multi-channel residual signature, and then performs anchor point correction or time mapping according to the anomaly type.

[0094] When a deviation occurs in the intermediate cumulative anchor point value, the deviation is recorded in the two anchoring intervals before and after the intermediate cumulative anchor point, causing the residuals of the two intervals to be in opposite directions and thus summing up and closing the gap. Time alignment anomalies change the time assignment of time-sharing data; using corresponding time mapping can reduce the interval residuals of the relevant channels. The anomaly diagnosis module distinguishes between anchor point value anomalies and time alignment anomalies based on this, and determines the subsequent processing path in conjunction with the physical feasibility status.

[0095] S301, the anomaly diagnosis module converts the accuracy levels of the metering equipment used to generate authoritative power increments and the metering equipment corresponding to independent time-division sources into relative error limits. The first component is obtained by multiplying the larger of the authoritative power increment and the reference power increment by the sum of the two relative error limits, and then adding the meter code resolution equivalent. The second component is obtained by multiplying the product of the allowable clock error and the maximum possible power of the channel. The third component is obtained by summing the three components to obtain the allowable deviation for the interval. Each parameter uses the registered values ​​from the metrological verification data and the metrological topology.

[0096] The interval residual is obtained by subtracting the reference electricity increment from the authoritative electricity increment. When the absolute value of the interval residual is not greater than the allowable deviation of the interval, the anchored interval is marked as a closable interval; when the absolute value of the interval residual is greater than the allowable deviation of the interval, proceed to S302. When an independent source is marked as unavailable, it is treated as a time-sharing data missing data; when there are undefinable continuous missing data, the already calculated allowable deviation of the interval is retained, and automatic anchor point correction is not performed on that interval.

[0097] S302, the anomaly diagnosis module uses a reverse method to identify reverse residual pairs between two adjacent anchoring intervals that share the same intermediate cumulative anchor point:

[0098] ;

[0099] In the formula: For any metering channel in the set of metering channels; The sequence number of the anchored interval; , Channels In the anchorage area Anchorage area The independent source availability identifier, with a value of 1 indicating availability and a value of 0 indicating unavailability; , These are the interval residuals for the two anchored intervals, respectively. , These are the allowable deviations between the two anchoring intervals; This is for absolute value operations.

[0100] When all conditions in the reverse approach are met, the intermediate cumulative anchor point is identified as a candidate anomaly anchor point, and the reverse residual direction status is recorded according to the sign of the interval residual of the previous anchoring interval. If the absolute values ​​of the residuals on both sides are greater than their respective interval allowable deviations but in the same direction, or if the absolute value of the sum of the residuals on both sides is greater than the sum of the two interval allowable deviations, a time mapping candidate is generated and tested. If only one independent source is available and its interval residual absolute value is greater than the interval allowable deviation, the relevant anchoring interval is merged into a pending confirmation phase, and recalculated after late data arrives. If the absolute values ​​of the interval residuals of multiple consecutive anchoring intervals are greater than their respective interval allowable deviations, all reverse residual pairs are identified and merged into a continuous anomaly anchor point group.

[0101] When directional constraints are not disabled, a physically feasible state is considered physically feasible if: the incremental power generation is not greater than the sum of the incremental power generation and the allowable deviation between the power generation channel and the power grid channel; the incremental power generation is not greater than the product of the rated AC output power and the anchoring interval duration plus the allowable deviation between the power generation channel and the grid channel; and the incremental power discharging is not less than the negative of the allowable deviation between the discharging channel and the grid channel. If any of these relationships are not met, the state is physically infeasible; and if necessary channel data is missing, the state is undeterminable. When directional constraints are disabled, the physically feasible state is determined based on energy storage metering constraints, energy storage charging capacity, and energy storage discharging capacity. The independent source availability identifiers, reverse residual directional states, and physically feasible states of each channel at the same time location are summarized to form a multi-channel residual signature.

[0102] Anomaly types are determined sequentially and mutually exclusively. If the physical feasible state is physically infeasible, the reverse formula is invalid, and the absolute value of the interval residual of each channel after all time mapping candidate trials is still greater than its respective interval allowable deviation, the anomaly type is physical relationship anomaly, and the relevant anchoring interval is set to settlement blocking and merged into pending confirmation.

[0103] If there is a unique time mapping result that reduces the sum of the absolute values ​​of the interval residuals and the absolute values ​​of the interval residuals of each channel after mapping are not greater than their respective interval allowable deviations, the anomaly type is time alignment anomaly and enters S305; if there are multiple optimal time mapping candidates, and the absolute value of the difference between the electricity consumption of the price segment formed by any two optimal time mapping candidates is not greater than the metering resolution, they are considered to have formed the same settlement result.

[0104] When at least one channel forms a reverse residual pair and is not a time alignment anomaly, the anomaly type is anchor value anomaly. A single candidate anomaly anchor point enters S303, and a group of consecutive anomaly anchor points enters S304. When an independent source is marked as unavailable or has undefined consecutive missing data, the anomaly type is time-sharing data missing, and the cumulative total is retained. When it corresponds to only one price version and its status is neither pending confirmation nor physical relationship anomaly, it enters S401, and the remaining settlement micro-elements are set to pending settlement.

[0105] Cases where the absolute value of the residual in other intervals exceeds the allowable deviation of the interval are identified as unidentified anomalies and included in the pending confirmation. In single-channel mode, the multi-channel residual signature consists of the independent source availability identifier of the corresponding channel, the reverse residual direction status, and the power determination result.

[0106] If both sides of an authoritative cumulative anchor point are closable intervals and the authoritative cumulative anchor point is not a candidate abnormal anchor point, or if both sides are closable intervals after S303 or S304 correction, the authoritative cumulative anchor point is determined as a reliable anchor point and is only used as a fixed boundary for S304 joint correction.

[0107] S303, the anomaly diagnosis module takes the anchor point correction weight as the reciprocal of the square of the interval allowable deviation; when the interval allowable deviation is zero, the table resolution equivalent is used to replace the interval allowable deviation. The anchor point correction amount is calculated using a correction formula:

[0108] ;

[0109] In the formula: For channel In the anchorage area Anchor point correction amount at the intermediate cumulative anchor point; , Channels In the anchorage area Anchorage area Anchor point correction weights; , These are the interval residuals for the corresponding anchored intervals.

[0110] The original value of the authoritative cumulative anchor point is added to the anchor point correction amount to obtain the corrected anchor point value and the corrected interval increment, and the total corrected amount is accumulated. An automatic correction result is formed and written to the correction view layer when the total electricity of two adjacent anchored intervals remains unchanged, the corrected cumulative sequence is not decreasing, the corrected interval increment is not less than zero, the power generation is not greater than the product of the rated AC output power and the anchored interval duration, a feasible solution exists for the multi-channel energy relationship, and the absolute values ​​of the residuals of the two recalculated intervals are not greater than their respective interval allowable deviations. If any of these conditions are not met, the automatic correction result is revoked, and the relevant authoritative cumulative anchor point is marked as an anchor point to be confirmed and merged into the pending confirmation category.

[0111] S304, the anomaly diagnosis and calibration module uses the reliable anchor points on both sides of the continuous anomaly anchor point group as fixed boundaries, and adopts a group calibration method to solve the anchor point correction amount:

[0112] ;

[0113] In the formula: For channel The set of anchor point correction values ​​within a group of consecutive abnormal anchor points; This refers to the cumulative anchor point number in the middle. , These are the sequence numbers of the first and last intermediate cumulative anchor points to be corrected, respectively. to It is a closed interval; , As reliable anchor points on both sides, , Take zero; , These are the anchor point correction values ​​for the current intermediate cumulative anchor point and the previous intermediate cumulative anchor point, respectively. Adjust the weights for the anchor points; For interval residuals; This represents the independent variable that minimizes the objective function; Indicates a closed interval to Sum.

[0114] Take the partial derivatives of the objective function with respect to the corrections at each anchor point and set the partial derivatives to zero to form a system of linear equations; the main diagonal coefficients are the sum of the correction weights of adjacent anchor points, the adjacent non-main diagonal coefficients are the negatives of the correction weights of the corresponding anchor points, and the corrections at the boundary anchor points are set to zero.

[0115] Solving the linear equations yields the anchor point correction values ​​within the continuous abnormal anchor point group. When the cumulative sequence is non-decreasing, the increment of the corrected interval is not less than zero, the power generation is not greater than the product of the rated AC output power and the anchor interval duration, and there is a feasible solution for the multi-channel energy relationship with all relevant intervals being closable, the automatic correction results are written to the correction view layer. If any credible anchor point is missing on either side or any condition is not met, the entire continuous abnormal anchor point group is marked as an anchor point to be confirmed and merged into the unconfirmed group.

[0116] S305, the anomaly diagnosis and calibration module generates zero offset mapping, receiving time substitution mapping, overall clock calibration difference shift mapping, relative timing rearrangement mapping within batch upload packets, and finite shift mapping based on the device's sampling time, server's receiving time, clock calibration records, batch upload flags, and communication recovery events.

[0117] The translation range of the finite translation mapping is no greater than the larger of the allowable clock error and the uncertainty window for communication recovery. The sampling period is used as the enumeration interval, and the two endpoints of the translation range, the zero offset, and the clock difference recorded in the clock calibration record are respectively included in the time mapping candidates. The relative timing within the batch upload packets is rearranged according to the packet sequence number or relative sampling time, maintaining timing monotonicity. The time mapping candidates only adjust the settlement micro-element to which the telemetry data belongs, and the adjustment range is limited to events between adjacent devices. The original cumulative electricity reading and price version boundary remain unchanged.

[0118] For each time mapping candidate, recalculate the interval residuals of the relevant channels, calculate the sum of squared interval residuals according to the anchor point correction weight, and select the mapping with the smallest calculation result as the optimal mapping. If the optimal mapping is unique, or if multiple optimal mappings exist and the absolute value of the difference between the electricity consumption of the price segment formed by any two optimal mappings is not greater than the metering resolution; if the absolute value of the interval residuals of each channel after mapping is not greater than their respective interval allowable deviation; and if the physical feasibility state after mapping is physically feasible, or if the physical feasibility state is undeterminable but the sum of the absolute values ​​of the interval residuals is not greater than before mapping, then generate a time mapping version and update the reference micro-element electricity consumption attribution.

[0119] When multiple optimal mappings meet the conditions, the candidate time mapping with the smallest absolute shift value is selected; when the absolute shift values ​​are the same, the time mapping version is determined according to the preset mapping selection order; the mapping selection order is as follows: zero offset mapping, receive time substitution mapping, clock calibration difference overall shift mapping, batch upload packet relative timing rearrangement mapping, and finite shift mapping.

[0120] When the optimal mapping is not unique and its impact on the electricity consumption of the price range is greater than the metering resolution, a non-unique mapping state is formed and merged into the pending confirmation state; when the absolute value of the interval residual of any related channel after mapping is still greater than the allowable interval deviation of that channel, the anomaly type determination of S302 is re-executed and processed according to the path corresponding to the anchor value anomaly, physical relationship anomaly, or unlocated anomaly.

[0121] See appendix Figure 6 In this embodiment, step S4 constructs a set of constrained feasible solutions, uses projection to determine the infinitesimal energy, uses distinguishable method to determine distinguishable markers and distinguishable energy, then calculates the energy absorbed within the boundary and generates settlement records for each layer.

[0122] The authoritative electricity increment is used to limit the cumulative total within the anchored interval, while the reference micro-element electricity is used to characterize the reference distribution of the cumulative total within each settlement micro-element. Closed projection does not change the correspondence between the cumulative total and the settlement micro-element and the price version; when there are multiple combinations of micro-element electricity in the constrained feasible solution set, the projection formula only generates temporary calculation results, and whether the electricity in the price segment can be settled is determined by the discriminative formula.

[0123] S401, the projection settlement module receives the closed interval, the anchor interval that has formed an automatic correction result, the anchor interval that has adopted the time mapping version, and the anchor interval that has an independent source available identifier that is unavailable but the cumulative total is valid, corresponds to only one price version and is not in a physical relationship abnormal state, or is pending confirmation or has an anchor point pending confirmation, and sets the micro-element electricity variable according to the metering channel and settlement micro-element.

[0124] When the anchorage boundaries of each channel are inconsistent, the anchorages that overlap and whose cumulative total constraints have not yet been closed, along with their settlement micro-elements, will be included in the same solution until all settlement micro-elements corresponding to each cumulative total constraint are fully included.

[0125] When an independent source is marked as available, the projection weight of the micro-element is taken as the reciprocal of the time-division error variance of the corresponding channel; when the time-division error variance is zero, the square of the error corresponding to the sampling resolution is used instead. When an independent source is marked as unavailable, the projection weight of the micro-element is zero, the reference micro-element charge is recorded as zero and used only for the numerical expression of the projection; since the projection weight of the micro-element is zero, the reference micro-element charge does not affect the objective function and the solution result, and the corresponding channel only participates in the price segment identification through the cumulative total constraint.

[0126] The feasible solution set consists of combinations of infinitesimal quantities that satisfy all constraints. The cumulative total constraint requires that the sum of the infinitesimal quantities of all settlement infinitesimal quantities within the same anchoring interval equals the authoritative quantity increment; if a corrected interval increment has already been formed, it equals the corrected interval increment. The power generation constraint requires that the power generation of each settlement infinitesimal quantity is not less than zero and not greater than the product of the rated AC output power and the settlement infinitesimal quantity's duration. When the directional constraint is not disabled, the on-grid quantity of each settlement infinitesimal quantity is not less than zero and not greater than the power generation of the corresponding settlement infinitesimal quantity, and the off-grid quantity is not less than zero.

[0127] When directional constraints are disabled, the energy storage charging and discharging power is allocated to the corresponding settlement micro-element according to the overlap between the data collection period and the settlement micro-element. Then, based on the energy storage metering constraints registered in the metering topology, the constraint relationship between the energy of each channel micro-element and the energy storage charging and discharging power is established. If the photovoltaic operation contract does not provide executable energy storage metering constraints, only the requirement that the cumulative total of each channel be closed and the energy of each settlement micro-element not be less than zero is retained, and the cross-channel time-sharing relationship of the corresponding anchoring interval is marked as not automatically confirmed.

[0128] When the feasible solution set is not empty, the projection method is used to solve for the infinitesimal charge element:

[0129] ;

[0130] In the formula: This represents the collection of micro-element electricity within each settlement micro-element for each channel; To constrain the combination of infinitesimal electric charges in the feasible solution set; For the set of feasible solutions to be constrained; This represents the summation of all settlement infinitesimals included in this solution; This represents the summation of all metering channels in the metering channel set; This is the serial number of the settlement micro-element; It is a collection of metering channels, including power generation channels, grid connection channels, and grid disconnection channels; For channel In settlement micro yuan The internal micro-elemental charge; For channel In settlement micro yuan Reference micro-element charge within; For channel In settlement micro yuan The projection weights of the infinitesimal elements within.

[0131] An equality constraint matrix is ​​established based on the cumulative total constraint, and an inequality constraint matrix is ​​established based on the power generation constraint, directional constraint, or energy storage metering constraint. A feasibility test is conducted to determine whether a combination of infinitesimal power elements that simultaneously satisfies all constraints exists. If it does, an initial feasible solution is obtained by first calculating a combination of infinitesimal power elements that satisfies all constraints using the equality and inequality constraint matrices. Then, the inequality constraints reaching the boundary are added to the effective constraint set, and the corresponding linear equations are solved. Iteration continues until all constraints are satisfied and there are no feasible adjustment directions that can further reduce the objective function, thus obtaining the infinitesimal power element.

[0132] Micro-element electricity quantities are only written to the temporary settlement layer. When there are multiple combinations of micro-element electricity quantities with the same minimum objective function value, one of these combinations is retained as a temporary calculation result, and the electricity quantity in the price segment is still determined based on the entire set of feasible solutions to the constraints. When the set of feasible solutions to the constraints is empty, the anomaly type is determined to be a physical relationship anomaly, the cumulative total is recorded as the total to be verified, and the relevant settlement micro-element is set to pending settlement.

[0133] S402, the projection settlement module determines whether the electricity consumption of each price segment is identifiable based on the set of feasible solutions, price segments, channel micro-element sets, and segment micro-element sets. When the set of feasible solutions is empty, the settlement micro-element to which the relevant price segment belongs is set to pending settlement; when the set of feasible solutions is not empty, the lower bound and upper bound of the segment electricity consumption are calculated using the identifiability formula.

[0134] ;

[0135] ;

[0136] In the formula: For channel In price range The lower limit of the section's electricity level; For channel In price range The upper limit of the section's electricity volume; and These are the minimum value operation and the maximum value operation, respectively. Price range The set of settlement micro-elements covered; For channel In settlement micro yuan The internal micro-elemental charge; Indicates price range Summing all settlement micro-elements covered.

[0137] The goal is to find the minimum and maximum values ​​of the sum of infinitesimal electricity quantities within each price segment, within the same set of feasible solutions with the same constraints. When a price segment spans multiple anchor intervals, if there is no cross-interval coupling relationship between the constraints of the anchor intervals, the lower and upper bounds of the segment electricity quantities for each intersection range are calculated and then added together. If there is a cross-interval coupling relationship, the lower and upper bounds of the segment electricity quantities for that price segment are directly calculated within the same set of feasible solutions covering all relevant anchor intervals.

[0138] When the difference between the upper and lower bounds of the segment's electricity volume is not greater than the metering resolution of the corresponding channel, the identifiable mark is marked as identifiable. The arithmetic mean of the two is taken and rounded according to the metering resolution to obtain the identifiable electricity volume. The metering resolution is determined based on the meter code resolution or the metering accuracy specified in the power grid settlement statement, with the specific value based on the metering topology registration value of the corresponding metering channel. When the difference between the upper and lower bounds of the segment's electricity volume is greater than the metering resolution, the relevant settlement micro-element is set to pending settlement, and the cumulative total is retained. It will be re-determined after late data arrives, a new authoritative cumulative anchor point is added, or confirmation results are obtained from an authorized entity. When the entire anchoring interval corresponds to only one price version, a unique identifiable electricity volume is still formed in the above manner.

[0139] S403, the projection settlement module calculates the electricity consumption within the boundary based on the micro-element electricity and identifiable markers, combined with the boundary consumption activation status, metering boundary, auxiliary power attribution, and line loss attribution. When the boundary consumption activation status is calculated by difference, the electricity generated within the same settlement micro-element is subtracted from the electricity fed into the grid; when the boundary consumption activation status is calculated by measurable consumption, the measurable consumption that has been registered and independently metered is taken.

[0140] When the metering period of the measurable consumption corresponds to only one settlement micro-element, the measurable consumption is included in the corresponding settlement micro-element; when the metering period of the measurable consumption spans multiple settlement micro-elements, if there are time-sharing consumption metering records, each time-sharing consumption metering record is included in the corresponding settlement micro-element according to the sampling time; if the time-sharing consumption metering records cross the boundary of the settlement micro-element, they are allocated according to the proportion of the overlap between the sampling period and each settlement micro-element; if there is only a cumulative value of the metering period, it is allocated according to the proportion of the overlap between the metering period and each settlement micro-element.

[0141] After allocation, the sum of the measurable consumption of each settlement micro-element should equal the measurable consumption before allocation. If the start and end times of the measurement period cannot be determined, or if the allocation result cannot satisfy the cumulative total closure, the relevant settlement micro-element is set to pending settlement. When the boundary consumption activation status is disabled, no settlement item corresponding to boundary consumption is generated.

[0142] The electricity consumed within the boundary refers to the amount consumed within the metering boundary as defined in the photovoltaic operation contract. When used for operation charges under the photovoltaic operation contract, the power generation metering confirmation status is "confirmed," the on-grid electricity is taken from the two-way settlement electricity meter, and the corresponding settlement micro-element meets the temporary settlement access requirements. When the boundary consumption activation status is "calculated by difference" and the measurable consumption has been registered, the measurable consumption allocated to the corresponding settlement micro-element is used as the upper limit for verification; if the electricity consumed within the boundary is less than zero or greater than this upper limit, no settlement account corresponding to the boundary consumption of that settlement micro-element is generated, and that settlement micro-element is set to pending settlement.

[0143] When the boundary consumption is enabled and calculated based on the measurable consumption amount, the measurable consumption amount allocated to the corresponding settlement micro-element and confirmed by independent metering is taken as the boundary consumption amount, and the upper limit of this value is no longer checked repeatedly.

[0144] S404, the projection settlement module performs transaction settlement based on the micro-element electricity, identifiable markers, identifiable electricity, and the electricity consumed within the boundary, combined with the price version and contract version. It reads the on-grid price, purchase price, and consumption settlement price from the price version to which the settlement micro-element belongs. When the identifiable marker for the corresponding price segment is identifiable, the settlement amount is calculated using the identifiable electricity of the corresponding channel and the electricity consumed within the boundary determined based on the identifiable electricity. When the anchored interval corresponds to only one price version and is a closable interval, the settlement amount is calculated using the total micro-element electricity within that price version. The micro-element electricity of each settlement micro-element is only counted in its respective price version.

[0145] The three amounts are summed according to the revenue and expenditure directions specified in the contract version to obtain the settlement amount, which is then recorded in the settlement status. Settlement within the boundary is only applicable to scenarios where a photovoltaic operation contract exists and the metering boundary is clearly defined.

[0146] Temporary settlement access is determined separately according to electricity volume category. The settlement accounts corresponding to the electricity sales settlement amount and the electricity purchase settlement amount should meet the aforementioned price identifiability condition, the status should not be settlement blocked, the status should not be pending confirmation, the pending confirmation anchor point or non-unique mapping status, and the cross-channel time-sharing relationship should not be non-automatically confirmable; in addition to meeting the above conditions, the settlement accounts corresponding to consumption within the boundary should also meet the conditions for activation of consumption within the boundary.

[0147] Settlement micro-element that meets the temporary settlement access criteria forms time-of-use accounts to be verified. The settlement amount, electricity sales settlement amount, electricity purchase settlement amount, and consumption settlement amount are written into the temporary settlement layer to form a temporary settlement record. Settlement micro-element that does not meet the temporary settlement access criteria is written into the pending settlement layer to form a pending settlement record.

[0148] If the settlement for consumption within the boundary does not meet the conditions for activation of consumption within the boundary, but the settlement amount for electricity sales or purchase meets the corresponding access conditions, only the settlement entries corresponding to consumption within the boundary will be written to the pending settlement layer. The settlement entries corresponding to electricity sales or purchase will still form time-of-use entries pending verification. The settlement entries in each layer of settlement records are marked with the electricity category.

[0149] See appendix Figure 7 In this embodiment, step S5 establishes settlement dependencies, determines the settlement impact domain based on data change events and performs local re-settlement, then performs grid single verification, and generates the final settlement record and settlement archive sheet.

[0150] Settlement dependencies are used to record the reference relationships between metering data, anomaly types, settlement micro-elements, and settlement items. When data changes, the settlement impact domain is determined based on the anomaly type and the electricity propagation relationship, and only the data within the settlement impact domain is reprocessed. Before performing grid bill reconciliation, the reconciliation scope between the grid settlement bill and the local reconciliation quantity is unified according to the account number, photovoltaic project identifier, metering point, metering direction, settlement cycle, and metering boundary.

[0151] In S501, the recalculation and verification module writes the original cumulative anchor point, correction view, time mapping version, settlement micro-element, independent time-sharing source, independent source availability identifier, multi-channel residual signature, anomaly type, identifiable marker, settlement status, settlement amount, electricity sales settlement amount, electricity purchase settlement amount, consumption settlement amount, and electricity category into the dependency information of the corresponding settlement item. The dependency information also includes the account number, photovoltaic project identifier, anchoring range, price version, contract version, processing mode, and the relationship between the original item and the replacement item.

[0152] Each settlement item is associated with the metering source according to the account number, photovoltaic project identifier, metering channel, settlement micro-element, and electricity category; the settlement basis according to the price version and contract version; and the processing result according to the anomaly type, identifiable marker, and settlement status. A substitute item is the entity used to replace the original item. The substitute record records the original item, the substitute item, the data change event, the correction amount, and the related relationships. Issued items, reversal records, and substitute records are used to determine whether the original item has been corrected. When a substitute record has already been generated for the same data change event, a substitute record will not be generated again for the same original item, thus forming a settlement dependency relationship.

[0153] S502, the recalculation and verification module receives data change events caused by anchor point correction, late anchor point insertion, time mapping change, price version change, contract change, or grid settlement statement verification discrepancies. Based on the settlement dependency relationship, it determines the set of affected micro-elements and the set of affected channels, and then filters and references the settlement accounts of the corresponding metering channels and settlement micro-elements according to the electricity category to form the settlement impact domain.

[0154] The propagation rules are executed according to the reference relationships between metering channels, settlement micro-elements, and electricity categories. Specifically, changes in power generation channels affect the accounts for power generation and electricity consumption within the boundary; changes in grid connection channels affect the accounts for grid-connected electricity and electricity consumption within the boundary; and changes in off-grid channels affect the accounts for off-grid electricity. When the time mapping changes, the set of affected micro-elements is the union of the old and new time mapping versions of the involved settlement micro-elements.

[0155] Price version changes affect settlement amount items within their effective scope, and recalculate electricity sales settlement amount, electricity purchase settlement amount, and electricity consumption settlement amount, while the basic metered electricity volume remains unchanged. Contract changes affect settlement rules and settlement items within the effective scope of the contract version. When a physical relationship anomaly is confirmed or resolved, all settlement items covered by the multi-channel residual signature enter the settlement impact domain. When a data change occurs in an anchored interval with a pending confirmation status, the settlement micro-element to which that anchored interval belongs is also included in the set of affected micro-elements.

[0156] S503, the recalculation and verification module limits the scope of local re-settlement based on the settlement impact domain and the type of data change event, and re-executes the authoritative cumulative anchor point extraction, settlement micro-element generation, interval residual calculation, anomaly type determination, closed projection, price segment identification determination and transaction settlement related to the data change event within the limited scope; price version changes only require re-execution of amount calculation.

[0157] When correcting a single anchor point, only the two anchor intervals before and after the authoritative cumulative anchor point are reprocessed; when a late anchor point is inserted, only the original anchor interval that was split is reprocessed, and the interval residuals of the anchor intervals that were originally in the pending confirmation state are recalculated; for a group of consecutive abnormal anchor points, the trusted anchor points on both sides are used as the processing boundary; when the time mapping changes, the union of the settlement micro-elements involved in the old time mapping version and the new time mapping version is processed.

[0158] When the price version changes, only the settlement amount, electricity sales settlement amount, electricity purchase settlement amount, and consumption settlement amount within the effective scope are recalculated. When the processing mode is downgraded, the automatic correction results that depend on the invalid independent time-of-use sources are revoked, the authoritative electricity increment corresponding to the original cumulative anchor point is restored in the correction view, the interval residuals, anomaly types, and identifiable markers within the settlement impact domain are recalculated, and time-of-use accounts that can no longer meet the temporary settlement access requirements are returned to the pending settlement layer.

[0159] Settlement items not yet issued are updated according to the recalculation results. Issued items are only corrected through replacement records; based on the relationship between the original and replacement items, combined with the issued items, reversal records, and replacement records, it is determined whether the original item has been reversed or replaced. When an item is replaced entirely, a reversal is performed; when an item does not need to be replaced entirely, the settlement amount of the replacement item is compared with that of the original item. If the settlement amount decreases, a refund is performed; if the settlement amount increases, a supplementary collection is performed. Replacement records are generated for all three. Metered electricity outside the settlement impact area remains unchanged from the settlement items, thus obtaining the resettled items.

[0160] S504, the recalculation and verification module aligns the grid settlement statement with the local verification quantity according to the household number, photovoltaic project identifier, power generation metering point, grid connection metering point or energy storage metering point, metering direction, settlement cycle, and metering boundary. If the electricity units are inconsistent, they are converted to the same units as specified in the grid settlement statement. If the verification scope cannot be consistent, the electricity quantities are not directly compared, and the differences in metering boundaries or settlement cycles are recorded.

[0161] When a corrected total exists for a channel, the local verification quantity is the corrected total confirmed by the authorized entity. If the corrected total has not been confirmed by the authorized entity, the corrected total or the original authoritative electricity increment will not be used to form the final local verification quantity for that verification range. The relevant final settlement remains pending and will be handled according to the cessation condition of the authorized entity's confirmation result being missing. When no corrected total exists for a channel, the local verification quantity is the authoritative electricity increment within the verification range. When there is a total to be verified in the anchored interval, the total to be verified is verified simultaneously, and compared with the temporary settlement record, the pending settlement record, the final settlement record, and the re-settled account items.

[0162] According to the preset grid statement reconciliation rules, the grid statement reconciliation result is consistent when the absolute value of the difference in electricity consumption between each channel is not greater than the metering resolution, or when the two values ​​are equal after processing according to the metering accuracy and rounding rules specified on the grid settlement statement. In other cases, the grid statement reconciliation result is inconsistent, and the involved metering channels, time ranges, electricity consumption categories, and reconciliation differences are recorded to form a reconciliation discrepancy. The preset grid statement reconciliation rules are set based on the metering accuracy, rounding digits, and local metering resolution specified on the grid settlement statement.

[0163] When the power grid single-verification result is consistent, the portion of the temporary settlement record that is identifiable and has no closed anomalies will be converted into the final settlement record. Unclosed anomalies are generated by the sum of pending confirmation, pending confirmation anchor points, settlement blockage, non-unique mapping status, and non-automatic confirmation status; if any of these statuses exist, the relevant settlement items will not be converted into the final settlement record.

[0164] When the power grid single-check results are inconsistent, the metering channel and time range are determined based on the check differences, resulting in an unsettled discrepancy statement, and the relevant final settlement remains pending. S301 and S302 are re-executed based on the unsettled discrepancy statement, and the subsequent processing path is determined according to the anomaly type: if the anomaly type is an anchor value anomaly and there is only a single candidate anomaly anchor point, proceed to S303; if a continuous group of anomaly anchor points is formed, proceed to S304; if the anomaly type is a time alignment anomaly, proceed to S305; if the anchoring interval meets the entry conditions of S401, proceed to S401.

[0165] When the anomaly type is physical relationship anomaly, unlocated anomaly, or when the anchoring interval is in a state of pending confirmation, pending confirmation anchor point, settlement blockage, non-unique mapping, or non-automatic confirmation, automatic correction, time mapping, or closed projection will not be performed. The corresponding anomaly type and state will be retained, and the relevant settlement items will remain pending.

[0166] After automatic correction, time mapping, or closed projection is completed, or when the anomaly type or settlement status changes, S502 determines the settlement impact domain, and S503 performs partial resettlement. When issued items need correction, a replacement record is generated; unissued settlement items are updated according to the recalculation results. If automatic processing can be completed and the local verification quantity changes, a new power grid single-item verification is performed. If the re-verification is still inconsistent, the unsettled difference sheet is updated, and the above processing is repeated. If automatic processing cannot be performed, the automatic processing does not produce a valid result, or no new data is obtained, the unsettled difference sheet is updated, and the automatic loop stop condition is determined; the same automatic correction, time mapping, or closed projection is not repeated.

[0167] The automatic loop determines whether to stop in the following order: It stops when the metering channel, time range, and difference corresponding to the discrepancy remain unchanged after two consecutive updates to the outstanding discrepancy statement; it stops when the grid settlement statement only contains summary data and cannot pinpoint the specific anchoring interval; it stops when the independent time-of-use source remains unavailable and no new authoritative cumulative anchor point is added; it stops when the confirmation result from the authorized entity is missing; and it stops when the local metering boundary differs from the metering boundary of the grid settlement statement. When a new authoritative cumulative anchor point is added, the independent time-of-use source becomes available again, or an authoritative entity confirmation result is obtained, the discrepancy is recalculated, and the number of consecutive updates is re-recorded.

[0168] After automatic cycling stops, the relevant final settlement remains pending, and unsettled discrepancies are transferred to the authorized entity for confirmation. Once the confirmation result from the authorized entity is obtained, the grid bill of exchange verification or settlement account is resumed according to the confirmation result. When the grid bill of exchange verification result is consistent and the anomaly is not closed (empty), a settlement archive is generated based on the final settlement record, the resettled account, the substitute record, the portion of the pending time-of-use account that has been transferred to the final settlement record, identifiable markers, settlement status, and the corrected total amount or authoritative electricity increment consistent with the grid bill of exchange verification. The final settlement record and the settlement archive are then written to the final settlement layer.

[0169] To further aid in understanding the technical solutions of the embodiments of the present invention, a specific application example is given below.

[0170] A residential photovoltaic (PV) project has a rated AC output power of 6kW, no energy storage, no other power sources, and its grid connection activation status is calculated based on the difference. The generation channel, grid connection channel, and off-grid connection all possess authoritative anchoring capabilities, and the generation channel has independent source capabilities. The processing mode is multi-channel mode. The time-of-use data sampling period is 15 minutes, the meter reading resolution and metering resolution are both 0.01kWh, the preset coverage threshold is 95%, and the allowable clock error is 60 seconds.

[0171] The processing period was selected from 09:00 to 11:00, with cumulative anchor points at 09:00, 10:00, 10:30, and 11:00, with 10:30 being the time of price version change. Independent time-sharing source integration of the power generation channel yielded reference electricity increments of 6.00 kWh, 3.80 kWh, and 4.20 kWh for the three anchoring intervals, respectively; the authoritative electricity increments calculated based on the original cumulative anchor points were also 6.80 kWh, 3.00 kWh, and 4.20 kWh, respectively. The allowable deviation for the first two anchoring intervals was 0.15 kWh, with corresponding interval residuals of 0.80 kWh and -0.80 kWh, respectively, satisfying the reverse formula. Therefore, the authoritative cumulative anchor point at 10:00 was identified as a candidate abnormal anchor point.

[0172] The anchor point correction weights for both anchoring intervals are the same, and the correction amount at 10:00 is -0.80 kWh obtained using the correction formula. The corrected interval increments are 6.00 kWh and 3.80 kWh respectively, and the total electricity generation for both anchoring intervals remains 9.80 kWh. The recalculated interval residuals are both 0. The corrected cumulative sequence is non-decreasing, the corrected interval increment is not less than zero, and the power generation does not exceed the product of the rated AC output power and the anchoring interval duration, thus forming an automatic correction result.

[0173] After the projection is closed, the power generation from 09:00 to 11:00 is 14.00 kWh, with 8.40 kWh fed into the grid and 2.20 kWh removed from the grid. Before and after 10:30, the identifiable grid-connected power is 5.80 kWh and 2.60 kWh, respectively, and the identifiable grid-removed power is 1.60 kWh and 0.60 kWh, respectively. The grid-connected electricity prices for the two pricing versions are RMB 0.38 / kWh and RMB 0.42 / kWh, respectively; the electricity purchase prices are RMB 0.56 / kWh and RMB 0.61 / kWh, respectively; and the settlement price for grid connection is RMB 0.12 / kWh for both.

[0174] The settlement amount for electricity sales is 5.80 × 0.38 + 2.60 × 0.42, which is 3.296 yuan; the settlement amount for electricity purchase is 1.60 × 0.56 + 0.60 × 0.61, which is 1.262 yuan. The electricity consumed within the boundary is 14.00 − 8.40, which is 5.60 kWh, and the settlement amount for consumption is 5.60 × 0.12, which is 0.672 yuan. Calculated according to the income minus expenditure direction, the settlement amount is 3.296 − 1.262 + 0.672, which is 2.706 yuan. The micro-element electricity volume of each settlement micro-element is only included in its respective price version.

[0175] To verify the identifiable price range determination and transaction settlement effectiveness, 60 complete metering periods were replayed for verification. Of these, 20 periods were injected with an intermediate cumulative anchor point deviation of 0.30kWh to 1.00kWh, 20 periods were injected with a time offset of 5 to 15 minutes, and 20 periods were set to have continuous time-sharing data missing. The metering results before the abnormal injection were used as the control value. The comparison method used cumulative meter code difference, and the electricity volume across price versions was allocated according to the settlement micro-element duration ratio.

[0176] In the 20 intermediate cumulative anchor point deviation periods, this method identifies 19 periods as single anchor point value anomalies, and the remaining period is included in the pending confirmation because the correction results do not meet all correction constraints; the comparison method records the deviations before and after the intermediate cumulative anchor points as interval anomalies respectively.

[0177] In 40 test periods involving price range boundaries, the median absolute error of the electricity consumption in the price range using this method was 0.02 kWh, with the first and third quartiles being 0.00 kWh and 0.05 kWh, respectively; while the median absolute error of the electricity consumption in the price range using the comparative method was 0.31 kWh, with the first and third quartiles being 0.18 kWh and 0.47 kWh, respectively.

[0178] The median absolute error of the settlement amount using this method is 0.01 yuan, with the first and third quartiles being 0.00 yuan and 0.02 yuan, respectively. The median absolute error of the settlement amount using the comparative method is 0.14 yuan, with the first and third quartiles being 0.08 yuan and 0.22 yuan, respectively. When time-of-use data is missing and the electricity consumption in a price segment does not meet the identifiable criteria, this method retains the cumulative total and sets the relevant settlement micro-element as pending settlement.

[0179] See appendix Figure 8 , Figure 8 In the box plot, the lower and upper boundaries of the box represent the first and third quartiles, respectively. The horizontal line inside the box represents the median, and the whiskers on both sides of the box represent the range of non-outlier data. The markers outside the whiskers represent outliers. Circular markers correspond to the median and outliers of this method, while square markers correspond to the median and outliers of the comparison method.

[0180] See appendix Figure 9 , Figure 9 The meaning of the box body, the horizontal lines inside the box body, the struts, and the markings outside the struts in the box molding. Figure 8 Consistent; circular markers correspond to the median and outliers of this method, while square markers correspond to the median and outliers of the comparison method.

[0181] Figure 8 and Figure 9 This indicates that the absolute error of electricity consumption in price ranges and the absolute error of settlement amount in the test data using this method are both lower than those of the comparative method.

[0182] While specific embodiments of the present invention have been described above, those skilled in the art should understand that these specific embodiments are merely illustrative. Those skilled in the art can omit, substitute, and modify the details of the above methods and systems in various ways without departing from the principles and essence of the present invention. For example, combining the above method steps to perform substantially the same function and achieve substantially the same result using substantially the same method falls within the scope of the present invention.

Claims

1. A method for statistical analysis and transaction settlement of household distributed photovoltaic power generation, characterized in that, Includes the following steps: Obtain the raw cumulative electricity readings, normalize and segment them into events to obtain normalized metering data; Authoritative cumulative anchor points are extracted from the standardized metering data to form anchor intervals and calculate authoritative power increments. Settlement micro-elements are generated. The authoritative cumulative anchor points are used as cumulative anchor points to be verified to verify independent time-sharing sources not derived from the cumulative anchor points to be verified, thereby obtaining reference power increments and available identifiers for independent sources. By combining the authoritative power increment, the reference power increment, and the independent source available identifier to calculate the interval residual, a reverse method is used to identify reverse residual pairs with opposite directions and closed sums on the anchored interval, forming a multi-channel residual signature and determining the anomaly type. When the anomaly type is an anchor value anomaly, the anchor correction amount is calculated to obtain the automatic correction result. When the anomaly type is a time alignment anomaly, a time mapping version is generated. Based on the authoritative power increment, the settlement micro-element, and the automatic correction result or time mapping version corresponding to the anomaly type, a set of constrained feasible solutions is constructed, the micro-element power is determined, and a distinguishable marker is identified. When the distinguishable marker is distinguishable, a temporary settlement record is generated. The settlement impact domain is determined based on the anomaly type and the preset propagation rules, and a partial re-settlement is performed. The temporary settlement records are then checked against the power grid single check to generate the final settlement record and settlement archive.

2. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 1, characterized in that, The steps to obtain normalized metering data by acquiring raw cumulative electricity readings, performing normalization and event segmentation include: Register the metering topology and obtain the metering resolution, collect raw data and obtain the raw cumulative electricity reading and the power grid settlement statement, perform normalization conversion on the raw cumulative electricity reading, divide the cumulative segment by event and calculate the effective increment to obtain the normalized metering data.

3. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 1, characterized in that, The steps to obtain the reference power increment and the available identifier of an independent source include: The authoritative cumulative anchor point is extracted from the standardized metering data to form the anchoring interval and obtain the authoritative electricity increment. The time boundary is then merged to generate the settlement micro-element. When the independent time-sharing source is not derived from the cumulative anchor point to be verified, calculate the reference power increment and the reference micro-element power, and mark the independent source as available. Otherwise, the independent source is marked as unavailable.

4. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 1, characterized in that, The steps for generating multi-channel residual signatures and determining anomaly types include: When the independent source is marked as available, the interval residual is obtained by subtracting the reference power increment from the authoritative power increment, and the interval allowable deviation is determined. When the absolute value of the interval residual is not greater than the allowable deviation of the interval, the anchoring interval is marked as a closable interval; The authoritative electricity increments of the power generation channel, grid connection channel, and grid disconnection channel are mapped to the power generation increment, grid connection increment, and grid disconnection increment, respectively. The authoritative cumulative anchor point shared by adjacent anchoring intervals is used as the intermediate cumulative anchor point. Using the reverse method, when the independent source available identifiers of the two intervals are both available, the interval residual directions are opposite, the absolute values ​​of each are greater than the allowable deviations of their respective intervals, and the absolute value of the sum of the two is not greater than the sum of the allowable deviations of the two intervals, the reverse residual pair is formed, the candidate abnormal anchor point is determined, and the reverse residual direction status is recorded. When the absolute values ​​of the residuals of both intervals are greater than their respective interval allowable deviations but in the same direction, or when the absolute value of the sum of the residuals of both intervals is greater than the sum of the allowable deviations of the two intervals, a time mapping candidate is generated and a trial calculation is performed. When the absolute value of the interval residual of multiple consecutive anchoring intervals is greater than the allowable deviation of each interval, they are merged into a group of consecutive abnormal anchor points. The physical feasibility status is determined based on the grid connection increment, power generation increment, and grid disconnection increment: when the grid connection increment is not greater than the power generation increment, the power generation increment is not greater than the product of the rated AC output power and the anchoring interval duration, and the grid disconnection increment is not less than zero, the physical feasibility status is physical feasibility; otherwise, it is physical infeasibility. The independent source availability identifiers, reverse residual direction status, and physical feasibility status of the power generation channel, grid connection channel, and grid disconnection channel are summarized to form the multi-channel residual signature; When the physically feasible state is physically infeasible and the reverse formula is not true, the exception type is physical relationship exception; When there is a unique time mapping candidate such that the absolute value of the interval residual of each channel after mapping is not greater than its respective interval allowable deviation, the anomaly type is time alignment anomaly. When the aforementioned reverse residual pair is formed and it does not belong to the time alignment anomaly, the anomaly type is anchor value anomaly; When the independent source is marked as unavailable, the exception type is time-sharing data missing; When the absolute value of the interval residual is greater than the allowable deviation of the interval and does not belong to the above-mentioned anomaly types, the anomaly type is unlocated anomaly.

5. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 4, characterized in that, The steps to obtain automatic correction results include: Authoritative cumulative anchor points whose previous and subsequent anchor intervals are both closable intervals and are not candidate abnormal anchor points are determined as reliable anchor points. When there is only a single candidate abnormal anchor point, the anchor point correction amount of the intermediate cumulative anchor point is calculated using a correction formula to form the automatic correction result and write it into the correction view. When a continuous abnormal anchor point group is formed, the anchor point correction amount of each intermediate cumulative anchor point is calculated using the two reliable anchor points as fixed boundaries, and the automatic correction result is generated and written into the correction view.

6. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 4, characterized in that, The steps to generate a time-mapped version include: Generate time mapping candidates and recalculate the interval residuals, then select the optimal mapping to obtain the time mapping version; When the optimal mapping is not unique and its impact on the electricity consumption of the price range is greater than the metering resolution, a non-unique mapping state is formed and marked as pending confirmation.

7. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 1, characterized in that, The steps for constructing a set of constrained feasible solutions and determining the infinitesimal quantity include: Construct a set of constrained feasible solutions, which must satisfy preset cumulative total constraints, power generation constraints, directional constraints, or energy storage metering constraints when the directional constraints are deactivated. When the independent source is marked as available, the reference micro-element power is obtained from the independent time-sharing source; When the set of feasible solutions to the constraints is not empty, the projective method is used to solve for the infinitesimal charge of each channel, so that the weighted squared difference between the infinitesimal charge of each channel and the reference infinitesimal charge is minimized.

8. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 7, characterized in that, The steps for generating a temporary settlement record when the identifiable marker is identifiable include: The price range is obtained from the settlement micro-element, and the lower bound and upper bound of the segment electricity are calculated using a distinguishable formula; When the difference between the upper limit of the segment electricity and the lower limit of the segment electricity is not greater than the metering resolution, the identifiable mark is identified, the identifiable electricity is obtained, and the temporary settlement record is generated. When the difference between the upper limit of the electricity consumption of the segment and the lower limit of the electricity consumption of the segment is greater than the metering resolution, the settlement micro-element to which the price segment belongs is set to pending settlement, and a pending settlement record is generated.

9. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 1, characterized in that, The steps of determining the settlement impact domain and performing local resettlement based on the anomaly type and preset propagation rules include: The settlement impact domain is constructed according to the propagation rules. The propagation rules are as follows: changes in power generation channels affect the accounts of power generation and consumption within the boundary; changes in grid connection channels affect the accounts of grid connection and consumption within the boundary; changes in off-grid channels affect the accounts of off-grid; changes in time mapping affect the union of settlement micro-elements involved in the old and new time mapping versions; and changes in price versions affect the accounts of monetary amounts while the basic metered electricity remains unchanged. When correcting a single anchor point, only the two anchoring intervals before and after are processed. For a group of consecutive abnormal anchor points, the reliable anchor points on both sides are used as the boundary. When the amount of issued accounts decreases, a refund is executed; when the amount increases, a supplementary collection is executed; when the entire account is replaced, a reversal is executed, generating a replacement record and obtaining the resettled account.

10. The method for statistical analysis and transaction settlement of household distributed photovoltaic power generation according to claim 9, characterized in that, The steps for generating the final settlement record and settlement archive include: Obtain the power grid settlement statement from the raw data, and then reconcile the power grid settlement statement with the local verification quantity. When a corrected total exists, the local verification quantity is the corrected total; otherwise, the authoritative power increment is used. When the verification results are consistent, the final settlement record is generated from the temporary settlement record, and the settlement archive is generated by combining the resettled accounts. When the verification results are inconsistent, an unsettled discrepancy form is generated. Based on the resettled items, the accounts are resettled locally according to the settlement impact domain and then verified again. When the verification results are consistent again, the final settlement record and the settlement archive form are generated.

Images