import math import json from datetime import datetime, timedelta from astral import LocationInfo from astral.sun import sun, elevation, azimuth import pytz class BifacialSolarSimulationSnow: def __init__(self): # Panel specifications (Hanersun N-TOPCon 700W) self.panel_power = 700 # watts self.panel_efficiency = 0.225 # 22.5% self.bifaciality_factor = 0.80 # Average of 75-85% self.panel_area = 2.384 * 1.303 # m² # Location: Lago Brown, Chile self.location = LocationInfo("Lago Brown", "Chile", "America/Santiago", -47.3919449, -72.3174973) self.timezone = pytz.timezone('America/Santiago') # Direct Normal Irradiance - constant when sun is shining (ideal conditions) self.dni_clear_sky = 900 # W/m² for clear sky conditions # Ground albedo (reflection factor) - SNOW ALBEDO for winter months self.ground_albedo = 0.8 # Snow albedo (high reflection) def get_solar_position(self, date, hour, minute=0): """Get precise solar position using astronomical calculations""" dt = datetime(date.year, date.month, date.day, hour, minute, 0) dt = self.timezone.localize(dt) # Get solar elevation and azimuth solar_elevation = elevation(self.location.observer, dt) solar_azimuth = azimuth(self.location.observer, dt) return solar_elevation, solar_azimuth def calculate_angle_of_incidence(self, solar_elevation, solar_azimuth, panel_tilt, panel_azimuth): """Calculate angle of incidence between sun and panel surface""" if solar_elevation <= 0: return 0 # Convert to radians se_rad = math.radians(solar_elevation) sa_rad = math.radians(solar_azimuth) pt_rad = math.radians(panel_tilt) pa_rad = math.radians(panel_azimuth) # Standard angle of incidence formula cos_incidence = ( math.sin(se_rad) * math.cos(pt_rad) + math.cos(se_rad) * math.sin(pt_rad) * math.cos(sa_rad - pa_rad) ) return cos_incidence def calculate_bifacial_irradiance(self, solar_elevation, solar_azimuth, panel_tilt, panel_azimuth): """Calculate irradiance with proper bifacial physics""" if solar_elevation <= 0: return 0, 0 # Standard angle of incidence for front side front_cos = self.calculate_angle_of_incidence(solar_elevation, solar_azimuth, panel_tilt, panel_azimuth) # Back side angle of incidence back_azimuth = (panel_azimuth + 180) % 360 back_cos = self.calculate_angle_of_incidence(solar_elevation, solar_azimuth, panel_tilt, back_azimuth) # Determine illumination scenario azimuth_to_front = abs(solar_azimuth - panel_azimuth) if azimuth_to_front > 180: azimuth_to_front = 360 - azimuth_to_front azimuth_to_back = abs(solar_azimuth - back_azimuth) if azimuth_to_back > 180: azimuth_to_back = 360 - azimuth_to_back # Determine primary illumination direction if azimuth_to_front < 90: # Front illumination (±90° from front) front_irr = self.dni_clear_sky * max(0, front_cos) back_irr = 0 elif azimuth_to_back < 90: # Back illumination - SPECIAL BIFACIAL PHYSICS front_irr = 0 if back_cos > 0: # For back illumination, vertical panels have significant advantage at low sun if solar_elevation < 30: # Low sun - bifacial geometry critical if panel_tilt == 90: # Vertical panel - best bifacial geometry bifacial_geometry_factor = 1.2 # 20% bonus for optimal geometry else: # Tilted panels lose efficiency for low sun back-illumination # More tilted = further from vertical = worse bifacial geometry vertical_deviation = abs(90 - panel_tilt) / 90.0 bifacial_geometry_factor = 1.0 - 0.4 * vertical_deviation # Stronger penalty else: bifacial_geometry_factor = 1.0 # High sun - less geometry effect back_irr = self.dni_clear_sky * back_cos * self.bifaciality_factor * bifacial_geometry_factor else: back_irr = 0 else: # Side illumination - minimal direct light front_irr = self.dni_clear_sky * max(0, front_cos) * 0.5 # Reduced efficiency back_irr = self.dni_clear_sky * max(0, back_cos) * self.bifaciality_factor * 0.5 return front_irr, back_irr def calculate_panel_irradiance(self, date, hour, minute, panel_tilt, panel_azimuth, panel_type='tilted'): """Calculate irradiance using simplified but physically accurate approach""" solar_elev, solar_az = self.get_solar_position(date, hour, minute) if solar_elev <= 0: return 0, 0 # No sun if panel_type == 'ew_vertical': # East-West vertical panel east_cos = self.calculate_angle_of_incidence(solar_elev, solar_az, panel_tilt, 90) west_cos = self.calculate_angle_of_incidence(solar_elev, solar_az, panel_tilt, 270) if east_cos > west_cos and east_cos > 0: front_direct = self.dni_clear_sky * east_cos back_direct = self.dni_clear_sky * west_cos * self.bifaciality_factor if west_cos > 0 else 0 else: front_direct = self.dni_clear_sky * west_cos if west_cos > 0 else 0 back_direct = self.dni_clear_sky * east_cos * self.bifaciality_factor if east_cos > 0 else 0 else: # North-South panel with proper bifacial physics front_direct, back_direct = self.calculate_bifacial_irradiance( solar_elev, solar_az, panel_tilt, panel_azimuth) # Diffuse irradiance (sky radiation) diffuse_factor = 0.15 sky_diffuse = self.dni_clear_sky * diffuse_factor front_diffuse = sky_diffuse * (1 + math.cos(math.radians(panel_tilt))) / 2 back_diffuse = sky_diffuse * (1 - math.cos(math.radians(panel_tilt))) / 2 # Ground reflected irradiance - SNOW ALBEDO (high reflection) ground_reflected = self.dni_clear_sky * math.sin(math.radians(solar_elev)) * self.ground_albedo front_ground = ground_reflected * (1 - math.cos(math.radians(panel_tilt))) / 2 back_ground = ground_reflected * (1 + math.cos(math.radians(panel_tilt))) / 2 # Apply bifacial efficiency to back side indirect components angle_from_vertical = abs(90 - panel_tilt) bifacial_efficiency = 1.0 - (angle_from_vertical / 90.0) back_indirect_factor = self.bifaciality_factor * (0.4 + 0.6 * bifacial_efficiency) # Total irradiance front_irradiance = front_direct + front_diffuse + front_ground back_irradiance = back_direct + (back_diffuse + back_ground) * back_indirect_factor return front_irradiance, back_irradiance def calculate_monthly_production(self, panel_tilt, panel_azimuth, panel_type='tilted'): """Calculate monthly energy production for given panel configuration""" monthly_production = {} for month in range(1, 13): # Representative date for month (15th day) date = datetime(2024, month, 15) daily_production = 0 half_hourly_data = [] # Calculate for each half hour of the day (48 intervals) for half_hour in range(48): hour = half_hour // 2 minute = 0 if half_hour % 2 == 0 else 30 front_irr, back_irr = self.calculate_panel_irradiance(date, hour, minute, panel_tilt, panel_azimuth, panel_type) if front_irr > 0 or back_irr > 0: total_irradiance = front_irr + back_irr # Power calculation: Irradiance (W/m²) × Panel area × Panel efficiency power = (total_irradiance / 1000) * self.panel_power # Energy for 30 minutes (0.5 hour) daily_production += power * 0.5 # Calculate bifacial contribution percentage bifacial_contribution = (back_irr / total_irradiance * 100) if total_irradiance > 0 else 0 time_str = f"{hour:02d}:{minute:02d}" half_hourly_data.append({ 'time': time_str, 'solar_elevation': round(self.get_solar_position(date, hour, minute)[0], 1), 'front_irradiance': round(front_irr, 2), 'back_irradiance': round(back_irr, 2), 'total_irradiance': round(total_irradiance, 2), 'power': round(power, 1), 'bifacial_contribution': round(bifacial_contribution, 1) }) else: time_str = f"{hour:02d}:{minute:02d}" half_hourly_data.append({ 'time': time_str, 'solar_elevation': 0, 'front_irradiance': 0, 'back_irradiance': 0, 'total_irradiance': 0, 'power': 0, 'bifacial_contribution': 0 }) # Days in month days_in_month = [31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31][month-1] monthly_energy = daily_production * days_in_month / 1000 # kWh monthly_production[month] = { 'energy_kwh': round(monthly_energy, 1), 'daily_production': round(daily_production, 1), 'half_hourly_data': half_hourly_data } return monthly_production def run_simulation(self): """Run simulation for all configurations""" configurations = [ {'name': 'Vertical North-South', 'tilt': 90, 'azimuth': 0, 'type': 'ns_vertical'}, {'name': 'Vertical East-West', 'tilt': 90, 'azimuth': 90, 'type': 'ew_vertical'}, {'name': 'Tilt 20°', 'tilt': 20, 'azimuth': 0, 'type': 'tilted'}, {'name': 'Tilt 30°', 'tilt': 30, 'azimuth': 0, 'type': 'tilted'}, {'name': 'Tilt 45°', 'tilt': 45, 'azimuth': 0, 'type': 'tilted'}, {'name': 'Tilt 60°', 'tilt': 60, 'azimuth': 0, 'type': 'tilted'}, {'name': 'Tilt 70°', 'tilt': 70, 'azimuth': 0, 'type': 'tilted'} ] results = {} for config in configurations: print(f" - {config['name']}") production = self.calculate_monthly_production(config['tilt'], config['azimuth'], config['type']) results[config['name']] = production return results def generate_html_report(self, results): """Generate HTML report with snow albedo data for May-July""" # Calculate annual totals annual_totals = {} for config_name, data in results.items(): annual_total = sum(data[month]['energy_kwh'] for month in range(1, 13)) annual_totals[config_name] = annual_total # Find best configuration best_config = max(annual_totals, key=annual_totals.get) best_value = annual_totals[best_config] months_czech = ['Leden', 'Únor', 'Březen', 'Duben', 'Květen', 'Červen', 'Červenec', 'Srpen', 'Září', 'Říjen', 'Listopad', 'Prosinec'] months_english = ['January', 'February', 'March', 'April', 'May', 'June', 'July', 'August', 'September', 'October', 'November', 'December'] html_content = f""" Analýza výkonu bifaciálního panelu - SNĚHOVÉ ALBEDO - Lago Brown, Chile

