SYSTEM, METHOD, AND INTERFACE FOR GOAL-ALLOCATION OF RESOURCES AND DYNAMIC MONITORING OF PROGRESS

Abstract

System, method, and interface for visualized resource allocation and algorithms for the reallocation of resources to achieve a goal. The system analyses an initial state of resource allocation, a cost function for undesirable resources, and a set of potential incremental improvements, each with an associated cost, and determines a step-wise path of applying the incremental improvements to achieve an ultimate resource-allocation goal in an economically feasible way. Simultaneously, a user interface depicts the state of the allocation at the beginning, at the end, and along the path, allowing an intuitive understanding of how the goal will be achieved.

Claims

1. One or more non-transitory computer-readable storage media storing computer-executable instructions that, when executed by at least one processor, perform a method of resource allocation of resources to achieve a goal and dynamic monitoring of progress, the method comprising: receiving source data indicative of flow types; receiving transformation data indicative of flow transformations, wherein the flow types and the flow transformations are classified by region; generating, by an interface, visualizations of the flow types based on the source data and the flow transformations based on the transformation data; monitoring the source data and the transformation data; and adjusting resource allocations to maintain an optimal allocation of the resources.

2. The media of claim 1, wherein the flow types are input flow, output flow, total flow, wasted flow, and remaining flow.

3. The media of claim 2, wherein the flow transformations are indicative of changes in public and private social infrastructure and transportation.

4. The media of claim 1, wherein the source data is indicative of social contributions to greenhouse gases, and wherein the resource allocations are adjusted to reduce greenhouse effects.

5. The media of claim 4, wherein the method further comprises adjusting coal, petroleum, and natural gas usage by residential, business, and transportation.

6. The media of claim 5, wherein the method further comprises: determining the region for the source data; calculating a greenhouse cost for the social contributions to the greenhouse gases; and calculating a resource cost for the resource allocation to reduce the greenhouse effects.

7. The media of claim 6, wherein the method further comprises displaying, for any region or data source within the region, a processing efficiency, emissions data, and region costs.

8. The media of claim 1, wherein the flow types comprise a benign inflow from a benign source and a harmful inflow from a harmful source, and wherein the method further comprises increasing the benign inflow by increments and incrementally reducing the harmful inflow.

9. The media of claim 8, wherein the method further comprises: tracking a cost of each of the benign inflow and the harmful inflow; and maximizing the benign inflow with a cost constraint.

10. The media of claim 1, wherein the method further comprises: monitoring the region and determining resource consumption patterns across the region; displaying, by the interface and across at least a present and a future timeline, the resource consumption patterns and a cost associated with the resource consumption patterns; receiving, by the interface, an input changing one or more parameters along the future timeline; and displaying, by the interface, changes to the cost and the resource allocation based on the one or more parameters.

11. The media of claim 10, wherein the method further comprises: evaluating the cost, a savings, and a potential cost across at least the present and the future timeline, wherein the potential cost is a difference between a maximum cost and the cost based on the one or more parameters; and further displaying the cost and the one or more parameters in a real-time cube view or a linear temporal view.

12. The media of claim 1, wherein the method further comprises: receiving the source data from a network device or a data acquisition machine, wherein the source data is monitored continuously or periodically based on an expected quantity of the flow types.

13. A method of resource allocation of resources to achieve a goal and dynamic monitoring of progress, the method comprising: receiving source data indicative of flow types; receiving transformation data indicative of flow transformations; generating, by an interface, visualizations of the flow types based on the source data and the flow transformations based on the transformation data; continuously monitoring the source data and the transformation data; and adjusting resource allocations to maintain an optimal allocation of the resources.

14. The method of claim 13, further comprising: determining a region for the source data; calculating a greenhouse cost for infrastructure contributions to a greenhouse gases in the region; calculating a resource cost for the resource allocations to reduce greenhouse effects; displaying, for any region or data source within the region, emissions data and region costs; and adjusting the resource allocations to reduce the greenhouse effects.

15. The method of claim 13, wherein the flow types are input flow, output flow, total flow, wasted flow, and remaining flow, and wherein the flow transformations are indicative of changes in public and private social infrastructure and transportation.

16. The method of claim 13, further comprising: receiving the source data from a network device or a data acquisition machine, wherein the source data is monitored continuously or periodically based on an expected quantity of the flow types.

17. One or more non-transitory computer-readable storage media storing computer-executable instructions that, when executed by at least one processor, perform a method of resource allocation of resources to achieve a goal and dynamic monitoring of progress, the method comprising: receiving source data indicative of flow types; receiving transformation data indicative of flow transformations; generating, by an interface, visualizations of the flow types based on the source data and the flow transformations based on the transformation data; continuously monitoring the source data and the transformation data; and adjusting resource allocations to maintain an optimal allocation of the resources.

18. The media of claim 17, wherein the method further comprises: receiving, by the interface, an input changing one or more parameters along a timeline; evaluating a cost, a savings, and a potential cost across the timeline; wherein the potential cost is a difference between a maximum cost and the cost based on the one or more parameters; and displaying, by the interface, the cost and the resource allocation based on the one or more parameters.

19. The media of claim 17, wherein the method further comprises: determining a region for the source data; calculating a greenhouse cost for estimated greenhouse gases emitted by infrastructure in the region; calculating a resource cost for the resource allocations to reduce greenhouse effects; and displaying, for any region or data source within the region, emissions data, and region costs associated with the resource allocations.

20. The media of claim 17, wherein the method further comprises: receiving the source data from a network device or a data acquisition machine, wherein the source data is monitored continuously or periodically based on an expected quantity of the flow types, wherein the flow types are input flow, output flow, total flow, wasted flow, and remaining flow, and wherein the flow transformations are indicative of changes in public and private social infrastructure and transportation.

Description

BRIEF DESCRIPTION OF THE DRAWING FIGURES

[0008] Embodiments of the invention are described in detail below with reference to the attached drawing figures, wherein:

[0009] FIG. 1 depicts an exemplary hardware platform for certain embodiments of the invention;

[0010] FIG. 2 depicts an exemplary display of a user interface presenting a three-dimensional data visualization in accordance with embodiments of the invention;

[0011] FIGS. 3A-3E depict a variety of additional views of the user interface of FIG. 2 in accordance with embodiments of the invention;

[0012] FIGS. 4-7 depict embodiments of relational queries;

[0013] FIG. 8 depicts one sector-type category of the user interface of FIG. 2 decomposed into inflows and transformations;

[0014] FIG. 9 depicts the same sector-type category as FIG. 8 of user interface of FIG. 2 decomposed instead into wasted energy and utilized energy;

[0015] FIG. 10 depicts a starting state of resource allocation in a system;

[0016] FIG. 11 depicts a final state of resource allocation in the system after reallocation;

[0017] FIGS. 12A-12B depict a visual illustration of the exemplary resource allocation process of FIGS. 10 and 11;

[0018] FIGS. 13A-13B depict the changes in the solution entities for the exemplary resource allocation process of FIGS. 10 and 11;

[0019] FIGS. 14A-14B depict the changes in the transformation waste ratios for the exemplary resource allocation process of FIGS. 10 and 11;

[0020] FIGS. 15A-15B depict visual illustrations of process towards meeting the goals for the exemplary resource allocation process of FIGS. 10 and 11;

[0021] FIG. 16 depicts a direct data entry method for accessing the data associated with the cube;

[0022] FIGS. 17-25 depict exemplary user interfaces for inputting data and displaying outputs of the system;

[0023] FIG. 26 depicts an exemplary topographical map for presenting data;

[0024] FIG. 27 depicts a hybrid visual/numeric user interface for a user's visualization of the goals and progress of the resource allocation process;

[0025] FIG. 28 depicts exemplary cubes and cubelets for visualizing energy source consumption;

[0026] FIG. 29 depicts an exemplary user interface providing selections to the user;

[0027] FIG. 30 depicts how an individual user's interface can be automatically updated to reflect new goals set by an authorized user;

[0028] FIG. 31 depicts an exemplary optimization flow; and

[0029] FIG. 32 depicts an exemplary flow diagram for resource allocation optimization and visualization.

[0030] The drawing figures do not limit the invention to the specific embodiments disclosed and described herein. The drawings are not necessarily to scale, emphasis instead being placed upon clearly illustrating the principles of the invention.

DETAILED DESCRIPTION

[0031] At a high level, embodiments of the invention provide for an interface for visualize resource allocation and algorithms for the reallocation of resources to achieve a goal. By visualizing the starting state, the ending state, and the transition between them, stakeholders can more easily grasp the path of a transition to a less ecologically harmful energy source balance. At the same time, the algorithm determines not merely the best final allocation of those energy sources, but an incremental path of the transition that is technologically and economically feasible.

[0032] The subject matter of embodiments of the invention is described in detail below to meet statutory requirements; however, the description itself is not intended to limit the scope of claims. Rather, the claimed subject matter might be embodied in other ways to include different steps or combinations of steps similar to the ones described in this document, in conjunction with other present or future technologies. Minor variations from the description below will be obvious to one skilled in the art, and are intended to be captured within the scope of the claimed invention. Terms should not be interpreted as implying any particular ordering of various steps described unless the order of individual steps is explicitly described.

[0033] The following detailed description of embodiments of the invention references the accompanying drawings that illustrate specific embodiments in which the invention can be practiced. The embodiments are intended to describe aspects of the invention in sufficient detail to enable those skilled in the art to practice the invention. Other embodiments can be utilized and changes can be made without departing from the scope of the invention. The following detailed description is, therefore, not to be taken in a limiting sense. The scope of embodiments of the invention is defined only by the appended claims, along with the full scope of equivalents to which such claims are entitled.

[0034] In this description, references to “one embodiment,” “an embodiment,” or “embodiments” mean that the feature or features being referred to are included in at least one embodiment of the technology. Separate reference to “one embodiment” “an embodiment”, or “embodiments” in this description do not necessarily refer to the same embodiment and are also not mutually exclusive unless so stated and/or except as will be readily apparent to those skilled in the art from the description. For example, a feature, structure, or act described in one embodiment may also be included in other embodiments but is not necessarily included. Thus, the technology can include a variety of combinations and/or integrations of the embodiments described herein.

[0035] Turning first to FIG. 1, an exemplary hardware platform for certain embodiments of the invention is depicted. Computer 102 can be a desktop computer, a laptop computer, a server computer, a mobile device such as a smartphone or tablet, or any other form factor of general- or special-purpose computing device. Depicted with computer 102 are several components, for illustrative purposes. In some embodiments, certain components may be arranged differently or absent. Additional components may also be present. Included in computer 102 is system bus 104, whereby other components of computer 102 can communicate with each other. In certain embodiments, there may be multiple busses or components may communicate with each other directly. Connected to system bus 104 is central processing unit (CPU) 106. Also attached to system bus 104 are one or more random-access memory (RAM) modules 108. Also attached to system bus 104 is graphics card 110. In some embodiments, graphics card 110 may not be a physically separate card, but rather may be integrated into the motherboard or the CPU 106. In some embodiments, graphics card 110 has a separate graphics-processing unit (GPU) 112, which can be used for graphics processing or for general purpose computing (GPGPU). Also on graphics card 110 is a processor 112 and GPU memory 114. Connected (directly or indirectly) to graphics card 110 is display 116 for user interaction. In some embodiments no display is present, while in others it is integrated into computer 102. Similarly, peripherals such as keyboard 118 and mouse 120 are connected to system bus 104. Like display 116, these peripherals may be integrated into computer 102 or absent. Also connected to system bus 104 is local storage 122, which may be any form of computer-readable media and may be internally installed in computer 102 or externally and removeably attached.

[0036] Computer-readable media include both volatile and nonvolatile media, removable and nonremovable media, and contemplate media readable by a database. For example, computer-readable media include (but are not limited to) RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile discs (DVD), holographic media or other optical disc storage, magnetic cassettes, magnetic tape, magnetic disk storage, and other magnetic storage devices. These technologies can store data temporarily or permanently. However, unless explicitly specified otherwise, the term “computer-readable media” should not be construed to include physical, but transitory, forms of signal transmission such as radio broadcasts, electrical signals through a wire, or light pulses through a fiber-optic cable. Examples of stored information include computer-useable instructions, data structures, program modules, and other data representations.

[0037] Finally, network interface card (NIC) 124 is also attached to system bus 104 and allows computer 102 to communicate over a network such as network 126. NIC 124 can be any form of network interface known in the art, such as Ethernet, ATM, fiber, Bluetooth, or Wi-Fi (i.e., the IEEE 802.11 family of standards). NIC 124 connects computer 102 to local network 126, which may also include one or more other computers, such as computer 128, and network storage, such as data store 130. Generally, a data store such as data store 130 may be any repository from which information can be stored and retrieved as needed. Examples of data stores include relational or object-oriented databases, spreadsheets, file systems, flat files, directory services such as LDAP and Active Directory, or email storage systems. A data store may be accessible via a complex API (such as, for example, Structured Query Language), a simple API providing only read, write and seek operations, or any level of complexity in between. Some data stores may additionally provide management functions for data sets stored therein such as backup or versioning. Data stores can be local to a single computer such as computer 128, accessible on a local network such as local network 126, or remotely accessible over Internet 132. Local network 126 is in turn connected to Internet 132, which connects many networks such as local network 126, remote network 134 or directly attached computers such as computer 136. In some embodiments, computer 102 can itself be directly connected to Internet 132.

