well-log-toolkit 0.1.146__py3-none-any.whl → 0.1.148__py3-none-any.whl

well_log_toolkit/manager.py

@@ -1618,6 +1618,420 @@ class _ManagerPropertyProxy:
             )
 
 
+class _ManagerMultiPropertyProxy:
+    """
+    Proxy for computing statistics across multiple properties on all wells.
+
+    Supports filter(), filter_intervals(), and sums_avg() methods.
+    Multi-property results nest property-specific stats under property names
+    while keeping common stats (depth_range, samples, thickness, etc.) at
+    the group level.
+    """
+
+    # Stats that are specific to each property (nested under property name)
+    PROPERTY_STATS = {'mean', 'median', 'mode', 'sum', 'std_dev', 'percentile', 'range'}
+
+    # Stats that are common across properties (stay at group level)
+    COMMON_STATS = {'depth_range', 'samples', 'thickness', 'gross_thickness', 'thickness_fraction', 'calculation'}
+
+    def __init__(
+        self,
+        manager: 'WellDataManager',
+        property_names: list[str],
+        filters: Optional[list[tuple]] = None,
+        custom_intervals: Optional[dict] = None
+    ):
+        self._manager = manager
+        self._property_names = property_names
+        self._filters = filters or []
+        self._custom_intervals = custom_intervals
+
+    def __getattr__(self, name: str) -> '_ManagerMultiPropertyProxy':
+        """
+        Attribute access as shorthand for filter().
+
+        Allows: manager.properties(['A', 'B']).Facies.sums_avg()
+        Same as: manager.properties(['A', 'B']).filter('Facies').sums_avg()
+        """
+        # Avoid recursion for private attributes
+        if name.startswith('_'):
+            raise AttributeError(f"'{type(self).__name__}' object has no attribute '{name}'")
+        # Treat as filter
+        return self.filter(name)
+
+    def filter(
+        self,
+        property_name: str,
+        insert_boundaries: Optional[bool] = None
+    ) -> '_ManagerMultiPropertyProxy':
+        """
+        Add a filter (discrete property) to group statistics by.
+
+        Parameters
+        ----------
+        property_name : str
+            Name of discrete property to group by
+        insert_boundaries : bool, optional
+            Whether to insert boundary values at filter transitions
+
+        Returns
+        -------
+        _ManagerMultiPropertyProxy
+            New proxy with filter added
+        """
+        new_filters = self._filters + [(property_name, insert_boundaries)]
+        return _ManagerMultiPropertyProxy(
+            self._manager, self._property_names, new_filters, self._custom_intervals
+        )
+
+    def filter_intervals(
+        self,
+        intervals: Union[str, list, dict],
+        name: str = "Custom_Intervals",
+        insert_boundaries: Optional[bool] = None,
+        save: Optional[str] = None
+    ) -> '_ManagerMultiPropertyProxy':
+        """
+        Filter by custom depth intervals.
+
+        Parameters
+        ----------
+        intervals : str, list, or dict
+            - str: Name of saved intervals to retrieve from each well
+            - list: List of interval dicts [{"name": "Zone_A", "top": 2500, "base": 2700}, ...]
+            - dict: Well-specific intervals {"well_name": [...], ...}
+        name : str, default "Custom_Intervals"
+            Name for the interval filter in results
+        insert_boundaries : bool, optional
+            Whether to insert boundary values at interval edges
+        save : str, optional
+            If provided, save intervals to wells with this name
+
+        Returns
+        -------
+        _ManagerMultiPropertyProxy
+            New proxy with custom intervals set
+        """
+        intervals_config = {
+            'intervals': intervals,
+            'name': name,
+            'insert_boundaries': insert_boundaries,
+            'save': save
+        }
+        return _ManagerMultiPropertyProxy(
+            self._manager, self._property_names, self._filters, intervals_config
+        )
+
+    def sums_avg(
+        self,
+        weighted: Optional[bool] = None,
+        arithmetic: Optional[bool] = None,
+        precision: int = 6
+    ) -> dict:
+        """
+        Compute statistics for multiple properties across all wells.
+
+        Multi-property results nest property-specific stats (mean, median, etc.)
+        under each property name, while common stats (depth_range, samples,
+        thickness, etc.) remain at the group level.
+
+        Parameters
+        ----------
+        weighted : bool, optional
+            Include depth-weighted statistics.
+            Default: True for continuous/discrete, False for sampled
+        arithmetic : bool, optional
+            Include arithmetic (unweighted) statistics.
+            Default: False for continuous/discrete, True for sampled
+        precision : int, default 6
+            Number of decimal places for rounding numeric results
+
+        Returns
+        -------
+        dict
+            Nested dictionary with structure:
+            {
+                "well_name": {
+                    "interval_name": {  # if using filter_intervals
+                        "filter_value": {
+                            "PropertyA": {"mean": ..., "median": ..., ...},
+                            "PropertyB": {"mean": ..., "median": ..., ...},
+                            "depth_range": {...},
+                            "samples": ...,
+                            "thickness": ...,
+                            ...
+                        }
+                    }
+                }
+            }
+
+        Examples
+        --------
+        >>> manager.properties(['PHIE', 'PERM']).filter('Facies').sums_avg()
+        >>> # Returns stats for both properties grouped by facies
+
+        >>> manager.properties(['PHIE', 'PERM']).filter_intervals("Zones").sums_avg()
+        >>> # Returns stats for both properties grouped by custom intervals
+
+        >>> # No filters - compute stats for full well
+        >>> manager.properties(['PHIE', 'PERM']).sums_avg()
+        """
+        result = {}
+
+        for well_name, well in self._manager._wells.items():
+            well_result = self._compute_sums_avg_for_well(
+                well, weighted, arithmetic, precision
+            )
+            if well_result is not None:
+                result[well_name] = well_result
+
+        return _sanitize_for_json(result)
+
+    def _compute_sums_avg_for_well(
+        self,
+        well,
+        weighted: Optional[bool],
+        arithmetic: Optional[bool],
+        precision: int
+    ):
+        """
+        Compute multi-property sums_avg for a single well.
+        """
+        # Check if this well has the required saved intervals (if using saved name)
+        if self._custom_intervals:
+            intervals = self._custom_intervals.get('intervals')
+            if isinstance(intervals, str):
+                # Saved filter name - check if this well has it
+                if intervals not in well._saved_filter_intervals:
+                    return None  # Skip wells that don't have this saved filter
+            elif isinstance(intervals, dict):
+                # Well-specific intervals - check if this well is in the dict
+                if well.name not in intervals and well.sanitized_name not in intervals:
+                    return None  # Skip wells not in the dict
+
+        # Collect results for each property
+        property_results = {}
+
+        for prop_name in self._property_names:
+            try:
+                prop = well.get_property(prop_name)
+
+                # Apply filter_intervals if set
+                if self._custom_intervals:
+                    prop = self._apply_filter_intervals(prop, well)
+                    if prop is None:
+                        continue  # Skip this property if intervals can't be applied
+
+                # Apply all filters
+                for filter_name, insert_boundaries in self._filters:
+                    if insert_boundaries is not None:
+                        prop = prop.filter(filter_name, insert_boundaries=insert_boundaries)
+                    else:
+                        prop = prop.filter(filter_name)
+
+                # Compute sums_avg
+                result = prop.sums_avg(
+                    weighted=weighted,
+                    arithmetic=arithmetic,
+                    precision=precision
+                )
+                property_results[prop_name] = result
+
+            except (PropertyNotFoundError, PropertyTypeError, AttributeError, KeyError):
+                # Property doesn't exist in this well or filter error, skip it
+                pass
+
+        if not property_results:
+            return None
+
+        # If no filters/intervals, return simple merged result (no grouping)
+        if not self._filters and not self._custom_intervals:
+            return self._merge_flat_results(property_results)
+
+        # Merge results: nest property-specific stats, keep common stats at group level
+        return self._merge_property_results(property_results)
+
+    def _apply_filter_intervals(self, prop, well):
+        """
+        Apply filter_intervals to a property if custom_intervals is set.
+
+        Returns None if the well doesn't have the required saved intervals.
+        """
+        if not self._custom_intervals:
+            return prop
+
+        intervals_config = self._custom_intervals
+        intervals = intervals_config['intervals']
+        name = intervals_config['name']
+        insert_boundaries = intervals_config['insert_boundaries']
+        save = intervals_config['save']
+
+        # Resolve intervals for this well
+        if isinstance(intervals, str):
+            # Saved filter name - check if this well has it
+            if intervals not in well._saved_filter_intervals:
+                return None  # Skip wells that don't have this saved filter
+            well_intervals = intervals
+        elif isinstance(intervals, dict):
+            # Well-specific intervals
+            well_intervals = None
+            if well.name in intervals:
+                well_intervals = intervals[well.name]
+            elif well.sanitized_name in intervals:
+                well_intervals = intervals[well.sanitized_name]
+            if well_intervals is None:
+                return None  # Skip wells not in the dict
+        elif isinstance(intervals, list):
+            # Direct list of intervals
+            well_intervals = intervals
+        else:
+            return None
+
+        # Apply filter_intervals
+        return prop.filter_intervals(
+            well_intervals,
+            name=name,
+            insert_boundaries=insert_boundaries,
+            save=save
+        )
+
+    def _merge_flat_results(self, property_results: dict) -> dict:
+        """
+        Merge results when no filters are applied (flat structure).
+
+        Returns a single dict with property-specific stats nested under property
+        names and common stats at the top level.
+
+        Parameters
+        ----------
+        property_results : dict
+            {property_name: sums_avg_result}
+
+        Returns
+        -------
+        dict
+            {
+                "PropertyA": {"mean": ..., "median": ..., ...},
+                "PropertyB": {"mean": ..., ...},
+                "depth_range": {...},
+                "samples": ...,
+                ...
+            }
+        """
+        if not property_results:
+            return {}
+
+        result = {}
+
+        # Add property-specific stats for each property
+        for prop_name, prop_result in property_results.items():
+            if isinstance(prop_result, dict):
+                # Extract property-specific stats
+                prop_stats = {
+                    k: v for k, v in prop_result.items()
+                    if k in self.PROPERTY_STATS
+                }
+                if prop_stats:
+                    result[prop_name] = prop_stats
+
+        # Add common stats from first property
+        first_result = next(iter(property_results.values()))
+        if isinstance(first_result, dict):
+            for k, v in first_result.items():
+                if k in self.COMMON_STATS:
+                    result[k] = v
+
+        return result
+
+    def _merge_property_results(self, property_results: dict) -> dict:
+        """
+        Merge results from multiple properties.
+
+        Nests property-specific stats under property names while keeping
+        common stats at the group level.
+
+        Parameters
+        ----------
+        property_results : dict
+            {property_name: sums_avg_result}
+
+        Returns
+        -------
+        dict
+            Merged result with structure:
+            {
+                "group_value": {
+                    "PropertyA": {"mean": ..., ...},
+                    "PropertyB": {"mean": ..., ...},
+                    "depth_range": {...},
+                    "samples": ...,
+                    ...
+                }
+            }
+        """
+        if not property_results:
+            return {}
+
+        # Use first property result as the structure template
+        first_prop = next(iter(property_results.keys()))
+        first_result = property_results[first_prop]
+
+        return self._merge_recursive(property_results, first_result)
+
+    def _merge_recursive(self, property_results: dict, template: dict) -> dict:
+        """
+        Recursively merge property results following the template structure.
+        """
+        result = {}
+
+        for key, value in template.items():
+            if isinstance(value, dict):
+                # Check if this is a stats dict (has property-specific keys)
+                if any(k in value for k in self.PROPERTY_STATS):
+                    # This is a leaf stats dict - merge property stats here
+                    merged = {}
+
+                    # Add property-specific stats for each property
+                    for prop_name, prop_result in property_results.items():
+                        # Navigate to the same key in this property's result
+                        prop_value = self._get_nested_value(prop_result, key)
+                        if prop_value and isinstance(prop_value, dict):
+                            # Extract property-specific stats
+                            prop_stats = {
+                                k: v for k, v in prop_value.items()
+                                if k in self.PROPERTY_STATS
+                            }
+                            if prop_stats:
+                                merged[prop_name] = prop_stats
+
+                    # Add common stats from the first property
+                    for k, v in value.items():
+                        if k in self.COMMON_STATS:
+                            merged[k] = v
+
+                    result[key] = merged
+                else:
+                    # This is an intermediate nesting level - recurse
+                    # Collect corresponding sub-dicts from all properties
+                    sub_property_results = {}
+                    for prop_name, prop_result in property_results.items():
+                        prop_value = self._get_nested_value(prop_result, key)
+                        if prop_value and isinstance(prop_value, dict):
+                            sub_property_results[prop_name] = prop_value
+
+                    if sub_property_results:
+                        result[key] = self._merge_recursive(sub_property_results, value)
+            else:
+                # Non-dict value, just copy from template
+                result[key] = value
+
+        return result
+
+    def _get_nested_value(self, d: dict, key: str):
+        """Get value from dict, returning None if key doesn't exist."""
+        return d.get(key) if isinstance(d, dict) else None
+
+
 class WellDataManager:
     """
     Global orchestrator for multi-well analysis.
@@ -1710,6 +2124,60 @@ class WellDataManager:
         # Return a proxy that can be used for operations across all wells
         return _ManagerPropertyProxy(self, name)
 
+    def properties(self, property_names: list[str]) -> _ManagerMultiPropertyProxy:
+        """
+        Create a multi-property proxy for computing statistics across multiple properties.
+
+        This allows computing statistics for multiple properties at once, with
+        property-specific stats (mean, median, etc.) nested under property names
+        and common stats (depth_range, samples, thickness, etc.) at the group level.
+
+        Parameters
+        ----------
+        property_names : list[str]
+            List of property names to include in statistics
+
+        Returns
+        -------
+        _ManagerMultiPropertyProxy
+            Proxy that supports filter(), filter_intervals(), and sums_avg()
+
+        Examples
+        --------
+        >>> # Compute stats for multiple properties grouped by facies
+        >>> manager.properties(['PHIE', 'PERM']).filter('Facies').sums_avg()
+        >>> # Returns:
+        >>> # {
+        >>> #     "well_A": {
+        >>> #         "Sand": {
+        >>> #             "PHIE": {"mean": 0.18, "median": 0.17, ...},
+        >>> #             "PERM": {"mean": 150, "median": 120, ...},
+        >>> #             "depth_range": {...},
+        >>> #             "samples": 387,
+        >>> #             "thickness": 29.4,
+        >>> #             ...
+        >>> #         }
+        >>> #     }
+        >>> # }
+
+        >>> # With custom intervals
+        >>> manager.properties(['PHIE', 'PERM']).filter('Facies').filter_intervals("Zones").sums_avg()
+        >>> # Returns:
+        >>> # {
+        >>> #     "well_A": {
+        >>> #         "Zone_1": {
+        >>> #             "Sand": {
+        >>> #                 "PHIE": {"mean": 0.18, ...},
+        >>> #                 "PERM": {"mean": 150, ...},
+        >>> #                 "depth_range": {...},
+        >>> #                 ...
+        >>> #             }
+        >>> #         }
+        >>> #     }
+        >>> # }
+        """
+        return _ManagerMultiPropertyProxy(self, property_names)
+
     def load_las(
         self,
         filepath: Union[str, Path, list[Union[str, Path]]],
@@ -2473,6 +2941,18 @@ class WellDataManager:
             # Delete sources marked for deletion
             well.delete_marked_sources(well_folder)
 
+            # Save filter intervals if any exist
+            if hasattr(well, '_saved_filter_intervals') and well._saved_filter_intervals:
+                import json
+                intervals_file = well_folder / "intervals.json"
+                with open(intervals_file, 'w') as f:
+                    json.dump(well._saved_filter_intervals, f, indent=2)
+            else:
+                # Remove intervals file if no intervals (in case they were deleted)
+                intervals_file = well_folder / "intervals.json"
+                if intervals_file.exists():
+                    intervals_file.unlink()
+
         # Save templates
         if self._templates:
             templates_folder = save_path / "templates"
@@ -2563,6 +3043,20 @@ class WellDataManager:
             for las_file in las_files:
                 self.load_las(las_file, silent=True)
 
+            # Load saved filter intervals if they exist
+            intervals_file = well_folder / "intervals.json"
+            if intervals_file.exists():
+                import json
+                try:
+                    with open(intervals_file, 'r') as f:
+                        saved_intervals = json.load(f)
+                    # Find the well for this folder and set its intervals
+                    well_key = well_folder.name  # e.g., "well_35_9_16_A"
+                    if well_key in self._wells:
+                        self._wells[well_key]._saved_filter_intervals = saved_intervals
+                except Exception as e:
+                    warnings.warn(f"Could not load intervals from {intervals_file}: {e}")
+
         return self
 
     def add_well(self, well_name: str) -> Well:

well_log_toolkit/property.py

@@ -1165,6 +1165,10 @@ class Property(PropertyOperationsMixin):
         new_prop._original_sample_count = len(self.depth)
         new_prop._boundary_samples_inserted = len(new_depth) - len(self.depth)
 
+        # Preserve custom intervals if they exist (from filter_intervals)
+        if hasattr(self, '_custom_intervals') and self._custom_intervals:
+            new_prop._custom_intervals = self._custom_intervals
+
         return new_prop
 
     def filter_intervals(

{well_log_toolkit-0.1.146.dist-info → well_log_toolkit-0.1.148.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: well-log-toolkit
-Version: 0.1.146
+Version: 0.1.148
 Summary: Fast LAS file processing with lazy loading and filtering for well log analysis
 Author-email: Kristian dF Kollsgård <kkollsg@gmail.com>
 License: MIT

{well_log_toolkit-0.1.146.dist-info → well_log_toolkit-0.1.148.dist-info}/RECORD

@@ -1,15 +1,15 @@
 well_log_toolkit/__init__.py,sha256=ilJAIIhh68pYfD9I3V53juTEJpoMN8oHpcpEFNpuXAQ,3793
 well_log_toolkit/exceptions.py,sha256=X_fzC7d4yaBFO9Vx74dEIB6xmI9Agi6_bTU3MPxn6ko,985
 well_log_toolkit/las_file.py,sha256=Tj0mRfX1aX2s6uug7BBlY1m_mu3G50EGxHGzD0eEedE,53876
-well_log_toolkit/manager.py,sha256=vT5-W21PS0XU5M4dSaueqpo2p8pK_tYe47qEVOgxWsQ,117930
+well_log_toolkit/manager.py,sha256=Jx0w1KpxRB2Oex6bC0FQFZOlJkoTKNexW3HMXxndL68,136713
 well_log_toolkit/operations.py,sha256=z8j8fGBOwoJGUQFy-Vawjq9nm3OD_dUt0oaNh8yuG7o,18515
-well_log_toolkit/property.py,sha256=O5Ti5ahWV3CTlBLGZ-ntEIed6GGyzsxnyO_EbYrNLP0,100752
+well_log_toolkit/property.py,sha256=GsiD9c4SfBw8ar7ZJXS0NNejPlpvFRHKck_eBR2lLmo,100965
 well_log_toolkit/regression.py,sha256=JDcRxaODJnFikAdPJyTq8eUV7iY0vCDmvnGufqlojxs,31625
 well_log_toolkit/statistics.py,sha256=_huPMbv2H3o9ezunjEM94mJknX5wPK8V4nDv2lIZZRw,16814
 well_log_toolkit/utils.py,sha256=O2KPq4htIoUlL74V2zKftdqqTjRfezU9M-568zPLme0,6866
 well_log_toolkit/visualization.py,sha256=nnpmFmbj44TbP0fsnLMR1GaKRkqKCEpI6Fd8Cp0oqBc,204716
 well_log_toolkit/well.py,sha256=n6XfaGSjGtyXCIaAr0ytslIK0DMUY_fSPQ_VCqj8jaU,106173
-well_log_toolkit-0.1.146.dist-info/METADATA,sha256=s_S72extRI5wIXZOZ4IIRFkmQPo-Gw5i2udBow2fRG4,63473
-well_log_toolkit-0.1.146.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
-well_log_toolkit-0.1.146.dist-info/top_level.txt,sha256=BMOo7OKLcZEnjo0wOLMclwzwTbYKYh31I8RGDOGSBdE,17
-well_log_toolkit-0.1.146.dist-info/RECORD,,
+well_log_toolkit-0.1.148.dist-info/METADATA,sha256=k7HnPCJ1iza1s8kvTFazQhfvYziyfVQ24JJYxWjrtvQ,63473
+well_log_toolkit-0.1.148.dist-info/WHEEL,sha256=wUyA8OaulRlbfwMtmQsvNngGrxQHAvkKcvRmdizlJi0,92
+well_log_toolkit-0.1.148.dist-info/top_level.txt,sha256=BMOo7OKLcZEnjo0wOLMclwzwTbYKYh31I8RGDOGSBdE,17
+well_log_toolkit-0.1.148.dist-info/RECORD,,