❄️ Analýza výkonu bifaciálního panelu - SNĚHOVÉ ALBEDO

Hanersun N-TOPCon 700W • Lago Brown, Chile • Květen-Červenec

⚠️ SPECIÁLNÍ SIMULACE - SNĚHOVÉ ALBEDO

Ground Albedo: 0.8 (sněhový povrch)

Tato simulace ukazuje výkon panelů s vysokým odrazem od sněhu místo standardního jílu (0.15)

Specifikace panelu a lokace

Panel Hanersun N-TOPCon

Výkon: 700 W
Účinnost: 22.5%
Bifacialita: 75-85%
Rozměry: 2384 × 1303 mm

Lokace Lago Brown

Zeměpisná šířka: -47.39°
Zeměpisná délka: -72.32°
Albedo terénu: 0.8 (sníh)
Klimatická oblast: Patagonská Chile

📅 Měsíční výroba energie - SNĚHOVÉ ALBEDO

""" # Monthly data table - only show May, June, July config_order = ['Vertical North-South', 'Vertical East-West', 'Tilt 20°', 'Tilt 30°', 'Tilt 45°', 'Tilt 60°', 'Tilt 70°'] for month in [5, 6, 7]: # Only May, June, July month_name = months_english[month-1] row_class = 'highlight-row' # Highlight these winter months html_content += f"" for config in config_order: value = results[config][month]['energy_kwh'] col_class = 'highlight-column' if config in ['Vertical North-South', 'Tilt 70°'] else '' if row_class and col_class: col_class += ' highlight-row' elif row_class: col_class = 'highlight-row' html_content += f"" html_content += "" html_content += """
MonthVertical North-SouthVertical East-WestTilt 20°Tilt 30°Tilt 45°Tilt 60°Tilt 70°
{month_name}{value}