[0038] Also depicted in FIG. 1 are a variety of data sources 138 used for classification, assignment of energy sources to an entity and/or region, timescale, etc. by applying appropriate unit conversions (e.g., between Watts, BTU, joules, kWh, volts, gallons, tons, ktoe (kilotonne of oil equivalent), therms, cubic feet, parts per million). Data may be stored in a database accessible over a network or locally. Real-time and/or live data can also be continuously received from a network device or data acquisition machine; similarly, historical or archival data can be accessed from a database, input via document ingestion, or directly entered by a user or another entity on user's behalf. In some exemplary embodiments, the data sources 138 may be energy sources generating energy such as through nuclear, coal, wind, and solar power. The data sources 138 may also be energy consumption through public and private social infrastructure and transportation. In some embodiments, the data sources 138 may be analyzed and displayed along with the analysis as described in embodiments below.

[0039] Turning now to FIG. 2, an exemplary display of a user interface in accordance with embodiments of the invention is depicted and referred to generally by reference numeral 200. As shown, the user interface 200 depicts a multidimensional data visualization 202 ideal for analyzing breakdowns of a whole from the data sources 138 along multiple axes simultaneously. As shown, for example, the whole of a city's greenhouse-gas-generating energy consumption is shown, broken down on a first axis, let's say J, by a consuming sector 204 of a carbon-based fuel (for example, “residential,” “business,” and “transportation”) and on a second axis, let's say K, by the particular type of carbon-based fuel 206 (for example, “coal,” “petroleum,” and “natural gas”). Thus, for example, it can readily be seen that little or no coal goes towards powering transportation. Instead, by far the largest share of carbon-based fuel used in transportation is petroleum-based. In some embodiments, the amount of greenhouse gasses generated by each sector-type category can be depicted instead of the amount of energy consumed.

[0040] Also shown in the user interface 200 is the amount of energy consumed 208 in each sector-type category that is wasted. The source of the wastage may differ in different categories. The proportional amount may be visualized by the height of the sector-type category along the third axis. For example, the wastage in the transportation-petroleum category may include energy consumed in transporting crude oil to refineries, energy consumed in refining the crude oil into gasoline and diesel fuel, and energy consumed in transporting the fuels to fueling stations. The transportation-petroleum category has relatively low impact from wind, biofuel, solar, geothermal, and nuclear energies. Evaluating the data sources 138 in this way, when determining how best to reduce greenhouse gas emissions, the indirect contributions of each sector-type can be taken into account as well as the direct contributions.

[0041] In some embodiments, the first and second axis of the visualization display the data source, or energy sources, coal, petroleum, and natural gas and the energy transitions or consumption processes: homes, business, and transportation, as displayed. In some embodiments the energy sources may be wind, solar, water, or any other method of generating energy and the consumption may be more detailed, or sub-categories, such as airplanes, automobiles, and boats. Further, the sub-categories may be provided in drill-down visualization methods described below.

[0042] Turning now to FIG. 3A, a second view of the user interface is depicted and referred to generally by reference numeral 300. FIG. 3A depicts a drill-down view 302 of a particular sector-type category. For example, the transportation-petroleum category might be further broken down into subcategories, X, along the transportation axis into “ships,” “cars,” “motorcycles,” “trains,” “planes,” “trucking,” and so forth, while the petroleum axis might be broken down into, for example, “residual oil,” “motor gasoline,” “jet fuel,” “diesel,” and so forth as Y.sub.k. Breaking down the axis into further sub-categories, or subsets, allows a user to visualized how each individual resource is consumed on the same data display.

[0043] In some embodiments, the breakdown of a particular axis is constant among the various categories and sub-categories of the other axis. Thus, for example, a “home heating oil” category would appear as a sub-category for petroleum even in the intersection with “transportation” where home heating oil is used little or not at all. In other embodiments, drill-down view 302 filters out inapplicable categories such that a “jet fuel” category is present in the drill-down view of the intersection of “transportation” and “petroleum” but not in the drill-down view of the intersection of “home” and “petroleum.” Similarly, in such an environment, “home heating oil” would be a subcategory of petroleum in the drill-down view of its intersection with “home” but not in the drill-down view of its intersection with “transportation.” Drill-downs are possible for a particular sector-type category, across an entire sector or category, or within a particular layer. This can be represented using relational algebra syntax as described below.

[0044] After classification of resource type Y.sub.k, time, t, usage location, p and entity S, the visual can be decomposed into more granular representations according to some standardized attribute (i.e. after conversion to the same or comparable physical units), such as distributions in terms of quantity, amount or density, as depicted in FIG. 3B. The standardization allows for relative comparison such that a user may easily visualize the usage and waste of the resources for each resource/energy consumption process. FIG. 3B depicts the cube visualization on the user interface 200 from above along with a drilled-down visualization 304. As shown, the top-view of the bottom layers of each individual layer can be converted to a two-dimensional x-y grid, where the values are laid across a spectrum 306. As shown in FIGS. 3B-E, the spectrum 306 may be provide texture. It should also be contemplated that the spectrum 306 may be a color spectrum, such as with lighter colors signifying lower values, and darker colors as higher values. The representation may also present textures, shading, lines, dots, or any other method of displaying various sections. For example, shades of color may be equivalent to height in the z-direction in the three-dimensional representation. The spectrum 306 may present a two-dimensional data visualization that may be easier to understand when many inputs are compared. Thus, the exemplary spectrum 306 provides a quick and easily understandable magnitude to the visualizations for a particular cross section of the graph.

[0045] Given three tables T.sub.1, T.sub.2, T.sub.3 corresponding to “residential home” (J=1), “commercial and industrial business” (J=2) and transportation (J=3) sectors, with each table containing information about consumption (in a given sector) of different types of energy source Y.sub.k (synonymous to A.sub.i . . . N and B.sub.i . . . M) distributed over geographic regions, p (for countries, states, cities, zipcodes) and time. t (for years, months, days, seconds), relational algebra operations can be used to represent and various input information from T.sub.1, T.sub.2, T.sub.3 for use in interface visualization, and/or optimization algorithm introduced later.

[0046] The system can aggregate J=1,2,3 sheets to make a treemap (302) where only environmentally harmful energy sources are shown. By eliminating unnecessary data columns using the operation of projection, filtering for this data can be performed:

Π.sub.K,Y,J=1(T.sub.1)

Π.sub.K,Y,J=2(T.sub.2)

Π.sub.K,Y,J=3(T.sub.3)   (1)

where T.sub.j is the table corresponding to the sheet J=1. Now, renaming the column “J=j” to “J” for all j=1,2,3:

ρ.sub.J/J=1(T.sub.1)

ρ.sub.J/J=2(T.sub.2)

ρ.sub.J/J=3(T.sub.3)   (2)

Finally, a new table is created with columns “K”, “Y”, “J”, etc. and whose rows are obtained by summing “J” values of corresponding rows of T.sub.1, T.sub.2, T.sub.3 (note that rows of T.sub.1, T.sub.2, T.sub.3 differ only in the “J”th column):

T=Σ.sub.J(T.sub.1, T.sub.2, T.sub.3).   (3)

Here Σ.sub.J(T.sub.1, T.sub.2, T.sub.3) denotes the row-wise summation table of T.sub.1, T.sub.2, T.sub.3.