⏰ Půlhodinová výroba - ZIMNÍ MĚSÍCE SE SNĚHEM

""" # Month tabs - only May, June, July winter_months = ['Květen', 'Červen', 'Červenec'] winter_month_nums = [5, 6, 7] for i, month_name in enumerate(winter_months): active_class = 'active' if i == 0 else '' html_content += f"" html_content += "
" # Half-hourly data for winter months only for i, month in enumerate(winter_month_nums): month_name = winter_months[i] month_name_en = months_english[month-1] active_class = 'active' if i == 0 else '' html_content += f"""

Půlhodinová výroba - {month_name_en} (15. den) - SNĚHOVÉ ALBEDO

""" # Find times with production half_hourly_data = {} for config in config_order: half_hourly_data[config] = results[config][month]['half_hourly_data'] # Get all daylight times all_times = [] date = datetime(2024, month, 15) for half_hour in range(48): hour = half_hour // 2 minute = 0 if half_hour % 2 == 0 else 30 solar_elev, solar_az = self.get_solar_position(date, hour, minute) if solar_elev > 0: time_str = f"{hour:02d}:{minute:02d}" all_times.append((time_str, solar_elev, solar_az)) if all_times: for time_str, solar_elev, solar_az in all_times: html_content += f"" for config in config_order: power = 0 bifacial_contrib = 0 for time_data in half_hourly_data[config]: if time_data['time'] == time_str: power = time_data['power'] bifacial_contrib = time_data['bifacial_contribution'] break if bifacial_contrib > 0: html_content += f"" else: html_content += f"" html_content += "" else: html_content += "" html_content += "
Time (Elevation/Azimuth)Vertical North-SouthVertical East-WestTilt 20°Tilt 30°Tilt 45°Tilt 60°Tilt 70°
{time_str}
({solar_elev:.1f}°/{solar_az:.0f}°)		{power:.1f} W
({bifacial_contrib:.1f}%)		{power:.1f} W
Žádná výroba v tomto měsíci
" html_content += """