[0047] The system can aggregate J=1,2,3 for a particular geographic region, p or a time interval, t. For this, T.sub.1, T.sub.2, T.sub.3 depend on an artifact R (considered earlier as dots on map (e.g., point diameter 316), or shape boundary (e.g., polygons 314) and time t, where t is given in one of the formats: “year”, “year”/“month”, “year”/“month”/“day” or “year”/“month”/“day”/“HH:MM:SS”. Therefore, the above mentioned aggregation may be performed at the level of each geographic region represented by the topographical map 308 and for every moment/interval of time.

[0048] The above-described aggregation may be performed for p=NYC and on the day t=2019/03/21. Steps (1) and (2) may be repeated for tables T.sub.1(R, t), T.sub.2(R, t), T.sub.3(R, t), where R runs through all artifacts whose geographic coordinates g is(x,y) are in p=NYC. Finally, the following summation operation is performed:

T(R,t):=Σ.sub.J(T.sub.1(R,t), T.sub.2(R, t), T.sub.3(R,t)),   (4)

and then one more summation over all artifacts R that are in p=NYC:

T(p, t):=Σ.sub.J({T(R, t): R is such that g is (R) ∈ p}).   (5)

[0049] In the last formula the summation operation is applied to tables T(R, t) along the column J. The resulting tables may be aggregated for p=NYC over a time interval, say from t.sub.1=2019/01/01 to t.sub.2=2019/06/30. The system then sums tables (5) (for appropriate t's, where t runs through all days between t.sub.1 and t.sub.2.):

T(p, [t.sub.2]):=Σ.sub.J({T(p, t): t.sub.1≤t≤t.sub.2}).   (6)

[0050] The system can search by extent or text identifier of a region, such as any given p (say New York City). In order to do this aggregation just at the level of one city or region p, the earlier steps can be repeated for T.sub.1(R), T.sub.2(R), T.sub.3(R) for all artifacts R whose coordinates g is(x, y) belong to p across different energy types, as shown in the drill down visualization 304. Here T.sub.j(R) is the table of the artifact R corresponding to the sector j=1,2,3; one can obtain it by summing T.sub.j(R, t) which was introduced earlier, over −∞<t<+∞ where t can vary from the first to the last time record.

T.sub.1(R)=Σ.sub.J=1, 1, . . . , 9({T.sub.1(R, t): −∞<t<∞}),

T.sub.2(R)=Σ.sub.J=2,1, . . . , 8,17({T.sub.2(R, t): −∞<t<∞}),

T.sub.3(R)=Σ.sub.J=3,1, . . . , 8({T.sub.3(R, t): −∞<t<∞}),   (7)

where Σ.sub.i.sub.1.sub., . . . , i.sub.d({T.sub.j(R, t): −∞<t<∞}) denotes the table obtained from tables T.sub.j(R, t) by summing them column-wise along columns i.sub.1, . . . , i.sub.d.

[0051] Applying steps (1) and (2) to T.sub.1(R), T.sub.2(R), T.sub.3(R) the system obtains a table T(R):=Σ.sub.J(T.sub.1(R), T.sub.2(R), T.sub.3(R)), which is the analogue of (3) for an individual artifact R, which can be vehicles, factories, power plants, buildings etc. Finally, performing the summation operation over all artifacts R whose coordinates g is (x,y) are in p:

T(p):=Σ.sub.J({T(R): R is such that g is(R) ∈ p})   (8)

[0052] To filter and conduct an aggregation from 1 for a given t (say t=2007), the system takes tables T.sub.J(R, t), j=1,2,3, and aggregates the result over all R, and not confined to just one p:

T(t):=Σ.sub.J({T(R, t): R}).   (9)

To aggregate for all regions p, for a time internal, such as from t.sub.1=2010 to t.sub.2=2020, it will be done similar to [t.sub.1, t.sub.2] and for a given city/region p, aggregated over all p:

T([t.sub.1, t.sub.2]):=Σ.sub.J({T(t): t.sub.1≤t≤t.sub.2}).   (10)

Within these, it is also possible filter only for a given sub-category, X=4 (cooking) and fuel type, Y=21 (K=B=2; kerosene) at J=2, but from t.sub.1=2010 to t.sub.2=2015. The row Y=21 (kerosene) of tables T.sub.2 (R, t) allows extraction of only this row using selection operation:

σ.sub.Y=21(T.sub.2(R, t)).   (11)

[0053] The table T.sub.2(R, t) has only one row (and same columns as before) and given the specific interest in the column X=4 (cooking) the system removes other columns via projection operation:

Π.sub.K,Y,X=4(T.sub.2(R, t)).   (12)

The only place where tables T.sub.2(R, t) differ (for different R and t) is the value in the column X=4 (in the unique row Y=21). To make aggregation over R:

T.sub.2(t):=Σ.sub.X=4({T.sub.2(R, t): R}),   (13)

where summation is along column X=4 (in (9) this was column J). To extract information over a time period [t.sub.1, t.sub.2], the system performs the following operation:

T.sub.2([t.sub.1, t.sub.2]):=Σ.sub.X=4({T.sub.2(t): t.sub.1≤t≤t.sub.2}),   (14)

where the summation is along column X=4.

[0054] As shown in FIG. 3C, the arbitrary positioning of the grid of FIG. 3B could be overlaid on more understandable geospatial grids, such as on topographical maps 308, or neighborhood blocks 310. Users can filter by type, such as flow types (input, output, total, wasted, or remaining), transformation flows (input, output, wasted, or consumed); or per-user flows of any type. For example, the country can be broken down into states, cities, counties, or even further. Each of these locations may be analyzed based on resource consumption in each category as described above. Using notation from relational algebra, we can consider tables which store these values in and conducting operations to extract necessary information.

[0055] Any identified entity and its data attributes can be further explained or visualized in multiple different ways. For example, the visualization may be provided in the grid distributions as described above or either as polygons of a predetermined shape, or as the diameter of a sphere, circle, or other shape which is centered at the centroid of a uniquely identified location or entity, as shown in FIG. 3D. The grid distributions, or spectrum 306 may be shown with a relative low-high indicator 312. In some embodiments, the low-high indicator may be colored to show the difference from low to high. Further, the grid may be broken down into more defined polygons 314. Further still, the relative magnitude of each region may be represented by a relative size of a shape, for example the size of the circles in the point diameter map 316. Any shapes and combination of shapes, textures, colors, and maps may be contemplated for visualization of data.

[0056] As depicted in FIG. 3E, embodiments of the invention may store polylines, features, feature-sets and feature-collections for each item, data point or element, searchable by extents. The user interface 200 in FIG. 3E depicts the drilled-down visualization 304 along with the topographical maps 308 and the more reduced geographical data set down to neighborhood blocks 310. This could be stored in multiple framework, datatypes, vectors, and languages such as JSON etc. This enables the user and other entities, such as regulators or solution providers, to more efficiently locate consumption patterns and evaluate deployment, or allocation of resources, to meet targets and goals, according to optimization principles discussed here, or other methodologies.

[0057] Exemplary queries are described below. FIGS. 4-7 depict exemplary embodiments of the queries described below. FIG. 4 depicts at least part of an interface where different J sectors are divided into G.sub.J subsectors, with options to switch to X.sub.j end-use types. These filters invoke the queries where different Y.sub.k energy types can be used to estimate emissions and their proportions (percentages) as shown in Query 3 and Query 4 below.

[0058] FIG. 5 depicts at least part of an interface where geospatial data details are visualized for a given city or location, p and given time t. The interface helps user select a set of R's for any given restrictions, such as Q (Query 1), η (Query 2) by using emissions factors (Query 6).

[0059] FIG. 6 depicts how different elements of the geospatial map can visualize different data entities for each R (shown in Query 5 and Query 8 below).

[0060] FIG. 7 depicts tabular formats of storing geospatial information at the R.sup.th level and aggregating to the regional p.sup.th level.

Query 1

[0061] The queries descried above and below, are shown in FIGS. 4-7 described above. The system can extract a set of artifacts R (buildings) from the table J1 (residential sector) that are built in the year Q1=Q1.sub.b (say Q1.sub.b=1977) and have square footage between Q2.sub.min and Q2.sub.max (say Q2.sub.min=1200 and Q2.sub.max=1300). By applying the operation of selection of rows with a specified year and square footage:

σ.sub.Q.sub.min.sub.≤Q2≤Q2.sub.max σ.sub.Q1=Q1.sub.b(J1)   (15)

[0062] The table J1 now contains only those rows with year Q1.sub.b and whose square footage is between Q2.sub.min and Q2.sub.max. If a user uses the system to extract only the ID numbers of those artifacts, this can be done by projecting everything on the corresponding column:

Π.sub.R(J1)   (16)

[0063] Given a table of unique identifiers (IDs) of these artifacts—the number of which is given by the number of rows in this new table of IDs, such as:

TABLE-US-00001 J1 R Q1 Q2 all other columns R1 1960 1494 . . . R2 1965 1076 . . . R3 1977 1243 . . . R4 1977 1802 . . .
After applying the first selection (artifacts with an appropriate year) we obtain:

TABLE-US-00002 σ.sub.Q1 = 1977(J1) R Q1 Q2 all other columns R3 1977 1243 . . . R4 1977 1802 . . .
Applying the second selection operation, we are left with one row only:

TABLE-US-00003 σ.sub.1200 ≤ Q2 ≤ 1300 σ.sub.Q1 = 1977(J1) R Q1 Q2 all other columns R3 1977 1243 . . .
Finally we project to the column of IDs:

TABLE-US-00004 Π.sub.R σ.sub.1200 ≤ Q2 ≤ 1300 σ.sub.Q1 = 1977(J1) R R3
There is only one artifact that meets the two requirements. Any set of parameters can have similar operations conducted, under different conditions and requirements, both on the interface as well as intermediate steps before visualization.

Query 2

[0064] Similarly to the above relational query 1, the system can also filter out artifacts (e.g. in the residential sector J1) whose aggregate efficiency factor η is at least 0.6 by applying the selection operation to the corresponding column:

σ.sub.η≥0.6(J1)   (17)

After this operation the table J1 contains only those artifacts whose efficiency factor satisfies η≥0.6, as shown:

TABLE-US-00005 J1 R η all other columns R1 0.34 . . . R2 0.91 . . . R3 0.45 . . . R4 0.65 . . .
only R2 and R4 satisfy the condition η≥0.6. Thus, applying further queries to this table, the resulting table is:

TABLE-US-00006 σ.sub.η ≥ 0.6(J1) R η all other columns R2 0.91 . . . R4 0.65 . . .

Query 3

[0065] The system can compute by summing emissions ϕ, for specific X.sub.j by extracting rows with X.sub.j=1 in a given J and return the total emission ϕ.sub.J,X j=1, and then represent this as a fraction of J or total Js. This is achieved by computing a fraction of ϕ.sub.J,X j=1 and ϕ.sub.J and a fraction of ϕ.sub.J,X j=1 and ϕ=Σ.sub.J ϕ.sub.j in the following ways: To each artifact R, the system can associate the corresponding aggregate emission ϕ(R) and storing it in the corresponding column. To compute the total emission of all artifacts in the residential sector J1 that have a non-zero consumption through sub-sector X.sub.j, the system selects those rows in the table J1 that have 1 in the column X.sub.j and computes the sum of values in the column Φ of the modified table:

SUM.sub.Φ(σ.sub.X j=1(J1)),   (18)

where SUM.sub.C(T) is the relational algebra operation of the summation of the values in the corresponding column C of the table T. To compute the ratio of emissions from artifacts in J1 that have nontrivial X.sub.j consumption to the total emission in J1, the system may perform the query:

[00001]$\begin{array}{cc}\frac{{\mathrm{SUM}}_{\Phi }\left({\sigma }_{X_j=1}\left(J1\right)\right)}{{\mathrm{SUM}}_{\Phi }\left(J1\right)}& \left(19\right)\end{array}$

as well as the ratio of the above numerator to the total emission across all sectors,

[00002]$\begin{array}{cc}\frac{{\mathrm{SUM}}_{\Phi }\left({\sigma }_{X_j=1}\left(J1\right)\right)}{{\mathrm{SUM}}_{\Phi }\left(J1\right)+{\mathrm{SUM}}_{\Phi }\left(J2\right)+{\mathrm{SUM}}_{\Phi }\left(J3\right)+{\mathrm{SUM}}_{\Phi }\left(J4\right)}.& \left(20\right)\end{array}$

[0066] Using a particular example, the system operator functions as follows:

TABLE-US-00007 J1 R X.sub.j Φ all other columns R1 0 230 . . . R2 0 120 . . . R3 1 430 . . . R4 1 279 . . .
430+279=709, since only R3 and R4 have X.sub.j=1, where the ratio equals 66%

[00003]$\begin{array}{cc}\frac{430+279}{230+120+430+279}=\frac{730}{1059}\approx 0.66.& \left(21\right)\end{array}$

Given total emissions in other sectors, the ratio is computed analogously.

Query 4

[0067] The system may compute the total emission of all artifacts in a specific subsector or subtype G.sub.j (say G.sub.j=Duplex,Single−family, etc.) of the residential sector J1. By first extracting all rows with G=G.sub.j using selection operation, the system can compute the sum of values in the column Φ of the new table:

SUM.sub.Φ(σ.sub.G=G.sub.j(J1)).   (22)

The system further computes the fraction of emissions of artifacts from the subsector G.sub.j relative to the total emission in J1,

[00004]$\begin{array}{cc}\frac{{\mathrm{SUM}}_{\Phi }\left({\sigma }_{G=G_j}\left(J1\right)\right)}{{\mathrm{SUM}}_{\Phi }\left(J1\right)},& \left(23\right)\end{array}$

and relative to the total emission across all sectors J1, J2, J3, J4,

[00005]$\begin{array}{cc}\frac{{\mathrm{SUM}}_{\Phi }\left({\sigma }_{G=G_j}\left(J1\right)\right)}{{\mathrm{SUM}}_{\Phi }\left(J1\right)+{\mathrm{SUM}}_{\Phi }\left(J2\right)+{\mathrm{SUM}}_{\Phi }\left(J3\right)+{\mathrm{SUM}}_{\Phi }\left(J4\right)}.& \left(24\right)\end{array}$

[0068] As an example:

TABLE-US-00008 J1 R G Φ all other columns R1 Single-family 130 . . . R2 Duplex 340 . . . R3 Duplex 503 . . . R4 Single-family 189 . . .
where G.sub.j is equal to Duplex, hence 340+503=843 since only R2 and R3 are Duplexes, while the ratio (23) equals 73%

[00006]$\begin{array}{cc}\frac{340+503}{130+340+503+189}=\frac{843}{1162}\approx 0.73.& \left(25\right)\end{array}$

Given the total emission in all 4 sectors, the ratio can be computed analogously.

Query 5

[0069] The system can further compute emissions of an artifact in each given energy subtype Y.sub.k and subsector X.sub.j. For each artifact R there is a table describing its consumptions at each fixed sublevels Y.sub.k and X.sub.k, where the emissions of all artifacts in the subsector X.sub.j (of the residential sector J1) that are consuming energy subtype Y.sub.k can be computed.

[0070] First, the system calculates (k, j)-emissions of a given artifact R and calls the table with Γ.sub.K which contains carbon contents and emissions factors of each energy subtypes of the energy type K, with two columns—energy subtypes and carbon content or emission factor per emission type. This process assumes one factor per type of emission, whereas many types of emissions might be present. These tables are used to compute (k, j)-emissions of R via the formula

Φ.sub.k,j(R):=SUM.sub.X.sub.j(σ.sub.F=Y.sub.k(J1(R))).Math.SUM.sub.Γ(σ.sub.F=Y.sub.k(Γ.sub.K)),

where SUM.sub.C simply means the system extracts the element (in the column C) of each of the two tables containing one row only. The computed (k, j)-emission of R is denoted by Φ.sub.k,j(R). Thus the total emission from consumption of the energy subtype Y.sub.k in the subsector X.sub.j is computed via the formula:

Φ.sub.k,j:=Σ.sub.R Φ.sub.k,j(R),   (26)

where the summation is performed over all artifacts R in the sector J1. The proportion of this quantity in the total emission in the sector J1:

[00007]$\begin{array}{cc}\frac{{\Phi }_{k,j}}{{\mathrm{SUM}}_{\Phi }\left(J1\right)}.& \left(27\right)\end{array}$

[0071] The system may filter only those artifacts R in the sector J1 that contribute to the (k, j)-emission (26) by first extracting the essential information from the tables J1(R) for all different artifacts R:

J1(R).sub.k,j:=.sub.R,X.sub.j (σ.sub.F=Y.sub.k(J1(R)))   (28)

[0072] The new table J1(R).sub.k,j contains only one row (corresponding to Y.sub.k) and two columns, R containing information about the artifact, and X.sub.j containing the value of the (k, j)-consumption, where combining all the resulting tables gives:

J1.sub.k,j:=U.sub.R J1(R).sub.k,j,   (29)

where the union is over the set of all artifacts in J1. J1.sub.k,j is a long table of (k, j)-consumptions of artifacts in J1, where applying the operation of natural join to J1 and J1.sub.k,j:

J1J1.sub.k,j.   (30)

Gives a table has an extra column X.sub.j containing information about (k, j)-consumptions, from where we can remove those artifacts in J1 that have trivial (k, j)-consumption,

σ.sub.X.sub.j>0(J.sup.1J1.sub.k,j),   (31)

and just extract the list of those artifacts in J1 with positive (k, j)-consumption,

Π.sub.R(σ.sub.X.sub.j>0(J1J1.sub.k,j))   (32)

[0073] In short, working with table J1 and selecting only those rows (R) for which the corresponding (k, j)th entry of the consumption table J1(R) is positive:

σ.sub.R: SUM .sub.X j.sub.(σ.sub.F=Y k.sub.(J1(R)))>0 (J1)   (33)

Query 6

[0074] In order to obtain the table of artifacts in a given sector (say J1) together with their localized (k, j)-emissions, the system may exploit the table J1.sub.k,j of (k, j)-consumptions and modify its second (X.sub.j) column by multiplying it with the carbon content:

Π.sub.R,X.sub.j.sub..Math.SUM.sub.Γ.sub.(σ.sub.F=Y k.sub.(Γ.sub.K.sub.)) (J1.sub.k,j).   (34)

If carbon contents are constant values and do not vary by sector, location or time, they may have fixed denotions, like Γ.sub.k for the carbon content or emissions factor in the energy subtype Y.sub.k, which can be simplified as:

ϕ.sub.k,j:=Π.sub.R,X.sub.j.sub..Math.Γ.sub.k(J1.sub.k,j).   (35)

Query 7

[0075] The system can further compute the fractions of emissions in a given sector (say J1) relative to the total emission in all 4 sectors:

[00008]$\begin{array}{cc}\frac{{\mathrm{SUM}}_{\Phi }\left(J1\right)}{{\mathrm{SUM}}_{\Phi }\left(J1\right)+{\mathrm{SUM}}_{\Phi }\left(J2\right)+{\mathrm{SUM}}_{\Phi }\left(J3\right)+{\mathrm{SUM}}_{\Phi }\left(J4\right)}.& \left(36\right)\end{array}$

If the values of emissions in sectors were SUM.sub.Φ(J1)=1290, SUM.sub.Φ(J2)=2120, SUM.sub.Φ(J3)=3400 and SUM.sub.Φ(J4)=1780, then the ratio for J1 would be 15%

[00009]$\begin{array}{cc}\frac{1290}{1290+2120+3400+1780}=\frac{1290}{8590}\approx 0.15& \left(37\right)\end{array}$

Query 8

[0076] Knowing for each artifact R the emissions ϕ as well as the energy consumption by X.sub.j and Y.sub.k, the system can aggregate based on the ownership of the artifacts, let's say S individual person (J=1) or individual company or organization (J=2). Getting S-level ϕ.sub.S by aggregating ϕ.sub.R in the selected p and t. The system will extract all R's belonging to each unique person or organization and compute the total emission of this entity and compare it to the total emission in this J and this city, p in any given time t.

[0077] The system can thus compute emissions of artifacts belonging to each entity. For simplicity, only artifacts confined to one city/region p is considered in this example, but this operation can be done across every city. Given that there is a table S encoding the relation between unique artifacts R and unique individuals S (ID number). If the system wants to solve the above computation for the individual with the ID number i, it filters all artifacts in S belonging to i:

S.sub.i:=σ.sub.I=i(S)   (38)

Then, the natural join operation of S.sub.i with J1:

S.sub.iϕ.sub.k,j,   (39)

where ϕ.sub.k,j is the table of emissions in J1. By joining rows from S.sub.i (corresponding to artifacts R) with the corresponding rows in ϕ.sub.k,j. The summed emissions in the resulting table to compute an entity's emissions in a given city/region, p over time, t is given by:

Φ.sub.k,j(i):=SUM.sub.X.sub.j(S.sub.iϕ.sub.k,j).   (40)

The system can also compare the computed entity's emissions to the total emissions in J1 sector by computing the respective fraction:

[00010]$\begin{array}{cc}\frac{{\Phi }_{k,j}\left(i\right)}{{\mathrm{SUM}}_{\Phi }\left(J1\right)}.& \left(41\right)\end{array}$

[0078] In some embodiments, the width of a particular category or subcategory on one axis may be proportional to the total amount of that category or subcategory relative to the other sub-categories. For example, if petroleum makes up 50% of fossil fuel usage, natural gas makes up 40% and coal makes up 10%, then the width of the petroleum bar on the type axis may be five times that of the coal bar, and so forth. The same bar scaling techniques can be applied on the sector axis, such that the bar for businesses is twice the depth of the bar for residential if business use accounts for twice as much fossil fuel consumption as residential use. In other embodiments, bars have equal widths and depths. In still other embodiments, bars of fixed and variable width and depth may be combined. For example, bars of fixed depth may be combined with bars of variable width, or bars of fixed with and depth at the top level may be combined with bars of variable width and depth at the drill-down level. Other combinations are also contemplated as being within the scope of the invention.

[0079] In some embodiments, user interface 200 is operable to track trends over time as well as evaluate the data at a particular point in time such as the present. For example, user interface 200 may include a slider to adjust the point in time for which the data is displayed. Alternatively or in addition, user interface 200 may animate to display trends in the data over time. This can be done via Queries 1-8 and FIGS. 4-7 described above. In some embodiments, the animation may display the current view (the full view or the drill-down view) as it changes over time. In other embodiment, the animation process resets the view to the full view. The user interface 200 may also display future projections. User interface 200 may also animate in other ways. For example, the viewing angle may rotate around the z-axis so as to better display, for example, sector-type categories or subcategories that would be obscured by other categories or subcategories. In still other embodiments, various greenhouse gas reduction strategies and models (as further discussed below) can be animated. For example, given a model for reducing the greenhouse gas emissions from a particular sector, the animation may show the projected future breakdown of the energy use in that category by non-carbon-based energy sources as time progresses.

[0080] In particular, user interface 200 may be used to generate or examine a greenhouse-gas mitigation strategy. In some embodiments, various sources of power may be colored differently. For example, carbon-based fuels may be dark grey or black, with wasted energy colored grey. By contrast, alternative, non-greenhouse-gas generating fuels may be green, or otherwise brightly colored. Thus, in such embodiments, as the energy mix becomes less carbon-dependent, user interface 200 will transition from black and grey to bright and multicolored.

[0081] The user interface 200 may show trends past, present, and future based inputs into the system and how the trends may react and changed based on the inputs. An algorithm for optimizing resource allocation and usage is described below.

Exemplary Implementation of the Algorithm

[0082] As mentioned previously, embodiments of the invention can be employed to generate or examine a greenhouse-gas mitigation strategy. Described below are two linear programming implementations of how such a model can be generated. The two programming implementations described below describe methods of receiving inputs of source data from benign and harmful energy sources, determining energy used and energy wasted by transfer processes. Further, the programming implementations provide methods of allocating the resources to minimize the harmful energy sources and maximize the benign energy sources. Further, the resource data from the inputs, the energy used and wasted in the transfer processes, and the optimization data and allocation data may be provided in the visualizations depicted in the figures.

[0083] Consider a system where input sources, Z(t) could be distributed into two main groups, either through discrete categorization (binary groups) or ranked according to a continuous value of ecological or environmental value creating or value destroying. Notation for describing the modelling process is given below.

[0084] A(t) is the total amount of resources coming from an ecologically benign source in any given time or period, and B(t) is the total amount of resources coming from an ecologically harmful source in any given time or period. Then Z(t)=A(t)+B(t). T broadly represents transformations of resources. T(t) is any transformation step that takes the intermediate source of resource before converting the resource to different form, TA(t) is the inflow of a given transformation step T(t) that takes the intermediate source of resource from an ecologically benign source A(t) before converting the resource to different form, and T.sup.B(t) is the inflow of a given transformation step T(t) that takes the intermediate source of resource from an ecologically harmful source B(t) before converting the resource to different form. Then, following from the above, T(t)=T.sup.A(t)+T.sup.B(t).

[0085] Next, J broadly describes the inflow of resources. J.sup.T.sup.(i) (t) is any inflow of resources transformed by step T(t) and then consumed by the i.sup.th end-user J. As an example, J.sup.T.sup.(1) (t) represents the amount of T transformed resources consumed by J.sub.1. J.sup.A.sup.(i) (t) and J.sup.B.sup.(i) (t) is any inflow of resources consumed by the i.sup.th end-user J, of form A, derived from A(t), and of form B, derived from B(t), respectively, such that J.sup.(i) (t)=J.sup.A.sup.(i) (t)+J.sup.B.sup.(i) (t)+J.sup.T.sup.(i) (t). FIG. 8 depicts one sector-type category 800 of user interface 200 decomposed into inflows 802 and transformations 804. As such, FIG. 8 depicts the relative relationship of the inflow from a benign source J.sup.A, the inflow from a harmful source J.sup.B, the transition of the inflow from the benign source T.sup.A, and the transition of the inflow from the harmful source T.sup.B.

[0086] J can broadly be decomposed into N and M. N.sup.(i) (t) is any inflow of resources that is unconsumed or wasted by the i.sup.th end-user J, and thus, not directly used. As an example, N.sup.(1) (t) represents the amount of resources unused by J.sup.(1). M.sup.(i) (t) is any inflow of resources, that represents total resource inputs used by the i.sup.th end-user J, and thus, directly used. As an example, M.sup.(1) (t) represents the amount of resources used by J.sup.(1).

[0087] Turning now to process efficiency, let θ.sup.(i) (t) be thus the ratio of wasted or unconsumed resource inflow (i.e., N.sup.(i) (t)) to the total inflow of J.sup.(i) (t) (i.e., N.sup.(i) (t)+M.sup.(i) (t)) by the i.sup.th end user: η.sup.(i)=N.sup.(i) (t)/(N.sup.(i) (t)+M.sup.(i) (t)). Let T.sup.C(t) is any outflow of T(t) transformation step that is wasted during the transformation process. Then μ.sub.T is thus the ratio of the wasted outflow T.sup.C(t) for the transformation step to the total outflow T(t):μ.sub.T=T.sup.C(t)/T(t). For notational convenience, let E(t) be the total outflow that represents total resources that are used or consumed during total process, including M.sup.(i) (t): E(t)=σ.sub.i.sup.n M.sup.(i) =M.sup.(1)(t)+M.sup.(2) (t); Where i=1, 2, . . . n. For one example, this description uses n=2; however, any value of n is contemplated. Similar to E, let C represent the total outflow that represents total resources that are wasted or unconsumed during total process (including N.sup.(i) (t) and T.sup.C (t): C(t)=Σ.sub.i.sup.n N.sup.(i)+T.sup.C (t). FIG. 9 depicts the same sector-type category 900 as FIG. 8 of user interface 200 decomposed instead into wasted energy C, utilized energy of the benign source E.sup.A, and utilized energy from the harmful source E.sup.B.

[0088] Now let T.sup.AC (t) be the proportion of any wasted outflow during the transformation process of step T(t) that is attributed to origin from an ecologically benign source, A(t) and let T.sup.BC (t) be the proportion of any wasted outflow during the transformation process of step T(t) that is attributed to origin from an environmentally harmful source, B(t). This gives T.sup.AN (t)+T.sup.BCN (t)=(1−η)*(p.sup.T(B)*J.sup.T)+(1−μ)*η*T.sup.A.Math.T.sup.AC.sup.(i) (t)) is the outflow of resources from the transformed step T(t) by J.sup.(i) attributed to environmentally benign source, A(t), so T.sup.AC (t)=p.sup.T (A)*T.sup.C, and T.sup.BC.sup.(i) (t) is the outflow of resources from the transformed step T(t) by J.sup.(i) attributed to environmentally harmful source, B(t), so T.sup.BC(t)=p.sup.T (B)*T.sup.C.

[0089] This leads to J.sup.AN.sup.(i) (t) being any inflow of resources originating from an environmentally benign resource, A(t), transformed by step T(t) and then consumed by the i.sup.th end-user J. As an example, J.sup.T,A.sup.(1) (t) represents the amount of T transformed resources of source A(t), consumed by J.sup.1: J.sup.T,A.sup.(i) (t)={T.sup.A(t)/T(t)}*J.sup.T.sup.(i) (t). Similarly, J.sup.BN.sup.(i) (t) is any inflow of resources originating from an environmentally harmful resource, B(t), transformed by step T(t) and then consumed by the i.sup.th end-user J. As an example, J.sup.T,B.sup.(i) (t) represents the amount of T transformed resources of source B(t), consumed by J.sup.(1): J.sup.T,B.sup.(i) (t)={T.sup.B (t)/T(t)}*J.sup.T.sup.(i) (t).

[0090] Broadly, Φ represents a quantification function for the environmental harmfulness of each unit of resource consumption and P represents economic value or cost associated with a resource source converted from the quantified value produced by Φ. Hence P(A(t)) represents the economic value or cost converted from the quantified factors Φ(A(t)) for the total resource inflow, A(t) and P(B(t)) represents the economic value or cost converted from the quantified factors Φ(B(t)) for the total resource inflow, B(t). P(E(t)) and P(C(t)) is thus the combined economic value or cost converted from the quantified source factors, P(A(t)) and P(B(t)). For Φ(E(t)) and Φ(C(t)), the algorithm can thus represent and compute in terms of sourced Φ(A(t)) and Φ(B(t) without regards to the economic value or cost function P. Now Φ(A(t)) represents the quantification factor for the environmental harmfulness of each unit of resource consumption of environmentally benign type A(t), such that Φ(A(t) can be approximated by 0. Φ(B(t)) represents the quantification factor for the environmental harmfulness of each unit of resource consumption of type B(t), such that Φ(B(t)>0. For simplification, some embodiments may assume a fixed value for Φ(B(t)) such as 2.

[0091] Turning to an exemplary embodiments depicted in FIG. 10, let a system have an overall quantity of throughput flow of Z (for concreteness, say 80 units), which can be of multiple types, (e.g., A and B) that collectively add up to Z (for example 60 units of B and 20 units of A). In particular, let p(A) and p(B) be the portion of Z which are made up of A and B respectively. Now, at least a portion of these flows from A and B go to a transformation step T, portions of which are represented by T.sup.A and T.sup.B that add up to T. Similar to p, let p.sup.T (A) and p.sup.T (B) be the proportions of T represented T.sup.A and T.sup.B.

[0092] The transformation step T has an efficiency factor, μ, where it separates the wasted portion, T.sup.C=μT, to C, from the remaining portion (J.sup.T) which is used and consumed along with direct A and B flow types, called J.sup.A and J.sup.B respectively (which add up to J). As above, let p.sup.J(A) and p.sup.J(B), proportions of J represented J.sup.A and J.sup.B. This, in turn, gives T.sup.C=T.sup.AC+T.sup.BC where T.sup.AC=p.sup.T (A)*T.sup.C and T.sup.BC=p.sup.T(B)*T.sup.C. T.sup.A, T.sup.B, J.sup.A and J.sup.B, collectively add up to total A and B, respectively, and thus to total system flow Z.

[0093] Finally, the J consumption stage has its own efficiency factor η, where it separates the actual used-portion to E, and the remaining wasted portion to C. In terms of T and J, first T.sup.A=T.sup.AN+T.sup.ACN+T.sup.AC and T.sup.B=T.sup.BN+T.sup.BCN+T.sup.BC, since T=T.sup.A+T.sup.B (where T.sup.ACN=(1−η)*(p.sup.T (A)*J.sup.T) and T.sup.BCN=(1−η)*(p.sup.T (B)*J.sup.T) and T.sup.AN=(1−μ)*η*T.sup.A and T.sup.BN=(1−μ)*η*T.sup.B)). Next, J.sup.A=J.sup.AN (η*J.sup.A)+J.sup.ACN (1−η*J.sup.A) and J.sup.B=J.sup.BN (η*J.sup.B)+J.sup.BCN (1−η*J.sup.B) (where J=J.sup.A+J.sup.B+J.sup.T); and finally J.sup.T=(1−μ)*T×J.sup.TA+J.sup.TB (where J.sup.TA=(1−μ)*T.sup.A and J.sup.TB=(1−μ)*T.sup.B).

[0094] The total Z could thus be visualized, and/or simplified in various reorganized terms, such as A, B, E, C or J and T—or a combination thereof. This described matrix system can then be enumerated within our example, with real values and terms. Then, the system can be optimized, where A and E should be maximized while B and C is minimized. After many iterations of the system with the above in mind, we can exit the optimization program when this above optimal has been achieved. Visually, over time t, the system should behave in this manner and user interface 200 should reflect the transition from ecologically harmful sources B to ecologically benign sources A based on the costing of quantification function Φ.

[0095] A concrete example of the above process is depicted in FIGS. 10-11. For the example of FIGS. 10-11, let Φ.sup.A be 0 and Φ.sup.B be 2 units of pollution, and assume that the costing function P gives $1.5 per unit of pollution. FIG. 10 depicts the starting state of the example above where B of 60 units, so Φ is 120 units and P is $180. This amount can be used to “purchase” items, v (representing reallocations and transformations of energy sources in different ways), thereby changes the system parameters in meeting the desired conditions (i.e., p(A) and p.sup.T (A) goes to 1, which implies p(B) and p.sup.T (B) goes to 0, and 72 goes to 1, such that μ goes to 0). We can introduce, a set of v where, for each i.sup.th v, we have a numeric condition that can satisfy required changes. For example, we might have v such that

v.sup.(i)=1; Δη=η.sup.(1)−η.sup.(0)=0.10; so η.sub.1=0.60

v.sup.(i)=4, Δη=η.sup.(1)−η.sup.(0)=0.15, so η.sub.1=0.65;

Δp.sup.T (A)=0.20, p.sup.T (A).sup.(1)−p.sup.T (A).sup.(0)=0.20 so p.sup.T (A).sup.(1)=0.40 implying p.sup.T(B).sup.(1)=0.60

v.sup.(i)=3, Δμ=μ.sup.(1)−μ.sup.(0)μ−0.20, so μ.sup.(1)=0.40; and

Δp(A)=0.15, p(A).sup.(1)−p(A).sup.(0)=0.15 so p(A).sup.(1)=0.40 implying p(B).sup.(1)=0.60.

[0096] This process can be iterated, computing the imputed pollution cost P(B(t)) at each iteration until p(B) such that P(B(t)), at each iteration, different techniques for selecting the v.sup.(i) to “purchase” using P(B(t)) are contemplated. For example, the v.sup.(i) with the least cost could be selected. As another example, the v.sup.(i) with the maximum reduction of Φ (i.e. maximized ϕ.sup.(1)−ϕ.sup.(0) or highest total pollution reduction) could be selected. Alternative, a hybrid of these methods could be used, such as maximizing (ϕ.sup.(1)−ϕ.sup.(0))/P(v). (i.e., maximizing pollution reduction per cost). After a number k of iterations, the system might result in the state depicted in FIG. 11. There, as can be seen, p(B) and μ have been minimized resulting in both B and C being approximately 0.

[0097] Thus, in summary, the exemplary embodiment of the system above can visualize (using user interface 200) how the initial conditions could evolve through iterations of parameterized changes and end at the optimal state, where due to step-wise increases of p(A), B and J.sup.B are eliminated and due to step-wise increases of p.sup.T (A), we eliminate T.sup.B, thereby achieving the goal of maximizing ecologically benign energy sources A while minimizing ecologically harmful energy sources B. Next, due to step-wise decreases of μ, we eliminate T.sup.C, one component of C, and due to step-wise increases of η, we eliminate the remaining component of C, thereby meeting our goal of minimizing wasted energy and maximizing energy utilization.

[0098] FIGS. 12A-15B depict an exemplary user interface 1200 of the above process visually. Initially, at FIG. 12A the goal of reducing the current pollution Φ of 500 units, by a target of ΔΦ of 100 units (to 400 total units) by 2030 depicted at evaluation 1202. This required ΔΦ is stored in the database for the specific user. Based on the generation of Φ units of pollution due to the consumption of M B's; a user can interactively enter change parameters (e.g., consumption of a particular B.sup.(i)) by moving the dials of the interface at input 1204, using a hardware device, like a computer mouse; or by using a touch-screen interface of a mobile device, tablet or kiosk monitor or via terminal line code and application programming interface syntax. From their initial positions, B (M=1) and A (N=1) are adjusted by a user or other entity, to direct changes in Δp.sup.T (A), Δp.sup.T (B) to the system's optimization model.

[0099] Following this parameter input, the system stores the changes (within a predetermined margin of error) into the database, and indicated to the user that the reduction in pollution was ΔΦ=20 units, meaning that 20% progress has been made if the indicated change is implemented, as depicted in FIG. 12B.

[0100] Next, as depicted in FIGS. 13A and 13B, is a solution evaluation visualization 1200 for each of the four solution entities, the corresponding η can be changed, such that, each Δη will now be stored in the database. After the change, which is displayed to the user, the progress metric is updated, communicating to the user or other entity, that an additional reduction in pollution of ΔΦ=15 units have been eliminated due to changes in the Δη matrix, as shown in FIG. 13B.

[0101] Further, as depicted in FIG. 14A, the user can also reduce Δμ, as shown, by 30% which is equated to a further reduction in pollution of ΔΦ=35 units. Thus, as depicted in other words, there remain only 30 units of Φ that to be eliminated, and the last Δμ step had reduced 35 units of Φ, as shown in FIG. 14B. As shown in the changes, these are stored as targets for each solution entity when they see their individual, separated cube interfaces and analytics, such that the total changes are also specific to the nuanced context of each solution entity or user.

[0102] As can be seen in FIGS. 15A and 15B, due to the above changes and storing of initial vs. final values stored in the database, a chart can be generated for the user to view the starting Φ value, and how the changes can lead to reductions in ϕ, knowing there is a goal of 100 units of Φ by a certain year or in a particular timeframe.

[0103] Each of these reductions in Φ is associated with or derived from changes in more granular parameters, such as Δp.sup.J (A), Δp.sup.J (B), Δp.sup.T (A), Δp.sup.T (B), Δp(A), Δp(B), Δη, and Δμ, etc., which is provided by different solutions. Each solution v may be associated with a cost, P(v), representing a monetary cost to implement individual solution v. This cost could be evaluated by the user by simply viewing the graph or chart, not only in terms of Φ, but in terms of P(Φ(t)) which can be a constant or a variable. Instead of doing $750 of damages by consuming 500 units, the user can instead invest the $656 required to purchase a number of solutions v.sup.(i); thereby making 75% progress to reaching its 2030 goal of 400 units.

[0104] Each user can thus evaluate their costs, total P(v) for all v.sup.(i), and their benefits, which are total Sv(t) savings over time, in addition to the $150 in damages ($750-$600) they would have incurred, had they not made the changes to reduce Φ by 100 units by adopting the chosen v.sup.(i). As shown in FIG. 15B, users can also move from the real-time cube view, to a more temporal progress view with real-time changes as interface dials are manipulated by the user or an entity on their behalf. As interface dials are changed, a set of v's can be automatically visualized on the interface, clicking which can show how the visuals and values will change depending on v-level parameters.

[0105] Similar to the user interface manipulations from FIGS. 12A-B, 13A-B and 14A-B, the analytical cube 1602 can be directly accessed via the user interface generally referred to by numeral 1600, as depicted in FIG. 16. An exemplary inflow and transition process is shown in the visualization 1204 similar to FIGS. 8 and 9. From the perspective of other entities, such as each v solution entity, the animated cube of any user (or group of users) can be viewed to determine the corresponding information such as parameters changed by that user, and of which type, and what is the corresponding cost (P(v)) to purchase, or savings created, Sv(t), or pollution/social and environmental costs reduces, in terms of P(B(t)) or Φ(B(t)), respectively.

[0106] This information may be stored in a database and made available to the user in the user interface as shown before. This interface may function like a bidding system, or marketplace, introducing different entities as buyers, sellers or intermediaries, like rating organizations, financial institutions, corporations or governmental agencies—programmatically asking user S to lower Φ or social costs, P(B(t)) by adopting certain v's which meet cost, P(v) or savings, Sv(t) thresholds. Such a user interface may additionally aggregate and visualize the incoming or stored data, set or suggest change targets, and load potential solutions or options the user has available to meet a particular target through any stored or accessed v-level change parameters Δp.sup.J (A), Δp.sup.J (B), Δp.sup.T (A), Δp.sup.T (B), Δp(A), Δp(B), Δη and Δμ, etc.

[0107] Details of each solution v.sup.(i) stored in the database may be depicted, such as the level of the change parameters (Δp.sup.J (A), Δp.sup.J (B), Δp.sup.T (A), Δp.sup.T (B), Δp(A), Δp(B), Δη and Δμ, etc.), the cost to purchase, P(v), the expected savings over time ΔSv(t), the estimated reduction in pollution, ΔΦ (or conversely, “progress made” or subsequently, the estimated reduction in economic harm or pollution cost P(t), by the a particular user in a given time period t).

[0108] FIG. 17 depicts a simple interface on how each of the v-solution data parameters may be inputted into the system by a user or set of users.

[0109] FIG. 18 depicts an interface on how each of the R-level data in a given region, p can be entered by S entities and encodes attributes such as Q, X.sub.j, G.sub.j, Y.sub.k, t etc. which can be used for filtering geospatially.

[0110] FIG. 19 depicts an interface to input other J sector R level data (such as J=transportation, R=passenger vehicle) for any given t.

[0111] FIG. 20 depicts an interface to input data on transformation T of an energy source (T.sup.A, T.sup.B) as well as the sector and location of the consumption, J.sup.T over time t. Same interface can be used compute emissions, as well as enter v solutions through “Act now” button.

[0112] FIG. 21 depicts a further interface on the act now button, which is similar to v-solution entry interface, listing Δp.sup.J (A), Δp.sup.J (B), Δp.sup.T (A), Δp.sup.T (B), Δp (A), Δp(B), Δη and Δμ as well as P(v), ΔΦ and ΔSv(t) for the algorithm inputs and storing outputs.

[0113] FIG. 22 depicts the input of each region p, denoting which plant or plants, R provide energy source, the cost and price of said energy and the emissions for that energy, for a given time t. Instead of p level regional data, S level entity data entry is also possible.

[0114] FIG. 23 depicts the input of plants R, denoting which location p, it provides energy to and the emissions as a result for a given time t.

[0115] FIG. 24 depicts the input of each plant R, by energy source, efficiency and emissions for any given region p, and time t.

[0116] FIG. 25 depicts v-solution entry interface, listing Δp.sup.J (A), Δp.sup.J (B), Δp.sup.T (A), Δp.sup.T (B), Δp(A), Δp(B), Δη and Δμ for the algorithm inputs, including solution purchase parameter like P(v), ΔΦ and ΔSv(t) etc.

[0117] FIG. 26 depicts the R.sup.th level plant, as well as energy source, efficiency and emissions data on can be shown on any topographical map, given known location in any given region p, and time t. Similar to earlier shown tabular (11) format of storing geospatial information at the R.sup.th level, where R is a powerplant, computing emissions from different energy sources and aggregating to the regional p.sup.th level in time t.

[0118] FIG. 27 presents a cost and savings user input 2700 that allows the user to input changes and dynamically visualize the past, present, and future changes to the allocation and costs. Visually, each of the changes to dials are reflected real-time on the interface via at least one of color, texture, and shape changes in the cube, (e.g. becoming more colorful as more progress is made). Additionally, the progress can be tracked numerically at numerical button 2702, and/or absolute Φ tracked at absolute button 2704. If the user is unable to buy solutions from the indicated v.sup.(i), they could also offset the Φ amount, such as communicating with a higher-order entity with the purchasing power to buy solutions that colors the user's cube, or reduces their Φ. For example, this could be portions of a municipal bond or climate/green bond, that increases Δp.sup.T (A) and reduces Δp.sup.T (B) that is only possible after pooling offset resources from multiple users, or securitized and transacted as a tradable instrument. Further, all progress may be tracked for all customizable conditions past, present, and future at selection block 2706. An embodiment of such an interface is succinctly shown in FIG. 27.

[0119] FIG. 28 depicts Any given user S can see the transformed visualization (200) of the cube and cubelets showing energy source consumption in each of the artifact, R's owned by the user, S, with S-level unique code identifier, with the ability to aggregate locationally—city P*, state P′, national q or global Q, based on relational queries identified earlier.

[0120] FIG. 29 depicts a user interface presenting data such that each user can understand the emissions Φ given their energy source visualization, and chose to either escrow the monetary amount, P(Φ) or P(B) in any given t in any currency to either (i) purchase relevant v solutions knowing the P(v) of each v solution, (ii) invest P(Φ) amount in community-level projects, either as a trackable donation or purchasing a security in an investment project, such as green bonds, (iii) make charitable donations of P(Φ) to other S entities, through unique alphanumeric IDs or (iv) buy carbon offsets or renewable energy credits in the amount of P(Φ) and get a discount code (another unique alphanumeric ID) from participating retailers.

[0121] FIG. 30 depicts an exemplary user interface generally referenced by the numeral 3000. The targets 3002 for each user could be edited or updated by the user or by other entities with the ability, ownership or legal/regulatory authority to change to new or different targets. For example, FIG. 30 shows what happens when a user or other authorized entity increases Δp(A) by 2 percentage points, from 14% to 16%; and decrease Δp(B) by 4 percentage points from 18% to 22%; the individual user's interface can be dynamically automatically updated to reflect new goals, and each solution table can be similarly automatically updated to match new goals. The recorded changes are presented in box 3006.

[0122] Turning now to FIG. 31, an alternative optimization algorithm is depicted and generally referenced by the numeral 3100. Starting from the previous optimization algorithm discussed above and ending at the starting point iteration of FIG. 10, then moving to an iterative process shown in FIG. 31 and described below, the system may then be generalized to a system such that A and B split into N and M sub-types A.sup.(1), . . . , A.sup.(N) and B.sup.(1), . . . , B.sup.(M), such that A=Σ.sub.i=1.sup.N A.sup.(i) and B=Σ.sub.i=1.sup.M B.sup.(i). Further,

[0123] S consumption stages J.sup.(1), . . . , J.sup.(S), such that J=Σ.sub.i=1j.sup.(i),

[0124] T.sup.C splits into S stages, say T.sup.C.sup.(1), . . . , T.sup.C.sup.(S) such that T.sup.C=Σ.sub.i=1.sup.S T.sup.C.sup.(i),

[0125] S actual-used portion E.sup.(1), . . . , E.sup.(S) such that E=Σ.sub.i=1.sup.S E.sup.(i),

[0126] S wasted portion C.sup.(1), . . . , C.sup.(S) such that C=Σ.sub.i=1.sup.S C.sup.(i).

[0127] The optimization problem comprises simultaneously maximizing the inflow from benign sources A.sup.(1), . . . , A.sup.(N) and the outflow of utilized energy E.sup.(1), . . . , E.sup.(N) and simultaneously minimizing the inflow from harmful sources B.sup.(1), . . . , B.sup.(M) and the outflow of wasted energy C.sup.(1), . . . , C.sup.(N) where A, E, B, and C are evaluated at each time (t). After k+1 iterations of the system, the program ends when the optimal solution is achieved. FIG. 31 depicts the flow process for achieving the optimal solution.

[0128] The inflow sources A and B depend on portion factors p(A) and p(B) such that p(A)+p(B)=1. Accordingly, p(A) comprises N components p(A.sup.(1)), . . . , p(A.sup.(N)) and p(B) comprises M components p(B.sup.(1)), . . . , p(B.sup.(M)) such that p(A)=Σ.sub.i=1.sup.N p(A.sup.(i)) and p(B)=Σ.sub.i=1.sup.M p(B.sup.(i)). Therefore, there exists N+M constants α.sup.(1), . . . , α.sup.(N), β.sup.(1), . . . , β.sup.(M) such that A.sup.(i)=Σ.sub.i=1.sup.N α.sup.(i)p(A.sup.(i))Z and B.sup.(i)=Σ.sub.i=1.sup.N β(B.sup.(i)p(B.sup.(i))Z. Therefore,

p(A)+p(B)=α.sup.(1)p(A) + . . . +α.sup.(N)p(A)+β.sup.(1)p(B) + . . . +β.sup.(M)p(B)=1

[0129] In order to maximize all quantities A.sup.(1), . . . , A.sup.(N) and simultaneously minimize all quantities B.sup.(1), . . . , B.sup.(M) each component α.sup.(i) p(A) should be maximized and each component β.sup.(i)p(B) should be minimized. Because Σ.sub.i=1.sup.N α.sup.(i)=1 and Σ.sub.i=1.sup.N β.sup.(i)=1 the total quantity p(A) should be maximized and the total quantity p(B) should be minimized.

[0130] As described above, the flows go to the transformation step T, portions of which are represented by T=Σ.sub.i=1.sup.N T.sup.A.sup.(i)+Σ.sub.i=1.sup.N T.sup.B.sup.(i). Accordingly, the proportion factors of the flow p.sup.T (A) and p.sup.T (B) split into p.sup.T (A.sup.(1)), . . . , p.sup.T (A.sup.(N)) and p.sup.T (B.sup.(1)), . . . , p.sup.T (B.sup.(M)) where p.sup.T(A.sup.(i))=ζ.sup.(i)p.sup.T (A) and p.sup.T (B.sup.(i))=ζ.sup.(i)p.sup.T (B) at every time (t) and for N+M proportion factors ζ.sup.(1), . . . , ζ.sup.(N) and ζ.sup.(1), . . . , ζ.sup.(N) such that Σ.sub.i=1.sup.N ζ.sup.(i)=1 and Σ.sub.i=1.sup.M ζ.sup.(i)=1. Therefore, T.sup.A=ζ.sup.(i)p.sup.T (A)T and T.sup.B=ζ.sup.(i)p.sup.T (B)T. Therefore,

p.sup.T (A)+p.sup.T (B)=ζ.sup.(1)p.sup.T (A) + . . . +ζ.sup.(N)p.sup.T (A)+ζ.sup.(1)p.sup.T (B) + . . . +ζ.sup.(M)p.sup.T (B)=1

[0131] Therefore, the transformation portions should be maximized and minimized accordingly. Specifically, the transformation portion of the inflow from benign sources p.sup.T (A) should be maximized and the transformation portion of the inflow from harmful sources p.sup.T (B) should be minimized.

[0132] At transformation step T at time (t) the inflow from benign sources and the inflow from harmful sources is mixed, then separated again into S components T.sup.(1), . . . , T.sup.(S), where total transition is T=Σ.sub.j=1.sup.S T.sup.(j). Therefore, each T.sup.(j) is a portion of T represented by proportion factors γ.sup.(1), . . . , γ.sup.(S). So, the total ideal transformation step may be represented by:

T=γ.sup.(j) (Σ.sub.i=1.sup.N T.sup.A.sup.(i)+Σ.sub.i=1.sup.N T.sup.B.sup.(i))=Σ.sub.i=1.sup.N γ.sup.(j)ζ.sup.(i)p.sup.T (A) T+Σ.sub.i=1.sup.M γ.sup.(j)ζ.sup.(i)p.sup.T (B) T

[0133] Further, the transformation step T has a transformation efficiency factor represented by μ(t), where it separates a wasted portion T.sup.C (t)=μ(t)T, to C. The used portion J.sup.T is consumed along with the direct flow A from benign sources and the direct flow B from harmful sources called J.sup.A and J.sup.B respectively which are dependent on the portion factors p.sup.J (A) and p.sup.J (B). The total used flow may be represented by J=J.sup.A+J.sup.B+J.sup.T.

[0134] There exists S efficiency components μ.sup.(1), . . . , μ.sup.(S) with corresponding proportion factors Ξ.sup.(1), . . . , Ξ.sup.(S) where Σ.sub.j=1.sup.S Ξ.sup.(j)=1, such that the total efficiency factor is represented by μ(t)=Σ.sub.j=1.sup.S Ξ.sup.(j) μ(t). For all j=1, . . . , S, the transformation loss is:

[00011]$T^{C^{\left(j\right)}}={\mu }^{\left(j\right)}T={\Xi }^{\left(j\right)}\mu \left(t\right){\gamma }^{\left(j\right)}\left({.Math.}_{i=1}^N{p}^T\left(A\right)T\right)+{\Xi }^{\left(j\right)}\mu \left(t\right){\gamma }^{\left(j\right)}\left({.Math.}_{i=1}^M{p}^T\left(B\right)T\right)$$\mathrm{and}$$J^{T^{\left(j\right)}}=\left(1-{\mu }^{\left(j\right)}\right)T=\left(1-{\Xi }^j\mu \left(t\right)\right){\gamma }^{\left(j\right)}\left({.Math.}_{i=1}^N{p}^T\left(A\right)T\right)+\left(1-{\Xi }^j\mu \left(t\right)\right){\gamma }^{\left(j\right)}\left({.Math.}_{i=1}^M{p}^T\left(B\right)T\right)$

In order to maximize the inflow from benign sources A.sup.(1), . . . , A.sup.(j) and the outflow of utilized energy E.sup.(1), . . . , E.sup.(j) and minimize the inflow from harmful resources B.sup.(1), . . . , B.sup.(j) and the outflow of wasted energy C.sup.(1), . . . , C.sup.(j), each parameters μ.sup.(1), . . . , μ.sup.(j) should be minimized.

[0135] Further, J.sup.T.sup.(j) may be written as,

T.sup.C.sup.(j)=T.sup.AC.sup.(j).sup.+T.sup.BC(j)=Σ.sub.i=1.sup.N μ.sup.(j)(t)γ.sup.(j)ζ.sup.(i)p.sup.T (A) T+Σ.sub.i=1.sup.M μ.sup.(j)(t)γ.sup.(j)ζ.sup.(i)p.sup.T (B) T

The proportion factors of J.sup.A and J.sup.B, p.sup.J (A) and p.sup.J (B) may be split into N and M components respectively, so there exists δ.sup.(1), . . . , δ.sup.(N) and ∈.sup.(1), . . . , ∈.sup.(M) such that Σ.sub.i=1.sup.N δ.sup.(i)=1 and Σ.sub.i=1.sup.M ∈.sup.(i)=1 and δ.sup.(i)p.sup.J (A)=p.sup.J (A.sup.(i)) for every i=1, . . . , N and ∈.sup.(i)p.sup.J (B)=p.sup.J (B.sup.(i)) for every i=1, . . . , M. Therefore,

p.sup.J(A)+p.sup.J(B)=δ.sup.(1)p.sup.J (A) + . . . +∈.sup.(N)p.sup.J (A)+∈.sup.(1)p.sup.T (B) + . . . +∈.sup.(M)p.sup.J (B)=1

Consequently,

[0136]
J.sup.A.sup.(i)=(Z−T)p.sup.J(A.sup.(i))=(Z−T)δ.sup.(i)p.sup.J (A) for all i=1, . . . , N.

J.sup.B.sup.(i)=(Z−T)p.sup.J (B.sup.(i))=(Z−T)∈.sup.(i)p.sup.J (B) for all i=1, . . . , M.

Because the consumption stage is represented by J=J.sup.A+J.sup.B+J.sup.T, then

J=(Σ.sub.i=1.sup.N(Z−T)δ.sup.(i)p.sup.J(A))+(Σ.sub.i=1.sup.M(Z−T)∈.sup.(i)p.sup.J(B))+(Σ.sub.i=1.sup.S γ.sup.(i)(1−μ.sup.(i)(t))T).

Further, J can be written as the sum of S components J.sup.(1), . . . , J.sup.(S), with proportion factors ϕ.sup.(1), . . . , ϕ.sup.(S). As such, any k-th component can be written as,

[00012]$J^{\left(k\right)}={\varphi }^{\left(k\right)}J={\varphi }^{\left(k\right)}\left[\left(\underset{i=1}{\overset{N}{.Math.}}\left(Z-T\right){\delta }^{\left(i\right)}{p}^J\left(A\right)\right)+\left(\underset{i=1}{\overset{M}{.Math.}}\left(Z-T\right){ϵ}^{\left(i\right)}{p}^J\left(B\right)\right)+\left(\underset{i=1}{\overset{S}{.Math.}}{\gamma }^{\left(i\right)}\left(1-{\mu }^{\left(i\right)}\left(t\right)\right)T\right)\right]$

Further still, the consumption stage J comprises a consumption efficiency factor η(t), where the actual used portion E is separated from the remaining wasted portion C. There are S portion factors θ.sup.(1), . . . , θ.sup.(S), such that Σ.sub.j=1.sup.S η.sup.(j) (t)=Σ.sub.j=1.sup.S θ.sup.(j)η(t). Throughout the example, S is used as an arbitrary index to denote the number of J's. In earlier derivations, the system considered 4 J's, such as (1) residential homes, (2) commercial, (3) industrial businesses, and (4) transportation, as shown in FIGS. 2 and 4, but S could also represent, beyond J's each constituent or disaggregated R's, across X.sub.j and G.sub.j, but also aggregated locations p's, to show the energy flow and subsequent optimization in any given time. J variable could thus be replaced by other such variables, do derive energy flow, and related maximization and minimization of artifacts, systems, or regions from microseconds to year-long time intervals.

[0137] In conclusion, the following series of equations are provided.

[00013]$Z=E+C$$E=\underset{k=1}{\overset{S}{.Math.}}{E}^{\left(k\right)}=\underset{k=1}{\overset{S}{.Math.}}{\eta }^{\left(k\right)}\left(t\right)J=\underset{k=1}{\overset{S}{.Math.}}{\theta }^{\left(k\right)}\eta \left(t\right)J$$C=\underset{k=1}{\overset{S}{.Math.}}{C}^{\left(k\right)}=\underset{k=1}{\overset{N}{.Math.}}\left(\left(1-{\eta }^{\left(k\right)}\left(t\right)\right)J+T^{C^{\left(k\right)}}\right)$

Therefore, to maximize the outflow of utilized energy E, consumption efficiency factor η.sup.(k)(t) should be maximized.

[0138] The optimization solution shown above may be visualized in any of the visualization methods provided in regards to the above-described optimization and visualization embodiments. Further, any combination of the visualizations of the optimization problems may be displayed. In some embodiments, the optimization methods may be displayed together.

[0139] FIG. 32 depicts a flow chart 3200 showing an exemplary process for embodiments of the invention. At step 3202, data from the data sources 138 are received. The received data may be indicative of the source of the generated energy and the amount of generated energy. The received data may also be indicative of an amount of energy consumed and wasted by each energy consumption method described above. For example, the energy source may be solar power and the energy consumption may be an office building or a residence. In another example, the energy source may be fossil fuels and the consumption may be transportation. Any energy source and energy consumption methods may be analyzed using the algorithms and displayed using the visualization described herein.

[0140] At step 3204, the energy that is used, wasted, and the sources and processes of using and wasting the energy are determined. In some embodiments, the energy used and wasted is determined for each energy consumption process such that the process efficiencies can be determined and compared based on each energy source. Further, process costs for both the economy and the environment may be determined. The environmental and economic cost of benign and harmful resources as well as wasted and used resources may be determined.

[0141] At step 3206, allocation of resources may be optimized based on the source of the energy and how the energy is being used and wasted. Determine the most optimal allocation of resources within a defined system to minimize the use of harmful resources and maximize the use of benign resources using the optimization algorithms described above. Further, the processes may be evaluated for find the most effective allocation of resources to limit waste and maximize benign resources. In some embodiments, the cost of the resources is taken into account and a weighted system process economic cost and environmental cost is determine. An allocation of resources maximizing benign resources and minimizing harmful resources may be determined. In some embodiments, wasted resources may be minimized and in some embodiments, monetary costs and savings may be minimized and maximized or otherwise weighted against along with the benign and harmful resources to provide environmental and economical solutions.

[0142] At step 3208, the visualization of the source data and optimized allocation of resources is displayed via the visualizations described above. The visualizations provide two- and three-dimensional data visualization via the visual three-dimensional cubes and two-dimensional maps described above. The user may customize any of the above-described visualizations to present customized past, present, and future visualizations based on the data sources, geospatial map information, transition processes, and the user input parameters.

[0143] At step 3210, the amount of P can be used to “purchase” items, v (representing reallocations and transformations of energy sources in different ways), thereby changing the system parameters in meeting the desired conditions (i.e., p(A) and p.sup.T (A) goes to 1, which implies p(B) and p.sup.T (B) goes to 0, and η goes to 1, such that μ goes to 0). The price to purchase may be maximized, the reduction of harmful inputs, may be minimized or any combination in between.

[0144] As described in reference to FIG. 29 above, each user can understand the emissions Φ given their energy source visualization, and chose to either escrow the monetary amount, P(Φ) or P(B) in any given t in any currency to either (i) purchase relevant v solutions knowing the P(v) of each v solution, (ii) invest P(Φ) amount in community-level projects, either as a trackable donation or purchasing a security in an investment project, such as green bonds, (iii) make charitable donations of P(Φ) to other S entities, through unique alphanumeric IDs or (iv) buy carbon offsets or renewable energy credits in the amount of P(Φ) and get a discount code (another unique alphanumeric ID) from participating retailers.

[0145] As shown in FIG. 29, the first option can be for an S user to offset $20 by using the system to match the $20 to a set of v's to “purchase” at least a portion, or a single v solution that their R, Q, X.sub.j matches that of the v. A v-card may be generated showing the monetary value, as well as savings Sv over time, including any δϕ.sub.v reduction; thereby optimizing the J.sup.A, J.sup.B and J.sup.T parameters of the user S. The v solutions provider may want to verify and give discount to the S consumers who have completed offsets or gone carbon neutral by checking if their unique alphanumeric code has been already spent from their account or not, to limit double spending. If already spent, a message may reiterate the note: “Thanks for your previous impact, we have already applied your discount for v purchase (<unique code>)” If not spent, says “Thanks for your impact, enjoy <discount (80%) on v purchase>”

[0146] The second option allows S user to use the system to match the $20 to a set of v's to “purchase” at least a portion, or a single v “project” that their R, p matches that of the v. A v-card may be generated showing the monetary value of this investment, as well as savings Sv over time or interest rate, such as bond coupons, including with it any δϕ.sub.v reduction; thereby optimizing the J.sup.A, J.sup.B and J.sup.T parameters of either the user or the region p. The v solutions provider may want to verify and give return on investments to the S consumers who have completed offsets or gone carbon neutral by checking if their unique alphanumeric code has been already spent from their account or not, to limit double spending. If already spent, a message may reiterate the note: “Thanks for your previous investment, we have already applied your impact funds to v project (<unique code>)” If not spent, says “Thanks for your impact, you now own minibond coupon <uniquecode>”

[0147] A third option can allow use of the system to match the $20 to a set of other S donations, to give away, at least a portion of their offset to another entity. A list may be generated showing the monetary value of this donation, as well as savings S.sub.v over time of the receiving entity, including with it any δϕ.sub.v reduction; thereby optimizing the J.sup.A, J.sup.B and J.sup.T parameters of the receiving user or entity. The donation receiver may want to verify through the system that the S consumers who have completed offsets or gone carbon neutral by checking if their unique alphanumeric code has been already spent from their account or not, to limit double spending in either other donations, v-purchases, bond investments or any other such system processes. If already spent, a message may reiterate the note: “Thanks for your previous donation, we have already applied your impact (<unique code>)” If not spent, says “Thanks for your impact <alphanumericode>”

[0148] A fourth option allows the S users to use the system to match the $20 to a combination of above processes, to get a discount code for any other entity who are not any of the above. A new alphanumeric code may be uniquely generated showing the monetary value of this donation, as well as savings S.sub.v through a discount percentage applied over time given by the retailer entity, including with it any δϕ.sub.v reduction; thereby optimizing the J.sup.A, J.sup.B and J.sup.T parameters of the receiving user or entity. The discount provider may want to verify through the system that the S donation receiver/consumers who have completed offsets or gone carbon neutral by checking if their unique alphanumeric code has been already spent from their account or not, to limit double spending in either other donations, v-purchases, bond investments or any other such system processes. If already spent, a message may reiterate the note: “Thanks for your previous impact, we have already applied your discount (<unique code>)” If not spent, says “Thanks for your impact, enjoy <discount (40%)>”. Many such transactions between different entities, subsystems, components and variables described above could be facilitated through the step-wise optimization process.

[0149] At step 3212, the source data is continuously received, optimized, and visualized in real-time. In some embodiments, the information is processed and updated continuously, or dynamically, or at any desired period. The optimization algorithms and the data visualization processes may continually receive input data from the data sources and provide current data visualizations and optimized resource allocations.

[0150] Many different arrangements of the various components depicted, as well as components not shown, are possible without departing from the scope of the claims below. Embodiments of the invention have been described with the intent to be illustrative rather than restrictive. Alternative embodiments will become apparent to readers of this disclosure after and because of reading it. Alternative means of implementing the aforementioned can be completed without departing from the scope of the claims below. Certain features and subcombinations are of utility and may be employed without reference to other features and subcombinations and are contemplated within the scope of the claims. Although the invention has been described with reference to the embodiments illustrated in the attached drawing figures, it is noted that equivalents may be employed and substitutions made herein without departing from the scope of the invention as contemplated.

[0151] The following function shows an example by receiving, as input, a generic set of parameters V.sub.i(t) and computes all the other quantities of the system corresponding to such V.sub.i(t). Every V.sub.i(t) comes with a received price, which for this example is generated randomly in the function.

TABLE-US-00009 Algorithm 1 V-EVALUATION Inputs 1: T, Z, p(A.sub.t−1.sup.(1)), . . . , p(A.sub.t−1.sup.(N)), p(B.sub.t−1.sup.(1)), . . . , p(B.sub.t−1.sup.(M)) 2: p.sup.T (A.sub.t−1.sup.(1)), . . . , p.sup.T (A.sub.t−1.sup.(N)), p.sup.T (B.sub.t−1.sup.(1)), . . . , p.sup.T (B.sub.t−1.sup.(M)) 3: p.sup.J (A.sub.t−1.sup.(1)), . . . , p.sup.J (A.sub.t−1.sup.(N)), p.sup.J (B.sub.t−1.sup.(1)), . . . , p.sup.J (B.sub.t−1.sup.(M)) 4: η.sup.(1) (t − 1), . . . , η.sup.(S) (t − 1), φ.sup.(1) (t − 1), . . . , φ.sup.(S) (t − 1) 5: γ.sup.(1) (t − 1), . . . , γ.sup.(S) (t − 1), μ.sup.(1) (t − 1), . . . , μ.sup.(S) (t − 1) 6: Δp(A.sub.t−1.sup.(1)), . . . , Δp(A.sub.t−1.sup.(N)), Δp(B.sub.t−1.sup.(1)), . . . , Δp(B.sub.t−1.sup.(M)) 7: Δp.sup.T (A.sub.t−1.sup.(1)), . . . , Δp.sup.T(A.sub.t−1.sup.(N)), Δp.sup.T (B.sub.t−1.sup.(1)), . . . , Δp.sup.T (B.sub.t−1.sup.(M)) 8: Δp.sup.J (A.sub.t−1.sup.(1)), . . . , Δp.sup.J (A.sub.t−1.sup.(N)), Δp.sup.J (B.sub.t−1.sup.(1)), . . . , Δp.sup.J (B.sub.t−1.sup.(M)) 9: Δη.sup.(1) (t − 1), . . . , Δη.sup.(S) (t − 1), Δφ.sup.(1) (t − 1), . . . , Δφ.sup.(S) (t − 1) 10: Δγ.sup.(1) (t − 1), . . . , Δγ.sup.(S) (t − 1), Δμ.sup.(1) (t − 1), . . . , Δμ.sup.(S) (t − 1) Steps 1: procedure V-EVALUATION (Inputs) 2: A.sub.t.sup.(i) ← p(A.sub.t.sup.(i)) * Z for i = 1, . . . , N 3: [00014]$A_t\leftarrow \underset{i}{\overset{N}{.Math.}}{A}_t^{\left(i\right)}$ 4: B.sub.t.sup.(i) ← p(B.sub.t.sup.(i) * Z for i = 1, . . . , M 5: [00015]$B_t\leftarrow \underset{i}{\overset{M}{.Math.}}{B}_t^{\left(i\right)}$ 6: p.sup.T (A.sub.t.sup.(i)) ← p.sup.T (A.sub.t−1.sup.(i)) + Δp.sup.T (A.sub.t.sup.(i)) for i = 1, . . . , N 7: p.sup.T (B.sub.t.sup.(i)) ← p.sup.T (B.sub.t−1.sup.(i)) + Δp.sup.T (B.sub.t.sup.(i)) for i = 1, . . . , M − 1 8: [00016]$p^T\left(B_t^{\left(M\right)}\right)\leftarrow 1-\left(\underset{i}{\overset{N}{.Math.}}{p}^T\left(A_t^{\left(i\right)}\right)+\underset{i}{\overset{M-1}{.Math.}}{p}^T\left(B_t^{\left(i\right)}\right)\right)$ 9: T.sub.t.sup.A.sup.(i) ← T * p.sup.T (A.sub.t.sup.(i)) for i = 1, . . . , N 10: T.sub.t.sup.B.sup.(i) ← T * p.sup.T (B.sub.t.sup.(i)) for i = 1, . . . , M − 1 11: T.sub.t.sup.B.sup.(M) ← T * p.sup.T (B.sub.t.sup.(M)) 12: p.sup.J (A.sub.t.sup.(i)) ← p.sup.J (A.sub.t−1.sup.(i)) + Δp.sup.J (A.sub.t.sup.(i)) for i = 1, . . . , N 13: p.sup.J (B.sub.t.sup.(i)) ← p.sup.J (B.sub.t−1.sup.(i)) + Δp.sup.J (B.sub.t.sup.(i)) for i = 1, . . . , M − 1 14: [00017]$p^J\left(B_t^{\left(M\right)}\right)\leftarrow 1-\left(\underset{i}{\overset{N}{.Math.}}{p}^J\left(A_t^{\left(i\right)}\right)+\underset{i}{\overset{M-1}{.Math.}}{p}^J\left(B_t^{\left(i\right)}\right)\right)$ 15: J.sup.A.sup.t.sup.A.sub.i ← (Z − T) * p.sup.J (A.sub.t.sup.(i)) for i = 1, . . . , N 16: J.sup.B.sup.t ← (Z − T) * p.sup.J (B.sub.t.sup.(i)) for i = 1, . . . , M 17: [00018]$J^{A_t}\leftarrow \underset{i}{\overset{N}{.Math.}}{J}^{A_t^{\left(\dot{\mathrm{.Math.}}\right)}}$ 18: [00019]$J^{B_t}\leftarrow \underset{i}{\overset{M}{.Math.}}{J}^{B_t^{\left(i\right)}}$ 19: γ.sup.(i) (t) ← γ.sup.(i) (t − 1) + Δγ.sup.(i) (t − 1) for i = 1, . . . , S 20: T.sub.t.sup.(i) ← T * γ.sup.(i) (t) for i = 1, . . . , S 21: μ.sup.(i) (t) ← μ.sup.(i) (t − 1) − Δμ.sup.(i) (t − 1) for i = 1, . . . , S 22: T.sub.t.sup.C.sup.(i) ← T.sub.t.sup.(i) * μ.sup.(i) (t) for i = 1, . . . , S 23: [00020]$T_t^C\leftarrow \underset{i}{.Math.}{T}_t^{C^{\left(\dot{\mathrm{.Math.}}\right)}}$ 24: J.sub.t.sup.T.sup.T ← T.sub.t.sup.(i) * (1 − μ.sup.(i) (t)) for i = 1, . . . , S 25: J.sub.t ← J.sup.A.sup.t + J.sup.B.sup.t + J.sub.t.sup.T 26: η.sup.(i) (t) ← η.sup.(i) (t − 1) + Δη.sup.(i) (t − 1) 27: [00021]$\eta \left(t\right)\leftarrow \underset{i}{.Math.}{\eta }^{\left(i\right)}\left(t\right)$ 28: E.sub.t.sup.(i) ← η.sup.(i) (t) * J.sub.t 29: C.sub.t.sup.(i) ← (1 − η.sup.(i) (t)) * J.sub.t.sup.(i) + T.sup.C.sup.(i) 30: T.sub.t.sup.BC ← T.sub.t.sup.C * p.sup.T (B.sub.t) 31: J.sup.BCN (t) ← (1 − η (t)) * J.sub.t.sup.B.sup.t 32: J.sup.BN (t) ← η (t) * J.sub.t.sup.B 33: Φ (B.sub.t) ← m * (T.sup.BC(t) + T.sup.BCN(t) + T.sup.BN(t) + J.sup.BCN(t) + J.sup.BN(t) where m = 2 34: Φ (B.sub.t) ← Σ Φ (B.sub.t) 35: P(ν.sub.t) ← Input pricing (cost, and savings) 36: P(B.sub.t) ← n* Φ (B.sub.t) where n = 1.5 37: P (ν.sub.i.sup.B (t)) ← P(B.sub.t) − P (ν.sub.i (t)) 38: T.sup.BC (t), J.sup.BCN (t), J.sup.BN (t), Φ (B.sub.t), P (ν.sub.i (t)), P(B.sub.t), P (ν.sub.i.sup.B (t))] 39: ν.sub.t ← List (p (A.sub.t.sup.(1)), . . . , p (A.sub.t.sup.(N)), p (B.sub.t.sup.(1)), . . . , p (B.sub.t.sup.(M)) 40: p.sup.T(A.sub.t.sup.(1)), . . . , p.sup.T(A.sub.t.sup.(N)), p.sup.T(B.sub.t.sup.(1)), . . . , p.sup.T(B.sub.t.sup.(M)), p.sup.J(A.sub.t.sup.(1)), . . . , p.sup.J(A.sub.t.sup.(N)), p.sup.J(B.sub.t−1.sup.(1)), . . . , p.sup.J(B.sub.t−1.sup.(M)), η.sup.(1) (t), . . . , η.sup.(S) (t), φ.sup.(1) (t), . . . , φ.sup.(S) (t), γ.sup.(1) (t), . . . , γ.sup.(S) (t), μ.sup.(1) (t), . . . , μ.sup.(S) (t)) 41: return v.sub.t 42: end procedure

[0152] At t=0 we have the above set v.sub.0 of parameters and we aim to produce the set of “all” possible candidates v.sub.t (1), . . . , v.sub.t (k) for t=1. We do so by trying “all” possible values for Δp(A), Δp.sup.T (A), Δp.sup.J (A), Δμ, Δη which we are going to add to the initial parameters obtaining the k new sets of parameters v.sub.t (1), . . . , v.sub.t (k). All the combinations of the 3(N+M−1)+4S parameters can be obtained with a 3(N+M−1)+4S-dimensional grid, that is computed with 3(N+M−1)+4S nested for-loops. The smaller each for-loop step is, the longer it will take to process the algorithm. The following function takes as input the set of parameters chosen for step t=n and returns a table with “all” the possible candidates for the following iteration t=n+1:

TABLE-US-00010 Algorithm 2 TABLE OF CANDIDATES v's Inputs  1: Parameters of the system at time t − 1, changes of parameters at time t, steps for the for-loops Steps  1: procedure candidates (Inputs)  2: block of N + M − 1 nested for-loops for parameters p.sub.i (A.sub.t.sup.(1)),. . . , p.sub.i (A.sub.t.sup.(N)), p.sub.i (B.sub.t.sup.(1)), ... , p.sub.i (B.sub.t.sup.(M−1)) do:  3: block of N + M − 1 nested for-loops for parameters p.sub.i.sup.J (A.sub.t.sup.(1)),. . . , p.sub.i.sup.J (A.sub.t.sup.(N)), p.sub.i.sup.J (B.sub.t.sup.(1)), ... , p.sub.i.sup.J (B.sub.t.sup.(M−1)) do:  4: block of N + M − 1 nested for-loops for parameters p.sub.i.sup.T (A.sub.t.sup.(1)),. . . , p.sub.i.sup.T (A.sub.t.sup.(N)), p.sub.i.sup.T (B.sub.t.sup.(1)), ... , p.sub.i.sup.T (B.sub.t.sup.(M−1)) do:  5: block of S nested for-loops for parameters η.sub.i.sup.(1) (t),. . . , η.sub.i.sup.(S)(t) do:  6: block of S nested for-loops for parameters φ.sub.i.sup.(1)(t),. . . , φ.sub.i.sup.(S)(t) do:  7: block of S nested for-loops for parameters γ.sub.i.sup.(1) (t),. . . , γ.sub.i.sup.(S)(t) do:  8: block of S nested for-loops for parameters μ.sub.i.sup.(1)(t),. . . , μ.sub.i.sup.(S)(t) do:  9: New Candidate ← v-evaluation (T, Z, p(A.sub.t−1.sup.(1)), p.sup.T (A.sub.t−1.sup.(1)), p.sup.J (A.sub.t−1.sup.(1)), η.sup.(1)(t − 1), μ.sup.(1)(t − 1), Δp(A.sub.t.sup.(1)), Δp.sup.T (A.sub.t.sup.(1)), Δp.sup.J (A.sub.t.sup.(1)), Δη.sup.(1)(t), Δμ.sup.(1)(t)) 10: end block 11: end block 12: end block 13: end block 14: end block 15: end block 16: end block 17: end procedure

[0153] At every iteration t, the choice of a candidate from the table can be made by placing some conditions on the values of the table in such a way that the system reaches the optimal condition as quickly as possible. For example, if for every iteration t we have the constrains:

P(B(t))≥P(v(t))

−|Φ(B(t)−Φ(B(t−1))| is maximized

[0154] We can also input a given ΔΦ(t), implying Φ to be reduced between two time periods, be it, the current time period and next, or two discrete time periods. Here we may generate the v's, and we want to be able to return a subset of them representing the lowest number of v's that can achieve the ΔΦ(t) reduction specified. Similarly, for the similar given ΔΦ reduction input, the system can also return the smallest summation of P(v) of those v's. The user could specify the lowest sum of P(v) in currency, OR the lowest number (i.e. count) of v's; and they may or may not be unique. Like ΔΦ, the input could be of any parameter, such as Δp(A), Δp.sup.T (A), Δp.sup.J (A), . . . Δμ, Δη, etc.

[0155] Given the initial set of parameters above, we can produce several possible sets of parameters. For example, for each v in the table, we can derive its time-based savings in currency, say ΔS.sub.v(t), or the payback period in time units, the return on investment (ROI) based on ΔS.sub.v(t) and P(v), the Internal rate of return (IRR) or calculate net-present value (NPV) etc. of each individual or collective sets of v's.

[0156] The next function deletes the candidates that do not satisfy the first condition, then sort the table by Φ, obtaining several good choices for v. Among the good candidates we choose the one with greatest E, and smallest B.

TABLE-US-00011 Algorithm 3 PICK CANDIDATE Inputs 1: Parameters of the system at time t − 1, changes of parameters at time t, steps for the for-loops Outputs 1: v.sub.i(t) from the table of possible v.sub.i(t)'s at time t that satisfies some conditions. Steps 1: procedure Pick candidate (Inputs) 2: Keep only candidates v.sub.i(t) of table such that P(v.sub.i.sup.B(t)) ≥ 0 3: Order the table by increasing Φ(v.sub.i(t)), then by decreasing E.sub.t, then by increasing B.sub.t 4: chosen ← first row of the table 5: return chosen 6: end procedure The next function puts together the collection of v.sub.i's chosen from the first to the last iteration.

TABLE-US-00012 Algorithm 4 OPTIMIZATION Inputs  1: Parameters of the system at time t − 1, changes of parameters at time t, initial parameters p(A.sub.0), p.sup.J (A.sub.0), p.sup.T (A.sub.0), η(0), μ(0). Outputs  1: Solution to the Optimization Problem Steps  1: procedure Solution (Inputs)  2: Final solution ← v evaluation (T, Z, p(A.sub.0), p.sup.T (A.sub.0),  p.sup.J (A.sub.0), η(0), μ(0), Δp (A.sub.t) = 0, Δp.sup.T (A.sub.t) = 0,  Δp.sup.J (A.sub.t) = 0, Δη(t) = 0, Δμ(t) = 0)  3: sol ← table of candidate v's (T, Z, p(A.sub.0), p.sup.T (A.sub.0), p.sup.J (A.sub.0), η(0), μ(0), Δp(A.sub.t) = 0, Δp.sup.T (A.sub.t) = 0, Δp.sup.J (A.sub.t) = 0, Δη(t) = 0, Δμ(t) = 0, steps for the for- loops)  4: while condition do  5: Decrease steps for the for-loops  6: sol ← Pick candidate (T, Z, p(A.sub.t−1.sup.(i)) , p.sup.T (A.sub.t−1.sup.(i)), , p.sup.J (A.sub.t−1.sup.(i)), η.sup.(i)(t − 1), μ.sup.(i) (t − 1), Δp(A.sub.t.sup.(i)), Δp.sup.T (A.sub.t.sup.(i)), Δp.sup.J (A.sub.t.sup.(i)), Δη.sup.(i)(t), Δμ.sup.(i)(t))  7: if Number of rows of sol = 0 then  8: condition = FALSE  9: else 10: Update parameters p(A.sub.t−1.sup.(i)), p.sup.T (A.sub.t−1.sup.(i)), p.sup.J (A.sub.t−1.sup.(i)), η.sup.(i)(t − 1), μ.sup.(i)(t − 1) 11: Final solution ← List (Final solution, sol) 12: end if 13: end while 14: return Final solution 15: end procedure