🎨 Legenda barevného označení:

9.7° - Elevace slunce (výška nad horizontem)
39° - Azimut slunce (směr od severu)
7.6% - Bifaciální příspěvek (% ze zadní strany panelu)

❄️ Simulace se sněhovým albedem (0.8) pro zimní měsíce květen-červenec
Datum: 12.08.2025
Poznámka: Vysoké albedo sněhu významně zvyšuje bifaciální příspěvek

""" with open('/home/www/claudev.cz/orb3/solar_panel_report_snow_albedo.html', 'w', encoding='utf-8') as f: f.write(html_content) def main(): simulation = BifacialSolarSimulationSnow() print("❄️ Spouštím simulaci se sněhovým albedem pro zimní měsíce...") print(f"Panel: Hanersun N-TOPCon 700W Bifacial") print(f"Lokace: {simulation.location.latitude}°, {simulation.location.longitude}°") print(f"Ground Albedo: {simulation.ground_albedo} (sníh)") print() # Run simulation print("Running snow albedo simulation...") results = simulation.run_simulation() # Generate HTML report simulation.generate_html_report(results) print() print("✅ HTML report se sněhovým albedem byl vytvořen: solar_panel_report_snow_albedo.html") # Show winter months data print() print("📊 Zimní měsíce se sněhovým albedem:") winter_months = ['Květen', 'Červen', 'Červenec'] for i, month_num in enumerate([5, 6, 7]): print(f" {winter_months[i]}:") for config in ['Vertical North-South', 'Vertical East-West', 'Tilt 45°', 'Tilt 70°']: energy = results[config][month_num]['energy_kwh'] print(f" {config}: {energy} kWh") print() if __name__ == "__main__": results = main()